Re: [petsc-users] Problem with MPI, MatAXPY and SAME_NONZERO_PATTERN

2015-03-30 Thread Klaus Kaiser
Hallo Metthew, at the weekend I updated petsc to the newest release version 3.5.2 and I still got this strange behavior. Best Klaus On 03/27/2015 05:07 PM, Klaus Kaiser wrote: Hallo Matthew, here is a short example I used for 8 cores: int d_nnz[3]; int o_nnz[3]; d_nnz[0] = 3;

Re: [petsc-users] Problem with MPI, MatAXPY and SAME_NONZERO_PATTERN

2015-03-30 Thread Klaus Kaiser
Hi Barry, I updated to the newest version but the behavior stays the same. Best Klaus On 03/27/2015 06:32 PM, Barry Smith wrote: On Mar 27, 2015, at 11:14 AM, Klaus Kaiser kai...@igpm.rwth-aachen.de wrote: Hallo Barry, I think you're right, but this is the setup I used for more than a

Re: [petsc-users] Problem with MPI, MatAXPY and SAME_NONZERO_PATTERN

2015-03-30 Thread Matthew Knepley
On Mon, Mar 30, 2015 at 5:18 AM, Klaus Kaiser kai...@igpm.rwth-aachen.de wrote: Hallo Metthew, at the weekend I updated petsc to the newest release version 3.5.2 and I still got this strange behavior. Okay, I will check it out. The first thing I am going to do is try with AIJ matrices

Re: [petsc-users] Understanding MatCreate bsize parameter

2015-03-30 Thread Matthew Knepley
On Mon, Mar 30, 2015 at 5:59 AM, Florian Lindner mailingli...@xgm.de wrote: Am Freitag, 27. März 2015, 07:34:56 schrieb Matthew Knepley: On Fri, Mar 27, 2015 at 7:31 AM, Florian Lindner mailingli...@xgm.de wrote: Am Freitag, 27. März 2015, 07:26:11 schrieb Matthew Knepley: On

Re: [petsc-users] Understanding MatCreate bsize parameter

2015-03-30 Thread Florian Lindner
Am Freitag, 27. März 2015, 07:34:56 schrieb Matthew Knepley: On Fri, Mar 27, 2015 at 7:31 AM, Florian Lindner mailingli...@xgm.de wrote: Am Freitag, 27. März 2015, 07:26:11 schrieb Matthew Knepley: On Fri, Mar 27, 2015 at 4:28 AM, Florian Lindner mailingli...@xgm.de wrote: Am

[petsc-users] Documentation for VTU format?

2015-03-30 Thread Justin Chang
Hi all, I have a couple questions: 1) Is there any documentation somewhere that talks specifically on how to output using the *.vtu format? Specifically for transient problems. When I look at TS ex11 it seems to me the user has to create several additional routines, but I was wondering if there

Re: [petsc-users] Documentation for VTU format?

2015-03-30 Thread Justin Chang
Jed thank you for the response. 1) In my case, I want to output a transient diffusion problem (with DMPlex). I see that in ex11 there's a TSMonitorSet(). Is this (and its calling function) be the place I have something along the lines that you have suggested? That is, if TSMonitorSet() calls

Re: [petsc-users] Documentation for VTU format?

2015-03-30 Thread Jed Brown
Justin Chang jchan...@uh.edu writes: Jed thank you for the response. 1) In my case, I want to output a transient diffusion problem (with DMPlex). I see that in ex11 there's a TSMonitorSet(). Is this (and its calling function) be the place I have something along the lines that you have

Re: [petsc-users] HermitianTranspose version of MatCreateTranspose.

2015-03-30 Thread Jed Brown
Andrew Spott ansp6...@colorado.edu writes: I submitted a pull request for this: https://bitbucket.org/petsc/petsc/pull-request/291/creates-a-matcreatehermitiantranspose/diff However, I forgot to pull master before submitting it, and someone added ex180 in the meantime.  Should I redo

Re: [petsc-users] HermitianTranspose version of MatCreateTranspose.

2015-03-30 Thread Andrew Spott
I submitted a pull request for this: https://bitbucket.org/petsc/petsc/pull-request/291/creates-a-matcreatehermitiantranspose/diff However, I forgot to pull master before submitting it, and someone added ex180 in the meantime.  Should I redo the pull request, or just let the merger fix it?

Re: [petsc-users] TimeStepper norm problems. EMIL Please read this

2015-03-30 Thread Andrew Spott
Emil:  If I want to fix this myself, do you have any idea where I should start looking?  I assume in the th-endpoint section of TSStep_Theta in http://www.mcs.anl.gov/petsc/petsc-current/src/ts/impls/implicit/theta/theta.c.html? Thanks, -Andrew On Tue, Mar 24, 2015 at 10:52 AM, Lisandro

[petsc-users] FEM 101

2015-03-30 Thread Justin Pogacnik
Hello, I'm beginning to write a finite element solid mechanics code in F90 and am having some problems getting started with finite element basics within Petsc. I'm using an unstructured grid and the DMPlex class and able to access most of the mesh information required for FE numerical

Re: [petsc-users] TimeStepper norm problems. EMIL Please read this

2015-03-30 Thread Andrew Spott
I’m sorry, I didn’t want to imply you were taking too much time.  I was just interested in trying to solve the problem myself, especially if it was a simpler fix. I don’t have enough of an understanding of the internals of PETSc though, if it depends on the compiler flags. Good luck!

Re: [petsc-users] TimeStepper norm problems. EMIL Please read this

2015-03-30 Thread Emil Constantinescu
Hi Andrew, Sorry for the delay; things are going slower with me over the spring break but we are working on it. The problem that you see may not be related to the integrator, which runs well (we get expected convergence rates and good results) for other problems, in particular, for problems

Re: [petsc-users] TimeStepper norm problems. EMIL Please read this

2015-03-30 Thread Emil Constantinescu
On 3/30/15 9:32 PM, Andrew Spott wrote: I’m sorry, I didn’t want to imply you were taking too much time. I was just interested in trying to solve the problem myself, especially if it was a simpler fix. Nothing of that sort. I just expected to figure it out sooner, especially when you were

Re: [petsc-users] FEM 101

2015-03-30 Thread Matthew Knepley
On Mon, Mar 30, 2015 at 4:39 PM, Justin Pogacnik j.pogac...@auckland.ac.nz wrote: Hello, I'm beginning to write a finite element solid mechanics code in F90 and am having some problems getting started with finite element basics within Petsc. I'm using an unstructured grid and the DMPlex

Re: [petsc-users] Documentation for VTU format?

2015-03-30 Thread Jed Brown
Justin Chang jchan...@uh.edu writes: Hi all, I have a couple questions: 1) Is there any documentation somewhere that talks specifically on how to output using the *.vtu format? Specifically for transient problems. When I look at TS ex11 it seems to me the user has to create several