[petsc-users] Advice on identifying boundary vertices

bel, but I need to discriminate among Face Sets of different strata. 2. what might be wrong in the aforementioned approach? Thanks, Aldo -- Dr. Aldo Bonfiglioli Associate professor of Fluid Machines Scuola di Ingegneria Universita' della Basilicata V.le dell'Ateneo lucano, 10 851

[petsc-users] markers

each side of the square domain has been given a different flag. How do I access the face sets information ? Thanks, Aldo -- Dr. Aldo Bonfiglioli Associate professor of Fluid Machines Scuola di Ingegneria Universita' della Basilicata V.le dell'Ateneo lucano, 10 85100 Potenza ITALY

[petsc-users] Advice on setting Dirichlet bcs in DMPlex

, the former does not. Regards, Aldo -- Dr. Aldo Bonfiglioli Associate professor of Fluid Machines Scuola di Ingegneria Universita' della Basilicata V.le dell'Ateneo lucano, 10 85100 Potenza ITALY tel:+39.0971.205203 fax:+39.0971.205215 web:https://urldefense.us/v3/__http://docen

[petsc-users] Advice on setting BCs in a PLEX

on or evaluating the rhs function I transfer data from the global to local representation, do the calculation, then transfer back. I am doing something wrong, though, which shows up in the attached log file. My simple code is also attached. Thank you for your advice. Aldo -- Dr. Aldo Bonfiglioli Asso

[petsc-users] makefile changes since release 12

. If I remove ${SOURCEF} from the ${PETSC_FCOMPILE} definition in $(PETSC_DIR)/lib/petsc/conf/variables things work, but I am not sure that this is the right thing to do. Thanks, Aldo PS My makefile is attached -- Dr. Aldo Bonfiglioli Associate professor of Fluid Machines Scuola di In

[petsc-users] More on BlockSize in release 3.3

, Aldo -- Dr. Aldo Bonfiglioli Associate professor of Fluid Flow Machinery Scuola di Ingegneria Universita' della Basilicata V.le dell'Ateneo lucano, 10 85100 Potenza ITALY tel:+39.0971.205203 fax:+39.0971.205215 Publications list <http://publicationslist.org/aldo.bonfiglioli> ---

[petsc-users] VecSetBlockSize with release 3.3

,DT,IFAIL) CALL VecSetType(DT,VECSEQ,IFAIL) CALL VecSetBlockSize(DT,NOFVAR,IFAIL) CALL VecSetSizes(DT,NPOIN*NOFVAR,PETSC_DECIDE,IFAIL) Is there a simpler (i.e. less library calls) way to achieve the same result? Regards, Aldo -- Dr. Aldo Bonfiglioli Associate professor of F

[petsc-users] PetscKernel_A_gets_inverse_A_

, but I need to change it as needed when a new PETSc version is released. Regards, Aldo -- Dr. Aldo Bonfiglioli Associate professor of Fluid Flow Machinery Scuola di Ingegneria Universita' della Basilicata V.le dell'Ateneo lucano, 10 85100 Potenza ITALY tel:+39.0971.205203 fax:+39.0

[petsc-users] Interpreting ksp_view

Dear all, I have a doubt concerning the interpretation of the output from -vec_view. I am using ILU(k) on a BAIJ matrix; the run is uni-processor. > KSP Object: > type: gmres > GMRES: restart=30, using Classical (unmodified) Gram-Schmidt > Orthogonalization with no iterative refinement >

[petsc-users] best PETSC solver for long nozzles

iteration. See, e.g.: http://dx.doi.org/10.1016/j.cam.2006.04.068 http://proceedings.fyper.com/eccomascfd2006/documents/398.pdf Aldo Dr. Aldo Bonfiglioli Dip.to di Ingegneria e Fisica dell'Ambiente (DIFA) Universita' della Basilicata V.le dell'Ateneo lucano, 10 85100 Potenza ITALY

Non repeatability issue

t tried yet. I know there exist smoother strategies than SER to rump the continuation parameter (I know this http://www.cs.kuleuven.ac.be/publicaties/rapporten/tw/TW304.ps.gz for instance) Aldo -- Dr. Aldo Bonfiglioli Dip.to di Ingegneria e Fisica dell'Ambiente (DIFA) Universita' del

Non repeatability issue

e ratio between the current and initial residuals of one of the conservation eqns (mass, tipically). In all previous plots, the fixed point iterations were not shown. The fixed point iteration was run only once to generate the initial solution used for all Newton runs. Aldo -- Dr. Aldo Bonfigliol

Non repeatability issue

ill trying to make sure there are no implementation bugs in our code. Aldo -- Dr. Aldo Bonfiglioli Dip.to di Ingegneria e Fisica dell'Ambiente (DIFA) Universita' della Basilicata V.le dell'Ateneo lucano, 10 85100 Potenza ITALY tel:+39.0971.205203 fax:+39.0971.205160

Non repeatability issue

changed the situation much. Out of 4 subsequent runs, 2 converge while 2 enter a limit cycle I had not seen previously (I mean without the -vecscatter_reproduce option). The initial solution is the same as that shown in the other plot. Aldo -- Dr. Aldo Bonfiglioli Dip.to di Ingegneria e Fisica d

Non repeatability issue

problem row scaling without seeing noticeable differences. But I guess this may be different from what you mention. Regards, Aldo -- Dr. Aldo Bonfiglioli Dip.to di Ingegneria e Fisica dell'Ambiente (DIFA) Universita' della Basilicata V.le dell'Ateneo lucano, 10 85100 Potenza ITALY tel:

Non repeatability issue

ns and see whether results > become more stable. Could you comment in particular on this latter? Is there a way we can modify PETSc's behaviour via command line options or otherwise? Regards, Aldo -- Dr. Aldo Bonfiglioli Dip.to di Ingegneria e Fisica dell'Ambiente (DIFA) Universita&

what's the difference between PetscViewerASCIIOpen() and PetscViewerBinaryOpen()?

in how. Regards, Aldo -- Dr. Aldo Bonfiglioli Dip.to di Ingegneria e Fisica dell'Ambiente (DIFA) Universita' della Basilicata V.le dell'Ateneo lucano, 10 85100 Potenza ITALY tel:+39.0971.205203 fax:+39.0971.205160 -- next part -- An embedded and charset

Non repeatability issue

Dear users, in the notes attached here we address an issue concerning nonrepeatability. We know this is a known issue in parallel floating point programs, but we would like to be sure we are interpreting things correctly a not mis-using PETSc. Comments are welcome. Regards, Aldo -- Dr. Aldo