mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
a18c15a4cf
- as per std::array and FixedList STYLE: rename Scalar.[CH] -> scalarImpl.[CH] (these are internal files) - added inclusion guards to protect against bad use
77 lines
1.4 KiB
C
77 lines
1.4 KiB
C
#include "vector.H"
|
|
#include "IOstreams.H"
|
|
#include "vectorTools.H"
|
|
#include "unitConversion.H"
|
|
|
|
using namespace Foam;
|
|
|
|
|
|
void test(const vector& a, const vector& b, const scalar tolerance)
|
|
{
|
|
Info<< "Vectors " << a << " and " << b
|
|
<< " are (to tolerance of " << tolerance << "): ";
|
|
|
|
if (vectorTools::areParallel(a, b, tolerance))
|
|
Info<< " parallel ";
|
|
|
|
if (vectorTools::areOrthogonal(a, b, tolerance))
|
|
Info<< " orthogonal ";
|
|
|
|
if (vectorTools::areAcute(a, b))
|
|
Info<< " acute ";
|
|
|
|
if (vectorTools::areObtuse(a, b))
|
|
Info<< " obtuse ";
|
|
|
|
Info<< ", angle = " << vectorTools::degAngleBetween(a, b);
|
|
|
|
Info<< endl;
|
|
}
|
|
|
|
|
|
int main()
|
|
{
|
|
vector a(vector::uniform(1));
|
|
vector b;
|
|
b.fill(2);
|
|
|
|
test(a, b, 0.0);
|
|
test(a, b, VSMALL);
|
|
test(a, b, SMALL);
|
|
test(a, b, 1e-3);
|
|
test(a, b, 1e-1);
|
|
|
|
a = vector(1,0,0);
|
|
b = vector(0,2,0);
|
|
|
|
test(a, b, 0.0);
|
|
test(a, b, VSMALL);
|
|
test(a, b, SMALL);
|
|
test(a, b, 1e-3);
|
|
test(a, b, 1e-1);
|
|
|
|
a = vector(1,0,0);
|
|
b = vector(-1,0,0);
|
|
|
|
test(a, b, 0.0);
|
|
test(a, b, VSMALL);
|
|
test(a, b, SMALL);
|
|
test(a, b, 1e-3);
|
|
test(a, b, 1e-1);
|
|
|
|
a = vector(1,0,0);
|
|
b = vector(-1,2,0);
|
|
|
|
test(a, b, 0.0);
|
|
test(a, b, VSMALL);
|
|
test(a, b, SMALL);
|
|
test(a, b, 1e-3);
|
|
test(a, b, 1e-1);
|
|
|
|
Info<< "\nEnd\n" << nl;
|
|
|
|
return 0;
|
|
}
|
|
|
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|