C++ (Cpp) NodeMatrixT 예제들

예제 #1

파일: Tetra3DLagrangeP1.cpp 프로젝트: Ivor23/coolfluid3

void Tetra3DLagrangeP1::mapped_coordinates(const CoordsT& coord, const NodeMatrixT& nodes, MappedCoordsT& mappedCoord)
{
  RealMatrix3 M;
  M.col(0) = nodes.row(1) - nodes.row(0);
  M.col(1) = nodes.row(2) - nodes.row(0);
  M.col(2) = nodes.row(3) - nodes.row(0);
  
  mappedCoord = M.inverse() * (coord - nodes.row(0).transpose());
}

예제 #2

파일: Triag2DLagrangeP2.cpp 프로젝트: Ivor23/coolfluid3

void Triag2DLagrangeP2::mapped_coordinates(const CoordsT& coord, const NodeMatrixT& nodes, MappedCoordsT& map_coord)
{
  throw Common::NotImplemented( FromHere(), "" );

  cf_assert(coord.size() == 2);
  cf_assert(map_coord.size() == 2);
  cf_assert(nodes.size() == 6);

  const Real invDet = 1. / jacobian_determinant(nodes);

  map_coord[KSI] = invDet * ((nodes(2, YY) - nodes(0, YY))*coord[XX] + (nodes(0, XX) - nodes(2, XX))*coord[YY] - nodes(0, XX)*nodes(2, YY) + nodes(2, XX)*nodes(0, YY));
  map_coord[ETA] = invDet * ((nodes(0, YY) - nodes(1, YY))*coord[XX] + (nodes(1, XX) - nodes(0, XX))*coord[YY] + nodes(0, XX)*nodes(1, YY) - nodes(1, XX)*nodes(0, YY));

}

예제 #3

Real Line2DLagrangeP2::area(const NodeMatrixT& nodes)
{
  throw Common::NotImplemented( FromHere(), "" );
  return (nodes.row(1)-nodes.row(0)).norm();
}

예제 #4

파일: Line3DLagrangeP1.cpp 프로젝트: Ivor23/coolfluid3

Real Line3DLagrangeP1::length(const NodeMatrixT& nodes)
{
  return (nodes.row(1)-nodes.row(0)).norm();
}