Function check_material_parameters

tools/Assembler/Materials/OrthotropicMaterial.h:13–30

namespace OrthotropicMaterialHelper {
    void check_material_parameters(size_t dim, const VectorF& young_modulus, 
            const VectorF& poisson_ratio, const VectorF& shear_modulus) {
        if (dim != 2 && dim != 3) {
            throw NotImplementedError("Only 2D and 3D materials are supported");
        }
        if (dim != young_modulus.size()) {
            throw RuntimeError("Wrong dimention for young's modulus");
        }
        if (dim*(dim-1)/2 != shear_modulus.size()) {
            throw RuntimeError("Wrong dimention for shear modulus");
        }
        if (young_modulus.minCoeff() < 0.0) {
            throw RuntimeError("Young's modulus must be positive!");
        }
        if (shear_modulus.minCoeff() < 0.0) {
            throw RuntimeError("Shear modulus must be positive!");
        }
    }

    bool is_symmetric(const MatrixF& M) {
        const Float TOL = 1e-6;