63   PetscInt        dim, Nv, Neq, c, v;  in DMPlexLabelToMetricConstraint()  local
74   Neq = (dim * (dim + 1)) / 2;  in DMPlexLabelToMetricConstraint()
75   PetscCall(PetscMalloc1(PetscSqr(Neq), &eqns));  in DMPlexLabelToMetricConstraint()
76   PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, Neq, Neq, eqns, &A));  in DMPlexLabelToMetricConstraint()
94           eqns[0 * Neq + r] = PetscSqr(e[0]);  in DMPlexLabelToMetricConstraint()
95           eqns[1 * Neq + r] = 2.0 * e[0] * e[1];  in DMPlexLabelToMetricConstraint()
96           eqns[2 * Neq + r] = PetscSqr(e[1]);  in DMPlexLabelToMetricConstraint()
99           eqns[0 * Neq + r] = PetscSqr(e[0]);  in DMPlexLabelToMetricConstraint()
100           eqns[1 * Neq + r] = 2.0 * e[0] * e[1];  in DMPlexLabelToMetricConstraint()
101           eqns[2 * Neq + r] = 2.0 * e[0] * e[2];  in DMPlexLabelToMetricConstraint()
102           eqns[3 * Neq + r] = PetscSqr(e[1]);  in DMPlexLabelToMetricConstraint()
103           eqns[4 * Neq + r] = 2.0 * e[1] * e[2];  in DMPlexLabelToMetricConstraint()
104           eqns[5 * Neq + r] = PetscSqr(e[2]);  in DMPlexLabelToMetricConstraint()