Searched refs:pmet (Results 1 – 1 of 1) sorted by relevance
144 PetscScalar *pmet; in ComputeMetric() local157 PetscCall(DMPlexPointLocalRef(dm, v, met, &pmet)); in ComputeMetric()161 if (i == 0) pmet[i * dim + j] = 1 / (h * h); in ComputeMetric()162 else pmet[i * dim + j] = lambda; in ComputeMetric()163 } else pmet[i * dim + j] = 0.0; in ComputeMetric()