Re: [petsc-users] Modified Taylor-Hood elements with piece-wise constant pressure for Stokes equation

2015-05-31 Thread Justin Chang
There are a few papers that discuss this modified/augmented Taylor-Hood elements for Stokes equations in detail (e.g., http://link.springer.com/article/10.1007%2Fs10915-011-9549-4). From what I have seem, it seems people primarily use this to ensure local mass conservation while attaining the desir

Re: [petsc-users] Modified Taylor-Hood elements with piece-wise constant pressure for Stokes equation

2015-05-31 Thread Jed Brown
Justin Chang writes: > I am referring to P2 / (P1 + P0) elements, I think this is the correct way > of expressing it. Some call it modified Taylor Hood, others call it > something else, but it's not Crouzeix-Raviart elements. Okay, thanks. This pressure space is not a disjoint union (the consta

Re: [petsc-users] Offloading linear solves in time stepper to GPU

2015-05-31 Thread Barry Smith
> On May 30, 2015, at 10:14 PM, Harshad Sahasrabudhe > wrote: > > Is your intent to solve a problem that matters in a way that makes sense for > a scientist or engineer > > I want to see if we can speed up the time stepper for a large system using > GPUs. For large systems with sparse matrix

Re: [petsc-users] MUMPS error and superLU error

2015-05-31 Thread Hong
venkatesh, As we discussed previously, even on smaller problems, both mumps and superlu_dist failed, although Mumps gave "OOM" error in numerical factorization. You acknowledged that B is singular, which may need additional reformulation for your eigenvalue problems. The option '-st_type sinvert'

Re: [petsc-users] NASM FormFunctionLocal

2015-05-31 Thread Rongliang Chen
Hi Barry, Thanks for your reply. I have no idea how to pass a vector into the FormFunctionLocal() right now. I will try to read the NASM source code and hope can find a way to do it. Best regards, Rongliang On 05/31/2015 10:41 AM, Barry Smith wrote: Rongliang, It appears the NASM

[petsc-users] MUMPS error and superLU error

2015-05-31 Thread Dave May
It failed due to a lack of memory. "OOM" stands for "out of memory". OOM killer terminated your job means you ran out of memory. On Sunday, 31 May 2015, venkatesh g > wrote: > Hi all, > > I tried to run my Generalized Eigenproblem in 120 x 24 = 2880 cores. > The matrix size of A = 20GB and B =

Re: [petsc-users] MUMPS error and superLU error

2015-05-31 Thread venkatesh g
Sorry the command was: aprun -n 2880 -N 24 ./ex7 -f1 a110t -f2 b110t -st_type sinvert -eps_nev 1 -st_ksp_type preonly -st_pc_type lu -st_pc_factor_mat_solver_package mumps -mat_mumps_cntl_1 1e-6 -mat_mumps_icntl_4 2 cheers, Venkatesh On Sun, May 31, 2015 at 6:07 PM, venkatesh g wrote: > Hi all,

Re: [petsc-users] MUMPS error and superLU error

2015-05-31 Thread venkatesh g
Hi all, I tried to run my Generalized Eigenproblem in 120 x 24 = 2880 cores. The matrix size of A = 20GB and B = 5GB. It got killed after 7 Hrs of run time. Please see the mumps error log. Why must it fail ? I gave the command: aprun -n 240 -N 24 ./ex7 -f1 a110t -f2 b110t -st_type sinvert -eps_n