nonlinear_solve.hpp

Go to the documentation of this file.

// Copyright (c) Lawrence Livermore National Security, LLC and
// other Smith Project Developers. See the top-level LICENSE file for
// details.
//
// SPDX-License-Identifier: (BSD-3-Clause)
 
#pragma once
 
#include <vector>
#include "smith/differentiable_numerics/field_state.hpp"
 
namespace smith {
 
class WeakForm;
class DifferentiableSolver;
class DifferentiableBlockSolver;
class BoundaryConditionManager;
class DirichletBoundaryConditions;
 
static constexpr size_t invalid_block_index = std::numeric_limits<size_t>::max() - 1;
 
FieldState solve(const WeakForm& residual_eval, const FieldState& shape_disp, const std::vector<FieldState>& states,
                 const std::vector<FieldState>& params, const TimeInfo& time_info, const DifferentiableSolver& solver,
                 const DirichletBoundaryConditions& bcs, size_t unknown_state_index = 0);
 
// with unknowns (with respect to the solver) being a, and b.
// r1 has unknowns a,b in the ‘slots’ 0, 1
// r2 has unknowns a,b, in the ‘slots’ 3,0
std::vector<FieldState> block_solve(const std::vector<WeakForm*>& residual_evals,
                                    const std::vector<std::vector<size_t>> block_indices, const FieldState& shape_disp,
                                    const std::vector<std::vector<FieldState>>& states,
                                    const std::vector<std::vector<FieldState>>& params, const TimeInfo& time_info,
                                    const DifferentiableBlockSolver* solver,
                                    const std::vector<const BoundaryConditionManager*>& bc_managers);
 
}  // namespace smith