397 #define MY_MEPS (PETSC_MACHINE_EPSILON * 100)  macro
401   return (va <= vb - MY_MEPS) ? 1 : (va > vb + MY_MEPS) ? -1 : 0; /* 0 for equal */  in gamg_hem_compare()
586       if (iter == 0 && max_e <= MY_MEPS) { // add BCs to fake aggregate  in MatCoarsenApply_HEM_private()
696 …PetscCheck(PetscRealPart(lid_max_ew[lid0]) >= e->weight - MY_MEPS, PETSC_COMM_SELF, PETSC_ERR_SUP,…  in MatCoarsenApply_HEM_private()
699         if (PetscRealPart(lid_max_ew[lid0]) > e->weight + MY_MEPS) {  in MatCoarsenApply_HEM_private()
707           if (PetscRealPart(lid_max_ew[lid1]) > e->weight + MY_MEPS) {  in MatCoarsenApply_HEM_private()
719           if (g_max_e1 > e->weight + MY_MEPS) {  in MatCoarsenApply_HEM_private()
722 …} else if (g_max_e1 >= e->weight - MY_MEPS && lghost_pe[ghost1_idx] > rank) { // is 'lghost_max_pe…  in MatCoarsenApply_HEM_private()
752 …if (ew >= PetscRealPart(lid_max_ew[lid0]) - MY_MEPS && lghost_max_pe[lidj] > rank) isOK = PETSC_FA…  in MatCoarsenApply_HEM_private()
753 …PetscCheck(ew <= max_e + MY_MEPS, PETSC_COMM_SELF, PETSC_ERR_SUP, "edge weight %e > max %e. ncols …  in MatCoarsenApply_HEM_private()
755 …MY_MEPS), ew >= PetscRealPart(lid_max_ew[lid0]) - MY_MEPS, lghost_pe[lidj] > rank, (double)ew, (do…  in MatCoarsenApply_HEM_private()
776 …if (ew >= PetscRealPart(lid_max_ew[lid1]) - MY_MEPS && lghost_max_pe[lidj] > rank) isOK = PETSC_FA…  in MatCoarsenApply_HEM_private()
777 …PetscCheck(ew <= max_e + MY_MEPS, PETSC_COMM_SELF, PETSC_ERR_SUP, "edge weight %e > max %e", (doub…  in MatCoarsenApply_HEM_private()
1027 …    if ((pe = lghost_pe[aj[jj]]) > max_pe && PetscRealPart(ap[jj]) >= max_e - MY_MEPS) max_pe = pe;  in MatCoarsenApply_HEM_private()