Re: [petsc-users] Some clarifications about TS ex3.c

2018-11-27 Thread Smith, Barry F. via petsc-users
PETSc uses reference counting to track when an object is no longer needed and can be freed. Thus XXXDestroy(&xxx) may not actually destroy an object it just decreases the reference count by 1 (and destroys the object if the reference count is 0). 1) TSSetIJacobian() increases the referenc

Re: [petsc-users] Fortran interface of some petsc routines seem to be missing

2018-11-27 Thread Smith, Barry F. via petsc-users
PETSc uses the Sowing packages bfort tool for automatically generating "Fortran stub functions" and interface definitions. There are certain C functions that bfort cannot handling including (at least) 1) functions with character string arguments 2) functions with function pointer argument

[petsc-users] Some clarifications about TS ex3.c

2018-11-27 Thread Sajid Ali via petsc-users
Hi, I wanted to ask a few questions about the TS example ex3.c when using the ifunc option 190: Mat J; 192: RHSMatrixHeat(ts,0.0,u,A,A,&appctx);193: MatDuplicate

Re: [petsc-users] Example 23 of ksp problems

2018-11-27 Thread Smith, Barry F. via petsc-users
This is an ad hoc way of checking that the error from the linear solve is "reasonable". Note that PETSC_MACHINE_EPSILON depends on the precision of the floating point used so for more precise floating point (in going from half precision to quad precision) we expect the error to be small

Re: [petsc-users] PetscPartitioner is missing for fortran

2018-11-27 Thread Jiaoyan Li via petsc-users
Matt, Thanks a lot for your reply, and I look forward to your further info. Best, Jiaoyan From: Matthew Knepley Date: Tuesday, November 27, 2018 at 09:06 To: Jiaoyan Li Cc: "Klaij, Christiaan" , PETSc Subject: Re: [petsc-users] PetscPartitioner is missing for fortran On Tue, Nov 27, 2018 at

Re: [petsc-users] Compile petsc using intel mpi

2018-11-27 Thread Balay, Satish via petsc-users
If you already have MPI from intel - use the MPI compiler wrappers instead of icc etc. ./configure --with-cc=mpiicc --with-cxx=mpiicpc --with-fc=mpiifort PETSC_ARCH=arch-intel-mpi Satish On Tue, 27 Nov 2018, Zhang, Junchao via petsc-users wrote: > > > On Tue, Nov 27, 2018 at 5:25 AM Edoard

Re: [petsc-users] Compile petsc using intel mpi

2018-11-27 Thread Zhang, Junchao via petsc-users
On Tue, Nov 27, 2018 at 5:25 AM Edoardo alinovi via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: Dear users, I have installed intel parallel studio on my workstation and thus I would like to take advantage of intel compiler. Before messing up my installation, have you got some guideli

Re: [petsc-users] PetscPartitioner is missing for fortran

2018-11-27 Thread Klaij, Christiaan via petsc-users
Personally, I never try to fix these things by myself, it's the job of the petsc developers, they know best and can make the fix available for all users. As a user, I just give feedback whenever I encounter a problem (which isn't often). Chris dr. ir. Christiaan Klaij | Senior Researcher | Rese

Re: [petsc-users] PetscPartitioner is missing for fortran

2018-11-27 Thread Jiaoyan Li via petsc-users
Thanks, Chris, for your kind reply. Yes, I believe the mailing list is the most convenient way for Petsc's to communicate with each other. Also, I am wondering if you may happen to use "PetscPartitioner" before? Have you ever tried to fix Petsc-Fortran problem by yourself before? I believe Pets

[petsc-users] IS Invert Not Permutation

2018-11-27 Thread Florian Lindner via petsc-users
Hello, I have a range of local input data indices that I want to use for row indexing, say { 3, 4, 6}. For that, I create a matrix with a local number of rows of 3 and map the indices {3, 4, 5} to these rows. I create an index set: ISCreateGeneral(comm, myIndizes.size(), myIndizes.data(), PET