[petsc-users] Use kspsolve repeatly

2013-07-31 Thread 丁老师
I need to use the kspsolve(A,b,x) repeatly in my code in the following style. but i have noticed from the system load monitor that during the code running, it allocate new memeory every step. I use the default setting for the PC. Could you please told me how to resolve this problem. for (

Re: [petsc-users] Is it possible to build PETSc project in Visual Studio?

2013-07-31 Thread Barry Smith
I have forwarded this to the pflotran developer who uses developer studio, hopefully he will be able to answer your questions. Barry On Jul 31, 2013, at 6:35 PM, Danyang Su wrote: > Hi All, > > I am currently developing codes under Windows utilizing intel visual fortran > + visual st

Re: [petsc-users] Segmentatoin fault

2013-07-31 Thread 丁老师
i get the infor from the system load monitor 在 2013-07-31 23:36:38,"Jed Brown" 写道: >丁老师 writes: > >> I write a FVM code for the convection diffusion problem on a 2D >> cartesian grid. the code can run correctly with a grid of >> 720*720. and it use only 26% of the system memory(8G) durin

Re: [petsc-users] PETSc problems with MVAPICH2

2013-07-31 Thread Barry Smith
http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind On Jul 31, 2013, at 10:56 AM, Jed Brown wrote: > wadud.m...@awe.co.uk writes: > >> Hello, >> >> I have built PETSc 3.4.2 with MVAPICH2 1.9b and my application code >> crashes when I spawn it with 16 processes; it works with 4 and

Re: [petsc-users] Segmentatoin fault

2013-07-31 Thread Barry Smith
http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind On Jul 31, 2013, at 7:46 AM, 丁老师 wrote: > I write a FVM code for the convection diffusion problem on a 2D cartesian > grid. > the code can run correctly with a grid of 720*720. and it use only 26% of the > system memory(8G) during

Re: [petsc-users] how to select values from certain positions of a vec from different process to process 0.

2013-07-31 Thread Barry Smith
On Jul 31, 2013, at 7:55 AM, 丁老师 wrote: > Does the function VecScatterCreateToZero resolve this problem? > What problem? I told you how to do it, there is no problem. VecScatterCreateToZero() moves ALL values to process 0 which for large problems is not scalable, so if you only want som

Re: [petsc-users] PETSc problems with MVAPICH2

2013-07-31 Thread Jed Brown
wadud.m...@awe.co.uk writes: > Hello, > > I have built PETSc 3.4.2 with MVAPICH2 1.9b and my application code > crashes when I spawn it with 16 processes; it works with 4 and 8 > processes. This occurs at the MatAssemblyEnd subroutine when the > non-local values are broadcast to other processes. H

[petsc-users] PETSc problems with MVAPICH2

2013-07-31 Thread Wadud.Miah
Hello, I have built PETSc 3.4.2 with MVAPICH2 1.9b and my application code crashes when I spawn it with 16 processes; it works with 4 and 8 processes. This occurs at the MatAssemblyEnd subroutine when the non-local values are broadcast to other processes. Has anyone else had issues with MVAPICH

Re: [petsc-users] FEM on 2D poisson equation

2013-07-31 Thread Jed Brown
Olivier Bonnefon writes: > Hello, > > You are right. I have to define the Jacobian function of the variational > formulation. I'm using snes and the petsFem struc (like in ex12). > > I need some information about the petscFem struct. I didn't find any > document about that, is there one ? > > T

Re: [petsc-users] Segmentatoin fault

2013-07-31 Thread Jed Brown
丁老师 writes: > I write a FVM code for the convection diffusion problem on a 2D > cartesian grid. the code can run correctly with a grid of > 720*720. and it use only 26% of the system memory(8G) during its > running. How are you testing the total application memory usage? > but if i increase

Re: [petsc-users] preallocation of mpiaij matrix

2013-07-31 Thread Jed Brown
Hoang Giang Bui writes: > Hi > > I want to consult the optimized way to allocate the mpiaij matrix for > FEM application. Typically, elements in each process contribute to > global matrix partially, therefore in each process I can define the > number of nonzeros in according row locally but no

Re: [petsc-users] FEM on 2D poisson equation

2013-07-31 Thread Olivier Bonnefon
Hello, You are right. I have to define the Jacobian function of the variational formulation. I'm using snes and the petsFem struc (like in ex12). I need some information about the petscFem struct. I didn't find any document about that, is there one ? The field f0Funcs is used for the term \

Re: [petsc-users] how to select values from certain positions of a vec from different process to process 0.

2013-07-31 Thread 丁老师
Does the function VecScatterCreateToZero resolve this problem? At 2013-07-31 09:08:26,"Barry Smith" wrote: > > > create an IS with the global indices you want on process 0 and a seq Vec of > the same size, > > create an IS of size 0 on all the other processes and a seq Vec of the size

[petsc-users] Segmentatoin fault

2013-07-31 Thread 丁老师
I write a FVM code for the convection diffusion problem on a 2D cartesian grid. the code can run correctly with a grid of 720*720. and it use only 26% of the system memory(8G) during its running. but if i increase the grid to 800*800, it gives me the following error? could you please give me some

Re: [petsc-users] question about special matrix-matrix-matrix product

2013-07-31 Thread Kubis, Tillmann C
Hi, Thanks a lot. The serial version is working fine now. It will take a while until I can try the parallel one. Thanks, Tillmann Kubis Tillmann Kubis, PhD Research Assistant Professor Network for Computational Nanotechnology 207 S Martin Jischke Drive Purdue

[petsc-users] preallocation of mpiaij matrix

2013-07-31 Thread Hoang Giang Bui
Hi I want to consult the optimized way to allocate the mpiaij matrix for FEM application. Typically, elements in each process contribute to global matrix partially, therefore in each process I can define the number of nonzeros in according row locally but not globally. Is there method in pets

Re: [petsc-users] MATMPIBAIJ matrix allocation

2013-07-31 Thread Jose E. Roman
El 31/07/2013, a las 12:03, Heikki Virtanen escribió: > Hi, Sorry to bother you, but I still have problems with MATMPIBAIJ matrices. > It is said here > > http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatCreateBAIJ.html > > that these matrices are created using block compress

Re: [petsc-users] MATMPIBAIJ matrix allocation

2013-07-31 Thread Heikki Virtanen
> El 29/07/2013, a las 14:54, Heikki Virtanen escribió: > > > Hi, I try to solve an eigenvalue problem (Ax = lambda B x) where A > > and B are complex matrices. Unfortunately, my data structures > > are not capable of handling directly complex numbers, yet. So, > > I have to use > > > >