int main(int /*arc*/, char** /*argv*/) {
double const PI = 3.1415926535897932;
const MyVector<double> radii(MV::fill, 3.,4.,5.);
double AnaVol = (4./3.)*PI*(radii[0]*radii[0]*radii[0] + \
radii[1]*radii[1]*radii[1] + \
radii[2]*radii[2]*radii[2]);
double CalcVol = TotalVolumeOfSpheresFromRadii(radii);
std::cout << "Analytical Volume Sum = " << AnaVol << std::endl;
std::cout << "Calculated Volume Sum = " << CalcVol << std::endl;
UtilsForTesting u;
IS_TRUE(CalcVol-AnaVol<1.0e-10,"Error in TotalVolumeOfSpheresFromRadii!! Calculated Volume =/= real volume!! You stupid bitch!!");
return u.NumberOfTestsFailed();
//return EXIT_SUCCESS;
}
int main(){
const MyVector<double> A(MV::fill,1.0,1.0,1.0,1.0);
double TotalV = TotalVolumeOfSpheresFromRadii(A);
std::cout<< "Total Volume = "<< TotalV <<std::endl;
}
