28 "Vector of Real values providing cut information");
39 const int cut_data_len = 9;
43 mooseError(
"Length of CircleCutUserObject cut_data must be 9");
52 _normal = rays.first.cross(rays.second);
55 std::pair<Real, Real> ray_radii =
58 if (
std::abs(ray_radii.first - ray_radii.second) > 1e-10)
59 mooseError(
"CircleCutUserObject only works for a circular cut");
61 _radius = 0.5 * (ray_radii.first + ray_radii.second);
62 _angle = std::acos((rays.first * rays.second) / (ray_radii.first * ray_radii.second));
74 const std::vector<Point>
77 std::vector<Point> crack_front_points(number_crack_front_points);
83 for (
unsigned int i = 0; i < number_crack_front_points; ++i)
86 crack_front_points[i] =
90 return crack_front_points;
93 const std::vector<RealVectorValue>
96 mooseError(
"getCrackPlaneNormals() is not implemented for this object.");
ADRealEigenVector< T, D, asd > sqrt(const ADRealEigenVector< T, D, asd > &)
virtual const std::vector< Point > getCrackFrontPoints(unsigned int num_crack_front_points) const override
get a set of points along a crack front from a XFEM GeometricCutUserObject
virtual bool isInsideCutPlane(Point p) const override
ADRealEigenVector< T, D, asd > abs(const ADRealEigenVector< T, D, asd > &)
registerMooseObject("XFEMApp", CircleCutUserObject)
static InputParameters validParams()
CircleCutUserObject(const InputParameters ¶meters)
std::vector< Real > _cut_data
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
std::vector< Point > _vertices
void mooseError(Args &&... args) const
void normalizePoint(Point &p)
virtual const std::vector< RealVectorValue > getCrackPlaneNormals(unsigned int num_crack_front_points) const override
get a set of normal vectors along a crack front from a XFEM GeometricCutUserObject ...
static InputParameters validParams()