[petsc-users] hdf5 viewer and zero sized local vectors

--- > [2]PETSC ERROR: VecView_MPI_HDF5() line 748 in > src/vec/vec/impls/mpi/pdvec.c > [2]PETSC ERROR: VecView_MPI() line 840 in src/vec/vec/impls/mpi/pdvec.c > [2]PETSC ERROR: VecView() line 710 in src/vec/vec/interface/vector.c > [2]PETSC ERROR: Output_grid_hdf5() line 2358 in "unknowndirectory/"Output.c > > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/5f152f55/attachment.htm>

[petsc-users] The ColPerm option for superlu

> > Thanks very much! > > Best regards, > > Rebecca > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/4e878a11/attachment.htm>

[petsc-users] hdf5 viewer and zero sized local vectors

n_local = N; } ierr = VecSetSizes(x,n_local, N); PETScErrAct(ierr); ierr = VecSetFromOptions(x); PETScErrAct(ierr); double *vcoord; ierr = VecGetArray(x, &vcoord); PETScErrAct(ierr); int i; for (i=0; ihttp://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/6263d8be/attachment-0001.htm>

[petsc-users] A few questions about hdf5 viewer

or granted before they begin their >> experiments is infinitely more interesting than any results to which their >> experiments lead. >> -- Norbert Wiener >> > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/ee1a6f18/attachment-0001.htm>

[petsc-users] A few questions about hdf5 viewer

;> >> >> If you want this special case, you have to manage it by hand. >> > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener ---

[petsc-users] A few questions about hdf5 viewer

mail/petsc-users/attachments/20120119/2c022bce/attachment.htm>

[petsc-users] MKL BLAS interface inconsistency

This (dirty) patch solves the problem: diff -r a3e9ca59ab58 src/vec/vec/impls/seq/bvec2.c --- a/src/vec/vec/impls/seq/bvec2.c Tue Jan 17 22:04:05 2012 -0600 +++ b/src/vec/vec/impls/seq/bvec2.c Thu Jan 19 17:28:39 2012 +0100 @@ -232,12 +232,13 @@ PetscErrorCodeierr; PetscInt

[petsc-users] A few questions about hdf5 viewer

art -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/04956938/attachment.htm>

[petsc-users] A few questions about hdf5 viewer

5 (Binary) the size is not reduced in the (*.h5) case. >>>>> In fact, it is slightly bigger than pure binary file!! >>>>> Is there any command we have to set in Petsc to tell HDF5 viewer to >>>>> use data compression? >>>>> >>>> >>>> We do not support it. We are happy to take patches that enable this. >>>> >>>> Thanks, >>>> >>>> Matt >>>> >>>> >>>>> Thanks for your patience, >>>>> Best, >>>>> Mohamad >>>>> >>>>> >>>>> >>>>> >>>>> >>>> >>>> >>>> -- >>>> What most experimenters take for granted before they begin their >>>> experiments is infinitely more interesting than any results to which their >>>> experiments lead. >>>> -- Norbert Wiener >>>> >>> >>> >> >> >> -- >> What most experimenters take for granted before they begin their >> experiments is infinitely more interesting than any results to which their >> experiments lead. >> -- Norbert Wiener >> > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/b8949b6e/attachment.htm>

[petsc-users] The ColPerm option for superlu

Hello, I use ksp_view to get the information about the lu pc as follow: PC Object: 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 matrix ordering: nd factor fill ratio given 0, needed 0 Factored matrix follows: Matrix O

[petsc-users] A few questions about hdf5 viewer

;>> >> >> We do not support it. We are happy to take patches that enable this. >> >> Thanks, >> >> Matt >> >> >>> Thanks for your patience, >>> Best, >>> Mohamad >>> >>> >>> >>> >>> >> >> >> -- >> What most experimenters take for granted before they begin their >> experiments is infinitely more interesting than any results to which their >> experiments lead. >> -- Norbert Wiener >> > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/06f4388f/attachment.htm>

[petsc-users] A few questions about hdf5 viewer

lization software to load the data using those >>>>>> coordinates. I am using orthogonal grid, that's why I don't need all the >>>>>> (x,y,z) coordinates for each cell. >>>>> >>>>> >>>>> If you want this specia

[petsc-users] A few questions about hdf5 viewer

gt;>> If you want this special case, you have to manage it by hand. >>> >> >> > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/d1dd7e10/attachment.htm>

[petsc-users] A few questions about hdf5 viewer

ll. > > > If you want this special case, you have to manage it by hand. > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/a28c714d/attachment.htm>

[petsc-users] A few questions about hdf5 viewer

using VecView() and HDF5 writer, what is the status of data >>>>>> compression? >>>>>> The reason that I am asking is that, I used the same example above >>>>>> and comparing two files saved via two different PetscViewers, i.e. (just) >>>>>> Binary and HDF5 (Binary) the size is not reduced in the (*.h5) case. >>>>>> In fact, it is slightly bigger than pure binary file!! >>>>>> Is there any command we have to set in Petsc to tell HDF5 viewer to >>>>>> use data compression? >>>>>> >>>>> >>>>> We do not support it. We are happy to take patches that enable this. >>>>> >>>>> Thanks, >>>>> >>>>> Matt >>>>> >>>>> >>>>>> Thanks for your patience, >>>>>> Best, >>>>>> Mohamad >>>>>> >>>>>> >>>>>> >>>>>> >>>>>> >>>>> >>>>> >>>>> -- >>>>> What most experimenters take for granted before they begin their >>>>> experiments is infinitely more interesting than any results to which their >>>>> experiments lead. >>>>> -- Norbert Wiener >>>>> >>>> >>>> >>> >>> >>> -- >>> What most experimenters take for granted before they begin their >>> experiments is infinitely more interesting than any results to which their >>> experiments lead. >>> -- Norbert Wiener >>> >> >> > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/c0b2a4da/attachment-0001.htm>

[petsc-users] MKL BLAS interface inconsistency

Dear petsc team, I've been struggling with this problem for a long time and seeking for your advice now. Let's take simple petsc program: #include int main(int argc,char **args) { Vecb; PetscReal norm; PetscInitialize(&argc,&args,(char *)0,NULL); VecCreateSeq(

[petsc-users] A few questions about hdf5 viewer

gt;> In fact, it is slightly bigger than pure binary file!! >>>> Is there any command we have to set in Petsc to tell HDF5 viewer to use >>>> data compression? >>>> >>> >>> We do not support it. We are happy to take patches that enable this. >>> >>> Thanks, >>> >>> Matt >>> >>> >>>> Thanks for your patience, >>>> Best, >>>> Mohamad >>>> >>>> >>>> >>>> >>>> >>> >>> >>> -- >>> What most experimenters take for granted before they begin their >>> experiments is infinitely more interesting than any results to which their >>> experiments lead. >>> -- Norbert Wiener >>> >> >> > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/131beeb3/attachment-0001.htm>

[petsc-users] VecView() and binary output

the size of the data. > Could someone tell me where that extra 8 byte is stored and how it is > going help? > > best > Mohamad > * > * > > > > -- > Regards, > Alexander > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/afc81b8e/attachment.htm>

[petsc-users] A few questions about hdf5 viewer

t;> Thanks for your patience, >> Best, >> Mohamad >> >> >> >> >> > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/5543660d/attachment.htm>

[petsc-users] [KSP] true norm residual

Hi, When I solve a linear system and monitor the true norm residual I got somethig like this (rtol is set to 1e-6) : Residual norms for stokes_ solve. 0 KSP preconditioned resid norm 1.270103385648e+00 true resid norm 2.319288792500e+06 ||r(i)||/||b|| 1.667494921925e-01 1 KSP preconditioned

[petsc-users] unstructured finite volume method matrix assembly for partitioned mesh

any results to which their > experiments lead. > -- Norbert Wiener > -- Zeus Numerix Pvt. Ltd. I2IT Campus, P-14 Rajiv Gandhi Infotech Park, Phase-1 Hinjewadi Pune -57 0ff. +91 64731511/9965 -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/750da5ae/attachment.htm>

[petsc-users] VecView() and binary output

gt; > best > > > Mohamad > > * > * -- Regards, Alexander -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/8907db0d/attachment.htm>

[petsc-users] A few questions about hdf5 viewer

What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/186d4d8f/attachment-0001.htm>

[petsc-users] access to matnest block (0,1) ?

hey share > the same partitioning but could have different stencils. How to use > the "split local space", DMComposite and MATNEST in this case? > If you order this way, then you don't need DMComposite or MatNest (although you can still make a MatNest that operates in this ordering, we just don't have a way to make it automatically). -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/9d2ff3aa/attachment.htm>

[petsc-users] access to matnest block (0,1) ?

> > I have two DMs which I add to a DMComposite and then use MATNEST when > > getting the corresponding matrix. This gives me block (0,0) and block > > (1,1). How do I set/get blocks (0,1) and (1,0)? Looking at ex28 I tried > > MatGetLocalSubMatrix but it gives a null arg... > > > > So the problem

[petsc-users] [KSP] true norm residual

> > I use GMRES with MG Preconditioning and GMRES+PCSHELL as smoother. The > PCSHELL at each level also implies a pre-solve scaling of the solution and > rhs vectors but these are also unscaled in a post-solve step. You cannot safely use GMRES inside GMRES so you should always run this cas

[petsc-users] A few questions about hdf5 viewer

ell HDF5 viewer to use data compression? Thanks for your patience, Best, Mohamad -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/3c7c2fc1/attachment-0001.htm>

[petsc-users] VecView() and binary output

is stored and how it is going help? best Mohamad * * -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120119/b6a9d41b/attachment.htm>