129   PetscInt             jj, tid, tt, idx, nselected_2;  in triangulateAndFormProl()  local
308               for (tt = 0; tt < 3; tt++) {  in triangulateAndFormProl()
309                 PetscInt cid2 = mid.trianglelist[3 * tid + tt];  in triangulateAndFormProl()
311                 AA[tt][0]     = coords[lid2];  in triangulateAndFormProl()
312                 AA[tt][1]     = coords[data_stride + lid2];  in triangulateAndFormProl()
313                 AA[tt][2]     = 1.0;  in triangulateAndFormProl()
314                 clids[tt]     = cid2; /* store for interp */  in triangulateAndFormProl()
317               for (tt = 0; tt < 3; tt++) alpha[tt] = (PetscScalar)fcoord[tt];  in triangulateAndFormProl()
324                 for (tt = 0, idx = 0; tt < 3; tt++) {  in triangulateAndFormProl()
325 … if (PetscRealPart(alpha[tt]) > (1.0 + EPS) || PetscRealPart(alpha[tt]) < -EPS) have = PETSC_FALSE;  in triangulateAndFormProl()
326                   if (PetscRealPart(alpha[tt]) < lowest) {  in triangulateAndFormProl()
327                     lowest = PetscRealPart(alpha[tt]);  in triangulateAndFormProl()
328                     idx    = tt;  in triangulateAndFormProl()
339                 for (tt = 0; tt < 3; tt++) {  in triangulateAndFormProl()
340                   PetscInt cid2 = mid.trianglelist[3 * tid + tt];  in triangulateAndFormProl()
342                   AA[tt][0]     = coords[lid2];  in triangulateAndFormProl()
343                   AA[tt][1]     = coords[data_stride + lid2];  in triangulateAndFormProl()
344                   AA[tt][2]     = 1.0;  in triangulateAndFormProl()
345                   clids[tt]     = cid2; /* store for interp */  in triangulateAndFormProl()
347                 for (tt = 0; tt < 3; tt++) alpha[tt] = fcoord[tt];  in triangulateAndFormProl()
353                   for (tt = 0; tt < 3 && have; tt++) {  in triangulateAndFormProl()
354 …   if (PetscRealPart(alpha[tt]) > 1.0 + EPS || PetscRealPart(alpha[tt]) < -EPS) have = PETSC_FALSE;  in triangulateAndFormProl()
355                     if ((v = PetscAbs(PetscRealPart(alpha[tt]) - 0.5)) > worst) worst = v;  in triangulateAndFormProl()
368               for (tt = 0; tt < 3; tt++) {  in triangulateAndFormProl()
369                 PetscInt cid2 = mid.trianglelist[3 * bestTID + tt];  in triangulateAndFormProl()
371                 AA[tt][0]     = coords[lid2];  in triangulateAndFormProl()
372                 AA[tt][1]     = coords[data_stride + lid2];  in triangulateAndFormProl()
373                 AA[tt][2]     = 1.0;  in triangulateAndFormProl()
374                 clids[tt]     = cid2; /* store for interp */  in triangulateAndFormProl()
376               for (tt = 0; tt < 3; tt++) alpha[tt] = fcoord[tt];  in triangulateAndFormProl()
387 …for (tt = 0; tt < bs; tt++, ii++, jj++) PetscCall(MatSetValues(a_Prol, 1, &ii, 1, &jj, &shp, INSER…  in triangulateAndFormProl()