CentralDifferencesJacobianAssembler.h

Go to the documentation of this file.

 
#pragma once
 
#include <memory>
#include "AbstractJacobianAssembler.h"
 
namespace BaseLib
{
class ConfigTree;
}  // BaseLib
 
namespace ProcessLib
{
class CentralDifferencesJacobianAssembler final
    : public AbstractJacobianAssembler
{
public:
    explicit CentralDifferencesJacobianAssembler(
        std::vector<double>&& absolute_epsilons);
 
    void assembleWithJacobian(LocalAssemblerInterface& local_assembler,
                              double const t, double const dt,
                              std::vector<double> const& local_x,
                              std::vector<double> const& local_xdot,
                              const double dxdot_dx, const double dx_dx,
                              std::vector<double>& local_M_data,
                              std::vector<double>& local_K_data,
                              std::vector<double>& local_b_data,
                              std::vector<double>& local_Jac_data) override;
 
private:
    std::vector<double> const _absolute_epsilons;
 
    // temporary data only stored here in order to avoid frequent memory
    // reallocations.
    std::vector<double> _local_M_data;
    std::vector<double> _local_K_data;
    std::vector<double> _local_b_data;
    std::vector<double> _local_x_perturbed_data;
    std::vector<double> _local_xdot_perturbed_data;
};
 
std::unique_ptr<CentralDifferencesJacobianAssembler>
createCentralDifferencesJacobianAssembler(BaseLib::ConfigTree const& config);
 
}  // namespace ProcessLib