Franz Georg Fuchs

When correcting it to: ierr = KSPSetOperators(ksp,M,P);CHKERRQ(ierr); it fails with ||--------------------------------------- || number of particles : 9261 || std of Gaussian (sigma) : 0.100000 || overlap ratio (h/sigma) : 1.000000...

the 2d code compiles after the same change (ierr = KSPSetOperators(ksp,M,P);CHKERRQ(ierr);)

actually, the 2d code also produces: ||--------------------------------------- || number of particles : 49729 || std of Gaussian (sigma) : 0.005000 || overlap ratio (h/sigma) : 0.900000 || non-overlapping subdomain :...

ping