On Thu, Jun 1, 2023 at 1:46 AM Duan Junming via petsc-users < [email protected]> wrote:
Dear all,
I have a simple demo code attached below, which gives a segmentation violation error.
Can you help me with this problem? I think the problem is due to the destroy function.
I am using version 3.19.2 with debugging.
Yes, the check is wrong. I will fix it. For now, you can just pass in the transpose argument as well. Thanks, Matt
#include <petsc.h>
static char help[] = "test.\n"; int main(int argc, char *argv[]) { PetscCall(PetscInitialize(&argc, &argv, 0, help)); PetscScalar *nodes; PetscScalar *weights; PetscScalar **diff; PetscInt n = 3; PetscCall(PetscMalloc2(n, &nodes, n, &weights)); PetscCall(PetscDTGaussLobattoLegendreQuadrature(n, PETSCGAUSSLOBATTOLEGENDRE_VIA_LINEAR_ALGEBRA, nodes, weights)); PetscCall(PetscGaussLobattoLegendreElementGradientCreate(n, nodes, weights, &diff, NULL)); PetscCall(PetscGaussLobattoLegendreElementGradientDestroy(n, nodes, weights, &diff, NULL)); PetscCall(PetscFree2(nodes, weights)); PetscCall(PetscFinalize()); return 0; }
Junming
-- 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 https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>