polyhedron tests

algoim/organizer/organizer.hpp

@@ -88,7 +88,8 @@ public:
 
             auto compositeTensorStdVector = xarray2StdVector(compositeTensor);
 
-            uvector<real, 3> testX(0.8, 0.8, 0.8);
+            // uvector<real, 3> testX(0.8, 0.8, 0.8);
+            uvector<real, 3> testX(0., 0., 0.);
             real             textEvalPower = evalPower(compositeTensor, testX);
 
             tensor3 phi(nullptr, 1 + faceCount);
@@ -102,6 +103,7 @@ public:
                 auto realX = x * range + xmin;
                 if (isInsidePowers(planeTensors, realX)) volume += w * integrand(realX);
             });
+            std::cout << "textEvalPower: " << textEvalPower << std::endl;
             std::cout << "quadraturePointCount: " << quadraturePointCount << std::endl;
         }
         volume *= pow(xmax - xmin, 3);