Re: [petsc-users] Inquiry about PetscDTSimplexQuadrature .

2023-06-24 Thread Matthew Knepley
On Fri, Jun 23, 2023 at 10:12 PM neil liu wrote: > Thanks, Matt. > It seems DMPlexComputeCellGeometryFEM works well with the quadrature > points to deliver Jacobian and inverse one. > Will it be a good choice ? > Yes, it is intended to compute these quantities for you. Let me know if it does

Re: [petsc-users] Inquiry about PetscDTSimplexQuadrature .

2023-06-23 Thread neil liu
Thanks, Matt. It seems DMPlexComputeCellGeometryFEM works well with the quadrature points to deliver Jacobian and inverse one. Will it be a good choice ? Have a good night. Thanks, On Fri, Jun 23, 2023 at 10:33 AM Matthew Knepley wrote: > On Fri, Jun 23, 2023 at 10:28 AM neil liu wrote: >

Re: [petsc-users] Inquiry about PetscDTSimplexQuadrature .

2023-06-23 Thread Matthew Knepley
On Fri, Jun 23, 2023 at 10:28 AM neil liu wrote: > Dear Petsc developers, > > I am learning *PetscDTSimplexQuadrature *and plan to use it. > I found that, in the biunit simplex (tetra), > (-1,-1,-1),(1,-1,-1),(-1,1,-1),(-1,-1,1), > degree 1: npoints 4, the sum of weights = 4/3(the volume of

[petsc-users] Inquiry about PetscDTSimplexQuadrature .

2023-06-23 Thread neil liu
Dear Petsc developers, I am learning *PetscDTSimplexQuadrature *and plan to use it. I found that, in the biunit simplex (tetra), (-1,-1,-1),(1,-1,-1),(-1,1,-1),(-1,-1,1), degree 1: npoints 4, the sum of weights = 4/3(the volume of this simplex) degree 2 : npoints 8; For my previous experience, I