void generate(size_t start, size_t end) const
{
float g[3][3];
float gTranpose[3][3];
float direction[3] = {0.0, 0.0, 0.0};
// Geneate all the Coordinates
for(size_t i = start; i < end; ++i)
{
FOrientArrayType eu(m_Eulers->getPointer(i * 3), 3);
FOrientArrayType om(9, 0.0);
OrientationTransforms<FOrientArrayType, float>::eu2om(eu, om);
om.toGMatrix(g);
MatrixMath::Transpose3x3(g, gTranpose);
// -----------------------------------------------------------------------------
// 001 Family
direction[0] = 0.0;
direction[1] = 0.0;
direction[2] = 1.0;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz001->getPointer(i * 6));
MatrixMath::Copy3x1(m_xyz001->getPointer(i * 6), m_xyz001->getPointer(i * 6 + 3));
MatrixMath::Multiply3x1withConstant(m_xyz001->getPointer(i * 6 + 3), -1);
// -----------------------------------------------------------------------------
// 011 Family
direction[0] = 1.0;
direction[1] = 0.0;
direction[2] = 0.0;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz011->getPointer(i * 6));
MatrixMath::Copy3x1(m_xyz011->getPointer(i * 6), m_xyz011->getPointer(i * 6 + 3));
MatrixMath::Multiply3x1withConstant(m_xyz011->getPointer(i * 6 + 3), -1);
// -----------------------------------------------------------------------------
// 111 Family
direction[0] = 0.0;
direction[1] = 1.0;
direction[2] = 0;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz111->getPointer(i * 6));
MatrixMath::Copy3x1(m_xyz111->getPointer(i * 6), m_xyz111->getPointer(i * 6 + 3));
MatrixMath::Multiply3x1withConstant(m_xyz111->getPointer(i * 6 + 3), -1);
}
}
void generate(size_t start, size_t end) const
{
float g[3][3];
float gTranpose[3][3];
float* currentEuler = NULL;
float direction[3] = {0.0, 0.0, 0.0};
for(size_t i = start; i < end; ++i)
{
currentEuler = eulers->GetPointer(i * 3);
OrientationMath::EulerToMat(currentEuler[0], currentEuler[1], currentEuler[2], g);
MatrixMath::Transpose3x3(g, gTranpose);
// -----------------------------------------------------------------------------
// 001 Family
direction[0] = 0.0;
direction[1] = 0.0;
direction[2] = 1.0;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz001->GetPointer(i * 6));
MatrixMath::Copy3x1(m_xyz001->GetPointer(i * 6), m_xyz001->GetPointer(i * 6 + 3));
MatrixMath::Multiply3x1withConstant(m_xyz001->GetPointer(i * 6 + 3), -1);
// -----------------------------------------------------------------------------
// 011 Family
direction[0] = 1.0;
direction[1] = 0.0;
direction[2] = 0.0;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz011->GetPointer(i * 6));
MatrixMath::Copy3x1(m_xyz011->GetPointer(i * 6), m_xyz011->GetPointer(i * 6 + 3));
MatrixMath::Multiply3x1withConstant(m_xyz011->GetPointer(i * 6 + 3), -1);
// -----------------------------------------------------------------------------
// 111 Family
direction[0] = 0.0;
direction[1] = 1.0;
direction[2] = 0;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz111->GetPointer(i * 6));
MatrixMath::Copy3x1(m_xyz111->GetPointer(i * 6), m_xyz111->GetPointer(i * 6 + 3));
MatrixMath::Multiply3x1withConstant(m_xyz111->GetPointer(i * 6 + 3), -1);
}
}
void generate(size_t start, size_t end) const
{
float g[3][3];
float gTranpose[3][3];
float direction[3] = {0.0, 0.0, 0.0};
for(size_t i = start; i < end; ++i)
{
FOrientArrayType eu(m_Eulers->getPointer(i * 3), 3);
FOrientArrayType om(9, 0.0);
OrientationTransforms<FOrientArrayType, float>::eu2om(eu, om);
om.toGMatrix(g);
MatrixMath::Transpose3x3(g, gTranpose);
// -----------------------------------------------------------------------------
// 001 Family
direction[0] = 1.0;
direction[1] = 0.0;
direction[2] = 0.0;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz001->getPointer(i * 18));
MatrixMath::Copy3x1(m_xyz001->getPointer(i * 18), m_xyz001->getPointer(i * 18 + 3));
MatrixMath::Multiply3x1withConstant(m_xyz001->getPointer(i * 18 + 3), -1);
direction[0] = 0.0;
direction[1] = 1.0;
direction[2] = 0.0;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz001->getPointer(i * 18 + 6));
MatrixMath::Copy3x1(m_xyz001->getPointer(i * 18 + 6), m_xyz001->getPointer(i * 18 + 9));
MatrixMath::Multiply3x1withConstant(m_xyz001->getPointer(i * 18 + 9), -1);
direction[0] = 0.0;
direction[1] = 0.0;
direction[2] = 1.0;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz001->getPointer(i * 18 + 12));
MatrixMath::Copy3x1(m_xyz001->getPointer(i * 18 + 12), m_xyz001->getPointer(i * 18 + 15));
MatrixMath::Multiply3x1withConstant(m_xyz001->getPointer(i * 18 + 15), -1);
// -----------------------------------------------------------------------------
// 011 Family
direction[0] = SIMPLib::Constants::k_1OverRoot2;
direction[1] = SIMPLib::Constants::k_1OverRoot2;
direction[2] = 0.0;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz011->getPointer(i * 36));
MatrixMath::Copy3x1(m_xyz011->getPointer(i * 36), m_xyz011->getPointer(i * 36 + 3));
MatrixMath::Multiply3x1withConstant(m_xyz011->getPointer(i * 36 + 3), -1);
direction[0] = SIMPLib::Constants::k_1OverRoot2;
direction[1] = 0.0;
direction[2] = SIMPLib::Constants::k_1OverRoot2;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz011->getPointer(i * 36 + 6));
MatrixMath::Copy3x1(m_xyz011->getPointer(i * 36 + 6), m_xyz011->getPointer(i * 36 + 9));
MatrixMath::Multiply3x1withConstant(m_xyz011->getPointer(i * 36 + 9), -1);
direction[0] = 0.0;
direction[1] = SIMPLib::Constants::k_1OverRoot2;
direction[2] = SIMPLib::Constants::k_1OverRoot2;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz011->getPointer(i * 36 + 12));
MatrixMath::Copy3x1(m_xyz011->getPointer(i * 36 + 12), m_xyz011->getPointer(i * 36 + 15));
MatrixMath::Multiply3x1withConstant(m_xyz011->getPointer(i * 36 + 15), -1);
direction[0] = -SIMPLib::Constants::k_1OverRoot2;
direction[1] = -SIMPLib::Constants::k_1OverRoot2;
direction[2] = 0.0;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz011->getPointer(i * 36 + 18));
MatrixMath::Copy3x1(m_xyz011->getPointer(i * 36 + 18), m_xyz011->getPointer(i * 36 + 21));
MatrixMath::Multiply3x1withConstant(m_xyz011->getPointer(i * 36 + 21), -1);
direction[0] = -SIMPLib::Constants::k_1OverRoot2;
direction[1] = 0.0;
direction[2] = SIMPLib::Constants::k_1OverRoot2;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz011->getPointer(i * 36 + 24));
MatrixMath::Copy3x1(m_xyz011->getPointer(i * 36 + 24), m_xyz011->getPointer(i * 36 + 27));
MatrixMath::Multiply3x1withConstant(m_xyz011->getPointer(i * 36 + 27), -1);
direction[0] = 0.0;
direction[1] = -SIMPLib::Constants::k_1OverRoot2;
direction[2] = SIMPLib::Constants::k_1OverRoot2;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz011->getPointer(i * 36 + 30));
MatrixMath::Copy3x1(m_xyz011->getPointer(i * 36 + 30), m_xyz011->getPointer(i * 36 + 33));
MatrixMath::Multiply3x1withConstant(m_xyz011->getPointer(i * 36 + 33), -1);
// -----------------------------------------------------------------------------
// 111 Family
direction[0] = SIMPLib::Constants::k_1OverRoot3;
direction[1] = SIMPLib::Constants::k_1OverRoot3;
direction[2] = SIMPLib::Constants::k_1OverRoot3;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz111->getPointer(i * 24));
MatrixMath::Copy3x1(m_xyz111->getPointer(i * 24), m_xyz111->getPointer(i * 24 + 3));
MatrixMath::Multiply3x1withConstant(m_xyz111->getPointer(i * 24 + 3), -1);
direction[0] = -SIMPLib::Constants::k_1OverRoot3;
direction[1] = SIMPLib::Constants::k_1OverRoot3;
direction[2] = SIMPLib::Constants::k_1OverRoot3;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz111->getPointer(i * 24 + 6));
MatrixMath::Copy3x1(m_xyz111->getPointer(i * 24 + 6), m_xyz111->getPointer(i * 24 + 9));
MatrixMath::Multiply3x1withConstant(m_xyz111->getPointer(i * 24 + 9), -1);
direction[0] = SIMPLib::Constants::k_1OverRoot3;
direction[1] = -SIMPLib::Constants::k_1OverRoot3;
direction[2] = SIMPLib::Constants::k_1OverRoot3;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz111->getPointer(i * 24 + 12));
MatrixMath::Copy3x1(m_xyz111->getPointer(i * 24 + 12), m_xyz111->getPointer(i * 24 + 15));
MatrixMath::Multiply3x1withConstant(m_xyz111->getPointer(i * 24 + 15), -1);
direction[0] = SIMPLib::Constants::k_1OverRoot3;
direction[1] = SIMPLib::Constants::k_1OverRoot3;
direction[2] = -SIMPLib::Constants::k_1OverRoot3;
MatrixMath::Multiply3x3with3x1(gTranpose, direction, m_xyz111->getPointer(i * 24 + 18));
MatrixMath::Copy3x1(m_xyz111->getPointer(i * 24 + 18), m_xyz111->getPointer(i * 24 + 21));
MatrixMath::Multiply3x1withConstant(m_xyz111->getPointer(i * 24 + 21), -1);
}
}
