Re: [petsc-users] Change the matrix row order

When and where do you want to change the row order? For PETSc and external factorization based solvers this is handled automatically, no need to do anything. The ordering of matrix rows corresponds to the ordering of the vectors; almost always one partitions and orders the vectors

Re: [petsc-users] PETSc 3.9 release

Randall : > Can you explain how VECNODE is different from normal vectors and where it > might be useful? > Vector uses on-node shared memory to store its entries. http://www.mcs.anl.gov/petsc/documentation/changes/39.html This is an experimental work for using mpi3 support of shared memory. The

Re: [petsc-users] PETSc version 3.9 with Fortran and DMPlex questions

Found my way to the answer of the first question: ksp/ksp/examples/tutorials/ex13f90.F90 With respect to question 2... I'm still confused about this behaviour of the DMPlex; or I am missing something. Best regards. Bernardo On Tue, Apr 10, 2018 at 3:19 PM, Bernardo Rocha <

[petsc-users] PETSc version 3.9 with Fortran and DMPlex questions

Hi everyone, I have two questions. (Q1) The first one is related to the new PETSc 3.9. I've been using Fortran and using the "userctx" to wrap some data for the function calls within my program. Just like : call func(...,userctx) and retrieving data as: PetscFortranAddr userctx(*) b =

Re: [petsc-users] Problems with DMDAVecGetArrayF90 + Intel

This should generate an SSE2 binary: 'COPTFLAGS=-g', 'FOPTFLAGS=-g', This should generate a KNL binary: 'COPTFLAGS=-g -xMIC-AVX512 -O3', 'FOPTFLAGS=-g -xMIC-AVX512 -O3', This should generate a SSE2 binary that also supports CORE-AVX2 dispatch. '--COPTFLAGS=-g

Re: [petsc-users] Problems with DMDAVecGetArrayF90 + Intel

On Tue, Apr 10, 2018 at 4:39 PM, Satish Balay wrote: > On Tue, 10 Apr 2018, Jeff Hammond wrote: > > > This should generate an SSE2 binary: > > > > 'COPTFLAGS=-g', > > 'FOPTFLAGS=-g', > > > > This should generate a KNL binary: > > > > 'COPTFLAGS=-g -xMIC-AVX512