Re: [petsc-users] Loading Laplacian as Module

I don't know much about modules so can't help, but PETSc variables are just like any other variables and so should behave in the same way. Barry > On Sep 23, 2016, at 2:31 PM, Manuel Valera wrote: > > Ok, i got this: > > RHS loaded, size: 213120 /

Re: [petsc-users] Loading Laplacian as Module

Ok, i got this: RHS loaded, size: 213120 / 213120 CSRMAt loaded, sizes: 213120 x 213120 8.43036175 s solveP pass: 1 RHS loaded, size: 213120 / 213120 [0]PETSC ERROR: MatGetSize() line 6295 in

Re: [petsc-users] Loading Laplacian as Module

Ok, so the problem is not memory corruption. > 0]PETSC ERROR: - Error Message > -- > [0]PETSC ERROR: Invalid argument > [0]PETSC ERROR: Wrong type of object: Parameter # 1 > [0]PETSC ERROR: #1 MatGetSize() line

Re: [petsc-users] Loading Laplacian as Module

Barry, that was awesome, all the valgrind error dissappeared after using the mpiexec from petsc folder, the more you know... Anyway this is my output from valgrind running with those options: Last Update: 9/23/201612: 5:12 ucmsSeamount Entering MAIN loop. RHS loaded, size: 213120 /

Re: [petsc-users] Loading Laplacian as Module

> On Sep 23, 2016, at 1:09 PM, Manuel Valera wrote: > > Thanks Barry, for the quick reply, > > I tried doing that once recently, not for this problem though, but it looks > like the model i'm working on isn't optimized at all for memory leaks, and > valgrind stopped

Re: [petsc-users] Loading Laplacian as Module

Thanks Barry, for the quick reply, I tried doing that once recently, not for this problem though, but it looks like the model i'm working on isn't optimized at all for memory leaks, and valgrind stopped with thousands of errors before reaching this part of the execution. Is there maybe an

Re: [petsc-users] Loading Laplacian as Module

Run with valgrind to find the exact location of the first memory corruption. http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind > On Sep 23, 2016, at 12:47 PM, Manuel Valera wrote: > > Hello all, > > I'm trying to load my laplacian matrix into a fortran

[petsc-users] Loading Laplacian as Module

Hello all, I'm trying to load my laplacian matrix into a fortran module, and i have implemented it and it works for the first iteration of laplacian solver, but when starts the second step the laplacian matrix object becomes corrupts and looks like it loses one of it's dimensions. Can you help

Re: [petsc-users] How to solve the pressure possion equations with four neuman bc

Any references on this topic except the users manual? I am interested in mathematics theory on this topic. Fande Kong, On Fri, Sep 23, 2016 at 9:58 AM, Jed Brown wrote: > 丁老师 writes: > > > Dear friends: > > In the projection method for the solution

Re: [petsc-users] How to solve the pressure possion equations with four neuman bc

"Kong, Fande" writes: > Any references on this topic except the users manual? I am interested in > mathematics theory on this topic. These are relevant. I've only skimmed them briefly, but they might suggest ways to improve PETSc's handling of singular systems. The

Re: [petsc-users] How to solve the pressure possion equations with four neuman bc

丁老师 writes: > Dear friends: > In the projection method for the solution of incompressible flow, a > pressure equation with four neuman bcs need to be solved. > but the pressure matrix is singular. it gives divergence solution. how to > solve this kind of equations.

Re: [petsc-users] Write binary to matrix

Florian: I pushed a fix in branch hzhang/fix_matview_mpisbaij (off petsc-maint) https://bitbucket.org/petsc/petsc/commits/d1654148bc9f02cde4d336bb9518a18cfb35148e After it is tested in our regression tests, it will be merged to petsc-maint and petsc-master. Thanks for reporting it! Hong On Fri,

Re: [petsc-users] Write binary to matrix

Florian: I can reproduce this error. This is a bug in PETSc library. I'll fix it and get back to you soon. Hong > > Am 22.09.2016 um 18:34 schrieb Hong: > > Florian: > > Would it work if replacing MATSBAIJ to MATAIJ or MATMPISBAIJ? > > MATAIJ works, but is not an option for my actual application.

Re: [petsc-users] Error during compiling my own code

On Fri, Sep 23, 2016 at 10:01 AM, Melanie Li Sing How wrote: > make[1]: Entering directory `/glade/u/home/lisingh/hippstr2.0/src' > make[2]: Entering directory `/glade/u/home/lisingh/hippstr2.0/src/library' > make[2]: `/glade/u/home/lisingh/hippstr2.0/lib/liblibrary.a' is up

Re: [petsc-users] Error during compiling my own code

On Fri, Sep 23, 2016 at 10:08 AM, Melanie Li Sing How wrote: > Ok Thank you. I only have include “mpif.h” in both modules that are > conflicting but thank you for your help. > Obviously from the message, you have a USE statement that eventually uses mpi. Matt > On Sep

Re: [petsc-users] Error during compiling my own code

make[1]: Entering directory `/glade/u/home/lisingh/hippstr2.0/src' make[2]: Entering directory `/glade/u/home/lisingh/hippstr2.0/src/library' make[2]: `/glade/u/home/lisingh/hippstr2.0/lib/liblibrary.a' is up to date. make[2]: Leaving directory `/glade/u/home/lisingh/hippstr2.0/src/library'

Re: [petsc-users] How to solve the pressure possion equations with four neuman bc

On Fri, Sep 23, 2016 at 9:38 AM, 丁老师 wrote: > Dear friends: > In the projection method for the solution of incompressible flow, a > pressure equation with four neuman bcs need to be solved. > but the pressure matrix is singular. it gives divergence solution. how to >

Re: [petsc-users] Error during compiling my own code

On Fri, Sep 23, 2016 at 9:53 AM, Melanie Li Sing How wrote: > Hi, > I happened to fall on your link (http://lists.mcs.anl.gov/ > pipermail/petsc-users/2010-April/006170.html) to the same problem I am > having. I am a yellowstone user and I could not find a contact detail for

Re: [petsc-users] Error during compiling my own code

Hi, I happened to fall on your link (http://lists.mcs.anl.gov/pipermail/petsc-users/2010-April/006170.html ) to the same problem I am having. I am a yellowstone user and I could not find a contact detail for the same error

[petsc-users] How to solve the pressure possion equations with four neuman bc

Dear friends: In the projection method for the solution of incompressible flow, a pressure equation with four neuman bcs need to be solved. but the pressure matrix is singular. it gives divergence solution. how to solve this kind of equations. Regards

Re: [petsc-users] DMPlex problem

On Fri, Sep 23, 2016 at 7:45 AM, Matthew Knepley wrote: > On Fri, Sep 23, 2016 at 3:48 AM, Morten Nobel-Jørgensen > wrote: > >> Dear PETSc developers >> >> Any update on this issue regarding DMPlex? Or is there any obvious >> workaround that we are unaware of? >>

Re: [petsc-users] DMPlex problem

On Fri, Sep 23, 2016 at 3:48 AM, Morten Nobel-Jørgensen wrote: > Dear PETSc developers > > Any update on this issue regarding DMPlex? Or is there any obvious > workaround that we are unaware of? > I have fixed this bug. It did not come up in nightly tests because we are not using

Re: [petsc-users] Write binary to matrix

Am 22.09.2016 um 18:34 schrieb Hong: > Florian: > Would it work if replacing MATSBAIJ to MATAIJ or MATMPISBAIJ? MATAIJ works, but is not an option for my actual application. MATMPISBAIJ does not work. Not very suprisingly, since afaik setting it to MATSBAIJ and executing it on multiple MPI

Re: [petsc-users] DMPlex problem

Dear PETSc developers Any update on this issue regarding DMPlex? Or is there any obvious workaround that we are unaware of? Also should we additionally register the issue on Bitbucket or is reporting the issue on the mailing list enough? Kind regards, Morten