Re: [petsc-users] Implementing TS routine

Thanks Zhang for your answer, I ended up getting a compiling and running TS routine... that does not give me the answers, and i am not sure what i am doing wrong, My TS code so far looks like this: (...initialization...) call TSCreate(PETSC_COMM_WORLD,ts,ierr) call TSSetProblemType(ts,TS_NONLIN

Re: [petsc-users] How to create a mapping of global to local indices?

Note that number ingof vertices/nodes whatever is just an arbitrary label; no one is better or more fundamental than another. In your code's interaction with PETSc vectors and matrices you must use the contiguous numbering, in your mesh management code you can do whatever you want. There ar

Re: [petsc-users] Issue with DMPlexRestoreCone from Fortran

Found the change... it was in Commit 6d563f2afc7bcefff9dc7869193a6debf1d639a5 Date: Wed Jun 26 12:32:01 2019 + remove remnants of non-existing DMPlexRestore{Cone,ConeOrientation,Support} The patch removed the fortran stubs but not the C functions. In this case I guess that the fortran

Re: [petsc-users] Communication during MatAssemblyEnd

Thank you Richard for your explanation. I first changed the way I was running the code. In the machine there's SLURM and I was using "srun -n ./my_program.x". I've seen more than 60% improvement in execution time by just running with "srun -n --ntasks-per-core=1 --cpu-bind=cores> ./my_program.x"

Re: [petsc-users] snes/ex19 issue with nvprof

My branch mark/gamg-fix-viennacl-rebased fixes this problem for me. Please give it a try. Thanks, Mark On Tue, Jul 2, 2019 at 10:41 PM Xiangdong via petsc-users < petsc-users@mcs.anl.gov> wrote: > Hello everyone, > > When I run the ex19 with cuda like this: > mpiexec -np 4 ./ex19 -da_refine 5 -s

Re: [petsc-users] How to create a mapping of global to local indices?

Thank you, Mark. My question is if I have some node indices like [0, 1, 3, 8, 9] handled by the current processor, is that necessary to renumber these indices to something like 0 to 4? or is there any alternative? Dongyu From: Mark Adams Sent: Wednesday, July

Re: [petsc-users] snes/ex19 issue with nvprof

I get this error on normal runs on SUMMIT when I have an empty process. It does not happen all of the time but it does look to be deterministic. I'm debugging it now. On Tue, Jul 2, 2019 at 10:41 PM Xiangdong via petsc-users < petsc-users@mcs.anl.gov> wrote: > Hello everyone, > > When I run the e

Re: [petsc-users] How to create a mapping of global to local indices?

PETSc matrices and vectors are created with a local size n or global size N and PETSC_DECIDE instead of n. The global PETSc indices are ordered from 0 to n_0 - 1 where n_0 is the number of equations on process 0. This numbering continues for all processes. You can use: PetscErrorCode

[petsc-users] How to create a mapping of global to local indices?

Hi, I am running a FEM program using PETSc. In the beginning, the mesh is partitioned in gmsh, and we read the partitioned mesh using our own reader. Now my question is: How can I get a global to local indices mapping? and do I need to renumber the indices after I read the partitioned mesh?

[petsc-users] Issue with DMPlexRestoreCone from Fortran

Dear Matt, Dear Petsc After my recent pull of petsc master I encounter test fails that use dmplex meshes. I attached a simple script to show an issue where DMPlexRestoreCone mangles the cone relationships in the DM. The behaviour happens with gcc 7.4.0 aswell as icc 19.0.3.199 @ current petsc/ma