4a61042f3f
- Can now retrieve or set a column/row of a tensor.
Either compile-time or run-time checks.
Get
t.col<1>(); t.col(1);
t.row<1>(); t.row(1);
Set
t.col<1>(vec); t.col(1,vec);
t.row<1>(vec); t.row(1,vec);
The templated versions are compile-time checked
t.col<3>();
t.col<3>(vec);
The parameter versions are run-time checked
t.col(3);
t.col(3,vec);
ENH: provide named access to tensor/tensor inner product as inner()
236 lines
6.3 KiB
C
236 lines
6.3 KiB
C
#include "tensor.H"
|
|
#include "triad.H"
|
|
#include "symmTensor.H"
|
|
#include "transform.H"
|
|
#include "stringList.H"
|
|
#include "IOstreams.H"
|
|
|
|
using namespace Foam;
|
|
|
|
int main()
|
|
{
|
|
tensor t1(1, 2, 3, 4, 5, 6, 7, 8, 9);
|
|
tensor t2(1, 2, 3, 1, 2, 3, 1, 2, 3);
|
|
|
|
tensor t3 = t1 + t2;
|
|
|
|
Info<< "tensor " << t1 << " + " << t2 << nl
|
|
<< " = " << t3 << nl << nl;
|
|
|
|
{
|
|
Info<< "rows:" << nl;
|
|
for (direction i=0; i < 3; ++i)
|
|
{
|
|
Info<< " (" << i << ") = " << t3.row(i) << nl;
|
|
}
|
|
}
|
|
{
|
|
Info<< "cols:" << nl;
|
|
for (direction i=0; i < 3; ++i)
|
|
{
|
|
Info<< " (" << i << ") = " << t3.col(i) << nl;
|
|
}
|
|
Info<< "col<0> = " << t3.col<0>() << nl;
|
|
Info<< "col<1> = " << t3.col<1>() << nl;
|
|
Info<< "col<2> = " << t3.col<2>() << nl;
|
|
// Compilation error: Info << "col<3> = " << t3.col<3>() << nl;
|
|
|
|
t3.col<1>({0, 2, 4});
|
|
Info<< "replaced col<1> = " << t3.col<1>() << nl;
|
|
Info<< "tensor " << t3 << nl;
|
|
|
|
t3.row<2>(Zero);
|
|
Info<< "replaced row<1> = " << t3.row<2>() << nl;
|
|
Info<< "tensor " << t3 << nl;
|
|
|
|
triad tr3(t3);
|
|
Info<< "triad " << tr3 << " :: T() " << tr3.T() << nl;
|
|
}
|
|
Info<< nl;
|
|
|
|
|
|
tensor t4(3,-2,1,-2,2,0,1, 0, 4);
|
|
|
|
Info<< inv(t4) << endl;
|
|
Info<< (inv(t4) & t4) << endl;
|
|
|
|
Info<< t1.x() << t1.y() << t1.z() << endl;
|
|
|
|
tensor t6(1,0,-4,0,5,4,-4,4,3);
|
|
//tensor t6(1,2,0,2,5,0,0,0,0);
|
|
|
|
Info<< "tensor " << t6 << endl;
|
|
vector e = eigenValues(t6);
|
|
Info<< "eigenvalues " << e << endl;
|
|
|
|
tensor ev = eigenVectors(t6);
|
|
Info<< "eigenvectors " << ev << endl;
|
|
|
|
Info<< "Check determinant " << e.x()*e.y()*e.z() << " " << det(t6) << endl;
|
|
|
|
Info<< "Check eigenvectors "
|
|
<< (eigenVectors(t6, e) & t6) << " "
|
|
<< (e.x() * eigenVectors(t6, e)).x()
|
|
<< (e.x() * eigenVectors(t6, e)).y()
|
|
<< (e.z() * eigenVectors(t6, e)).z()
|
|
<< endl;
|
|
|
|
Info<< "Check eigenvalues for symmTensor "
|
|
<< eigenValues(symm(t6)) - eigenValues(tensor(symm(t6))) << endl;
|
|
|
|
Info<< "Check eigenvectors for symmTensor "
|
|
<< eigenVectors(symm(t6)) - eigenVectors(tensor(symm(t6))) << endl;
|
|
|
|
tensor t7(1, 2, 3, 2, 4, 5, 3, 5, 6);
|
|
|
|
Info<< "Check transformation "
|
|
<< (t1 & t7 & t1.T()) << " " << transform(t1, t7) << endl;
|
|
|
|
symmTensor st1(1, 2, 3, 4, 5, 6);
|
|
symmTensor st2(7, 8, 9, 10, 11, 12);
|
|
|
|
Info<< "Check symmetric transformation "
|
|
<< transform(t1, st1) << endl;
|
|
|
|
Info<< "Check dot product of symmetric tensors "
|
|
<< (st1 & st2) << endl;
|
|
|
|
Info<< "Check inner sqr of a symmetric tensor "
|
|
<< innerSqr(st1) << " " << innerSqr(st1) - (st1 & st1) << endl;
|
|
|
|
Info<< "Check symmetric part of dot product of symmetric tensors "
|
|
<< twoSymm(st1&st2) - ((st1&st2) + (st2&st1)) << endl;
|
|
|
|
tensor sk1 = skew(t6);
|
|
Info<< "Check dot product of symmetric and skew tensors "
|
|
<< twoSymm(st1&sk1) - ((st1&sk1) - (sk1&st1)) << endl;
|
|
|
|
vector v1(1, 2, 3);
|
|
|
|
Info<< sqr(v1) << endl;
|
|
Info<< symm(t7) << endl;
|
|
Info<< twoSymm(t7) << endl;
|
|
Info<< magSqr(st1) << endl;
|
|
Info<< mag(st1) << endl;
|
|
|
|
Info<< (symm(t7) && t7) - (0.5*(t7 + t7.T()) && t7) << endl;
|
|
Info<< (t7 && symm(t7)) - (t7 && 0.5*(t7 + t7.T())) << endl;
|
|
|
|
|
|
/*
|
|
// Lots of awkward eigenvector tests ...
|
|
|
|
tensor T_rand_real
|
|
(
|
|
0.9999996423721313, 0.3330855667591095, 0.6646450161933899,
|
|
0.9745196104049683, 0.0369445420801640, 0.0846728682518005,
|
|
0.6474838852882385, 0.1617118716239929, 0.2041363865137100
|
|
);
|
|
DebugVar(T_rand_real);
|
|
vector L_rand_real(eigenValues(T_rand_real));
|
|
DebugVar(L_rand_real);
|
|
tensor U_rand_real(eigenVectors(T_rand_real));
|
|
DebugVar(U_rand_real);
|
|
|
|
Info << endl << endl;
|
|
|
|
tensor T_rand_imag
|
|
(
|
|
0.8668024539947510, 0.1664607226848602, 0.8925783634185791,
|
|
0.9126510620117188, 0.7408077120780945, 0.1499115079641342,
|
|
0.0936608463525772, 0.7615650296211243, 0.8953040242195129
|
|
);
|
|
DebugVar(T_rand_imag);
|
|
vector L_rand_imag(eigenValues(T_rand_imag));
|
|
DebugVar(L_rand_imag);
|
|
tensor U_rand_imag(eigenVectors(T_rand_imag));
|
|
DebugVar(U_rand_imag);
|
|
|
|
Info << endl << endl;
|
|
|
|
tensor T_rand_symm
|
|
(
|
|
1.9999992847442627, 1.3076051771640778, 1.3121289014816284,
|
|
1.3076051771640778, 0.0738890841603279, 0.2463847398757935,
|
|
1.3121289014816284, 0.2463847398757935, 0.4082727730274200
|
|
);
|
|
DebugVar(T_rand_symm);
|
|
vector L_rand_symm(eigenValues(T_rand_symm));
|
|
DebugVar(L_rand_symm);
|
|
tensor U_rand_symm(eigenVectors(T_rand_symm));
|
|
DebugVar(U_rand_symm);
|
|
|
|
Info << endl << endl;
|
|
|
|
symmTensor T_rand_Symm
|
|
(
|
|
1.9999992847442627, 1.3076051771640778, 1.3121289014816284,
|
|
0.0738890841603279, 0.2463847398757935,
|
|
0.4082727730274200
|
|
);
|
|
DebugVar(T_rand_Symm);
|
|
vector L_rand_Symm(eigenValues(T_rand_Symm));
|
|
DebugVar(L_rand_Symm);
|
|
tensor U_rand_Symm(eigenVectors(T_rand_Symm));
|
|
DebugVar(U_rand_Symm);
|
|
|
|
Info << endl << endl;
|
|
|
|
tensor T_rand_diag
|
|
(
|
|
0.8668024539947510, 0, 0,
|
|
0, 0.7408077120780945, 0,
|
|
0, 0, 0.8953040242195129
|
|
);
|
|
DebugVar(T_rand_diag);
|
|
vector L_rand_diag(eigenValues(T_rand_diag));
|
|
DebugVar(L_rand_diag);
|
|
tensor U_rand_diag(eigenVectors(T_rand_diag));
|
|
DebugVar(U_rand_diag);
|
|
|
|
Info << endl << endl;
|
|
|
|
tensor T_repeated
|
|
(
|
|
0, 1, 1,
|
|
1, 0, 1,
|
|
1, 1, 0
|
|
);
|
|
DebugVar(T_repeated);
|
|
vector L_repeated(eigenValues(T_repeated));
|
|
DebugVar(L_repeated);
|
|
tensor U_repeated(eigenVectors(T_repeated));
|
|
DebugVar(U_repeated);
|
|
|
|
Info << endl << endl;
|
|
|
|
tensor T_repeated_zero
|
|
(
|
|
1, 1, 1,
|
|
1, 1, 1,
|
|
1, 1, 1
|
|
);
|
|
DebugVar(T_repeated_zero);
|
|
vector L_repeated_zero(eigenValues(T_repeated_zero));
|
|
DebugVar(L_repeated_zero);
|
|
tensor U_repeated_zero(eigenVectors(T_repeated_zero));
|
|
DebugVar(U_repeated_zero);
|
|
|
|
Info << endl << endl;
|
|
|
|
tensor T_triple
|
|
(
|
|
2, 0, 0,
|
|
0, 2, 0,
|
|
0, 0, 2
|
|
);
|
|
DebugVar(T_triple);
|
|
vector L_triple(eigenValues(T_triple));
|
|
DebugVar(L_triple);
|
|
tensor U_triple(eigenVectors(T_triple));
|
|
DebugVar(U_triple);
|
|
*/
|
|
|
|
return 0;
|
|
}
|