[petsc-users] Vector field ordering question

Dear development team, I'm trying to write a code to deal with a multi-field problem. Currently I find that the vector ordering for the field is (a1,a2,a3,b1,b2,b3), here assume a1,a2,a3 belongs to one field and bs for another field. Can I get a point-block ordering for the unknowns as (a1,b1,a

Re: [petsc-users] R: How to use Intel OneApi mpi wrappers on Linux

Hi, fwiw, I tried to compile with ipcx too, without mpi wrappers... However, I had other problems... check here: https://gitlab.com/petsc/petsc/-/issues/1255 Anyone have compiled PETSc with the latest Intel OneAPI release? Can you give a working configure line? Thanks, Eric On 2022-10-03

Re: [petsc-users] PetscLogView produces nan's instead of timing data when using GPUs

It prints Nan to indicate that the time for that event is not known accurately. But the times for the larger events that contain these events are known. So for example the time for KSPSolve is know but not the time for VecNorm. The other numbers in the events, like number of times called etc

[petsc-users] PetscLogView produces nan's instead of timing data when using GPUs

Hi PETSc-developers, I'm having trouble with getting performance logs from an application that uses PETSc. There are no issues when I run it on a CPU, but every time a GPU is used there is no timing data and almost all times are replaced by times that are just `nan` (on two different clusters). I

[petsc-users] About Q3 tensor product Hermite element

Dear all, I need to use Q3 tensor product Hermite element in 2D (point value, gradient, and mixed derivative at 4 vertices in a cell as unknowns). Is it available in PETSc FEM module now? I found that only Lagrange element is available. If not, what is the correct path to implement Q3 tensor

Re: [petsc-users] clarification on extreme eigenvalues from KSPComputeEigenvalues

On Tue, Oct 4, 2022 at 5:20 PM feng wang wrote: > Hi Mark, > > Thanks for your reply. Below is the output if I call KSPComputeEigenvalues > > 0.330475 -0.0485014 > 0.521211 0.417409 > 0.684726 -0.377126 > 0.885941 0.354342 > 0.957845 -0.0508471 > 0.964676 -0.241642 > 1.05921 0.0742963 > 1.82065 -

Re: [petsc-users] code with TS throws error at the end

Can you try running with valgrind? https://petsc.org/release/faq/?highlight=valgrind#what-does-corrupt-argument-or-caught-signal-or-segv-or-segmentation-violation-or-bus-error-mean-can-i-use-valgrind-or-cuda-memcheck-to-debug-memory-corruption-issues

Re: [petsc-users] Is the results of `DMAdaptLabel` as expected in `src/dm/impls/plex/tests/ex20.c`

On Wed, Oct 5, 2022 at 7:40 AM Zongze Yang wrote: > Matthew Knepley 于2022年10月5日周三 00:33写道： > >> On Tue, Oct 4, 2022 at 3:19 PM Zongze Yang wrote: >> >>> Hi everyone, >>> >>> I am learning how to use the `DMAdaptLabel` for `DMPlex`, and found the >>> example `src/dm/impls/plex/tests/ex20.c` whic