OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Added relative errors to tests/TestContactAngle

Browse Source
This commit is contained in:
James McClure 2014-12-30 10:45:17 -05:00
parent 68c6240e7f
commit 60539cb0f9
1 changed files with 9 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
tests/TestContactAngle.cpp
View File

@ -200,29 +200,29 @@ int main (int argc, char *argv[])
printf("-------------------------------- \n");
printf("NWP volume = %f \n", nwp_volume);
analytical = 2*PI*(CAPHEIGHT*CAPHEIGHT+RADIUS*RADIUS);
RelError = (awn-analytical)/analytical;
RelError = fabs(awn-analytical)/analytical;
printf("Area wn = %f, Analytical = %f, Rel. Error = %f \n", awn,analytical,RelError);
analytical = 2*PI*RADIUS*(N-2)-4*PI*RADIUS*(CAPRAD-CAPHEIGHT);
RelError = (ans-analytical)/analytical;
RelError = fabs(ans-analytical)/analytical;
printf("Area ns = %f, Analytical = %f, Rel. Error = %f \n", ans, analytical,RelError);
analytical = 2*PI*RADIUS*(CAPRAD-CAPHEIGHT);
RelError = (aws-analytical)/analytical;
RelError = fabs(aws-analytical)/analytical;
printf("Area ws = %f, Analytical = %f, Rel. Error = %f \n", aws, analytical,RelError);
analytical = 2*PI*RADIUS*(N-2);
RelError = (As-analytical)/analytical;
RelError = fabs(As-analytical)/analytical;
printf("Area s = %f, Analytical = %f, Rel. Error = %f \n", As, analytical,RelError);
analytical = 4*PI*RADIUS;
RelError = (lwns-analytical)/analytical;
RelError = fabs(lwns-analytical)/analytical;
printf("Length wns = %f, Analytical = %f, Rel. Error = %f \n", lwns, analytical,RelError);
analytical = 1.0*RADIUS/CAPRAD;
RelError = (efawns-analytical)/analytical;
RelError = fabs(efawns-analytical)/analytical;
printf("Cos(theta_wns) = %f, Analytical = %f, Rel. Error = %f \n",efawns,analytical,RelError);
analytical = 2.0/RADIUS;
RelError = (Jwn-analytical)/analytical;
analytical = 2.0/CAPRAD;
RelError = fabs(Jwn-analytical)/analytical;
printf("Mean curvature (wn) = %f, Analytical = %f, Rel. Error = %f \n", Jwn, analytical,RelError);
printf("Geodesic curvature (wns) = %f, Analytical, Rel. Error = %f = %f \n", KGwns, 0.0, 0.0);
analytical = 1.0/RADIUS;
RelError = (KNwns-analytical)/analytical;
RelError = fabs(KNwns-analytical)/analytical;
printf("Normal curvature (wns) = %f, Analytical = %f, Rel. Error = %f \n", KNwns, analytical,RelError);
printf("-------------------------------- \n");
//.........................................................................