[skip ci] And we add another solution...

test/sumfact/poisson/facedir-facemod-variation/sumfact_poisson_dg_3d_facedir_facemod_variation_driver.cc

@@ -253,6 +253,8 @@ int main(int argc, char** argv){
 
       std::vector<RangeType> solution_2 {-0.22216923038681022, -0.003662938655271696, -0.23743729880973413, -0.0030364756689728193, -0.29395296456019571, -0.0018314693276358439, -0.29973173153779054, -0.0015182378344864129, 0.0036629386552717003, 0.2309838664360746, 0.0030364756689728349, 0.24202149140895873, 0.001831469327635845, 0.28737146661901281, 0.001518237834486417, 0.29291440083048437};
 
+      std::vector<RangeType> solution_3 {-0.21436796383145873, -0.23383260730195662, -0.28383755411465533, -0.29526590985481704, -0.0045507639649666034, -0.0034803883238202818, -0.0036071199470257055, -0.0024060631441813391, 0.22702984288940109, 0.0045507639649666042, 0.27429663847448937, 0.0036071199470257059, 0.23989650875067281, 0.0034803883238202818, 0.28608104498832443, 0.0024060631441813391};
+
       std::cout.precision(17);
       std::vector<RangeType> residual(16);
       for (std::size_t i=0; i<16; ++i){
@@ -268,7 +270,9 @@ int main(int argc, char** argv){
                 << fuzzy_is_permutation(solution_1, residual) << std::endl;
       std::cout << "palpo is_permuation to solution_2: "
                 << fuzzy_is_permutation(solution_2, residual) << std::endl;
-      if (!fuzzy_is_permutation(solution_1, residual) and !fuzzy_is_permutation(solution_2, residual) )
+      std::cout << "palpo is_permuation to solution_3: "
+                << fuzzy_is_permutation(solution_3, residual) << std::endl;
+      if (!fuzzy_is_permutation(solution_1, residual) and !fuzzy_is_permutation(solution_2, residual) and !fuzzy_is_permutation(solution_3, residual))
         return 1;
     }
     if (initree.get<bool>("printmatrix", false)) {