JSONUtils.cpp

Go to the documentation of this file.

#include "JSONUtils.hpp"
 
#include <polyfem/utils/StringUtils.hpp>
#include <polyfem/utils/Logger.hpp>
 
#include <fstream>
#include <filesystem>
 
#include <jse/jse.h>
 
#include <Eigen/Geometry>
 
namespace polyfem
{
    namespace utils
    {
        void apply_common_params(json &args)
        {
            if (!args.contains("common"))
                return;
 
            const std::string common_params_path = resolve_path(args["common"], args["root_path"]);
 
            if (common_params_path.empty())
                return;
 
            std::ifstream file(common_params_path);
            if (!file.is_open())
                log_and_throw_error("Unable to open common params {} file", common_params_path);
 
            json common_params;
            file >> common_params;
            file.close();
 
            // Recursively apply common params
            const bool has_root_path = common_params.contains("root_path");
            if (has_root_path)
                common_params["root_path"] = resolve_path(common_params["root_path"], common_params_path);
            else
                common_params["root_path"] = common_params_path;
            apply_common_params(common_params);
 
            // If there is a root path in the common params, it overrides the one in the current params.
            // This is somewhat backwards as normally current params override common params, but this is
            // an easy way to make sure that the relative paths in common are correct.
            if (has_root_path)
                args["root_path"] = common_params["root_path"];
 
            json patch;
            if (args.contains("patch"))
            {
                patch = args["patch"];
                args.erase("patch");
            }
 
            common_params.merge_patch(args);
            if (!patch.empty())
                common_params = common_params.patch(patch);
            args = common_params;
 
            args.erase("common"); // Remove common params from the final json
        }
 
        void expand_bc_sidecars(json &args, const json &rules)
        {
            if (!args.contains("boundary_conditions"))
                return;
 
            json &bcs = args["boundary_conditions"];
            if (!bcs.contains("dirichlet_boundary") || !bcs["dirichlet_boundary"].is_array())
                return;
 
            const json root_path = args.contains("root_path") ? args["root_path"] : json("");
 
            json expanded = json::array();
            for (const auto &entry : bcs["dirichlet_boundary"])
            {
                if (entry.is_string())
                {
                    const std::string path = resolve_path(entry, root_path);
                    if (std::filesystem::path(path).extension() == ".json")
                    {
                        std::ifstream file(path);
                        if (!file.is_open())
                            log_and_throw_error("Unable to open dirichlet_boundary {} file", path);
 
                        json sidecar;
                        try
                        {
                            file >> sidecar;
                        }
                        catch (const std::exception &e)
                        {
                            log_and_throw_error("Invalid JSON in dirichlet_boundary file {}: {}", path, e.what());
                        }
                        file.close();
 
                        if (!sidecar.is_array())
                            log_and_throw_error("dirichlet_boundary file {} must contain an array", path);
 
                        jse::JSE jse;
                        for (const auto &e : sidecar)
                        {
                            json filled = e;
                            expanded.push_back(jse.inject_defaults(filled, rules));
                        }
                        continue;
                    }
                }
                expanded.push_back(entry);
            }
            bcs["dirichlet_boundary"] = expanded;
        }
 
        Eigen::Matrix3d to_rotation_matrix(const json &jr, std::string mode)
        {
            std::transform(mode.begin(), mode.end(), mode.begin(), ::tolower);
 
            if (jr.is_array() && jr.empty())
            {
                return Eigen::Matrix3d::Identity(3, 3);
            }
 
            Eigen::VectorXd r;
            if (jr.is_number())
            {
                r.setZero(3);
                assert(mode.size() == 1); // must be either "x", "y", or "z"
                int i = mode[0] - 'x';
                assert(i >= 0 && i < 3);
                r[i] = jr.get<double>();
            }
            else
            {
                assert(jr.is_array());
                r = jr;
            }
 
            if (mode == "axis_angle")
            {
                assert(r.size() == 4);
                double angle = deg2rad(r[0]); // NOTE: assumes input angle is in degrees
                Eigen::Vector3d axis = r.tail<3>().normalized();
                return Eigen::AngleAxisd(angle, axis).toRotationMatrix();
            }
 
            if (mode == "quaternion")
            {
                assert(r.size() == 4);
                Eigen::Vector4d q = r.normalized();
                return Eigen::Quaterniond(q).toRotationMatrix();
            }
 
            // The following expect the input is given in degrees
            r = deg2rad(r);
 
            if (mode == "rotation_vector")
            {
                assert(r.size() == 3);
                double angle = r.norm();
                if (angle != 0)
                {
                    return Eigen::AngleAxisd(angle, r / angle).toRotationMatrix();
                }
                else
                {
                    return Eigen::Matrix3d::Identity();
                }
            }
 
            Eigen::Matrix3d R = Eigen::Matrix3d::Identity();
 
            for (int i = 0; i < mode.size(); i++)
            {
                int j = mode[i] - 'x';
                assert(j >= 0 && j < 3);
                Eigen::Vector3d axis = Eigen::Vector3d::Zero();
                axis[j] = 1;
                R = Eigen::AngleAxisd(r[j], axis).toRotationMatrix() * R;
            }
 
            return R;
        }
 
        bool is_param_valid(const json &params, const std::string &key)
        {
            return !params.is_null() && params.contains(key) && !params[key].is_null();
        }
    } // namespace utils
} // namespace polyfem