On Mon, Jul 13, 2015 at 7:12 AM, Alejandro D Otero aot...@fi.uba.ar wrote:
Hi, I seems to be having problems with PETSc vector VTK viewer. The output
file I get has some values which are different from what I actually have.
I have a scalar field represented by a vector. If I print it on
Thank you, Lisandro.
On Monday July 13 2015 at 10:09 Lisandro Dalcin wrote:
On 12 July 2015 at 13:52, Francesco Caimmi francesco.cai...@polimi.it
wrote:
These functions are not currently exposed in petsc4py. They should be
easy to add, but right now I don't have the time right now. Please
Do we have any data with respect to the number of PETSc users? Better yet
number of GAMG users?
Mark
You can save a (or multiple) linear system(s) with additional options
-ksp_view_mat binary:filename -ksp_view_rhs binary:filename
and then simply email the file named filename to petsc-ma...@mcs.anl.gov with
number of processes you are using and solver options and we can run the same
Greg, I am forwarding this to the PETSc mailing list.
Please send the entire output from this run. As I recall you were getting
a message that all values were not the same on all processors in GMRES. I
have seen this when I get NaNs in the system.
While you are doing this you should use a
On Jul 13, 2015, at 2:05 PM, Anthony Haas a...@email.arizona.edu wrote:
Hi,
I ran my program under Valgrind with 2 processors using the following:
${PETSC_DIR}/${PETSC_ARCH}/bin/mpiexec -n 2 valgrind --tool=memcheck
--leak-check=full --show-leak-kinds=all --num-callers=20
Run under valgrind first, see if it gives any more details about the memory
issue http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind
Can you send the code that produces this problem?
On Jul 13, 2015, at 10:56 AM, Florian Lindner mailingli...@xgm.de wrote:
Hello,
our petsc
You gnu compiler is ancient and supports little of C++11. Likely there is a
more recent one on the system (like 4.8) that you can switch to by using the
right module.
Barry
On Jul 13, 2015, at 10:35 AM, venkatesh g venkateshg...@gmail.com wrote:
sorry, I have attached the
Hi all, is it posible to build a rectangular matrix that relates two fields
with (posible) different number of DoF per point defining the nonzero
structure from 2 petscsections representing one field each?
As an example, I want to build an operator that given a 2D velocity
computes the
Hello,
our petsc application suffers from a memory error (double free or corruption).
Situation is a like that:
A KSP is private member of a C++ class. In its constructor I call KSPCreate.
Inbetween it may haben that I call KSPREset. In the class' destructor I call
KSPDestroy. That's where
Thanks for the suggestions. Suppose though that I wanted to avoid using a
self-consistent approach and just attack the full nonlinear problem. Is there
a “smart” way to handle my data where the Vec q, which stores both the function
on the mesh, and the eigenvalue parameter. Suppose the
Ok I got it thanks. The vector sol_seq that I used in
VecScatterCreateToZero (see below) was also created before with a
VecCreate() but only destroyed once. I didn't realize that the vector
was automatically created in the context of VecScatterCreateToZero. I
modified the code and the bugs
It lists the line numbers in __baseflow_MOD_poiseuille_3d_petsc
(module_baseflows.F90:182) __baseflow_MOD_poiseuille_3d_petsc
(module_baseflows.F90:356) where apparently two vectors are created but never
get destroyed. Perhaps they get created twice and you only destroy them once?
Or you
Please respond to all, to get the petsc mailing list.
The out file here looks fine. It would help of you run in debug mode (and
get -g) so that gdb can give line numbers.
If the error is in the dot product it looks like you are getting an Inf or
Nan someplace. If the first evaluation of the
Hi Barry,
If you look in module_baseflows.F90 line 467-473, I have:
call KSPDestroy(ksp,ierr)
call VecDestroy(frhs,ierr)
call VecDestroy(frhs2,ierr)
call VecDestroy(sol,ierr)
call VecDestroy(sol_seq,ierr)
call MatDestroy(DXX,ierr)
call MatDestroy(DYY,ierr)
I checked
You can use a DMCOMPOSITE with a DMDA and a DMREDUNDANT see
src/snes/examples/tutorials/ex21.c It is intended for this type of situation.
Barry
On Jul 13, 2015, at 4:05 PM, Gideon Simpson gideon.simp...@gmail.com wrote:
Thanks for the suggestions. Suppose though that I wanted to
Hi, when I use --with-cxx-dialect=C++11, I get the error
Could not determine compiler flag for with-cxx-dialect=C++11, use CXXFLAGS
Kindly let me know the syntax or how to enable this.
w.r.t --with-mpi-dir on Crays, I will keep it in mind and will not pass
them as you have told Barry, they
17 matches
Mail list logo