Revision 6cfd76397657f949473c3ee5f3402e22d8cfbc5a authored by Dmitry Yu. Naumov on 23 January 2021, 11:56:43 UTC, committed by Dmitry Yu. Naumov on 23 January 2021, 11:56:43 UTC
Fixed a bug that skips output at a accept time step with rejections before. See merge request ogs/ogs!3416
CreateJacobianAssembler.cpp
/**
* \file
* \copyright
* Copyright (c) 2012-2021, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*
*/
#include "CreateJacobianAssembler.h"
#include "BaseLib/Error.h"
#include "AnalyticalJacobianAssembler.h"
#include "CentralDifferencesJacobianAssembler.h"
#include "CompareJacobiansJacobianAssembler.h"
namespace ProcessLib
{
std::unique_ptr<AbstractJacobianAssembler> createJacobianAssembler(
boost::optional<BaseLib::ConfigTree> const& config)
{
if (!config)
{
return std::make_unique<AnalyticalJacobianAssembler>();
}
//! \ogs_file_param{prj__processes__process__jacobian_assembler__type}
auto const type = config->peekConfigParameter<std::string>("type");
if (type == "Analytical") {
config->ignoreConfigParameter("type");
return std::make_unique<AnalyticalJacobianAssembler>();
}
if (type == "CentralDifferences")
{
return createCentralDifferencesJacobianAssembler(*config);
}
if (type == "CompareJacobians")
{
return createCompareJacobiansJacobianAssembler(*config);
}
OGS_FATAL("Unknown Jacobian assembler type: `{:s}'.", type);
}
} // namespace ProcessLib
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...