test_teqpcpp.cpp

Go to the documentation of this file.

 
#include <iostream>
#include "teqp/cpp/teqpcpp.hpp"
#include "teqp/cpp/derivs.hpp"
 
int main() {
    nlohmann::json j = { 
        {"kind", "multifluid"}, 
        {"model", {
            {"components", {"../mycp/dev/fluids/Methane.json","../mycp/dev/fluids/Ethane.json"}},
            {"BIP", "../mycp/dev/mixtures/mixture_binary_pairs.json"},
            {"departure", "../mycp/dev/mixtures/mixture_departure_functions.json"}
        }
    }};
    //std::cout << j.dump(2);
    auto am = teqp::cppinterface::make_model(j);
 
    auto z = (Eigen::ArrayXd(2) << 0.5, 0.5).finished();
    double Ar01 = am->get_Arxy(0, 1, 300, 3, z);
    std::cout << Ar01 << std::endl;
    
    auto fugcoeff = am->get_fugacity_coefficients(300.0, z*300.0);
    std::cout << fugcoeff << std::endl;
    
    auto vhat = am->get_partial_molar_volumes(300.0, z*300.0);
    std::cout << vhat << std::endl;
    
    auto mat = am->get_deriv_mat2(300.0, 300.0, z);
    std::cout << mat << std::endl;
    
    const std::vector<char> vars = {'T','D','P','S'};
    auto im = teqp::cppinterface::build_iteration_Jv(vars, mat, mat, 8.3144, 300.0, 300.0, z);
    std::cout << im.J << std::endl;
    
    const auto& p = std::get<teqp::multifluid_t>(am->get_model());
    std::cout << p.get_meta() << std::endl;
    
}