IShape3D.cpp File Reference

Description

Implements default methods of interface IShape3D.

Homepage:\n http://www.bornagainproject.org

License:\n GNU General Public License v3 or higher (see COPYING)

Copyright

Forschungszentrum Jülich GmbH 2018

Authors

Scientific Computing Group at MLZ (see CITATION, AUTHORS)

Definition in file IShape3D.cpp.

Include dependency graph for IShape3D.cpp:

Go to the source code of this file.

Functions
std::vector< R3 >	EllipseVerticesXtrunc (double x, double r_y, double r_z, double z_b, double z_t)
	Generate vertices of centered vertical ellipse with given semi-axes at position x and truncation heights z_b (bottom) and z_t (top) More...
std::vector< R3 >	EllipseVerticesZ (double r_x, double r_y, double z)
	Generate vertices of centered horizontal ellipse with given semi-axes at height z. More...
std::vector< R3 >	RectangleVertices (double length, double width, double z)
	Helper functions to construct lists of vertices. More...

Function Documentation

EllipseVerticesXtrunc()

std::vector<R3> EllipseVerticesXtrunc	(	double	x,
		double	r_y,
		double	r_z,
		double	z_b,
		double	z_t
	)

Generate vertices of centered vertical ellipse with given semi-axes at position x and truncation heights z_b (bottom) and z_t (top)

Definition at line 50 of file IShape3D.cpp.

{
    static constexpr double delta_angle = 2.0 * M_PI / IShape3D::N_Circle;
    // for full ellipse
    //    std::vector<R3> result(IShape3D::N_Circle);
    //    for (size_t i = 0; i < IShape3D::N_Circle; ++i) {
    //        double angle = i * delta_angle;
    //        double y = r_y * std::cos(angle);
    //        double z = r_z * std::sin(angle) + r_z; // bottom is at 0, top at 2*r_z
    //        result[i] = R3(x, y, z);
    //    }
 
    // for truncated ellipse
    std::vector<R3> result;
    result.reserve(IShape3D::N_Circle + 4);
    for (size_t i = 0; i < IShape3D::N_Circle; ++i) {
        double angle = i * delta_angle;
        double y = r_y * std::cos(angle);
        double z = r_z * std::sin(angle) - z_b; // bottom is at 0, top at z_t-z_b
 
        if (0 <= z && z <= (z_t - z_b)) {
            result.push_back(R3(x, y, z));
        }
    }
    double alpha_t = atan(z_t / sqrt(r_z * r_z - z_t * z_t) * r_z / r_y);
    double alpha_b = atan(z_b / sqrt(r_z * r_z - z_b * z_b) * r_z / r_y);
 
    double y_t = r_y * std::cos(alpha_t);
    double y_b = r_y * std::cos(alpha_b);
 
    result.push_back(R3(x, y_t, z_t - z_b));
    result.push_back(R3(x, -y_t, z_t - z_b));
    result.push_back(R3(x, y_b, 0));
    result.push_back(R3(x, -y_b, 0));
 
    return result;
}

References M_PI, and IShape3D::N_Circle.

Referenced by DoubleEllipseX::DoubleEllipseX().

EllipseVerticesZ()

std::vector<R3> EllipseVerticesZ	(	double	r_x,
		double	r_y,
		double	z
	)

Generate vertices of centered horizontal ellipse with given semi-axes at height z.

Definition at line 37 of file IShape3D.cpp.

{
    static constexpr double delta_angle = 2.0 * M_PI / IShape3D::N_Circle;
    std::vector<R3> result(IShape3D::N_Circle);
    for (size_t i = 0; i < IShape3D::N_Circle; ++i) {
        double angle = i * delta_angle;
        double x = r_x * std::cos(angle);
        double y = r_y * std::sin(angle);
        result[i] = R3(x, y, z);
    }
    return result;
}

References M_PI, and IShape3D::N_Circle.

Referenced by DoubleEllipseZ::DoubleEllipseZ(), and TruncatedEllipsoidNet::TruncatedEllipsoidNet().

RectangleVertices()

std::vector<R3> RectangleVertices	(	double	length,
		double	width,
		double	z
	)

Helper functions to construct lists of vertices.

Generate vertices of centered rectangle at height z

Definition at line 28 of file IShape3D.cpp.

{
    std::vector<R3> result = {{length / 2.0, width / 2.0, z},
                              {-length / 2.0, width / 2.0, z},
                              {-length / 2.0, -width / 2.0, z},
                              {length / 2.0, -width / 2.0, z}};
    return result;
}

Referenced by BoxNet::BoxNet().