[petsc-users] Functions for judging the properties of the matrix

f you want estimates for the original operator. http://www.mcs.anl.gov/petsc/documentation/faq.html#conditionnumber -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/ffec87ae/attachment.html>

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

n F90Array4dAccess 1>libpetsc.lib(f90_cwrap.o) : error LNK2019: unresolved external symbol F90ARRAY4dDESTROYSCALAR referenced in function F90Array4dDestroy 1>c:\obj_tmp\dm_test2d\Debug/dm_test2d.exe : fatal error LNK1120: 6 unresolved externals 1> 1>Build log written to "file://c:\obj_tmp\dm_test2d\Debug\BuildLog.htm"/ Any other solutions? -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/54007e5e/attachment.html>

[petsc-users] Draw '2D' Vector (FFT)

Hello list, what is the easiest way to 'draw' (PETSC_VIEWER_DRAW_WORLD) a 2D plot of a petsc vector which was created with MatGetVecsFFTW()? The employed FFT matrix has DIM=2 and therefore the mentioned vector has N1*N2 components. As far as I know VecView(x, PETSC_VIEWER_DRAW_WORLD) only creates

[petsc-users] Functions for judging the properties of the matrix

. Jim -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/0d26346a/attachment.html>

[petsc-users] Does this mean the matrix is ill-conditioned?

Hi, all, The following is the output from the solution of a Poisson equation from Darcy's law. To compute the condition number of matrix, I did not use PC and use GMRES KSP to do the test. It seems like that the condition number keep increasing during the iterative solution. Does this mean the m

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

t;>>> symbol >>>>>> F90ARRAY4dCREATESCALAR referenced in function F90Array4dCreate >>>>>> 1>libpetsc.lib(f90_cwrap.o) : error LNK2019: unresolved external >>>>>> symbol >>>>>> F90ARRAY4dACCESSFORTRANADDR referenced in function F90Array4dAccess >>>>>> 1>libpetsc.lib(f90_cwrap.o) : error LNK2019: unresolved external >>>>>> symbol >>>>>> F90ARRAY4dACCESSINT referenced in function F90Array4dAccess >>>>>> 1>libpetsc.lib(f90_cwrap.o) : error LNK2019: unresolved external >>>>>> symbol >>>>>> F90ARRAY4dACCESSREAL referenced in function F90Array4dAccess >>>>>> 1>libpetsc.lib(f90_cwrap.o) : error LNK2019: unresolved external >>>>>> symbol >>>>>> F90ARRAY4dACCESSSCALAR referenced in function F90Array4dAccess >>>>>> 1>libpetsc.lib(f90_cwrap.o) : error LNK2019: unresolved external >>>>>> symbol >>>>>> F90ARRAY4dDESTROYSCALAR referenced in function F90Array4dDestroy >>>>>> 1>c:\obj_tmp\ibm2d_high_Re_staggered_old\Debug/ibm2d_high_Re_staggered.exe >>>>>> : >>>>>> fatal error LNK1120: 6 unresolved externals/ >>>>>> >>>>>> I wonder if it's fixed in the new petsc-dev. >>>>>> >>>>>> Thanks >>>>>>> On Jul 6, 2012, at 6:30 AM, Blaise Bourdin wrote: >>>>>>> >>>>>>>> Hi, >>>>>>>> >>>>>>>> I have added the creation, destruction and accessor functions for >>>>>>>> 4d >>>>>>>> vectors in F90. The accessor was missing and needed for >>>>>>>> DMDAVecGetArrayF90 >>>>>>>> with a 3d DMDA and>1 dof. As far as I can test, ex11f90 in DM >>>>>>>> should >>>>>>>> now >>>>>>>> completely work with the intel compilers. >>>>>>>> >>>>>>>> Some of the functions are probably not used (F90Array4dAccessReal, >>>>>>>> F90Array4dAccessInt, F90Array4dAccessFortranAddr, for instance), >>>>>>>> but I >>>>>>>> added them anyway. Let me know if you want me to submit a patch >>>>>>>> without >>>>>>>> them. >>>>>>>> >>>>>>>> Regards, >>>>>>>> >>>>>>>> Blaise >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> -- >>>>>>>> Department of Mathematics and Center for Computation& Technology >>>>>>>> Louisiana State University, Baton Rouge, LA 70803, USA >>>>>>>> Tel. +1 (225) 578 1612, Fax +1 (225) 578 4276 >>>>>>>> http://www.math.lsu.edu/~bourdin >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/8634863b/attachment.html>

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

On 11/7/2012 5:12 PM, Satish Balay wrote: > You would revert the old one - and apply the new one. i.e > > patch -Np1 -R< DMDAVecRestoreArrayF90-2.patch > patch -Np1< DMDAVecRestoreArrayF90-2.1.patch > > Satish Hi, Just tried but I still get the error: $ !494 patch -p1 < DMDAVecRestoreArrayF90

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

f the functions are probably not used (F90Array4dAccessReal, >>>> F90Array4dAccessInt, F90Array4dAccessFortranAddr, for instance), but I >>>> added them anyway. Let me know if you want me to submit a patch without >>>> them. >>>> >>>> Regards, >>>> >>>> Blaise >>>> >>>> >>>> >>>> -- >>>> Department of Mathematics and Center for Computation& Technology >>>> Louisiana State University, Baton Rouge, LA 70803, USA >>>> Tel. +1 (225) 578 1612, Fax +1 (225) 578 4276 >>>> http://www.math.lsu.edu/~bourdin >>>> >>>> >>>> >>>> >>>> >>>> >>>> >> >> -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/e0b3ff2c/attachment-0001.html>

[petsc-users] boomeramg memory usage

317 49 32 45, I www.marin.nl > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/409c35d9/attachment.html>

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

I added the patch to both 3.3 and petsc-dev repos. Hopefully thats easier to use. Satish On Wed, 11 Jul 2012, Blaise Bourdin wrote: > Rectification: clone petsc-3.3, not petsc-dev, of course... > B > > > There is something fishy in your tree or your build: the patch that Satish > > sent addre

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

ath.lsu.edu/~bourdin -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/aa046b61/attachment.html>

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

state: - rename your current petsc-dev tree - clone petsc-dev from the official tree (BTW, do you really need petsc-dev, and not 3.3?) - apply the patch Satish sent you or, quoting Satich >> edit src/sys/f90-src/f90_cwrap.c and replace all occurrences of 'F90ARRAY4d' >> with 'F90ARRAY4D' - run configure with your favorite options, then make This should fix your problem. Blaise -- Department of Mathematics and Center for Computation & Technology Louisiana State University, Baton Rouge, LA 70803, USA Tel. +1 (225) 578 1612, Fax +1 (225) 578 4276 http://www.math.lsu.edu/~bourdin -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/cde280b1/attachment.html>

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

Matt -- 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/20120711/4183179c/attachment-0001.html>

[petsc-users] Draw '2D' Vector (FFT)

or has > N1*N2 components. As far as I know VecView(x, PETSC_VIEWER_DRAW_WORLD) > only creates a lineplot. > Can the Vec be associated with a DM? If so, it will do a 2D plot. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/piper

[petsc-users] Does this mean the matrix is ill-conditioned?

ults to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/be6c890b/attachment.html>

[petsc-users] SLEPc mantain eigenvalue/eigenvector ordering compatible with problem matrix

El 11/07/2012, a las 13:05, Alessio Cardillo escribi?: > Hello to all, > > I am using SLEPc library to solve an eigenvalue problem using the lapack > solver. > > the solver works fine (both eigenvalue and eigenvectors are correct) > but the order in which solutions are given is different from

[petsc-users] boomeramg memory usage

11, F +31 317 49 32 45, I www.marin.nl > -- Ir. B. Sanderse Centrum Wiskunde en Informatica Science Park 123 1098 XG Amsterdam t: +31 20 592 4161 e: sanderse at cwi.nl -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/a7ed92ba/attachment.html>

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

On Wed, 11 Jul 2012, TAY wee-beng wrote: > > I just tried on a different system and the same error occurs. The 1st command > : > / > / > > /patch -Np1 -R < DMDAVecRestoreArrayF90-2.patch > > $ patch -Np1 -R < DMDAVecRestoreArrayF90-2.patch > (Stripping trailing CRs from patch.) > patching file

[petsc-users] SLEPc mantain eigenvalue/eigenvector ordering compatible with problem matrix

Hello to all, I am using SLEPc library to solve an eigenvalue problem using the lapack solver. the solver works fine (both eigenvalue and eigenvectors are correct) but the order in which solutions are given is different from the original one. My problem is: I need to retrieve the pair (eigenval

[petsc-users] boomeramg memory usage

1 317 49 39 11, F +31 317 49 32 45, I www.marin.nl > -- Ir. B. Sanderse Centrum Wiskunde en Informatica Science Park 123 1098 XG Amsterdam t: +31 20 592 4161 e: sanderse at cwi.nl -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/61ebfd3e/attachment.html>

[petsc-users] About DIVERGED_ITS

Jim -- >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/20120711/7200307f/attachment-0001.html>

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

;> Department of Mathematics and Center for Computation & Technology >> Louisiana State University, Baton Rouge, LA 70803, USA >> Tel. +1 (225) 578 1612, Fax +1 (225) 578 4276 >> http://www.math.lsu.edu/~bourdin >> >> >> >> >> >> >> -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/2b38f398/attachment.html>

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

Hm - if you have a mercurial clone of petsc-dev - you don't need to apply/revert DMDAVecRestoreArrayF90-2.patch. You can just apply DMDAVecRestoreArrayF90-2.1.patch. But I'm not sure why you are getting errors. >> sbalay at ps3 ~/petsc-dev $ hg revert -a sbalay at ps3 ~/petsc-dev $ h

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

You would revert the old one - and apply the new one. i.e patch -Np1 -R < DMDAVecRestoreArrayF90-2.patch patch -Np1 < DMDAVecRestoreArrayF90-2.1.patch Satish On Wed, 11 Jul 2012, TAY wee-beng wrote: > > > On 11/7/2012 3:40 PM, Satish Balay wrote: > > The second patch has typos. Attaching the

[petsc-users] KSPSetOperators/sub-preconditioner type

1958d66 > Awesome, thanks! -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/f20c47a6/attachment.html>

[petsc-users] Functions for judging the properties of the matrix

he original operator. http://www.mcs.anl.gov/petsc/documentation/faq.html#conditionnumber -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/2731f34e/attachment.html>

[petsc-users] boomeramg memory usage

I'm solving a 3D Poisson equation on a 200x200x100 grid using CG and algebraic multigrid as preconditioner. I noticed that with boomeramg, the memory increases as the number of procs increases: 1 proc:2.8G 2 procs: 2.8G 4 procs: 5.2G 8 procs: >12G (max reached, swapping) The memory usag

[petsc-users] KSPSetOperators/sub-preconditioner type

-pc_type bjacobi -sub_pc_type ilu > -ksp_view -ksp_monitor -sub_pc_factor_mat_solver_package superlu > -sub_pc_factor_levels 1 > > Is the correct solution just to call KSPSetFromOptions every time? Will > that be noticeable at all? I hope this was less vague than my original > message. Thanks again. > > Best, > > Paul > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120711/81bf2eff/attachment.html>

[petsc-users] [petsc-dev] second patch for DMDAVecGetArrayF90

The second patch has typos. Attaching the modified patch. If it works for you - I'll add it to petsc-3.3/petsc-dev thanks, Satish On Wed, 11 Jul 2012, TAY wee-beng wrote: > On 6/7/2012 5:48 PM, Barry Smith wrote: > > Blaise, > > > >Thanks. > > > > Satish, > > > >If th