Re: [petsc-users] [petsc-maint] Fwd: DMPlex global to natural problem using DmPlexGetVertexNumbering or DMPlexGlobalToNatural

On 18-11-29 06:13 PM, Matthew Knepley wrote: On Thu, Nov 29, 2018 at 7:40 PM Danyang Su via petsc-maint mailto:petsc-ma...@mcs.anl.gov>> wrote: Dear PETSc developers & users, Sorry to bother you again. I just encounter some difficulties in DMPlex global to natural order. This

Re: [petsc-users] Problem with large grid size

No > On Nov 29, 2018, at 8:54 PM, Fazlul Huq wrote: > > Sorry! I made mistake in running the code. > > It's actually working now until matrix size of 999! > > But when I try to extend to one order more it give me error message: >

Re: [petsc-users] Problem with large grid size

> On Nov 29, 2018, at 8:55 PM, Fazlul Huq wrote: > > So, if I try to run the same code in a super computer then I can run with > even larger matrix? Within reason. But every time you increase n by 10 it will require more than 100 times as much memory so even on the biggest machine you

Re: [petsc-users] Problem with large grid size

Sorry! I made mistake in running the code. It's actually working now until matrix size of 999! But when I try to extend to one order more it give me error message: ** Out

Re: [petsc-users] Problem with large grid size

You need to run it on more processors, this one processor doesn't have enough memory to fit the vectors (which by the way are huge 1,614,907,707,076) Barry > On Nov 29, 2018, at 8:26 PM, Fazlul Huq wrote: > > Thanks. > > I have configured with 64-bit and then when I run, I got the

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

> On Nov 29, 2018, at 4:50 PM, Sajid Ali > wrote: > > Another question I want to ask is about dttol that is used in the Monitor > function as follows : > > 417: dttol = .0001; > 418: PetscOptionsGetReal(NULL,NULL,"-dttol",,NULL); > 419: if (dt < dttol) { > 420: dt *= .999; >

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

Another question I want to ask is about dttol that is used in the Monitor function as follows : 417: dttol = .0001; 418: PetscOptionsGetReal(NULL,NULL,"-dttol",,NULL); 419: if (dt < dttol) { 420: dt *= .999; 421: TSSetTimeStep(ts,dt); 422: } Is this decreasing dt if it's less

Re: [petsc-users] Problem with large grid size

> On Nov 29, 2018, at 1:46 PM, Fazlul Huq via petsc-users > wrote: > > Hello PETSc Developers, > > I am trying to run the code (attached herewith) with the following command > and it works until the size of the matrix is 9X9. But when I try to > run with 99X99 then I got

[petsc-users] Problem with large grid size

Hello PETSc Developers, I am trying to run the code (attached herewith) with the following command and it works until the size of the matrix is 9X9. But when I try to run with 99X99 then I got weird result. The command is: ./poisson_m -n 99 -pc_type hypre -pc_hypre_type

[petsc-users] SPAM on mailing lists

With the increase in SPAM on the mailing lists (petsc-dev and petsc-users) - I've disabled posting by non-subscribers [which was previously enabled for user convince]. My current plan is to use a white-list for valid non-subscribers. This way - valid non-subscribers can continue participating in

Re: [petsc-users] Solving complex linear sparse matrix in parallel + external library

It seems to be working now, $ mpiexec -n 4 ./solveSystem -ksp_type gmres -pc_type lu -pc_factor_mat_solver_type mumps -ksp_max_it 1 -ksp_monitor_true_residual 0 KSP preconditioned resid norm 2.536556771280e+03 true resid norm 4.215007800243e+05 ||r(i)||/||b|| 1.e+00 1 KSP

Re: [petsc-users] A question regarding a potential use case for DMNetwork

A simple example of DMNetwork, an electric circuit: petsc/src/ksp/ksp/examples/tutorials/network/ex1.c Hong On Wed, Nov 28, 2018 at 5:08 PM Matthew Knepley via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: On Wed, Nov 28, 2018 at 3:51 PM Markus Lohmayer via petsc-users

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

Hi Edoardo, I also have the Intel Parallel Studio XE compilers and MPI installed, and I use it to build PETSc as follows. # Either add these to your .bashrc or run them on the command line before beginning the PETSc installation source /opt/intel/parallel_studio_xe_2018/bin/psxevars.sh intel64

Re: [petsc-users] IS Invert Not Permutation

Am 28.11.18 um 18:12 schrieb Matthew Knepley: > On Wed, Nov 28, 2018 at 12:06 PM Florian Lindner > wrote: > > Hey, > > thanks for your quick reply! > > As far as I understand how MatSetLocalToGlobalMapping works, there is no > way to create such a