Revision 1f037930ff0f8bb09bb223b084fcd5f17decaceb authored by FZill on 17 December 2021, 10:30:00 UTC, committed by FZill on 20 December 2021, 08:08:45 UTC
- THM: added computeElasticTangentStiffness
- THM: added ElasticTangentStiffness to getBulkModulus
- T/THM: updating BGRaCreepAndInitialStressAtIP_AREHS

the solid_material fix slightly changed the test result,
the relative differences to the data sets of the previous result
are about 1e-7
1 parent a4cabe1
Raw File
AutoCheckGenerators.cpp
/**
 * \file
 * \copyright
 * Copyright (c) 2012-2021, OpenGeoSys Community (http://www.opengeosys.org)
 *            Distributed under a Modified BSD License.
 *              See accompanying file LICENSE.txt or
 *              http://www.opengeosys.org/project/license
 */

#include "Tests/GeoLib/AutoCheckGenerators.h"

namespace autocheck
{
// reflect point p on the point c in x-y plane
MathLib::Point3d reflect(MathLib::Point3d const& c, MathLib::Point3d const& p)
{
    return MathLib::Point3d(
        std::array<double, 3>{{2 * c[0] - p[0], 2 * c[1] - p[1], 0.0}});
}

GeoLib::LineSegment translate(Eigen::Vector3d const& translation,
                              GeoLib::LineSegment const& line_seg)
{
    auto a = std::make_unique<GeoLib::Point>(line_seg.getBeginPoint());
    auto b = std::make_unique<GeoLib::Point>(line_seg.getEndPoint());
    for (std::size_t k(0); k < 3; ++k)
    {
        (*a)[k] += translation[k];
    }
    for (std::size_t k(0); k < 3; ++k)
    {
        (*b)[k] += translation[k];
    }
    return GeoLib::LineSegment{a.release(), b.release(), true};
}

}  // namespace autocheck
back to top