Re: [petsc-users] solve problem with pastix

2019-11-05 Thread Smith, Barry F. via petsc-users
Google finds this https://gforge.inria.fr/forum/forum.php?thread_id=32824_id=599_id=186 > On Nov 5, 2019, at 7:01 PM, Matthew Knepley via petsc-users > wrote: > > I have no idea. That is a good question for the PasTix list. > > Thanks, > > Matt > > On Tue, Nov 5, 2019 at 5:32

Re: [petsc-users] solve problem with pastix

2019-11-05 Thread Matthew Knepley via petsc-users
I have no idea. That is a good question for the PasTix list. Thanks, Matt On Tue, Nov 5, 2019 at 5:32 PM hg wrote: > Should thread affinity be invoked? I set -mat_pastix_threadnbr 1 and also > OMP_NUM_THREADS to 1 > > Giang > > > On Tue, Nov 5, 2019 at 10:50 PM Matthew Knepley wrote:

Re: [petsc-users] VecDuplicate for FFTW-Vec causes VecDestroy to fail conditionally on VecLoad

2019-11-05 Thread Zhang, Junchao via petsc-users
Fixed in https://gitlab.com/petsc/petsc/merge_requests/2262 --Junchao Zhang On Fri, Nov 1, 2019 at 6:51 PM Sajid Ali mailto:sajidsyed2...@u.northwestern.edu>> wrote: Hi Junchao/Barry, It doesn't really matter what the h5 file contains, so I'm attaching a lightly edited script of

Re: [petsc-users] Select a preconditioner for SLEPc eigenvalue solver Jacobi-Davidson

2019-11-05 Thread Jose E. Roman via petsc-users
JD sets STPRECOND at EPSSetUp(), if it was not set before. So I guess you need to add -st_type precond on the command line, so that it is set at EPSSetFromOptions(). Jose > El 5 nov 2019, a las 18:13, Fande Kong escribió: > > How about I want to determine the ST type on runtime? > >

Re: [petsc-users] Select a preconditioner for SLEPc eigenvalue solver Jacobi-Davidson

2019-11-05 Thread Fande Kong via petsc-users
How about I want to determine the ST type on runtime? mpirun -n 1 ./ex3 -eps_type jd -st_ksp_type gmres -st_pc_type none -eps_view -eps_target 0 -eps_monitor -st_ksp_monitor ST is indeed STPrecond, but the passed preconditioning matrix is still ignored. EPS Object: 1 MPI processes type:

Re: [petsc-users] Select a preconditioner for SLEPc eigenvalue solver Jacobi-Davidson

2019-11-05 Thread Jose E. Roman via petsc-users
Currently, the function that passes the preconditioner matrix is specific of STPRECOND, so you have to add ierr = STSetType(st,STPRECOND);CHKERRQ(ierr); before ierr = STPrecondSetMatForPC(st,B);CHKERRQ(ierr); otherwise this latter call is ignored. We may be changing a little bit the way in

Re: [petsc-users] SLEPC no speedup in parallel

2019-11-05 Thread Perceval Desforges via petsc-users
Hello, After carefully looking at the example in the tutorial suggested in section 3.4.5 of the manual, I managed to determine that the problem was caused by me calling EPSSetFromOptions(eps) after EPSKrylovSchurSetPartitions(eps,size) and not before. It now works fine. Sorry! Best regards,