ForwardDifferencesJacobianAssembler.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include "NumericalJacobianAssembler.h"
 
namespace ProcessLib
{
class ForwardDifferencesJacobianAssembler final
    : public NumericalJacobianAssembler
{
public:
    explicit ForwardDifferencesJacobianAssembler(
        std::vector<double>&& absolute_epsilons);

    void assembleWithJacobian(std::size_t const mesh_item_id,
                              LocalAssemblerInterface& local_assembler,
                              double const t, double const dt,
                              std::vector<double> const& local_x_data,
                              std::vector<double> const& local_x_prev_data,
                              std::vector<double>& local_b_data,
                              std::vector<double>& local_Jac_data) override;
 
    std::unique_ptr<AbstractJacobianAssembler> copy() const override;
 
private:
    // 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;
};
 
}  // namespace ProcessLib