Re: [petsc-dev] KSPComputeExtremeSingularValues_GMRES floating point exception

wrote: > > On Jul 23, 2014, at 9:54 AM, John Mousel wrote: > > > > > > > > > > > I'm getting an FPE in LAPACKgesvd of > KSPComputeExtremeSingularValues_GMRES. > > ierr = PetscFPTrapPush(PETSC_FP_TRAP_OFF);CHKERRQ(ierr); > #if !defined(PETSC_USE_

[petsc-dev] KSPComputeExtremeSingularValues_GMRES floating point exception

I'm getting an FPE in LAPACKgesvd of KSPComputeExtremeSingularValues_GMRES. I'm not sure where this routine is required during the algorithm as I'm using BiCG + richardson/SOR on all multigrid levels. Also, even though I've set mg_coarse_ksp_type richardson, KSPView is showing that preonly is used

[petsc-dev] Floating Point Error in VecSet

I'm getting a floating point error in VecSet. I see a similar error reported in a nightly log here: http://ftp.mcs.anl.gov/pub/petsc/nightlylogs/archive/2014/06/18/examples_master_arch-opensolaris-cmplx-pkgs-dbg_n-gage.log I've attached the error message below. Any suggestions? Thanks, John [

[petsc-dev] ParMETIS assertion with pc_gamg_repartition

I'm getting a ParMETIS assertion: [53] ***ASSERTION failed on line 176 of file /nics/b/home/jmousel/NumericalLibraries/petsc/intel-debug/externalpackages/parmetis-4.0.3-p1/libparmetis/comm.c: j == nnbrs [110] ***ASSERTION failed on line 176 of file /nics/b/home/jmousel/NumericalLibraries/petsc/int

Re: [petsc-dev] Installation crash on test

Satish, The script that pulled petsc-dev uses: export PETSC_DIR=/home/augspurger/NumericalLibraries/petsc git clone https://bitbucket.org/petsc/petsc cd petsc export PETSC_ARCH=gnu-opt ./config/configure.py --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-openmpi --download-metis --down

Re: [petsc-dev] Error in PetscMallocValidate

-pc_gamg_agg_nsmooths 1 -pc_gamg_sym_graph true -ksp_monitor_true_residual John On Tue, Sep 24, 2013 at 9:58 AM, Jed Brown wrote: > John Mousel writes: > > > I cloned a fresh repo and built from scratch. I'm seeing the same error I > > previously reported with both

Re: [petsc-dev] Error in PetscMallocValidate

ction John On Mon, Sep 23, 2013 at 6:50 PM, John Mousel wrote: > Make test ran correctly for both debug and optimized builds. It runs fine > with plain bJacobi/ilu preconditioning. > > > On Monday, September 23, 2013, Matthew Knepley wrote: > >> On Mon, Sep 23, 2013

Re: [petsc-dev] Error in PetscMallocValidate

Make test ran correctly for both debug and optimized builds. It runs fine with plain bJacobi/ilu preconditioning. On Monday, September 23, 2013, Matthew Knepley wrote: > On Mon, Sep 23, 2013 at 1:14 PM, John Mousel > > > wrote: > >> I just pulled from petsc-dev an

Re: [petsc-dev] Valgrind invalid read in GAMG

Matt, Here is the output of ksp_view_binary. John binaryoutput<https://docs.google.com/file/d/0B3DmQgVAlCBUVENZSVhtSkE4eTA/edit?usp=drive_web> On Fri, Aug 9, 2013 at 9:14 AM, Matthew Knepley wrote: > On Fri, Aug 9, 2013 at 8:49 AM, John Mousel wrote: > >> I'

[petsc-dev] Valgrind invalid read in GAMG

I'm getting the following invalid read when I use GAMG. I've included the output of KSPView at the bottom. ==18312== Memcheck, a memory error detector ==18312== Copyright (C) 2002-2012, and GNU GPL'd, by Julian Seward et al. ==18312== Using Valgrind-3.8.1 and LibVEX; rerun with -h for copyright in

[petsc-dev] GAMG MatMatMult error on 1 core

, 2013 at 10:31 AM, John Mousel gmail.com>wrote: > >> So this >> >> Ap => proj%A >> CALL MatCreate(solver%comm,Ap,ierr) >> CALL MatSetSizes(Ap,nlp,nlp,ngp,ngp,ierr) >> CALL MatSetType(Ap,MATMPIAIJ,ierr) >> CALL MatMPIAIJSetPr

[petsc-dev] GAMG MatMatMult error on 1 core

CALL MatSetSizes(Ap,nlp,nlp,ngp,ngp,ierr) CALL MatSetType(Ap,MATAIJ,ierr) CALL MatXAIJSetPreallocation(Ap,0,dnnz,0,onnz,ierr) ? On Fri, Mar 29, 2013 at 10:25 AM, Jed Brown wrote: > On Fri, Mar 29, 2013 at 10:22 AM, John Mousel wrote: > >> I get the following error when I try to use

[petsc-dev] GAMG MatMatMult error on 1 core

I get the following error when I try to use GAMG on 1 core. I just pulled dev and it still exists. John [0]PETSC ERROR: - Error Message [0]PETSC ERROR: Arguments are incompatible! [0]PETSC ERROR: MatMatMult requires A, mpiaij, to be compati

[petsc-dev] Vector operations on ghosted vectors

standing is occuring. I do appreciate your help. John On Sat, Mar 2, 2013 at 7:06 AM, Matthew Knepley wrote: > On Fri, Mar 1, 2013 at 12:40 PM, John Mousel wrote: > >> So am I right in reading this as you suggesting to: >> >> 1. Compute VecAXPY on all the u(n),u(n-1)...

[petsc-dev] Vector operations on ghosted vectors

is way. > > Matt > > >> John >> >> J. L. Guermond, P. Minev, and J. Shen. An overview of projection methods >> for incompressible flows. Computer Methods in Applied Mechanics and >> Engineering, 195(44-47):6011?6045, 2006. >> >> >>

[petsc-dev] Vector operations on ghosted vectors

ible flows. Computer Methods in Applied Mechanics and Engineering, 195(44-47):6011?6045, 2006. On Fri, Mar 1, 2013 at 9:28 AM, Matthew Knepley wrote: > On Fri, Mar 1, 2013 at 10:24 AM, John Mousel wrote: > >> I use the data in the vectors for things besides pure vector operations. >

[petsc-dev] Vector operations on ghosted vectors

portion of the vector in the u(k)->u(k-1) operation, but the ghosts as well. Doesn't copying the entire ghosted vector make sense in this context? John On Fri, Mar 1, 2013 at 8:09 AM, Jed Brown wrote: > On Wed, Feb 27, 2013 at 4:56 PM, John Mousel wrote: > >> Is there a pos

[petsc-dev] Vector operations on ghosted vectors

Let me rephrase. Is a functionality like this something that is worth adding to the PETSc source, or should I just write and use my own wrapper. On Wed, Feb 27, 2013 at 4:51 PM, Matthew Knepley wrote: > On Wed, Feb 27, 2013 at 4:56 PM, John Mousel wrote: > >> Is there a possibility

[petsc-dev] Vector operations on ghosted vectors

Is there a possibility of adding a wrapper function around a few basic vector operations such as VecCopy, VecAXPY, VECAXPYPZ... to operate on ghosted vectors? I perform a lot of vector operations including the ghost region to avoid communication. It really clutters code to see stuff like CALL

[petsc-dev] corrupted ta->count with GAMG

des in the same aggregate, this time... *ML*WRN* (file ./Coarsen/ml_agg_METIS.c, line 966) I'm about to build, and I didn't know if I'm missing an option that I need. John On Thu, Oct 25, 2012 at 10:24 AM, Mark F. Adams wrote: > > On Oct 25, 2012, at 7:30 AM, John Mouse

[petsc-dev] corrupted ta->count with GAMG

t-petscsf-without-one-sided > > > On Thu, Oct 25, 2012 at 9:23 AM, Mark F. Adams columbia.edu>wrote: > >> The problem is detected in MatTranspose and Jed was messing with this a >> few weeks ago ... >> >> >> On Oct 25, 2012, at 6:48 AM, John Mousel wrot

[petsc-dev] corrupted ta->count with GAMG

ately that does me no good! John On Thu, Oct 25, 2012 at 9:23 AM, Mark F. Adams wrote: > The problem is detected in MatTranspose and Jed was messing with this a > few weeks ago ... > > > On Oct 25, 2012, at 6:48 AM, John Mousel wrote: > > Mark, > > Yes, this worked be

[petsc-dev] corrupted ta->count with GAMG

Mark, Yes, this worked before pulling. I hadn't pulled on that machine in a few weeks though. John On Thu, Oct 25, 2012 at 1:18 AM, Mark F. Adams wrote: > Did this work before the pull? > > Mark > > On Oct 24, 2012, at 8:04 PM, John Mousel wrote: > > > I just pul

[petsc-dev] corrupted ta->count with GAMG

I just pulled and I'm getting the error below when running with GAMG. [20]PETSC ERROR: - Error Message [20]PETSC ERROR: Corrupted Petsc object! [20]PETSC ERROR: corrupted ta->count! [20]PETSC ERROR: --

[petsc-dev] PetscMallocValidate Error

, 2012 at 11:52 AM, John Mousel wrote: > >> I just pulled and updated, and I'm getting an error in >> PetscMallocValidate. >> I've attached the error output. >> >> Thanks, >> John >> > > -- next part -- An HTML at

[petsc-dev] PetscMallocValidate Error

I just pulled and updated, and I'm getting an error in PetscMallocValidate. I've attached the error output. Thanks, John -- next part -- An HTML attachment was scrubbed... URL: ---

[petsc-dev] GAMG MatSetValues error

nneeded,2591 used [1] MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 15 [1] MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 13 John On Tue, Sep 25, 2012 at 5:41 AM, Matthew Knepley wrote: > On Mon, Sep 24, 2012 at 9:48 PM, John Mousel wrote: > >> I'm having a problem

[petsc-dev] GAMG MatSetValues error

I'm having a problem with GAMG causing a MatSetValues() is [some non-zero number] after [0] PCSetUp(): Setting up new PC I don't get this with either ILU or HYPRE/BoomerAMG. I've attached the output of -info. I see there is a recent change to some pre-allocation for GAMG here: http://petsc.cs.i

[petsc-dev] [petsc-maint #119133] petsc-dev configure crash

Thanks Jed. Somehow I forgot to put the post-pull in .hg/hgrc for that build. On Wed, Jun 6, 2012 at 10:03 AM, Jed Brown wrote: > You need to pull BuildSystem. > > > On Wed, Jun 6, 2012 at 10:00 AM, John Mousel wrote: > >> I'm getting the following after

[petsc-dev] make test error

I just pulled petsc-dev and am getting the following error in the tests. John Running test examples to verify correct installation Using PETSC_DIR=/home/vchivukula/NumericalLibraries/petsc-dev and PETSC_ARCH=linux-intel Possible error running C/C++ src/snes/examples/tutorials/ex19 with 1 MPI pro

[petsc-dev] ParMetis4.0

ParMetis4.0 has just been released as a stable version. It supposedly has a lot more error checking. Is there any way this could be added to the --download option? Thanks, John -- next part -- An HTML attachment was scrubbed... URL:

[petsc-dev] PCREDISTRIBUTE

Could someone throw a #define PCREDISTRIBUTE 'redistribute' in /include/finclude/petscpcdef.h? Thanks, John -- next part -- An HTML attachment was scrubbed... URL: