#include <transient_rb_construction.h>

Inheritance diagram for libMesh::TransientRBConstruction:

[legend]

Public Types
typedef TransientRBConstruction	sys_type

typedef TransientSystem< RBConstruction >	Parent

typedef std::map< std::string, SparseMatrix< Number > * >::iterator	matrices_iterator

typedef std::map< std::string, SparseMatrix< Number > * >::const_iterator	const_matrices_iterator

typedef std::map< std::string, NumericVector< Number > * >::iterator	vectors_iterator

typedef std::map< std::string, NumericVector< Number > * >::const_iterator	const_vectors_iterator

Public Member Functions
	TransientRBConstruction (EquationSystems &es, const std::string &name, const unsigned int number)

virtual	~TransientRBConstruction ()

virtual void	clear () libmesh_override

virtual void	initialize_rb_construction (bool skip_matrix_assembly=false, bool skip_vector_assembly=false) libmesh_override

virtual Real	truth_solve (int write_interval) libmesh_override

virtual Real	train_reduced_basis (const bool resize_rb_eval_data=true) libmesh_override

virtual void	process_parameters_file (const std::string &parameters_filename) libmesh_override

virtual void	print_info () libmesh_override

virtual bool	greedy_termination_test (Real abs_greedy_error, Real initial_greedy_error, int count) libmesh_override

virtual void	assemble_all_affine_operators () libmesh_override

virtual void	assemble_misc_matrices () libmesh_override

void	assemble_L2_matrix (SparseMatrix< Number > *input_matrix, bool apply_dirichlet_bc=true)

void	assemble_mass_matrix (SparseMatrix< Number > *input_matrix)

void	add_scaled_mass_matrix (Number scalar, SparseMatrix< Number > *input_matrix)

void	mass_matrix_scaled_matvec (Number scalar, NumericVector< Number > &dest, NumericVector< Number > &arg)

void	set_L2_assembly (ElemAssembly &L2_assembly_in)

ElemAssembly &	get_L2_assembly ()

void	assemble_Mq_matrix (unsigned int q, SparseMatrix< Number > *input_matrix, bool apply_dirichlet_bc=true)

SparseMatrix< Number > *	get_M_q (unsigned int q)

SparseMatrix< Number > *	get_non_dirichlet_M_q (unsigned int q)

virtual void	get_all_matrices (std::map< std::string, SparseMatrix< Number > * > &all_matrices) libmesh_override

virtual void	truth_assembly () libmesh_override

int	get_max_truth_solves () const

void	set_max_truth_solves (int max_truth_solves_in)

Real	get_POD_tol () const

void	set_POD_tol (const Real POD_tol_in)

void	set_delta_N (const unsigned int new_delta_N)

virtual void	load_rb_solution () libmesh_override

const NumericVector< Number > &	get_error_temporal_data ()

void	update_RB_initial_condition_all_N ()

virtual void	write_riesz_representors_to_files (const std::string &riesz_representors_dir, const bool write_binary_residual_representors) libmesh_override

virtual void	read_riesz_representors_from_files (const std::string &riesz_representors_dir, const bool write_binary_residual_representors) libmesh_override

sys_type &	system ()

virtual void	reinit () libmesh_override

virtual std::string	system_type () const libmesh_override

Number	old_solution (const dof_id_type global_dof_number) const

Number	older_solution (const dof_id_type global_dof_number) const

virtual void	solve_for_matrix_and_rhs (LinearSolver< Number > &input_solver, SparseMatrix< Number > &input_matrix, NumericVector< Number > &input_rhs)

void	set_rb_evaluation (RBEvaluation &rb_eval_in)

RBEvaluation &	get_rb_evaluation ()

bool	is_rb_eval_initialized () const

RBThetaExpansion &	get_rb_theta_expansion ()

void	set_rb_assembly_expansion (RBAssemblyExpansion &rb_assembly_expansion_in)

RBAssemblyExpansion &	get_rb_assembly_expansion ()

virtual Real	compute_max_error_bound ()

const RBParameters &	get_greedy_parameter (unsigned int i)

void	set_rel_training_tolerance (Real new_training_tolerance)

Real	get_rel_training_tolerance ()

void	set_abs_training_tolerance (Real new_training_tolerance)

Real	get_abs_training_tolerance ()

void	set_normalize_rb_bound_in_greedy (bool normalize_rb_bound_in_greedy)

bool	get_normalize_rb_bound_in_greedy ()

unsigned int	get_Nmax () const

virtual void	set_Nmax (unsigned int Nmax)

void	set_quiet_mode (bool quiet_mode_in)

bool	is_quiet () const

virtual void	load_basis_function (unsigned int i)

SparseMatrix< Number > *	get_inner_product_matrix ()

SparseMatrix< Number > *	get_Aq (unsigned int q)

SparseMatrix< Number > *	get_non_dirichlet_Aq (unsigned int q)

NumericVector< Number > *	get_Fq (unsigned int q)

NumericVector< Number > *	get_non_dirichlet_Fq (unsigned int q)

NumericVector< Number > *	get_output_vector (unsigned int n, unsigned int q_l)

NumericVector< Number > *	get_non_dirichlet_output_vector (unsigned int n, unsigned int q_l)

virtual void	get_all_vectors (std::map< std::string, NumericVector< Number > * > &all_vectors)

virtual void	get_output_vectors (std::map< std::string, NumericVector< Number > * > &all_vectors)

void	assemble_inner_product_matrix (SparseMatrix< Number > *input_matrix, bool apply_dof_constraints=true)

void	assemble_Aq_matrix (unsigned int q, SparseMatrix< Number > *input_matrix, bool apply_dof_constraints=true)

void	assemble_Fq_vector (unsigned int q, NumericVector< Number > *input_vector, bool apply_dof_constraints=true)

void	add_scaled_Aq (Number scalar, unsigned int q_a, SparseMatrix< Number > *input_matrix, bool symmetrize)

virtual void	recompute_all_residual_terms (const bool compute_inner_products=true)

void	set_rb_construction_parameters (unsigned int n_training_samples_in, bool deterministic_training_in, unsigned int training_parameters_random_seed_in, bool quiet_mode_in, unsigned int Nmax_in, Real rel_training_tolerance_in, Real abs_training_tolerance_in, bool normalize_rb_error_bound_in_greedy_in, RBParameters mu_min_in, RBParameters mu_max_in, std::map< std::string, std::vector< Real > > discrete_parameter_values_in, std::map< std::string, bool > log_scaling)

void	print_basis_function_orthogonality ()

unsigned int	get_delta_N () const

void	set_inner_product_assembly (ElemAssembly &inner_product_assembly_in)

ElemAssembly &	get_inner_product_assembly ()

void	zero_constrained_dofs_on_vector (NumericVector< Number > &vector)

void	set_convergence_assertion_flag (bool flag)

numeric_index_type	get_n_training_samples () const

numeric_index_type	get_local_n_training_samples () const

numeric_index_type	get_first_local_training_index () const

numeric_index_type	get_last_local_training_index () const

virtual void	initialize_training_parameters (const RBParameters &mu_min, const RBParameters &mu_max, unsigned int n_training_parameters, std::map< std::string, bool > log_param_scale, bool deterministic=true)

virtual void	load_training_set (std::map< std::string, std::vector< Number > > &new_training_set)

void	broadcast_parameters (unsigned int proc_id)

void	set_training_random_seed (unsigned int seed)

void	set_deterministic_training_parameter_name (const std::string &name)

const std::string &	get_deterministic_training_parameter_name () const

void	set_deterministic_training_parameter_repeats (unsigned int repeats)

unsigned int	get_deterministic_training_parameter_repeats () const

virtual void	assemble () libmesh_override

virtual void	restrict_solve_to (const SystemSubset *subset, const SubsetSolveMode subset_solve_mode=SUBSET_ZERO) libmesh_override

virtual void	solve () libmesh_override

virtual LinearSolver< Number > *	get_linear_solver () const libmesh_override

virtual void	release_linear_solver (LinearSolver< Number > *) const libmesh_override

virtual void	assembly (bool get_residual, bool get_jacobian, bool apply_heterogeneous_constraints=false, bool apply_no_constraints=false) libmesh_override

unsigned int	n_linear_iterations () const

Real	final_linear_residual () const

void	attach_shell_matrix (ShellMatrix< Number > *shell_matrix)

void	detach_shell_matrix ()

ShellMatrix< Number > *	get_shell_matrix ()

virtual void	disable_cache () libmesh_override

virtual std::pair< unsigned int, Real >	get_linear_solve_parameters () const

virtual void	assemble_residual_derivatives (const ParameterVector &parameters) libmesh_override

virtual std::pair< unsigned int, Real >	sensitivity_solve (const ParameterVector &parameters) libmesh_override

virtual std::pair< unsigned int, Real >	weighted_sensitivity_solve (const ParameterVector &parameters, const ParameterVector &weights) libmesh_override

virtual std::pair< unsigned int, Real >	adjoint_solve (const QoISet &qoi_indices=QoISet()) libmesh_override

virtual std::pair< unsigned int, Real >	weighted_sensitivity_adjoint_solve (const ParameterVector &parameters, const ParameterVector &weights, const QoISet &qoi_indices=QoISet()) libmesh_override

virtual void	adjoint_qoi_parameter_sensitivity (const QoISet &qoi_indices, const ParameterVector &parameters, SensitivityData &sensitivities) libmesh_override

virtual void	forward_qoi_parameter_sensitivity (const QoISet &qoi_indices, const ParameterVector &parameters, SensitivityData &sensitivities) libmesh_override

virtual void	qoi_parameter_hessian (const QoISet &qoi_indices, const ParameterVector &parameters, SensitivityData &hessian) libmesh_override

virtual void	qoi_parameter_hessian_vector_product (const QoISet &qoi_indices, const ParameterVector &parameters, const ParameterVector &vector, SensitivityData &product) libmesh_override

SparseMatrix< Number > &	add_matrix (const std::string &mat_name)

bool	have_matrix (const std::string &mat_name) const

const SparseMatrix< Number > *	request_matrix (const std::string &mat_name) const

SparseMatrix< Number > *	request_matrix (const std::string &mat_name)

const SparseMatrix< Number > &	get_matrix (const std::string &mat_name) const

SparseMatrix< Number > &	get_matrix (const std::string &mat_name)

virtual unsigned int	n_matrices () const libmesh_override

virtual void	assemble_qoi (const QoISet &qoi_indices=QoISet()) libmesh_override

virtual void	assemble_qoi_derivative (const QoISet &qoi_indices=QoISet(), bool include_liftfunc=true, bool apply_constraints=true) libmesh_override

void	init ()

virtual void	reinit_constraints ()

bool	is_initialized ()

virtual void	update ()

bool	is_adjoint_already_solved () const

void	set_adjoint_already_solved (bool setting)

virtual void	qoi_parameter_sensitivity (const QoISet &qoi_indices, const ParameterVector &parameters, SensitivityData &sensitivities)

virtual bool	compare (const System &other_system, const Real threshold, const bool verbose) const

const std::string &	name () const

void	project_solution (FunctionBase< Number > f, FunctionBase< Gradient > g=libmesh_nullptr) const

void	project_solution (FEMFunctionBase< Number > f, FEMFunctionBase< Gradient > g=libmesh_nullptr) const

void	project_solution (Number fptr(const Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name), Gradient gptr(const Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name), const Parameters &parameters) const

void	project_vector (NumericVector< Number > &new_vector, FunctionBase< Number > f, FunctionBase< Gradient > g=libmesh_nullptr, int is_adjoint=-1) const

void	project_vector (NumericVector< Number > &new_vector, FEMFunctionBase< Number > f, FEMFunctionBase< Gradient > g=libmesh_nullptr, int is_adjoint=-1) const

void	project_vector (Number fptr(const Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name), Gradient gptr(const Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name), const Parameters &parameters, NumericVector< Number > &new_vector, int is_adjoint=-1) const

void	boundary_project_solution (const std::set< boundary_id_type > &b, const std::vector< unsigned int > &variables, FunctionBase< Number > f, FunctionBase< Gradient > g=libmesh_nullptr)

void	boundary_project_solution (const std::set< boundary_id_type > &b, const std::vector< unsigned int > &variables, Number fptr(const Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name), Gradient gptr(const Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name), const Parameters &parameters)

void	boundary_project_vector (const std::set< boundary_id_type > &b, const std::vector< unsigned int > &variables, NumericVector< Number > &new_vector, FunctionBase< Number > f, FunctionBase< Gradient > g=libmesh_nullptr, int is_adjoint=-1) const

void	boundary_project_vector (const std::set< boundary_id_type > &b, const std::vector< unsigned int > &variables, Number fptr(const Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name), Gradient gptr(const Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name), const Parameters &parameters, NumericVector< Number > &new_vector, int is_adjoint=-1) const

unsigned int	number () const

void	update_global_solution (std::vector< Number > &global_soln) const

void	update_global_solution (std::vector< Number > &global_soln, const processor_id_type dest_proc) const

const MeshBase &	get_mesh () const

MeshBase &	get_mesh ()

const DofMap &	get_dof_map () const

DofMap &	get_dof_map ()

const EquationSystems &	get_equation_systems () const

EquationSystems &	get_equation_systems ()

bool	active () const

void	activate ()

void	deactivate ()

void	set_basic_system_only ()

vectors_iterator	vectors_begin ()

const_vectors_iterator	vectors_begin () const

vectors_iterator	vectors_end ()

const_vectors_iterator	vectors_end () const

NumericVector< Number > &	add_vector (const std::string &vec_name, const bool projections=true, const ParallelType type=PARALLEL)

void	remove_vector (const std::string &vec_name)

bool &	project_solution_on_reinit (void)

bool	have_vector (const std::string &vec_name) const

const NumericVector< Number > *	request_vector (const std::string &vec_name) const

NumericVector< Number > *	request_vector (const std::string &vec_name)

const NumericVector< Number > *	request_vector (const unsigned int vec_num) const

NumericVector< Number > *	request_vector (const unsigned int vec_num)

const NumericVector< Number > &	get_vector (const std::string &vec_name) const

NumericVector< Number > &	get_vector (const std::string &vec_name)

const NumericVector< Number > &	get_vector (const unsigned int vec_num) const

NumericVector< Number > &	get_vector (const unsigned int vec_num)

const std::string &	vector_name (const unsigned int vec_num) const

const std::string &	vector_name (const NumericVector< Number > &vec_reference) const

void	set_vector_as_adjoint (const std::string &vec_name, int qoi_num)

int	vector_is_adjoint (const std::string &vec_name) const

void	set_vector_preservation (const std::string &vec_name, bool preserve)

bool	vector_preservation (const std::string &vec_name) const

NumericVector< Number > &	add_adjoint_solution (unsigned int i=0)

NumericVector< Number > &	get_adjoint_solution (unsigned int i=0)

const NumericVector< Number > &	get_adjoint_solution (unsigned int i=0) const

NumericVector< Number > &	add_sensitivity_solution (unsigned int i=0)

NumericVector< Number > &	get_sensitivity_solution (unsigned int i=0)

const NumericVector< Number > &	get_sensitivity_solution (unsigned int i=0) const

NumericVector< Number > &	add_weighted_sensitivity_adjoint_solution (unsigned int i=0)

NumericVector< Number > &	get_weighted_sensitivity_adjoint_solution (unsigned int i=0)

const NumericVector< Number > &	get_weighted_sensitivity_adjoint_solution (unsigned int i=0) const

NumericVector< Number > &	add_weighted_sensitivity_solution ()

NumericVector< Number > &	get_weighted_sensitivity_solution ()

const NumericVector< Number > &	get_weighted_sensitivity_solution () const

NumericVector< Number > &	add_adjoint_rhs (unsigned int i=0)

NumericVector< Number > &	get_adjoint_rhs (unsigned int i=0)

const NumericVector< Number > &	get_adjoint_rhs (unsigned int i=0) const

NumericVector< Number > &	add_sensitivity_rhs (unsigned int i=0)

NumericVector< Number > &	get_sensitivity_rhs (unsigned int i=0)

const NumericVector< Number > &	get_sensitivity_rhs (unsigned int i=0) const

unsigned int	n_vectors () const

unsigned int	n_vars () const

unsigned int	n_variable_groups () const

unsigned int	n_components () const

dof_id_type	n_dofs () const

dof_id_type	n_active_dofs () const

dof_id_type	n_constrained_dofs () const

dof_id_type	n_local_constrained_dofs () const

dof_id_type	n_local_dofs () const

unsigned int	add_variable (const std::string &var, const FEType &type, const std::set< subdomain_id_type > *const active_subdomains=libmesh_nullptr)

unsigned int	add_variable (const std::string &var, const Order order=FIRST, const FEFamily=LAGRANGE, const std::set< subdomain_id_type > *const active_subdomains=libmesh_nullptr)

unsigned int	add_variables (const std::vector< std::string > &vars, const FEType &type, const std::set< subdomain_id_type > *const active_subdomains=libmesh_nullptr)

unsigned int	add_variables (const std::vector< std::string > &vars, const Order order=FIRST, const FEFamily=LAGRANGE, const std::set< subdomain_id_type > *const active_subdomains=libmesh_nullptr)

const Variable &	variable (unsigned int var) const

const VariableGroup &	variable_group (unsigned int vg) const

bool	has_variable (const std::string &var) const

const std::string &	variable_name (const unsigned int i) const

unsigned short int	variable_number (const std::string &var) const

void	get_all_variable_numbers (std::vector< unsigned int > &all_variable_numbers) const

unsigned int	variable_scalar_number (const std::string &var, unsigned int component) const

unsigned int	variable_scalar_number (unsigned int var_num, unsigned int component) const

const FEType &	variable_type (const unsigned int i) const

const FEType &	variable_type (const std::string &var) const

bool	identify_variable_groups () const

void	identify_variable_groups (const bool)

Real	calculate_norm (const NumericVector< Number > &v, unsigned int var, FEMNormType norm_type, std::set< unsigned int > *skip_dimensions=libmesh_nullptr) const

Real	calculate_norm (const NumericVector< Number > &v, const SystemNorm &norm, std::set< unsigned int > *skip_dimensions=libmesh_nullptr) const

void	read_header (Xdr &io, const std::string &version, const bool read_header=true, const bool read_additional_data=true, const bool read_legacy_format=false)

void	read_legacy_data (Xdr &io, const bool read_additional_data=true)

template<typename ValType >
void	read_serialized_data (Xdr &io, const bool read_additional_data=true)

void	read_serialized_data (Xdr &io, const bool read_additional_data=true)

template<typename InValType >
std::size_t	read_serialized_vectors (Xdr &io, const std::vector< NumericVector< Number > * > &vectors) const

std::size_t	read_serialized_vectors (Xdr &io, const std::vector< NumericVector< Number > * > &vectors) const

template<typename InValType >
void	read_parallel_data (Xdr &io, const bool read_additional_data)

void	read_parallel_data (Xdr &io, const bool read_additional_data)

void	write_header (Xdr &io, const std::string &version, const bool write_additional_data) const

void	write_serialized_data (Xdr &io, const bool write_additional_data=true) const

std::size_t	write_serialized_vectors (Xdr &io, const std::vector< const NumericVector< Number > * > &vectors) const

void	write_parallel_data (Xdr &io, const bool write_additional_data) const

std::string	get_info () const

void	attach_init_function (void fptr(EquationSystems &es, const std::string &name))

void	attach_init_object (Initialization &init)

void	attach_assemble_function (void fptr(EquationSystems &es, const std::string &name))

void	attach_assemble_object (Assembly &assemble)

void	attach_constraint_function (void fptr(EquationSystems &es, const std::string &name))

void	attach_constraint_object (Constraint &constrain)

void	attach_QOI_function (void fptr(EquationSystems &es, const std::string &name, const QoISet &qoi_indices))

void	attach_QOI_object (QOI &qoi)

void	attach_QOI_derivative (void fptr(EquationSystems &es, const std::string &name, const QoISet &qoi_indices, bool include_liftfunc, bool apply_constraints))

void	attach_QOI_derivative_object (QOIDerivative &qoi_derivative)

virtual void	user_initialization ()

virtual void	user_assembly ()

virtual void	user_constrain ()

virtual void	user_QOI (const QoISet &qoi_indices)

virtual void	user_QOI_derivative (const QoISet &qoi_indices=QoISet(), bool include_liftfunc=true, bool apply_constraints=true)

virtual void	restrict_vectors ()

virtual void	prolong_vectors ()

Number	current_solution (const dof_id_type global_dof_number) const

Number	point_value (unsigned int var, const Point &p, const bool insist_on_success=true) const

Number	point_value (unsigned int var, const Point &p, const Elem &e) const

Number	point_value (unsigned int var, const Point &p, const Elem *e) const

Gradient	point_gradient (unsigned int var, const Point &p, const bool insist_on_success=true) const

Gradient	point_gradient (unsigned int var, const Point &p, const Elem &e) const

Gradient	point_gradient (unsigned int var, const Point &p, const Elem *e) const

Tensor	point_hessian (unsigned int var, const Point &p, const bool insist_on_success=true) const

Tensor	point_hessian (unsigned int var, const Point &p, const Elem &e) const

Tensor	point_hessian (unsigned int var, const Point &p, const Elem *e) const

void	local_dof_indices (const unsigned int var, std::set< dof_id_type > &var_indices) const

void	zero_variable (NumericVector< Number > &v, unsigned int var_num) const

bool &	hide_output ()

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

void	initialize_parameters (const RBParameters &mu_min_in, const RBParameters &mu_max_in, const std::map< std::string, std::vector< Real > > &discrete_parameter_values)

void	initialize_parameters (const RBParametrized &rb_parametrized)

unsigned int	get_n_params () const

unsigned int	get_n_continuous_params () const

unsigned int	get_n_discrete_params () const

std::set< std::string >	get_parameter_names () const

const RBParameters &	get_parameters () const

void	set_parameters (const RBParameters &params)

const RBParameters &	get_parameters_min () const

const RBParameters &	get_parameters_max () const

Real	get_parameter_min (const std::string &param_name) const

Real	get_parameter_max (const std::string &param_name) const

void	print_parameters () const

void	write_parameter_data_to_files (const std::string &continuous_param_file_name, const std::string &discrete_param_file_name, const bool write_binary_data)

void	read_parameter_data_from_files (const std::string &continuous_param_file_name, const std::string &discrete_param_file_name, const bool read_binary_data)

bool	is_discrete_parameter (const std::string &mu_name) const

const std::map< std::string, std::vector< Real > > &	get_discrete_parameter_values () const

void	print_discrete_parameter_values () const

Real	get_delta_t () const

void	set_delta_t (const Real delta_t_in)

Real	get_euler_theta () const

void	set_euler_theta (const Real euler_theta_in)

unsigned int	get_time_step () const

void	set_time_step (const unsigned int k)

unsigned int	get_n_time_steps () const

void	set_n_time_steps (const unsigned int K)

Real	get_control (const unsigned int k) const

void	set_control (const std::vector< Real > &control)

void	process_temporal_parameters_file (const std::string &parameters_filename)

void	pull_temporal_discretization_data (RBTemporalDiscretization &other)

Static Public Member Functions
static void	print_info (std::ostream &out=libMesh::out)

static void	print_info (std::ostream &out=libMesh::out)

static UniquePtr< DirichletBoundary >	build_zero_dirichlet_boundary_object ()

static std::string	get_info ()

static std::string	get_info ()

static unsigned int	n_objects ()

static unsigned int	n_objects ()

static void	enable_print_counter_info ()

static void	enable_print_counter_info ()

static void	disable_print_counter_info ()

static void	disable_print_counter_info ()

static Real	get_closest_value (Real value, const std::vector< Real > &list_of_values)

Public Attributes
UniquePtr< SparseMatrix< Number > >	L2_matrix

UniquePtr< SparseMatrix< Number > >	non_dirichlet_L2_matrix

std::vector< SparseMatrix< Number > * >	M_q_vector

std::vector< SparseMatrix< Number > * >	non_dirichlet_M_q_vector

std::vector< std::vector< Number > >	truth_outputs_all_k

bool	nonzero_initialization

bool	compute_truth_projection_error

std::string	init_filename

UniquePtr< NumericVector< Number > >	old_local_solution

UniquePtr< NumericVector< Number > >	older_local_solution

std::vector< Real >	training_error_bounds

UniquePtr< LinearSolver< Number > >	inner_product_solver

LinearSolver< Number > *	extra_linear_solver

UniquePtr< SparseMatrix< Number > >	inner_product_matrix

std::vector< Number >	truth_outputs

std::vector< std::vector< Number > >	output_dual_innerprods

std::vector< NumericVector< Number > * >	Fq_representor

std::vector< Number >	Fq_representor_innerprods

bool	exit_on_repeated_greedy_parameters

bool	impose_internal_fluxes

bool	compute_RB_inner_product

bool	store_non_dirichlet_operators

bool	use_empty_rb_solve_in_greedy

bool	Fq_representor_innerprods_computed

UniquePtr< LinearSolver< Number > >	linear_solver

SparseMatrix< Number > *	matrix

bool	zero_out_matrix_and_rhs

NumericVector< Number > *	rhs

bool	assemble_before_solve

bool	use_fixed_solution

int	extra_quadrature_order

UniquePtr< NumericVector< Number > >	solution

UniquePtr< NumericVector< Number > >	current_local_solution

Real	time

std::vector< Number >	qoi

bool	verbose_mode

Protected Types
typedef std::map< std::string, std::pair< unsigned int, unsigned int > >	Counts

typedef std::map< std::string, std::pair< unsigned int, unsigned int > >	Counts

Protected Member Functions
virtual void	allocate_data_structures () libmesh_override

virtual void	assemble_affine_expansion (bool skip_matrix_assembly, bool skip_vector_assembly) libmesh_override

virtual void	initialize_truth ()

virtual SparseMatrix< Number > &	get_matrix_for_output_dual_solves () libmesh_override

void	add_IC_to_RB_space ()

virtual void	enrich_RB_space () libmesh_override

virtual void	update_system () libmesh_override

virtual void	update_RB_system_matrices () libmesh_override

virtual void	update_residual_terms (bool compute_inner_products) libmesh_override

Number	set_error_temporal_data ()

virtual void	init_data () libmesh_override

virtual void	re_update () libmesh_override

virtual UniquePtr< DGFEMContext >	build_context ()

void	update_greedy_param_list ()

void	add_scaled_matrix_and_vector (Number scalar, ElemAssembly elem_assembly, SparseMatrix< Number > input_matrix, NumericVector< Number > *input_vector, bool symmetrize=false, bool apply_dof_constraints=true)

virtual void	set_context_solution_vec (NumericVector< Number > &vec)

virtual void	assemble_all_affine_vectors ()

virtual void	assemble_all_output_vectors ()

virtual void	compute_output_dual_innerprods ()

virtual void	compute_Fq_representor_innerprods (bool compute_inner_products=true)

virtual Real	get_RB_error_bound ()

virtual void	init_context (FEMContext &)

bool	get_convergence_assertion_flag () const

void	check_convergence (LinearSolver< Number > &input_solver)

RBParameters	get_params_from_training_set (unsigned int index)

void	set_params_from_training_set (unsigned int index)

virtual void	set_params_from_training_set_and_broadcast (unsigned int index)

virtual void	init_matrices ()

void	project_vector (NumericVector< Number > &, int is_adjoint=-1) const

void	project_vector (const NumericVector< Number > &, NumericVector< Number > &, int is_adjoint=-1) const

void	increment_constructor_count (const std::string &name)

void	increment_constructor_count (const std::string &name)

void	increment_destructor_count (const std::string &name)

void	increment_destructor_count (const std::string &name)

Static Protected Member Functions
static void	get_global_max_error_pair (const Parallel::Communicator &communicator, std::pair< numeric_index_type, Real > &error_pair)

static void	generate_training_parameters_random (const Parallel::Communicator &communicator, std::map< std::string, bool > log_param_scale, std::map< std::string, NumericVector< Number > * > &training_parameters_in, unsigned int n_training_samples_in, const RBParameters &min_parameters, const RBParameters &max_parameters, int training_parameters_random_seed=-1, bool serial_training_set=false)

static void	generate_training_parameters_deterministic (const Parallel::Communicator &communicator, std::map< std::string, bool > log_param_scale, std::map< std::string, NumericVector< Number > * > &training_parameters_in, unsigned int n_training_samples_in, const RBParameters &min_parameters, const RBParameters &max_parameters, bool serial_training_set=false)

Protected Attributes
Real	POD_tol

int	max_truth_solves

ElemAssembly *	L2_assembly

DenseVector< Number >	RB_ic_proj_rhs_all_N

unsigned int	Nmax

unsigned int	delta_N

bool	quiet_mode

bool	output_dual_innerprods_computed

bool	assert_convergence

bool	serial_training_set

UniquePtr< NumericVector< Number > >	inner_product_storage_vector

unsigned int	_n_linear_iterations

Real	_final_linear_residual

ShellMatrix< Number > *	_shell_matrix

const SystemSubset *	_subset

SubsetSolveMode	_subset_solve_mode

const Parallel::Communicator &	_communicator

Static Protected Attributes
static Counts	_counts

static Counts	_counts

static Threads::atomic< unsigned int >	_n_objects

static Threads::atomic< unsigned int >	_n_objects

static Threads::spin_mutex	_mutex

static Threads::spin_mutex	_mutex

static bool	_enable_print_counter = true

static bool	_enable_print_counter = true

Private Attributes
std::vector< NumericVector< Number > * >	temporal_data

Detailed Description

This class is part of the rbOOmit framework.

TransientRBConstruction extends RBConstruction to add functionality relevant in the time-dependent case.

We can handle time controls on the RHS as h(t)*f(x, $\mu$ ). See Martin Grepl's thesis for more details.

Author: David J. Knezevic

Date: 2009

Definition at line 48 of file transient_rb_construction.h.

Member Typedef Documentation

typedef std::map<std::string, SparseMatrix<Number> *>::const_iterator libMesh::ImplicitSystem::const_matrices_iterator

inherited

Definition at line 283 of file implicit_system.h.

typedef std::map<std::string, NumericVector<Number> *>::const_iterator libMesh::System::const_vectors_iterator

inherited

Definition at line 749 of file system.h.

typedef std::map<std::string, std::pair<unsigned int, unsigned int> > libMesh::ReferenceCounter::Counts

protectedinherited

Data structure to log the information. The log is identified by the class name.

Definition at line 119 of file reference_counter.h.

typedef std::map<std::string, std::pair<unsigned int, unsigned int> > libMesh::ReferenceCounter::Counts

protectedinherited

Data structure to log the information. The log is identified by the class name.

Definition at line 119 of file reference_counter.h.

typedef std::map<std::string, SparseMatrix<Number> *>::iterator libMesh::ImplicitSystem::matrices_iterator

inherited

Matrix iterator typedefs.

Definition at line 282 of file implicit_system.h.

typedef TransientSystem<RBConstruction> libMesh::TransientRBConstruction::Parent

The type of the parent.

Definition at line 73 of file transient_rb_construction.h.

typedef TransientRBConstruction libMesh::TransientRBConstruction::sys_type

The type of system.

Definition at line 68 of file transient_rb_construction.h.

typedef std::map<std::string, NumericVector<Number> *>::iterator libMesh::System::vectors_iterator

inherited

Vector iterator typedefs.

Definition at line 748 of file system.h.

Constructor & Destructor Documentation

libMesh::TransientRBConstruction::TransientRBConstruction	(	EquationSystems &	es,
		const std::string &	name,
		const unsigned int	number
	)

Constructor. Optionally initializes required data structures.

virtual libMesh::TransientRBConstruction::~TransientRBConstruction ( )

virtual

Destructor.

Member Function Documentation

void libMesh::System::activate ( )

inlineinherited

Activates the system. Only active systems are solved.

Definition at line 2052 of file system.h.

References libMesh::System::_active.

Referenced by libMesh::System::get_equation_systems().

 {
   _active = true;
 }

bool libMesh::System::active ( ) const

inlineinherited

Returns: true if the system is active, false otherwise. An active system will be solved.

Definition at line 2044 of file system.h.

References libMesh::System::_active.

Referenced by libMesh::System::get_equation_systems().

 {
   return _active;
 }

NumericVector< Number > & libMesh::System::add_adjoint_rhs ( unsigned int i = 0 )

inherited

Returns: A reference to one of the system's adjoint rhs vectors, by default the one corresponding to the first qoi. Creates the vector if it doesn't already exist.

Definition at line 1035 of file system.C.

References libMesh::System::add_vector().

Referenced by libMesh::ExplicitSystem::assemble_qoi_derivative(), libMesh::FEMSystem::assemble_qoi_derivative(), and libMesh::System::project_solution_on_reinit().

 {
   std::ostringstream adjoint_rhs_name;
   adjoint_rhs_name << "adjoint_rhs" << i;
 
   return this->add_vector(adjoint_rhs_name.str(), false);
 }

NumericVector< Number > & libMesh::System::add_adjoint_solution ( unsigned int i = 0 )

inherited

Returns: A reference to one of the system's adjoint solution vectors, by default the one corresponding to the first qoi. Creates the vector if it doesn't already exist.

Definition at line 971 of file system.C.

References libMesh::System::add_vector(), and libMesh::System::set_vector_as_adjoint().

Referenced by libMesh::ImplicitSystem::adjoint_solve(), and libMesh::System::project_solution_on_reinit().

 {
   std::ostringstream adjoint_name;
   adjoint_name << "adjoint_solution" << i;
 
   NumericVector<Number> & returnval = this->add_vector(adjoint_name.str());
   this->set_vector_as_adjoint(adjoint_name.str(), i);
   return returnval;
 }

void libMesh::TransientRBConstruction::add_IC_to_RB_space ( )

protected

Initialize RB space by adding the truth initial condition as the first RB basis function.

SparseMatrix< Number > & libMesh::ImplicitSystem::add_matrix ( const std::string & mat_name )

inherited

Adds the additional matrix mat_name to this system. Only allowed prior to assemble(). All additional matrices have the same sparsity pattern as the matrix used during solution. When not System but the user wants to initialize the mayor matrix, then all the additional matrices, if existent, have to be initialized by the user, too.

Definition at line 209 of file implicit_system.C.

References libMesh::ImplicitSystem::_can_add_matrices, libMesh::ImplicitSystem::_matrices, libMesh::SparseMatrix< T >::build(), libMesh::ParallelObject::comm(), and libMesh::ImplicitSystem::have_matrix().

Referenced by libMesh::ImplicitSystem::add_system_matrix(), libMesh::EigenTimeSolver::init(), and libMesh::NewmarkSystem::NewmarkSystem().

 {
   // only add matrices before initializing...
   if (!_can_add_matrices)
     libmesh_error_msg("ERROR: Too late.  Cannot add matrices to the system after initialization"
                       << "\n any more.  You should have done this earlier.");
 
   // Return the matrix if it is already there.
   if (this->have_matrix(mat_name))
     return *(_matrices[mat_name]);
 
   // Otherwise build the matrix and return it.
   SparseMatrix<Number> * buf = SparseMatrix<Number>::build(this->comm()).release();
   _matrices.insert (std::make_pair (mat_name, buf));
 
   return *buf;
 }

void libMesh::RBConstruction::add_scaled_Aq	(	Number	scalar,
		unsigned int	q_a,
		SparseMatrix< Number > *	input_matrix,
		bool	symmetrize
	)

inherited

Add the scaled q^th affine matrix to input_matrix. If symmetrize==true, then we symmetrize Aq before adding it.

Referenced by libMesh::RBConstruction::is_quiet().

void libMesh::TransientRBConstruction::add_scaled_mass_matrix	(	Number	scalar,
		SparseMatrix< Number > *	input_matrix
	)

Add the scaled mass matrix (assembled for the current parameter) to input_matrix.

void libMesh::RBConstruction::add_scaled_matrix_and_vector	(	Number	scalar,
		ElemAssembly *	elem_assembly,
		SparseMatrix< Number > *	input_matrix,
		NumericVector< Number > *	input_vector,
		bool	symmetrize = `false`,
		bool	apply_dof_constraints = `true`
	)

protectedinherited

This function loops over the mesh and applies the specified interior and/or boundary assembly routines, then adds the scaled result to input_matrix and/or input_vector. If symmetrize==true then we assemble the symmetric part of the matrix, 0.5*(A + A^T)

NumericVector< Number > & libMesh::System::add_sensitivity_rhs ( unsigned int i = 0 )

inherited

Returns: A reference to one of the system's sensitivity rhs vectors, by default the one corresponding to the first parameter. Creates the vector if it doesn't already exist.

Definition at line 1065 of file system.C.

References libMesh::System::add_vector().

Referenced by libMesh::ImplicitSystem::assemble_residual_derivatives(), and libMesh::System::project_solution_on_reinit().

 {
   std::ostringstream sensitivity_rhs_name;
   sensitivity_rhs_name << "sensitivity_rhs" << i;
 
   return this->add_vector(sensitivity_rhs_name.str(), false);
 }

NumericVector< Number > & libMesh::System::add_sensitivity_solution ( unsigned int i = 0 )

inherited

Returns: A reference to one of the system's solution sensitivity vectors, by default the one corresponding to the first parameter. Creates the vector if it doesn't already exist.

Definition at line 920 of file system.C.

References libMesh::System::add_vector().

Referenced by libMesh::System::project_solution_on_reinit(), and libMesh::ImplicitSystem::sensitivity_solve().

 {
   std::ostringstream sensitivity_name;
   sensitivity_name << "sensitivity_solution" << i;
 
   return this->add_vector(sensitivity_name.str());
 }

unsigned int libMesh::System::add_variable	(	const std::string &	var,
		const FEType &	type,
		const std::set< subdomain_id_type > *const	active_subdomains = `libmesh_nullptr`
	)

inherited

Adds the variable var to the list of variables for this system.

Returns: The index number for the new variable.

Definition at line 1095 of file system.C.

References libMesh::System::_variable_groups, libMesh::System::_variable_numbers, libMesh::System::_variables, libMesh::System::add_variables(), libMesh::VariableGroup::append(), libMesh::System::identify_variable_groups(), libMesh::System::is_initialized(), libMesh::libmesh_assert(), libmesh_nullptr, libMesh::System::n_variable_groups(), libMesh::System::n_vars(), libMesh::System::number(), libMesh::System::variable_name(), and libMesh::System::variable_type().

Referenced by libMesh::DifferentiableSystem::add_second_order_dot_vars(), libMesh::System::add_variable(), libMesh::ErrorVector::plot_error(), and libMesh::System::project_solution_on_reinit().

 {
   libmesh_assert(!this->is_initialized());
 
   // Make sure the variable isn't there already
   // or if it is, that it's the type we want
   for (unsigned int v=0; v<this->n_vars(); v++)
     if (this->variable_name(v) == var)
       {
         if (this->variable_type(v) == type)
           return _variables[v].number();
 
         libmesh_error_msg("ERROR: incompatible variable " << var << " has already been added for this system!");
       }
 
   // Optimize for VariableGroups here - if the user is adding multiple
   // variables of the same FEType and subdomain restriction, catch
   // that here and add them as members of the same VariableGroup.
   //
   // start by setting this flag to whatever the user has requested
   // and then consider the conditions which should negate it.
   bool should_be_in_vg = this->identify_variable_groups();
 
   // No variable groups, nothing to add to
   if (!this->n_variable_groups())
     should_be_in_vg = false;
 
   else
     {
       VariableGroup & vg(_variable_groups.back());
 
       // get a pointer to their subdomain restriction, if any.
       const std::set<subdomain_id_type> * const
         their_active_subdomains (vg.implicitly_active() ?
                                  libmesh_nullptr : &vg.active_subdomains());
 
       // Different types?
       if (vg.type() != type)
         should_be_in_vg = false;
 
       // they are restricted, we aren't?
       if (their_active_subdomains && !active_subdomains)
         should_be_in_vg = false;
 
       // they aren't restriced, we are?
       if (!their_active_subdomains && active_subdomains)
         should_be_in_vg = false;
 
       if (their_active_subdomains && active_subdomains)
         // restricted to different sets?
         if (*their_active_subdomains != *active_subdomains)
           should_be_in_vg = false;
 
       // OK, after all that, append the variable to the vg if none of the conditions
       // were violated
       if (should_be_in_vg)
         {
           const unsigned short curr_n_vars = cast_int<unsigned short>
             (this->n_vars());
 
           vg.append (var);
 
           _variables.push_back(vg(vg.n_variables()-1));
           _variable_numbers[var] = curr_n_vars;
           return curr_n_vars;
         }
     }
 
   // otherwise, fall back to adding a single variable group
   return this->add_variables (std::vector<std::string>(1, var),
                               type,
                               active_subdomains);
 }

unsigned int libMesh::System::add_variable	(	const std::string &	var,
		const Order	order = `FIRST`,
		const FEFamily	family = `LAGRANGE`,
		const std::set< subdomain_id_type > *const	active_subdomains = `libmesh_nullptr`
	)

inherited

Adds the variable var to the list of variables for this system. Same as before, but assumes LAGRANGE as default value for FEType.family.

Definition at line 1173 of file system.C.

References libMesh::System::add_variable().

 {
   return this->add_variable(var,
                             FEType(order, family),
                             active_subdomains);
 }

unsigned int libMesh::System::add_variables	(	const std::vector< std::string > &	vars,
		const FEType &	type,
		const std::set< subdomain_id_type > *const	active_subdomains = `libmesh_nullptr`
	)

inherited

Adds the variable var to the list of variables for this system.

Returns: The index number for the new variable.

Definition at line 1185 of file system.C.

References libMesh::System::_variable_groups, libMesh::System::_variable_numbers, libMesh::System::_variables, libMesh::System::is_initialized(), libMesh::libmesh_assert(), libmesh_nullptr, libMesh::System::n_components(), libMesh::System::n_vars(), libMesh::System::number(), libMesh::System::variable_name(), and libMesh::System::variable_type().

Referenced by libMesh::System::add_variable(), libMesh::System::add_variables(), and libMesh::System::project_solution_on_reinit().

 {
   libmesh_assert(!this->is_initialized());
 
   // Make sure the variable isn't there already
   // or if it is, that it's the type we want
   for (std::size_t ov=0; ov<vars.size(); ov++)
     for (unsigned int v=0; v<this->n_vars(); v++)
       if (this->variable_name(v) == vars[ov])
         {
           if (this->variable_type(v) == type)
             return _variables[v].number();
 
           libmesh_error_msg("ERROR: incompatible variable " << vars[ov] << " has already been added for this system!");
         }
 
   const unsigned short curr_n_vars = cast_int<unsigned short>
     (this->n_vars());
 
   const unsigned int next_first_component = this->n_components();
 
   // Add the variable group to the list
   _variable_groups.push_back((active_subdomains == libmesh_nullptr) ?
                              VariableGroup(this, vars, curr_n_vars,
                                            next_first_component, type) :
                              VariableGroup(this, vars, curr_n_vars,
                                            next_first_component, type, *active_subdomains));
 
   const VariableGroup & vg (_variable_groups.back());
 
   // Add each component of the group individually
   for (std::size_t v=0; v<vars.size(); v++)
     {
       _variables.push_back (vg(v));
       _variable_numbers[vars[v]] = cast_int<unsigned short>
         (curr_n_vars+v);
     }
 
   libmesh_assert_equal_to ((curr_n_vars+vars.size()), this->n_vars());
 
   // BSK - Defer this now to System::init_data() so we can detect
   // VariableGroups 12/28/2012
   // // Add the variable group to the _dof_map
   // _dof_map->add_variable_group (vg);
 
   // Return the number of the new variable
   return cast_int<unsigned int>(curr_n_vars+vars.size()-1);
 }

unsigned int libMesh::System::add_variables	(	const std::vector< std::string > &	vars,
		const Order	order = `FIRST`,
		const FEFamily	family = `LAGRANGE`,
		const std::set< subdomain_id_type > *const	active_subdomains = `libmesh_nullptr`
	)

inherited

Adds the variable var to the list of variables for this system. Same as before, but assumes LAGRANGE as default value for FEType.family.

Definition at line 1238 of file system.C.

References libMesh::System::add_variables().

 {
   return this->add_variables(vars,
                              FEType(order, family),
                              active_subdomains);
 }

NumericVector< Number > & libMesh::System::add_vector	(	const std::string &	vec_name,
		const bool	projections = `true`,
		const ParallelType	type = `PARALLEL`
	)

inherited

Adds the additional vector vec_name to this system. All the additional vectors are similarly distributed, like the solution, and inititialized to zero.

By default vectors added by add_vector are projected to changed grids by reinit(). To zero them instead (more efficient), pass "false" as the second argument

Definition at line 676 of file system.C.

References libMesh::System::_dof_map, libMesh::System::_is_initialized, libMesh::System::_vector_is_adjoint, libMesh::System::_vector_projections, libMesh::System::_vector_types, libMesh::System::_vectors, libMesh::NumericVector< T >::build(), libMesh::ParallelObject::comm(), libMesh::GHOSTED, libMesh::System::have_vector(), libMesh::NumericVector< T >::init(), libMesh::System::n_dofs(), and libMesh::System::n_local_dofs().

Referenced by libMesh::System::add_adjoint_rhs(), libMesh::System::add_adjoint_solution(), libMesh::System::add_sensitivity_rhs(), libMesh::System::add_sensitivity_solution(), libMesh::ExplicitSystem::add_system_rhs(), libMesh::System::add_weighted_sensitivity_adjoint_solution(), libMesh::System::add_weighted_sensitivity_solution(), libMesh::AdjointRefinementEstimator::estimate_error(), libMesh::SecondOrderUnsteadySolver::init(), libMesh::UnsteadySolver::init(), libMesh::OptimizationSystem::init_data(), libMesh::ContinuationSystem::init_data(), libMesh::NewmarkSystem::NewmarkSystem(), libMesh::FrequencySystem::set_frequencies(), libMesh::FrequencySystem::set_frequencies_by_range(), and libMesh::FrequencySystem::set_frequencies_by_steps().

 {
   // Return the vector if it is already there.
   if (this->have_vector(vec_name))
     return *(_vectors[vec_name]);
 
   // Otherwise build the vector
   NumericVector<Number> * buf = NumericVector<Number>::build(this->comm()).release();
   _vectors.insert (std::make_pair (vec_name, buf));
   _vector_projections.insert (std::make_pair (vec_name, projections));
 
   _vector_types.insert (std::make_pair (vec_name, type));
 
   // Vectors are primal by default
   _vector_is_adjoint.insert (std::make_pair (vec_name, -1));
 
   // Initialize it if necessary
   if (_is_initialized)
     {
       if (type == GHOSTED)
         {
 #ifdef LIBMESH_ENABLE_GHOSTED
           buf->init (this->n_dofs(), this->n_local_dofs(),
                      _dof_map->get_send_list(), false,
                      GHOSTED);
 #else
           libmesh_error_msg("Cannot initialize ghosted vectors when they are not enabled.");
 #endif
         }
       else
         buf->init (this->n_dofs(), this->n_local_dofs(), false, type);
     }
 
   return *buf;
 }

NumericVector< Number > & libMesh::System::add_weighted_sensitivity_adjoint_solution ( unsigned int i = 0 )

inherited

Returns: A reference to one of the system's weighted sensitivity adjoint solution vectors, by default the one corresponding to the first qoi. Creates the vector if it doesn't already exist.

Definition at line 1003 of file system.C.

References libMesh::System::add_vector(), and libMesh::System::set_vector_as_adjoint().

Referenced by libMesh::System::project_solution_on_reinit(), and libMesh::ImplicitSystem::weighted_sensitivity_adjoint_solve().

 {
   std::ostringstream adjoint_name;
   adjoint_name << "weighted_sensitivity_adjoint_solution" << i;
 
   NumericVector<Number> & returnval = this->add_vector(adjoint_name.str());
   this->set_vector_as_adjoint(adjoint_name.str(), i);
   return returnval;
 }

NumericVector< Number > & libMesh::System::add_weighted_sensitivity_solution ( )

inherited

Returns: A reference to the solution of the last weighted sensitivity solve Creates the vector if it doesn't already exist.

Definition at line 950 of file system.C.

References libMesh::System::add_vector().

Referenced by libMesh::System::project_solution_on_reinit(), and libMesh::ImplicitSystem::weighted_sensitivity_solve().

 {
   return this->add_vector("weighted_sensitivity_solution");
 }

void libMesh::ImplicitSystem::adjoint_qoi_parameter_sensitivity	(	const QoISet &	qoi_indices,
		const ParameterVector &	parameters,
		SensitivityData &	sensitivities
	)

virtualinherited

Solves for the derivative of each of the system's quantities of interest q in qoi[qoi_indices] with respect to each parameter in parameters, placing the result for qoi i and parameter j into sensitivities[i][j].

Uses adjoint_solve() and the adjoint sensitivity method.

Currently uses finite differenced derivatives (partial q / partial p) and (partial R / partial p).

Reimplemented from libMesh::System.

Definition at line 693 of file implicit_system.C.

References std::abs(), libMesh::ImplicitSystem::adjoint_solve(), libMesh::SensitivityData::allocate_data(), libMesh::ExplicitSystem::assemble_qoi(), libMesh::ImplicitSystem::assemble_residual_derivatives(), libMesh::NumericVector< T >::dot(), libMesh::System::get_sensitivity_rhs(), libMesh::QoISet::has_index(), libMesh::System::is_adjoint_already_solved(), std::max(), libMesh::System::qoi, libMesh::Real, libMesh::ParameterVector::size(), and libMesh::TOLERANCE.

Referenced by libMesh::ImplicitSystem::assembly().

 {
   ParameterVector & parameters =
     const_cast<ParameterVector &>(parameters_in);
 
   const unsigned int Np = cast_int<unsigned int>
     (parameters.size());
   const unsigned int Nq = cast_int<unsigned int>
     (qoi.size());
 
   // An introduction to the problem:
   //
   // Residual R(u(p),p) = 0
   // partial R / partial u = J = system matrix
   //
   // This implies that:
   // d/dp(R) = 0
   // (partial R / partial p) +
   // (partial R / partial u) * (partial u / partial p) = 0
 
   // We first do an adjoint solve:
   // J^T * z = (partial q / partial u)
   // if we havent already or dont have an initial condition for the adjoint
   if (!this->is_adjoint_already_solved())
     {
       this->adjoint_solve(qoi_indices);
     }
 
   this->assemble_residual_derivatives(parameters_in);
 
   // Get ready to fill in senstivities:
   sensitivities.allocate_data(qoi_indices, *this, parameters);
 
   // We use the identities:
   // dq/dp = (partial q / partial p) + (partial q / partial u) *
   //         (partial u / partial p)
   // dq/dp = (partial q / partial p) + (J^T * z) *
   //         (partial u / partial p)
   // dq/dp = (partial q / partial p) + z * J *
   //         (partial u / partial p)
 
   // Leading to our final formula:
   // dq/dp = (partial q / partial p) - z * (partial R / partial p)
 
   // In the case of adjoints with heterogenous Dirichlet boundary
   // function phi, where
   // q :=  S(u) - R(u,phi)
   // the final formula works out to:
   // dq/dp = (partial S / partial p) - z * (partial R / partial p)
   // Because we currently have no direct access to
   // (partial S / partial p), we use the identity
   // (partial S / partial p) = (partial q / partial p) +
   //                           phi * (partial R / partial p)
   // to derive an equivalent equation:
   // dq/dp = (partial q / partial p) - (z-phi) * (partial R / partial p)
 
   // Since z-phi degrees of freedom are zero for constrained indices,
   // we can use the same constrained -(partial R / partial p) that we
   // use for forward sensitivity solves, taking into account the
   // differing sign convention.
   //
   // Since that vector is constrained, its constrained indices are
   // zero, so its product with phi is zero, so we can neglect the
   // evaluation of phi terms.
 
   for (unsigned int j=0; j != Np; ++j)
     {
       // We currently get partial derivatives via central differencing
 
       // (partial q / partial p) ~= (q(p+dp)-q(p-dp))/(2*dp)
       // (partial R / partial p) ~= (rhs(p+dp) - rhs(p-dp))/(2*dp)
 
       Number old_parameter = *parameters[j];
 
       const Real delta_p =
         TOLERANCE * std::max(std::abs(old_parameter), 1e-3);
 
       *parameters[j] = old_parameter - delta_p;
       this->assemble_qoi(qoi_indices);
       std::vector<Number> qoi_minus = this->qoi;
 
       NumericVector<Number> & neg_partialR_partialp = this->get_sensitivity_rhs(j);
 
       *parameters[j] = old_parameter + delta_p;
       this->assemble_qoi(qoi_indices);
       std::vector<Number> & qoi_plus = this->qoi;
 
       std::vector<Number> partialq_partialp(Nq, 0);
       for (unsigned int i=0; i != Nq; ++i)
         if (qoi_indices.has_index(i))
           partialq_partialp[i] = (qoi_plus[i] - qoi_minus[i]) / (2.*delta_p);
 
       // Don't leave the parameter changed
       *parameters[j] = old_parameter;
 
       for (unsigned int i=0; i != Nq; ++i)
         if (qoi_indices.has_index(i))
           sensitivities[i][j] = partialq_partialp[i] +
             neg_partialR_partialp.dot(this->get_adjoint_solution(i));
     }
 
   // All parameters have been reset.
   // Reset the original qoi.
 
   this->assemble_qoi(qoi_indices);
 }

std::pair< unsigned int, Real > libMesh::ImplicitSystem::adjoint_solve ( const QoISet & qoi_indices = QoISet() )

virtualinherited

Assembles & solves the linear system (dR/du)^T*z = dq/du, for those quantities of interest q specified by qoi_indices.

Leave qoi_indices empty to solve all adjoint problems.

Returns: A pair with the total number of linear iterations performed and the (sum of the) final residual norms

Reimplemented from libMesh::System.

Reimplemented in libMesh::DifferentiableSystem.

Definition at line 364 of file implicit_system.C.

References libMesh::System::add_adjoint_solution(), libMesh::LinearSolver< T >::adjoint_solve(), libMesh::System::assemble_before_solve, libMesh::ExplicitSystem::assemble_qoi_derivative(), libMesh::ImplicitSystem::assembly(), libMesh::DofMap::enforce_adjoint_constraints_exactly(), libMesh::System::get_adjoint_rhs(), libMesh::System::get_adjoint_solution(), libMesh::System::get_dof_map(), libMesh::ImplicitSystem::get_linear_solve_parameters(), libMesh::ImplicitSystem::get_linear_solver(), libMesh::QoISet::has_index(), libMesh::ImplicitSystem::matrix, libMesh::System::qoi, and libMesh::ImplicitSystem::release_linear_solver().

Referenced by libMesh::ImplicitSystem::adjoint_qoi_parameter_sensitivity(), libMesh::DifferentiableSystem::adjoint_solve(), libMesh::ImplicitSystem::assembly(), libMesh::ImplicitSystem::qoi_parameter_hessian(), and libMesh::ImplicitSystem::qoi_parameter_hessian_vector_product().

 {
   // Log how long the linear solve takes.
   LOG_SCOPE("adjoint_solve()", "ImplicitSystem");
 
   if (this->assemble_before_solve)
     // Assemble the linear system
     this->assembly (/* get_residual = */ false,
                     /* get_jacobian = */ true);
 
   // The adjoint problem is linear
   LinearSolver<Number> * linear_solver = this->get_linear_solver();
 
   // Reset and build the RHS from the QOI derivative
   this->assemble_qoi_derivative(qoi_indices,
                                 /* include_liftfunc = */ false,
                                 /* apply_constraints = */ true);
 
   // Our iteration counts and residuals will be sums of the individual
   // results
   std::pair<unsigned int, Real> solver_params =
     this->get_linear_solve_parameters();
   std::pair<unsigned int, Real> totalrval = std::make_pair(0,0.0);
 
   for (std::size_t i=0; i != this->qoi.size(); ++i)
     if (qoi_indices.has_index(i))
       {
         const std::pair<unsigned int, Real> rval =
           linear_solver->adjoint_solve (*matrix, this->add_adjoint_solution(i),
                                         this->get_adjoint_rhs(i),
                                         solver_params.second,
                                         solver_params.first);
 
         totalrval.first  += rval.first;
         totalrval.second += rval.second;
       }
 
   this->release_linear_solver(linear_solver);
 
   // The linear solver may not have fit our constraints exactly
 #ifdef LIBMESH_ENABLE_CONSTRAINTS
   for (std::size_t i=0; i != this->qoi.size(); ++i)
     if (qoi_indices.has_index(i))
       this->get_dof_map().enforce_adjoint_constraints_exactly
         (this->get_adjoint_solution(i), i);
 #endif
 
   return totalrval;
 }

virtual void libMesh::TransientRBConstruction::allocate_data_structures ( )

protectedvirtual

Helper function that actually allocates all the data structures required by this class.

Reimplemented from libMesh::RBConstruction.

virtual void libMesh::LinearImplicitSystem::assemble ( )

inlinevirtualinherited

Prepares matrix and _dof_map for matrix assembly. Does not actually assemble anything. For matrix assembly, use the assemble() in derived classes. Should be overridden in derived classes.

Reimplemented from libMesh::ImplicitSystem.

Reimplemented in libMesh::FrequencySystem, and libMesh::NewmarkSystem.

Definition at line 104 of file linear_implicit_system.h.

References libMesh::ImplicitSystem::assemble(), libMesh::LinearImplicitSystem::assembly(), libMesh::LinearImplicitSystem::get_linear_solver(), libMesh::LinearImplicitSystem::release_linear_solver(), libMesh::LinearImplicitSystem::restrict_solve_to(), libMesh::LinearImplicitSystem::solve(), and libMesh::SUBSET_ZERO.

Referenced by libMesh::NewmarkSystem::assemble(), libMesh::FrequencySystem::assemble(), libMesh::LinearImplicitSystem::assembly(), and libMesh::LinearImplicitSystem::solve().

104 { ImplicitSystem::assemble(); }

libMesh::ImplicitSystem::assemble

virtual void assemble() libmesh_override

Definition: implicit_system.C:190

virtual void libMesh::TransientRBConstruction::assemble_affine_expansion	(	bool	skip_matrix_assembly,
		bool	skip_vector_assembly
	)

protectedvirtual

Override assemble_affine_expansion to also initialize RB_ic_proj_rhs_all_N, if necessary.

Reimplemented from libMesh::RBConstruction.

virtual void libMesh::TransientRBConstruction::assemble_all_affine_operators ( )

virtual

Assemble and store all the affine operators. Override to assemble the mass matrix operators.

Reimplemented from libMesh::RBConstruction.

virtual void libMesh::RBConstruction::assemble_all_affine_vectors ( )

protectedvirtualinherited

Assemble and store the affine RHS vectors.

virtual void libMesh::RBConstruction::assemble_all_output_vectors ( )

protectedvirtualinherited

Assemble and store the output vectors.

void libMesh::RBConstruction::assemble_Aq_matrix	(	unsigned int	q,
		SparseMatrix< Number > *	input_matrix,
		bool	apply_dof_constraints = `true`
	)

inherited

Assemble the q^th affine matrix and store it in input_matrix.

Referenced by libMesh::RBConstruction::is_quiet().

void libMesh::RBConstruction::assemble_Fq_vector	(	unsigned int	q,
		NumericVector< Number > *	input_vector,
		bool	apply_dof_constraints = `true`
	)

inherited

Assemble the q^th affine vector and store it in input_matrix.

Referenced by libMesh::RBConstruction::is_quiet().

void libMesh::RBConstruction::assemble_inner_product_matrix	(	SparseMatrix< Number > *	input_matrix,
		bool	apply_dof_constraints = `true`
	)

inherited

Assemble the inner product matrix and store it in input_matrix.

Referenced by libMesh::RBConstruction::is_quiet().

void libMesh::TransientRBConstruction::assemble_L2_matrix	(	SparseMatrix< Number > *	input_matrix,
		bool	apply_dirichlet_bc = `true`
	)

Assemble the L2 matrix.

void libMesh::TransientRBConstruction::assemble_mass_matrix ( SparseMatrix< Number > * input_matrix )

Assemble the mass matrix at the current parameter and store it in input_matrix.

virtual void libMesh::TransientRBConstruction::assemble_misc_matrices ( )

virtual

Override to assemble the L2 matrix as well.

Reimplemented from libMesh::RBConstruction.

void libMesh::TransientRBConstruction::assemble_Mq_matrix	(	unsigned int	q,
		SparseMatrix< Number > *	input_matrix,
		bool	apply_dirichlet_bc = `true`
	)

Assemble the q^th affine term of the mass matrix and store it in input_matrix.

void libMesh::ExplicitSystem::assemble_qoi ( const QoISet & qoi_indices = QoISet() )

virtualinherited

Prepares qoi for quantity of interest assembly, then calls user qoi function. Can be overridden in derived classes.

Reimplemented from libMesh::System.

Reimplemented in libMesh::FEMSystem.

Definition at line 89 of file explicit_system.C.

References libMesh::System::assemble_qoi(), libMesh::QoISet::has_index(), and libMesh::System::qoi.

Referenced by libMesh::ImplicitSystem::adjoint_qoi_parameter_sensitivity(), libMesh::ImplicitSystem::forward_qoi_parameter_sensitivity(), libMesh::ImplicitSystem::qoi_parameter_hessian(), libMesh::ImplicitSystem::qoi_parameter_hessian_vector_product(), and libMesh::ExplicitSystem::system().

 {
   // The user quantity of interest assembly gets to expect to
   // accumulate on initially zero values
   for (std::size_t i=0; i != qoi.size(); ++i)
     if (qoi_indices.has_index(i))
       qoi[i] = 0;
 
   Parent::assemble_qoi (qoi_indices);
 }

void libMesh::ExplicitSystem::assemble_qoi_derivative	(	const QoISet &	qoi_indices = `QoISet()`,
		bool	include_liftfunc = `true`,
		bool	apply_constraints = `true`
	)

virtualinherited

Prepares adjoint_rhs for quantity of interest derivative assembly, then calls user qoi derivative function. Can be overridden in derived classes.

Reimplemented from libMesh::System.

Reimplemented in libMesh::FEMSystem.

Definition at line 102 of file explicit_system.C.

References libMesh::System::add_adjoint_rhs(), libMesh::System::assemble_qoi_derivative(), libMesh::QoISet::has_index(), libMesh::System::qoi, and libMesh::NumericVector< T >::zero().

Referenced by libMesh::ImplicitSystem::adjoint_solve(), libMesh::ImplicitSystem::forward_qoi_parameter_sensitivity(), libMesh::ImplicitSystem::qoi_parameter_hessian(), libMesh::ImplicitSystem::qoi_parameter_hessian_vector_product(), libMesh::ExplicitSystem::system(), and libMesh::ImplicitSystem::weighted_sensitivity_adjoint_solve().

 {
   // The user quantity of interest derivative assembly gets to expect
   // to accumulate on initially zero vectors
   for (std::size_t i=0; i != qoi.size(); ++i)
     if (qoi_indices.has_index(i))
       this->add_adjoint_rhs(i).zero();
 
   Parent::assemble_qoi_derivative (qoi_indices, include_liftfunc,
                                    apply_constraints);
 }

void libMesh::ImplicitSystem::assemble_residual_derivatives ( const ParameterVector & parameters )

virtualinherited

Residual parameter derivative function.

Uses finite differences by default.

This will assemble the sensitivity rhs vectors to hold -(partial R / partial p_i), making them ready to solve the forward sensitivity equation.

Can be overridden in derived classes.

Reimplemented from libMesh::System.

Definition at line 650 of file implicit_system.C.

References std::abs(), libMesh::System::add_sensitivity_rhs(), libMesh::ImplicitSystem::assembly(), libMesh::NumericVector< T >::close(), std::max(), libMesh::Real, libMesh::ExplicitSystem::rhs, libMesh::ParameterVector::size(), and libMesh::TOLERANCE.

Referenced by libMesh::ImplicitSystem::adjoint_qoi_parameter_sensitivity(), libMesh::ImplicitSystem::assembly(), and libMesh::ImplicitSystem::sensitivity_solve().

 {
   ParameterVector & parameters =
     const_cast<ParameterVector &>(parameters_in);
 
   const unsigned int Np = cast_int<unsigned int>
     (parameters.size());
 
   for (unsigned int p=0; p != Np; ++p)
     {
       NumericVector<Number> & sensitivity_rhs = this->add_sensitivity_rhs(p);
 
       // Approximate -(partial R / partial p) by
       // (R(p-dp) - R(p+dp)) / (2*dp)
 
       Number old_parameter = *parameters[p];
 
       const Real delta_p =
         TOLERANCE * std::max(std::abs(old_parameter), 1e-3);
 
       *parameters[p] -= delta_p;
 
       //      this->assembly(true, false, true);
       this->assembly(true, false, false);
       this->rhs->close();
       sensitivity_rhs = *this->rhs;
 
       *parameters[p] = old_parameter + delta_p;
 
       //      this->assembly(true, false, true);
       this->assembly(true, false, false);
       this->rhs->close();
 
       sensitivity_rhs -= *this->rhs;
       sensitivity_rhs /= (2*delta_p);
       sensitivity_rhs.close();
 
       *parameters[p] = old_parameter;
     }
 }

void libMesh::LinearImplicitSystem::assembly	(	bool	get_residual,
		bool	get_jacobian,
		bool	apply_heterogeneous_constraints = `false`,
		bool	apply_no_constraints = `false`
	)

virtualinherited

Assembles a residual in rhs and/or a jacobian in matrix, as requested.

Reimplemented from libMesh::ImplicitSystem.

Definition at line 366 of file linear_implicit_system.C.

References libMesh::NumericVector< T >::add_vector(), libMesh::LinearImplicitSystem::assemble(), libMesh::NumericVector< T >::close(), libMesh::SparseMatrix< T >::close(), libMesh::ImplicitSystem::matrix, libMesh::ExplicitSystem::rhs, and libMesh::System::solution.

Referenced by libMesh::LinearImplicitSystem::assemble().

 {
   // Residual R(u(p),p) := A(p)*u(p) - b(p)
   // partial R / partial u = A
 
   this->assemble();
   this->rhs->close();
   this->matrix->close();
 
   *(this->rhs) *= -1.0;
   this->rhs->add_vector(*this->solution, *this->matrix);
 }

void libMesh::System::attach_assemble_function ( void fptrEquationSystems &es, const std::string &name )

inherited

Register a user function to use in assembling the system matrix and RHS.

Definition at line 1801 of file system.C.

References libMesh::System::_assemble_system_function, libMesh::System::_assemble_system_object, libMesh::libmesh_assert(), libmesh_nullptr, and libMesh::out.

Referenced by libMesh::System::read_parallel_data().

 {
   libmesh_assert(fptr);
 
   if (_assemble_system_object != libmesh_nullptr)
     {
       libmesh_here();
       libMesh::out << "WARNING:  Cannot specify both assembly function and object!"
                    << std::endl;
 
       _assemble_system_object = libmesh_nullptr;
     }
 
   _assemble_system_function = fptr;
 }

void libMesh::System::attach_assemble_object ( System::Assembly & assemble_in )

inherited

Register a user object to use in assembling the system matrix and RHS.

Definition at line 1820 of file system.C.

References libMesh::System::_assemble_system_function, libMesh::System::_assemble_system_object, libmesh_nullptr, and libMesh::out.

Referenced by libMesh::System::read_parallel_data().

 {
   if (_assemble_system_function != libmesh_nullptr)
     {
       libmesh_here();
       libMesh::out << "WARNING:  Cannot specify both assembly object and function!"
                    << std::endl;
 
       _assemble_system_function = libmesh_nullptr;
     }
 
   _assemble_system_object = &assemble_in;
 }

void libMesh::System::attach_constraint_function ( void fptrEquationSystems &es, const std::string &name )

inherited

Register a user function for imposing constraints.

Definition at line 1836 of file system.C.

References libMesh::System::_constrain_system_function, libMesh::System::_constrain_system_object, libMesh::libmesh_assert(), libmesh_nullptr, and libMesh::out.

Referenced by libMesh::System::read_parallel_data().

 {
   libmesh_assert(fptr);
 
   if (_constrain_system_object != libmesh_nullptr)
     {
       libmesh_here();
       libMesh::out << "WARNING:  Cannot specify both constraint function and object!"
                    << std::endl;
 
       _constrain_system_object = libmesh_nullptr;
     }
 
   _constrain_system_function = fptr;
 }

void libMesh::System::attach_constraint_object ( System::Constraint & constrain )

inherited

Register a user object for imposing constraints.

Definition at line 1855 of file system.C.

References libMesh::System::_constrain_system_function, libMesh::System::_constrain_system_object, libmesh_nullptr, and libMesh::out.

Referenced by libMesh::System::read_parallel_data().

 {
   if (_constrain_system_function != libmesh_nullptr)
     {
       libmesh_here();
       libMesh::out << "WARNING:  Cannot specify both constraint object and function!"
                    << std::endl;
 
       _constrain_system_function = libmesh_nullptr;
     }
 
   _constrain_system_object = &constrain;
 }

void libMesh::System::attach_init_function ( void fptrEquationSystems &es, const std::string &name )

inherited

Register a user function to use in initializing the system.

Definition at line 1766 of file system.C.

References libMesh::System::_init_system_function, libMesh::System::_init_system_object, libMesh::libmesh_assert(), libmesh_nullptr, and libMesh::out.

Referenced by libMesh::System::read_parallel_data().

 {
   libmesh_assert(fptr);
 
   if (_init_system_object != libmesh_nullptr)
     {
       libmesh_here();
       libMesh::out << "WARNING:  Cannot specify both initialization function and object!"
                    << std::endl;
 
       _init_system_object = libmesh_nullptr;
     }
 
   _init_system_function = fptr;
 }

void libMesh::System::attach_init_object ( System::Initialization & init_in )

inherited

Register a user class to use to initialize the system.

Note: This is exclusive with the attach_init_function.

Definition at line 1785 of file system.C.

References libMesh::System::_init_system_function, libMesh::System::_init_system_object, libmesh_nullptr, and libMesh::out.

Referenced by libMesh::System::read_parallel_data().

 {
   if (_init_system_function != libmesh_nullptr)
     {
       libmesh_here();
       libMesh::out << "WARNING:  Cannot specify both initialization object and function!"
                    << std::endl;
 
       _init_system_function = libmesh_nullptr;
     }
 
   _init_system_object = &init_in;
 }

void libMesh::System::attach_QOI_derivative ( void fptrEquationSystems &es, const std::string &name, const QoISet &qoi_indices, bool include_liftfunc, bool apply_constraints )

inherited

Register a user function for evaluating derivatives of a quantity of interest with respect to test functions, whose values should be placed in System::rhs

Definition at line 1907 of file system.C.

References libMesh::System::_qoi_evaluate_derivative_function, libMesh::System::_qoi_evaluate_derivative_object, libMesh::libmesh_assert(), libmesh_nullptr, and libMesh::out.

Referenced by libMesh::System::read_parallel_data().

 {
   libmesh_assert(fptr);
 
   if (_qoi_evaluate_derivative_object != libmesh_nullptr)
     {
       libmesh_here();
       libMesh::out << "WARNING:  Cannot specify both QOI derivative function and object!"
                    << std::endl;
 
       _qoi_evaluate_derivative_object = libmesh_nullptr;
     }
 
   _qoi_evaluate_derivative_function = fptr;
 }

void libMesh::System::attach_QOI_derivative_object ( QOIDerivative & qoi_derivative )

inherited

Register a user object for evaluating derivatives of a quantity of interest with respect to test functions, whose values should be placed in System::rhs

Definition at line 1926 of file system.C.

References libMesh::System::_qoi_evaluate_derivative_function, libMesh::System::_qoi_evaluate_derivative_object, libmesh_nullptr, and libMesh::out.

Referenced by libMesh::System::read_parallel_data().

 {
   if (_qoi_evaluate_derivative_function != libmesh_nullptr)
     {
       libmesh_here();
       libMesh::out << "WARNING:  Cannot specify both QOI derivative object and function!"
                    << std::endl;
 
       _qoi_evaluate_derivative_function = libmesh_nullptr;
     }
 
   _qoi_evaluate_derivative_object = &qoi_derivative;
 }

void libMesh::System::attach_QOI_function ( void fptrEquationSystems &es, const std::string &name, const QoISet &qoi_indices )

inherited

Register a user function for evaluating the quantities of interest, whose values should be placed in System::qoi

Definition at line 1871 of file system.C.

References libMesh::System::_qoi_evaluate_function, libMesh::System::_qoi_evaluate_object, libMesh::libmesh_assert(), libmesh_nullptr, and libMesh::out.

Referenced by libMesh::System::read_parallel_data().

 {
   libmesh_assert(fptr);
 
   if (_qoi_evaluate_object != libmesh_nullptr)
     {
       libmesh_here();
       libMesh::out << "WARNING:  Cannot specify both QOI function and object!"
                    << std::endl;
 
       _qoi_evaluate_object = libmesh_nullptr;
     }
 
   _qoi_evaluate_function = fptr;
 }

void libMesh::System::attach_QOI_object ( QOI & qoi )

inherited

Register a user object for evaluating the quantities of interest, whose values should be placed in System::qoi

Definition at line 1891 of file system.C.

References libMesh::System::_qoi_evaluate_function, libMesh::System::_qoi_evaluate_object, libmesh_nullptr, and libMesh::out.

Referenced by libMesh::System::read_parallel_data().

 {
   if (_qoi_evaluate_function != libmesh_nullptr)
     {
       libmesh_here();
       libMesh::out << "WARNING:  Cannot specify both QOI object and function!"
                    << std::endl;
 
       _qoi_evaluate_function = libmesh_nullptr;
     }
 
   _qoi_evaluate_object = &qoi_in;
 }

void libMesh::LinearImplicitSystem::attach_shell_matrix ( ShellMatrix< Number > * shell_matrix )

inherited

This function enables the user to provide a shell matrix, i.e. a matrix that is not stored element-wise, but as a function. When you register your shell matrix using this function, calling solve() will no longer use the matrix member but the registered shell matrix instead. You can reset this behaviour to its original state by supplying a NULL pointer to this function.

Definition at line 158 of file linear_implicit_system.C.

References libMesh::LinearImplicitSystem::_shell_matrix.

Referenced by libMesh::LinearImplicitSystem::detach_shell_matrix(), and libMesh::LinearImplicitSystem::final_linear_residual().

 {
   _shell_matrix = shell_matrix;
 }

void libMesh::System::boundary_project_solution	(	const std::set< boundary_id_type > &	b,
		const std::vector< unsigned int > &	variables,
		FunctionBase< Number > *	f,
		FunctionBase< Gradient > *	g = `libmesh_nullptr`
	)

inherited

Projects arbitrary boundary functions onto a vector of degree of freedom values for the current system. Only degrees of freedom which affect the function's trace on a boundary in the set b are affected. Only degrees of freedom associated with the variables listed in the vector variables are projected. The function value f and its gradient g are user-provided cloneable functors. A gradient g is only required/used for projecting onto finite element spaces with continuous derivatives. If non-default Parameters are to be used, they can be provided in the parameters argument.

This method projects an arbitary boundary function onto the solution via L2 projections and nodal interpolations on each element.

Definition at line 999 of file system_projection.C.

References libMesh::System::boundary_project_vector().

Referenced by libMesh::System::project_vector(), and libMesh::System::system_type().

 {
   this->boundary_project_vector(b, variables, *solution, f, g);
 
   solution->localize(*current_local_solution);
 }

void libMesh::System::boundary_project_solution	(	const std::set< boundary_id_type > &	b,
		const std::vector< unsigned int > &	variables,
		Number	fptrconst Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name,
		Gradient	gptrconst Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name,
		const Parameters &	parameters
	)

inherited

Projects arbitrary boundary functions onto a vector of degree of freedom values for the current system. Only degrees of freedom which affect the function's trace on a boundary in the set b are affected. Only degrees of freedom associated with the variables listed in the vector variables are projected. The function value fptr and its gradient gptr are represented by function pointers. A gradient gptr is only required/used for projecting onto finite element spaces with continuous derivatives.

void libMesh::System::boundary_project_vector	(	const std::set< boundary_id_type > &	b,
		const std::vector< unsigned int > &	variables,
		NumericVector< Number > &	new_vector,
		FunctionBase< Number > *	f,
		FunctionBase< Gradient > *	g = `libmesh_nullptr`,
		int	is_adjoint = `-1`
	)		const

inherited

Projects arbitrary boundary functions onto a vector of degree of freedom values for the current system. Only degrees of freedom which affect the function's trace on a boundary in the set b are affected. Only degrees of freedom associated with the variables listed in the vector variables are projected. The function value f and its gradient g are user-provided cloneable functors. A gradient g is only required/used for projecting onto finite element spaces with continuous derivatives. If non-default Parameters are to be used, they can be provided in the parameters argument.

Constrain the new vector using the requested adjoint rather than primal constraints if is_adjoint is non-negative.

This method projects an arbitrary function via L2 projections and nodal interpolations on each element.

Definition at line 1041 of file system_projection.C.

References libMesh::NumericVector< T >::close(), libMesh::get_dof_map(), and libMesh::Threads::parallel_for().

Referenced by libMesh::System::boundary_project_solution(), and libMesh::System::system_type().

 {
   LOG_SCOPE ("boundary_project_vector()", "System");
 
   Threads::parallel_for
     (ConstElemRange (this->get_mesh().active_local_elements_begin(),
                      this->get_mesh().active_local_elements_end() ),
      BoundaryProjectSolution(b, variables, *this, f, g,
                              this->get_equation_systems().parameters,
                              new_vector)
      );
 
   // We don't do SCALAR dofs when just projecting the boundary, so
   // we're done here.
 
   new_vector.close();
 
 #ifdef LIBMESH_ENABLE_CONSTRAINTS
   if (is_adjoint == -1)
     this->get_dof_map().enforce_constraints_exactly(*this, &new_vector);
   else if (is_adjoint >= 0)
     this->get_dof_map().enforce_adjoint_constraints_exactly(new_vector,
                                                             is_adjoint);
 #endif
 }

void libMesh::System::boundary_project_vector	(	const std::set< boundary_id_type > &	b,
		const std::vector< unsigned int > &	variables,
		Number	fptrconst Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name,
		Gradient	gptrconst Point &p, const Parameters &parameters, const std::string &sys_name, const std::string &unknown_name,
		const Parameters &	parameters,
		NumericVector< Number > &	new_vector,
		int	is_adjoint = `-1`
	)		const

inherited

Projects arbitrary boundary functions onto a vector of degree of freedom values for the current system. Only degrees of freedom which affect the function's trace on a boundary in the set b are affected. Only degrees of freedom associated with the variables listed in the vector variables are projected. The function value fptr and its gradient gptr are represented by function pointers. A gradient gptr is only required/used for projecting onto finite element spaces with continuous derivatives.

Constrain the new vector using the requested adjoint rather than primal constraints if is_adjoint is non-negative.

void libMesh::RBConstructionBase< LinearImplicitSystem >::broadcast_parameters ( unsigned int proc_id )

inherited

Broadcasts parameters on processor proc_id to all processors.

virtual UniquePtr<DGFEMContext> libMesh::RBConstruction::build_context ( )

protectedvirtualinherited

Builds a DGFEMContext object with enough information to do evaluations on each element. We use DGFEMContext since it allows for both DG and continuous Galerkin formulations.

static UniquePtr<DirichletBoundary> libMesh::RBConstruction::build_zero_dirichlet_boundary_object ( )

staticinherited

It's helpful to be able to generate a DirichletBoundary that stores a ZeroFunction in order to impose Dirichlet boundary conditions.

Referenced by libMesh::RBConstruction::get_delta_N().

Real libMesh::System::calculate_norm	(	const NumericVector< Number > &	v,
		unsigned int	var,
		FEMNormType	norm_type,
		std::set< unsigned int > *	skip_dimensions = `libmesh_nullptr`
	)		const

inherited

Returns: A norm of variable var in the vector v, in the specified norm (e.g. L2, L_INF, H1)

Definition at line 1393 of file system.C.

References libMesh::DISCRETE_L1, libMesh::DISCRETE_L2, libMesh::DISCRETE_L_INF, libMesh::System::discrete_var_norm(), libMesh::L2, libMesh::System::n_vars(), and libMesh::Real.

Referenced by libMesh::AdaptiveTimeSolver::calculate_norm(), libMesh::UnsteadySolver::du(), and libMesh::System::project_solution_on_reinit().

 {
   //short circuit to save time
   if (norm_type == DISCRETE_L1 ||
       norm_type == DISCRETE_L2 ||
       norm_type == DISCRETE_L_INF)
     return discrete_var_norm(v,var,norm_type);
 
   // Not a discrete norm
   std::vector<FEMNormType> norms(this->n_vars(), L2);
   std::vector<Real> weights(this->n_vars(), 0.0);
   norms[var] = norm_type;
   weights[var] = 1.0;
   Real val = this->calculate_norm(v, SystemNorm(norms, weights), skip_dimensions);
   return val;
 }

Real libMesh::System::calculate_norm	(	const NumericVector< Number > &	v,
		const SystemNorm &	norm,
		std::set< unsigned int > *	skip_dimensions = `libmesh_nullptr`
	)		const

inherited

Returns: A norm of the vector v, using component_norm and component_scale to choose and weight the norms of each variable.

Definition at line 1415 of file system.C.

References libMesh::System::_dof_map, libMesh::System::_mesh, std::abs(), libMesh::MeshBase::active_local_elements_begin(), libMesh::MeshBase::active_local_elements_end(), libMesh::TypeVector< T >::add_scaled(), libMesh::TypeTensor< T >::add_scaled(), libMesh::FEGenericBase< OutputType >::build(), libMesh::NumericVector< T >::build(), libMesh::ParallelObject::comm(), libMesh::FEType::default_quadrature_rule(), libMesh::Elem::dim(), libMesh::DISCRETE_L1, libMesh::DISCRETE_L2, libMesh::DISCRETE_L_INF, libMesh::System::discrete_var_norm(), libMesh::DofMap::dof_indices(), libMesh::MeshBase::elem_dimensions(), libMesh::FEGenericBase< OutputType >::get_d2phi(), libMesh::System::get_dof_map(), libMesh::FEGenericBase< OutputType >::get_dphi(), libMesh::FEAbstract::get_JxW(), libMesh::System::get_mesh(), libMesh::FEGenericBase< OutputType >::get_phi(), libMesh::H1, libMesh::H1_SEMINORM, libMesh::H2, libMesh::H2_SEMINORM, libMesh::SystemNorm::is_discrete(), libMesh::L1, libMesh::NumericVector< T >::l1_norm(), libMesh::L2, libMesh::NumericVector< T >::l2_norm(), libMesh::L_INF, libMesh::libmesh_assert(), libmesh_nullptr, libMesh::NumericVector< T >::linfty_norm(), libMesh::NumericVector< T >::localize(), std::max(), libMesh::Parallel::Communicator::max(), libMesh::QBase::n_points(), libMesh::System::n_vars(), libMesh::TypeVector< T >::norm(), libMesh::TypeTensor< T >::norm(), libMesh::TensorTools::norm_sq(), libMesh::TypeVector< T >::norm_sq(), libMesh::TypeTensor< T >::norm_sq(), libMesh::Real, libMesh::FEAbstract::reinit(), libMesh::SERIAL, libMesh::NumericVector< T >::size(), libMesh::Parallel::Communicator::sum(), libMesh::SystemNorm::type(), libMesh::DofMap::variable_type(), libMesh::W1_INF_SEMINORM, libMesh::W2_INF_SEMINORM, libMesh::SystemNorm::weight(), and libMesh::SystemNorm::weight_sq().

 {
   // This function must be run on all processors at once
   parallel_object_only();
 
   LOG_SCOPE ("calculate_norm()", "System");
 
   // Zero the norm before summation
   Real v_norm = 0.;
 
   if (norm.is_discrete())
     {
       //Check to see if all weights are 1.0 and all types are equal
       FEMNormType norm_type0 = norm.type(0);
       unsigned int check_var = 0;
       for (; check_var != this->n_vars(); ++check_var)
         if ((norm.weight(check_var) != 1.0) || (norm.type(check_var) != norm_type0))
           break;
 
       //All weights were 1.0 so just do the full vector discrete norm
       if (check_var == this->n_vars())
         {
           if (norm_type0 == DISCRETE_L1)
             return v.l1_norm();
           if (norm_type0 == DISCRETE_L2)
             return v.l2_norm();
           if (norm_type0 == DISCRETE_L_INF)
             return v.linfty_norm();
           else
             libmesh_error_msg("Invalid norm_type0 = " << norm_type0);
         }
 
       for (unsigned int var=0; var != this->n_vars(); ++var)
         {
           // Skip any variables we don't need to integrate
           if (norm.weight(var) == 0.0)
             continue;
 
           v_norm += norm.weight(var) * discrete_var_norm(v, var, norm.type(var));
         }
 
       return v_norm;
     }
 
   // Localize the potentially parallel vector
   UniquePtr<NumericVector<Number> > local_v = NumericVector<Number>::build(this->comm());
   local_v->init(v.size(), true, SERIAL);
   v.localize (*local_v, _dof_map->get_send_list());
 
   // I'm not sure how best to mix Hilbert norms on some variables (for
   // which we'll want to square then sum then square root) with norms
   // like L_inf (for which we'll just want to take an absolute value
   // and then sum).
   bool using_hilbert_norm = true,
     using_nonhilbert_norm = true;
 
   // Loop over all variables
   for (unsigned int var=0; var != this->n_vars(); ++var)
     {
       // Skip any variables we don't need to integrate
       Real norm_weight_sq = norm.weight_sq(var);
       if (norm_weight_sq == 0.0)
         continue;
       Real norm_weight = norm.weight(var);
 
       // Check for unimplemented norms (rather than just returning 0).
       FEMNormType norm_type = norm.type(var);
       if ((norm_type==H1) ||
           (norm_type==H2) ||
           (norm_type==L2) ||
           (norm_type==H1_SEMINORM) ||
           (norm_type==H2_SEMINORM))
         {
           if (!using_hilbert_norm)
             libmesh_not_implemented();
           using_nonhilbert_norm = false;
         }
       else if ((norm_type==L1) ||
                (norm_type==L_INF) ||
                (norm_type==W1_INF_SEMINORM) ||
                (norm_type==W2_INF_SEMINORM))
         {
           if (!using_nonhilbert_norm)
             libmesh_not_implemented();
           using_hilbert_norm = false;
         }
       else
         libmesh_not_implemented();
 
       const FEType & fe_type = this->get_dof_map().variable_type(var);
 
       // Allow space for dims 0-3, even if we don't use them all
       std::vector<FEBase *> fe_ptrs(4,libmesh_nullptr);
       std::vector<QBase *> q_rules(4,libmesh_nullptr);
 
       const std::set<unsigned char> & elem_dims = _mesh.elem_dimensions();
 
       // Prepare finite elements for each dimension present in the mesh
       for (std::set<unsigned char>::const_iterator d_it = elem_dims.begin();
            d_it != elem_dims.end(); ++d_it)
         {
           if (skip_dimensions && skip_dimensions->find(*d_it) != skip_dimensions->end())
             continue;
 
           q_rules[*d_it] =
             fe_type.default_quadrature_rule (*d_it).release();
 
           // Construct finite element object
 
           fe_ptrs[*d_it] = FEBase::build(*d_it, fe_type).release();
 
           // Attach quadrature rule to FE object
           fe_ptrs[*d_it]->attach_quadrature_rule (q_rules[*d_it]);
         }
 
       std::vector<dof_id_type> dof_indices;
 
       // Begin the loop over the elements
       MeshBase::const_element_iterator       el     =
         this->get_mesh().active_local_elements_begin();
       const MeshBase::const_element_iterator end_el =
         this->get_mesh().active_local_elements_end();
 
       for ( ; el != end_el; ++el)
         {
           const Elem * elem = *el;
           const unsigned int dim = elem->dim();
 
           if (skip_dimensions && skip_dimensions->find(dim) != skip_dimensions->end())
             continue;
 
           FEBase * fe = fe_ptrs[dim];
           QBase * qrule = q_rules[dim];
           libmesh_assert(fe);
           libmesh_assert(qrule);
 
           const std::vector<Real> &               JxW = fe->get_JxW();
           const std::vector<std::vector<Real> > * phi = libmesh_nullptr;
           if (norm_type == H1 ||
               norm_type == H2 ||
               norm_type == L2 ||
               norm_type == L1 ||
               norm_type == L_INF)
             phi = &(fe->get_phi());
 
           const std::vector<std::vector<RealGradient> > * dphi = libmesh_nullptr;
           if (norm_type == H1 ||
               norm_type == H2 ||
               norm_type == H1_SEMINORM ||
               norm_type == W1_INF_SEMINORM)
             dphi = &(fe->get_dphi());
 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
           const std::vector<std::vector<RealTensor> > *   d2phi = libmesh_nullptr;
           if (norm_type == H2 ||
               norm_type == H2_SEMINORM ||
               norm_type == W2_INF_SEMINORM)
             d2phi = &(fe->get_d2phi());
 #endif
 
           fe->reinit (elem);
 
           this->get_dof_map().dof_indices (elem, dof_indices, var);
 
           const unsigned int n_qp = qrule->n_points();
 
           const unsigned int n_sf = cast_int<unsigned int>
             (dof_indices.size());
 
           // Begin the loop over the Quadrature points.
           for (unsigned int qp=0; qp<n_qp; qp++)
             {
               if (norm_type == L1)
                 {
                   Number u_h = 0.;
                   for (unsigned int i=0; i != n_sf; ++i)
                     u_h += (*phi)[i][qp] * (*local_v)(dof_indices[i]);
                   v_norm += norm_weight *
                     JxW[qp] * std::abs(u_h);
                 }
 
               if (norm_type == L_INF)
                 {
                   Number u_h = 0.;
                   for (unsigned int i=0; i != n_sf; ++i)
                     u_h += (*phi)[i][qp] * (*local_v)(dof_indices[i]);
                   v_norm = std::max(v_norm, norm_weight * std::abs(u_h));
                 }
 
               if (norm_type == H1 ||
                   norm_type == H2 ||
                   norm_type == L2)
                 {
                   Number u_h = 0.;
                   for (unsigned int i=0; i != n_sf; ++i)
                     u_h += (*phi)[i][qp] * (*local_v)(dof_indices[i]);
                   v_norm += norm_weight_sq *
                     JxW[qp] * TensorTools::norm_sq(u_h);
                 }
 
               if (norm_type == H1 ||
                   norm_type == H2 ||
                   norm_type == H1_SEMINORM)
                 {
                   Gradient grad_u_h;
                   for (unsigned int i=0; i != n_sf; ++i)
                     grad_u_h.add_scaled((*dphi)[i][qp], (*local_v)(dof_indices[i]));
                   v_norm += norm_weight_sq *
                     JxW[qp] * grad_u_h.norm_sq();
                 }
 
               if (norm_type == W1_INF_SEMINORM)
                 {
                   Gradient grad_u_h;
                   for (unsigned int i=0; i != n_sf; ++i)
                     grad_u_h.add_scaled((*dphi)[i][qp], (*local_v)(dof_indices[i]));
                   v_norm = std::max(v_norm, norm_weight * grad_u_h.norm());
                 }
 
 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
               if (norm_type == H2 ||
                   norm_type == H2_SEMINORM)
                 {
                   Tensor hess_u_h;
                   for (unsigned int i=0; i != n_sf; ++i)
                     hess_u_h.add_scaled((*d2phi)[i][qp], (*local_v)(dof_indices[i]));
                   v_norm += norm_weight_sq *
                     JxW[qp] * hess_u_h.norm_sq();
                 }
 
               if (norm_type == W2_INF_SEMINORM)
                 {
                   Tensor hess_u_h;
                   for (unsigned int i=0; i != n_sf; ++i)
                     hess_u_h.add_scaled((*d2phi)[i][qp], (*local_v)(dof_indices[i]));
                   v_norm = std::max(v_norm, norm_weight * hess_u_h.norm());
                 }
 #endif
             }
         }
 
       // Need to delete the FE and quadrature objects to prevent a memory leak
       for (std::size_t i=0; i<fe_ptrs.size(); i++)
         if (fe_ptrs[i])
           delete fe_ptrs[i];
 
       for (std::size_t i=0; i<q_rules.size(); i++)
         if (q_rules[i])
           delete q_rules[i];
     }
 
   if (using_hilbert_norm)
     {
       this->comm().sum(v_norm);
       v_norm = std::sqrt(v_norm);
     }
   else
     {
       this->comm().max(v_norm);
     }
 
   return v_norm;
 }

void libMesh::RBConstruction::check_convergence ( LinearSolver< Number > & input_solver )

protectedinherited

Check if the linear solver reports convergence. Throw an error when that is not the case.

Referenced by libMesh::RBConstruction::init_context().

virtual void libMesh::TransientRBConstruction::clear ( )

virtual

Clear all the data structures associated with the system.

Reimplemented from libMesh::TransientSystem< RBConstruction >.

const Parallel::Communicator& libMesh::ParallelObject::comm ( ) const

inlineinherited

Returns: A reference to the Parallel::Communicator object used by this mesh.

Definition at line 87 of file parallel_object.h.

References libMesh::ParallelObject::_communicator.

88 { return _communicator; }

libMesh::ParallelObject::_communicator

const Parallel::Communicator & _communicator

Definition: parallel_object.h:105

bool libMesh::System::compare	(	const System &	other_system,
		const Real	threshold,
		const bool	verbose
	)		const

virtualinherited

Returns: true when the other system contains identical data, up to the given threshold. Outputs some diagnostic info when verbose is set.

Definition at line 526 of file system.C.

References libMesh::System::_is_initialized, libMesh::System::_sys_name, libMesh::System::_vectors, libMesh::System::get_vector(), libMesh::libmesh_assert(), libMesh::System::n_vectors(), libMesh::System::name(), libMesh::out, and libMesh::System::solution.

Referenced by libMesh::EquationSystems::compare(), and libMesh::System::set_adjoint_already_solved().

 {
   // we do not care for matrices, but for vectors
   libmesh_assert (_is_initialized);
   libmesh_assert (other_system._is_initialized);
 
   if (verbose)
     {
       libMesh::out << "  Systems \"" << _sys_name << "\"" << std::endl;
       libMesh::out << "   comparing matrices not supported." << std::endl;
       libMesh::out << "   comparing names...";
     }
 
   // compare the name: 0 means identical
   const int name_result = _sys_name.compare(other_system.name());
   if (verbose)
     {
       if (name_result == 0)
         libMesh::out << " identical." << std::endl;
       else
         libMesh::out << "  names not identical." << std::endl;
       libMesh::out << "   comparing solution vector...";
     }
 
 
   // compare the solution: -1 means identical
   const int solu_result = solution->compare (*other_system.solution.get(),
                                              threshold);
 
   if (verbose)
     {
       if (solu_result == -1)
         libMesh::out << " identical up to threshold." << std::endl;
       else
         libMesh::out << "  first difference occured at index = "
                      << solu_result << "." << std::endl;
     }
 
 
   // safety check, whether we handle at least the same number
   // of vectors
   std::vector<int> ov_result;
 
   if (this->n_vectors() != other_system.n_vectors())
     {
       if (verbose)
         {
           libMesh::out << "   Fatal difference. This system handles "
                        << this->n_vectors() << " add'l vectors," << std::endl
                        << "   while the other system handles "
                        << other_system.n_vectors()
                        << " add'l vectors." << std::endl
                        << "   Aborting comparison." << std::endl;
         }
       return false;
     }
   else if (this->n_vectors() == 0)
     {
       // there are no additional vectors...
       ov_result.clear ();
     }
   else
     {
       // compare other vectors
       for (const_vectors_iterator pos = _vectors.begin();
            pos != _vectors.end(); ++pos)
         {
           if (verbose)
             libMesh::out << "   comparing vector \""
                          << pos->first << "\" ...";
 
           // assume they have the same name
           const NumericVector<Number> & other_system_vector =
             other_system.get_vector(pos->first);
 
           ov_result.push_back(pos->second->compare (other_system_vector,
                                                     threshold));
 
           if (verbose)
             {
               if (ov_result[ov_result.size()-1] == -1)
                 libMesh::out << " identical up to threshold." << std::endl;
               else
                 libMesh::out << " first difference occured at" << std::endl
                              << "   index = " << ov_result[ov_result.size()-1] << "." << std::endl;
             }
 
         }
 
     } // finished comparing additional vectors
 
 
   bool overall_result;
 
   // sum up the results
   if ((name_result==0) && (solu_result==-1))
     {
       if (ov_result.size()==0)
         overall_result = true;
       else
         {
           bool ov_identical;
           unsigned int n    = 0;
           do
             {
               ov_identical = (ov_result[n]==-1);
               n++;
             }
           while (ov_identical && n<ov_result.size());
           overall_result = ov_identical;
         }
     }
   else
     overall_result = false;
 
   if (verbose)
     {
       libMesh::out << "   finished comparisons, ";
       if (overall_result)
         libMesh::out << "found no differences." << std::endl << std::endl;
       else
         libMesh::out << "found differences." << std::endl << std::endl;
     }
 
   return overall_result;
 }

virtual void libMesh::RBConstruction::compute_Fq_representor_innerprods ( bool compute_inner_products = true )

protectedvirtualinherited

Compute the terms that are combined `online' to determine the dual norm of the residual. Here we compute the terms associated with the right-hand side. These terms are basis independent, hence we separate them from the rest of the calculations that are done in update_residual_terms. By default, inner product terms are also computed, but you can turn this feature off e.g. if you are already reading in that data from files.

virtual Real libMesh::RBConstruction::compute_max_error_bound ( )

virtualinherited

(i) Compute the a posteriori error bound for each set of parameters in the training set, (ii) set current_parameters to the parameters that maximize the error bound, and (iii) return the maximum error bound.

virtual void libMesh::RBConstruction::compute_output_dual_innerprods ( )

protectedvirtualinherited

Compute and store the dual norm of each output functional.

Number libMesh::System::current_solution ( const dof_id_type global_dof_number ) const

inherited

Returns: The current solution for the specified global DOF.

Definition at line 187 of file system.C.

References libMesh::System::_dof_map, and libMesh::System::current_local_solution.

Referenced by libMesh::ExactSolution::_compute_error(), libMesh::UniformRefinementEstimator::_estimate_error(), libMesh::HPCoarsenTest::add_projection(), libMesh::ExactErrorEstimator::estimate_error(), libMesh::System::point_gradient(), libMesh::System::point_hessian(), libMesh::System::point_value(), libMesh::HPCoarsenTest::select_refinement(), libMesh::EnsightIO::write_scalar_ascii(), and libMesh::EnsightIO::write_vector_ascii().

 {
   // Check the sizes
   libmesh_assert_less (global_dof_number, _dof_map->n_dofs());
   libmesh_assert_less (global_dof_number, current_local_solution->size());
 
   return (*current_local_solution)(global_dof_number);
 }

void libMesh::System::deactivate ( )

inlineinherited

Deactivates the system. Only active systems are solved.

Definition at line 2060 of file system.h.

References libMesh::System::_active.

Referenced by libMesh::System::get_equation_systems().

 {
   _active = false;
 }

void libMesh::LinearImplicitSystem::detach_shell_matrix ( )

inlineinherited

Detaches a shell matrix. Same as attach_shell_matrix(libmesh_nullptr).

Definition at line 179 of file linear_implicit_system.h.

References libMesh::LinearImplicitSystem::attach_shell_matrix(), and libmesh_nullptr.

179 { attach_shell_matrix(libmesh_nullptr); }

libMesh::LinearImplicitSystem::attach_shell_matrix

void attach_shell_matrix(ShellMatrix< Number > *shell_matrix)

Definition: linear_implicit_system.C:158

libmesh_nullptr

const class libmesh_nullptr_t libmesh_nullptr

void libMesh::ImplicitSystem::disable_cache ( )

virtualinherited

Avoids use of any cached data that might affect any solve result. Should be overridden in derived systems.

Reimplemented from libMesh::System.

Definition at line 298 of file implicit_system.C.

References libMesh::System::assemble_before_solve, libMesh::ImplicitSystem::get_linear_solver(), and libMesh::LinearSolver< T >::reuse_preconditioner().

                                     {
   this->assemble_before_solve = true;
   this->get_linear_solver()->reuse_preconditioner(false);
 }

void libMesh::ReferenceCounter::disable_print_counter_info ( )

staticinherited

Definition at line 107 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

Referenced by libMesh::LibMeshInit::LibMeshInit(), and libMesh::ReferenceCounter::n_objects().

 {
   _enable_print_counter = false;
   return;
 }

void libMesh::ReferenceCounter::disable_print_counter_info ( )

staticinherited

Definition at line 107 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

Referenced by libMesh::LibMeshInit::LibMeshInit(), and libMesh::ReferenceCounter::n_objects().

 {
   _enable_print_counter = false;
   return;
 }

void libMesh::ReferenceCounter::enable_print_counter_info ( )

staticinherited

Methods to enable/disable the reference counter output from print_info()

Definition at line 101 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

Referenced by libMesh::ReferenceCounter::n_objects().

 {
   _enable_print_counter = true;
   return;
 }

void libMesh::ReferenceCounter::enable_print_counter_info ( )

staticinherited

Methods to enable/disable the reference counter output from print_info()

Definition at line 101 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

Referenced by libMesh::ReferenceCounter::n_objects().

 {
   _enable_print_counter = true;
   return;
 }

virtual void libMesh::TransientRBConstruction::enrich_RB_space ( )

protectedvirtual

Add a new basis functions to the RB space. In the transient case we first perform a POD of the time-dependent "truth" and then add a certain number of POD modes to the reduced basis.

Reimplemented from libMesh::RBConstruction.

Real libMesh::LinearImplicitSystem::final_linear_residual ( ) const

inlineinherited

Returns: The final residual for the linear system solve.

Definition at line 163 of file linear_implicit_system.h.

References libMesh::LinearImplicitSystem::_final_linear_residual, and libMesh::LinearImplicitSystem::attach_shell_matrix().

163 { return _final_linear_residual; }

libMesh::LinearImplicitSystem::_final_linear_residual

Real _final_linear_residual

Definition: linear_implicit_system.h:198

void libMesh::ImplicitSystem::forward_qoi_parameter_sensitivity	(	const QoISet &	qoi_indices,
		const ParameterVector &	parameters,
		SensitivityData &	sensitivities
	)

virtualinherited

Solves for the derivative of each of the system's quantities of interest q in qoi[qoi_indices] with respect to each parameter in parameters, placing the result for qoi i and parameter j into sensitivities[i][j].

Uses the forward sensitivity method.

Currently uses finite differenced derivatives (partial q / partial p) and (partial R / partial p).

Reimplemented from libMesh::System.

Definition at line 804 of file implicit_system.C.

References std::abs(), libMesh::SensitivityData::allocate_data(), libMesh::ExplicitSystem::assemble_qoi(), libMesh::ExplicitSystem::assemble_qoi_derivative(), libMesh::ImplicitSystem::assembly(), libMesh::NumericVector< T >::close(), libMesh::SparseMatrix< T >::close(), libMesh::NumericVector< T >::dot(), libMesh::System::get_adjoint_rhs(), libMesh::System::get_sensitivity_solution(), libMesh::QoISet::has_index(), libMesh::ImplicitSystem::matrix, std::max(), libMesh::System::qoi, libMesh::Real, libMesh::ExplicitSystem::rhs, libMesh::ImplicitSystem::sensitivity_solve(), libMesh::ParameterVector::size(), and libMesh::TOLERANCE.

Referenced by libMesh::ImplicitSystem::assembly().

 {
   ParameterVector & parameters =
     const_cast<ParameterVector &>(parameters_in);
 
   const unsigned int Np = cast_int<unsigned int>
     (parameters.size());
   const unsigned int Nq = cast_int<unsigned int>
     (qoi.size());
 
   // An introduction to the problem:
   //
   // Residual R(u(p),p) = 0
   // partial R / partial u = J = system matrix
   //
   // This implies that:
   // d/dp(R) = 0
   // (partial R / partial p) +
   // (partial R / partial u) * (partial u / partial p) = 0
 
   // We first solve for (partial u / partial p) for each parameter:
   // J * (partial u / partial p) = - (partial R / partial p)
 
   this->sensitivity_solve(parameters);
 
   // Get ready to fill in senstivities:
   sensitivities.allocate_data(qoi_indices, *this, parameters);
 
   // We use the identity:
   // dq/dp = (partial q / partial p) + (partial q / partial u) *
   //         (partial u / partial p)
 
   // We get (partial q / partial u) from the user
   this->assemble_qoi_derivative(qoi_indices,
                                 /* include_liftfunc = */ true,
                                 /* apply_constraints = */ false);
 
   // We don't need these to be closed() in this function, but libMesh
   // standard practice is to have them closed() by the time the
   // function exits
   for (std::size_t i=0; i != this->qoi.size(); ++i)
     if (qoi_indices.has_index(i))
       this->get_adjoint_rhs(i).close();
 
   for (unsigned int j=0; j != Np; ++j)
     {
       // We currently get partial derivatives via central differencing
 
       // (partial q / partial p) ~= (q(p+dp)-q(p-dp))/(2*dp)
 
       Number old_parameter = *parameters[j];
 
       const Real delta_p =
         TOLERANCE * std::max(std::abs(old_parameter), 1e-3);
 
       *parameters[j] = old_parameter - delta_p;
       this->assemble_qoi(qoi_indices);
       std::vector<Number> qoi_minus = this->qoi;
 
       *parameters[j] = old_parameter + delta_p;
       this->assemble_qoi(qoi_indices);
       std::vector<Number> & qoi_plus = this->qoi;
 
       std::vector<Number> partialq_partialp(Nq, 0);
       for (unsigned int i=0; i != Nq; ++i)
         if (qoi_indices.has_index(i))
           partialq_partialp[i] = (qoi_plus[i] - qoi_minus[i]) / (2.*delta_p);
 
       // Don't leave the parameter changed
       *parameters[j] = old_parameter;
 
       for (unsigned int i=0; i != Nq; ++i)
         if (qoi_indices.has_index(i))
           sensitivities[i][j] = partialq_partialp[i] +
             this->get_adjoint_rhs(i).dot(this->get_sensitivity_solution(j));
     }
 
   // All parameters have been reset.
   // We didn't cache the original rhs or matrix for memory reasons,
   // but we can restore them to a state consistent solution -
   // principle of least surprise.
   this->assembly(true, true);
   this->rhs->close();
   this->matrix->close();
   this->assemble_qoi(qoi_indices);
 }

static void libMesh::RBConstructionBase< LinearImplicitSystem >::generate_training_parameters_deterministic	(	const Parallel::Communicator &	communicator,
		std::map< std::string, bool >	log_param_scale,
		std::map< std::string, NumericVector< Number > * > &	training_parameters_in,
		unsigned int	n_training_samples_in,
		const RBParameters &	min_parameters,
		const RBParameters &	max_parameters,
		bool	serial_training_set = `false`
	)

staticprotectedinherited

Static helper function for generating a deterministic set of parameters. Only works with 1 or 2 parameters (as defined by the lengths of min/max parameters vectors), otherwise throws an error.

static void libMesh::RBConstructionBase< LinearImplicitSystem >::generate_training_parameters_random	(	const Parallel::Communicator &	communicator,
		std::map< std::string, bool >	log_param_scale,
		std::map< std::string, NumericVector< Number > * > &	training_parameters_in,
		unsigned int	n_training_samples_in,
		const RBParameters &	min_parameters,
		const RBParameters &	max_parameters,
		int	training_parameters_random_seed = `-1`,
		bool	serial_training_set = `false`
	)

staticprotectedinherited

Static helper function for generating a randomized set of parameters.

Real libMesh::RBConstruction::get_abs_training_tolerance ( )

inlineinherited

Definition at line 186 of file rb_construction.h.

References libMesh::RBConstruction::abs_training_tolerance.

186 { return abs_training_tolerance; }

libMesh::RBConstruction::abs_training_tolerance

Real abs_training_tolerance

Definition: rb_construction.h:772

NumericVector< Number > & libMesh::System::get_adjoint_rhs ( unsigned int i = 0 )

inherited

Returns: A reference to one of the system's adjoint rhs vectors, by default the one corresponding to the first qoi. This what the user's QoI derivative code should assemble when setting up an adjoint problem

Definition at line 1045 of file system.C.

References libMesh::System::get_vector().

Referenced by libMesh::ImplicitSystem::adjoint_solve(), libMesh::ImplicitSystem::forward_qoi_parameter_sensitivity(), libMesh::System::project_solution_on_reinit(), libMesh::ImplicitSystem::qoi_parameter_hessian(), libMesh::ImplicitSystem::qoi_parameter_hessian_vector_product(), and libMesh::ImplicitSystem::weighted_sensitivity_adjoint_solve().

 {
   std::ostringstream adjoint_rhs_name;
   adjoint_rhs_name << "adjoint_rhs" << i;
 
   return this->get_vector(adjoint_rhs_name.str());
 }

const NumericVector< Number > & libMesh::System::get_adjoint_rhs ( unsigned int i = 0 ) const

inherited

Returns: A reference to one of the system's adjoint rhs vectors, by default the one corresponding to the first qoi.

Definition at line 1055 of file system.C.

References libMesh::System::get_vector().

 {
   std::ostringstream adjoint_rhs_name;
   adjoint_rhs_name << "adjoint_rhs" << i;
 
   return this->get_vector(adjoint_rhs_name.str());
 }

NumericVector< Number > & libMesh::System::get_adjoint_solution ( unsigned int i = 0 )

inherited

Returns: A reference to one of the system's adjoint solution vectors, by default the one corresponding to the first qoi.

Definition at line 983 of file system.C.

References libMesh::System::get_vector().

Referenced by libMesh::UniformRefinementEstimator::_estimate_error(), libMesh::ImplicitSystem::adjoint_solve(), libMesh::AdjointRefinementEstimator::estimate_error(), libMesh::AdjointResidualErrorEstimator::estimate_error(), libMesh::System::project_solution_on_reinit(), libMesh::ImplicitSystem::qoi_parameter_hessian(), libMesh::ImplicitSystem::qoi_parameter_hessian_vector_product(), and libMesh::ImplicitSystem::weighted_sensitivity_adjoint_solve().

 {
   std::ostringstream adjoint_name;
   adjoint_name << "adjoint_solution" << i;
 
   return this->get_vector(adjoint_name.str());
 }

const NumericVector< Number > & libMesh::System::get_adjoint_solution ( unsigned int i = 0 ) const

inherited

Returns: A reference to one of the system's adjoint solution vectors, by default the one corresponding to the first qoi.

Definition at line 993 of file system.C.

References libMesh::System::get_vector().

 {
   std::ostringstream adjoint_name;
   adjoint_name << "adjoint_solution" << i;
 
   return this->get_vector(adjoint_name.str());
 }

virtual void libMesh::TransientRBConstruction::get_all_matrices ( std::map< std::string, SparseMatrix< Number > * > & all_matrices )

virtual

Get a map that stores pointers to all of the matrices.

Reimplemented from libMesh::RBConstruction.

void libMesh::System::get_all_variable_numbers ( std::vector< unsigned int > & all_variable_numbers ) const

inherited

Fills all_variable_numbers with all the variable numbers for the variables that have been added to this system.

Definition at line 1272 of file system.C.

References libMesh::System::_variable_numbers, and libMesh::System::n_vars().

Referenced by libMesh::System::project_solution_on_reinit().

 {
   all_variable_numbers.resize(n_vars());
 
   // Make sure the variable exists
   std::map<std::string, unsigned short int>::const_iterator
     it = _variable_numbers.begin();
   std::map<std::string, unsigned short int>::const_iterator
     it_end = _variable_numbers.end();
 
   unsigned int count = 0;
   for ( ; it != it_end; ++it)
     {
       all_variable_numbers[count] = it->second;
       count++;
     }
 }

virtual void libMesh::RBConstruction::get_all_vectors ( std::map< std::string, NumericVector< Number > * > & all_vectors )

virtualinherited

Get a map that stores pointers to all of the vectors.

Referenced by libMesh::RBConstruction::is_quiet().

SparseMatrix<Number>* libMesh::RBConstruction::get_Aq ( unsigned int q )

inherited

Get a pointer to Aq.

Referenced by libMesh::RBConstruction::is_quiet().

static Real libMesh::RBParametrized::get_closest_value	(	Real	value,
		const std::vector< Real > &	list_of_values
	)

staticinherited

Returns: The closest entry to value from list_of_values.

Real libMesh::RBTemporalDiscretization::get_control ( const unsigned int k ) const

inherited

Get/set the RHS control.

bool libMesh::RBConstruction::get_convergence_assertion_flag ( ) const

protectedinherited

Getter for the flag determining if convergence should be checked after each solve.

Referenced by libMesh::RBConstruction::init_context().

unsigned int libMesh::RBConstruction::get_delta_N ( ) const

inlineinherited

Get delta_N, the number of basis functions we add to the RB space per iteration of the greedy algorithm. For steady-state systems, this should be 1, but can be more than 1 for time-dependent systems.

Definition at line 393 of file rb_construction.h.

References libMesh::RBConstruction::build_zero_dirichlet_boundary_object(), libMesh::RBConstruction::delta_N, libMesh::RBConstruction::get_inner_product_assembly(), libMesh::RBConstruction::set_convergence_assertion_flag(), libMesh::RBConstruction::set_inner_product_assembly(), and libMesh::RBConstruction::zero_constrained_dofs_on_vector().

393 { return delta_N; }

libMesh::RBConstruction::delta_N

unsigned int delta_N

Definition: rb_construction.h:698

Real libMesh::RBTemporalDiscretization::get_delta_t ( ) const

inherited

Get/set delta_t, the time-step size.

const std::string& libMesh::RBConstructionBase< LinearImplicitSystem >::get_deterministic_training_parameter_name ( ) const

inherited

Get the name of the parameter that we will generate deterministic training parameters for.

unsigned int libMesh::RBConstructionBase< LinearImplicitSystem >::get_deterministic_training_parameter_repeats ( ) const

inherited

Get the number of times each sample of the deterministic training parameter is repeated.

const std::map< std::string, std::vector<Real> >& libMesh::RBParametrized::get_discrete_parameter_values ( ) const

inherited

Get a const reference to the discrete parameter values.

const DofMap & libMesh::System::get_dof_map ( ) const

inlineinherited

Returns: A constant reference to this system's _dof_map.

Definition at line 2028 of file system.h.

References libMesh::System::_dof_map.

 {
   return *_dof_map;
 }

DofMap & libMesh::System::get_dof_map ( )

inlineinherited

Returns: A writable reference to this system's _dof_map.

Definition at line 2036 of file system.h.

References libMesh::System::_dof_map.

 {
   return *_dof_map;
 }

const EquationSystems& libMesh::System::get_equation_systems ( ) const

inlineinherited

Returns: A constant reference to this system's parent EquationSystems object.

Definition at line 712 of file system.h.

References libMesh::System::_equation_systems.

712 { return _equation_systems; }

libMesh::System::_equation_systems

EquationSystems & _equation_systems

Definition: system.h:1869

EquationSystems& libMesh::System::get_equation_systems ( )

inlineinherited

Returns: A reference to this system's parent EquationSystems object.

Definition at line 717 of file system.h.

References libMesh::System::_equation_systems, libMesh::System::activate(), libMesh::System::active(), libMesh::System::deactivate(), and libMesh::System::set_basic_system_only().

717 { return _equation_systems; }

libMesh::System::_equation_systems

EquationSystems & _equation_systems

Definition: system.h:1869

const NumericVector<Number>& libMesh::TransientRBConstruction::get_error_temporal_data ( )

Get the column of temporal_data corresponding to the current time level. This gives access to the truth projection error data. If the RB basis is empty, then this corresponds to the truth solution data itself.

Referenced by set_delta_N().

Real libMesh::RBTemporalDiscretization::get_euler_theta ( ) const

inherited

Get/set euler_theta, parameter that determines the temporal discretization.

numeric_index_type libMesh::RBConstructionBase< LinearImplicitSystem >::get_first_local_training_index ( ) const

inherited

Get the first local index of the training parameters.

NumericVector<Number>* libMesh::RBConstruction::get_Fq ( unsigned int q )

inherited

Get a pointer to Fq.

Referenced by libMesh::RBConstruction::is_quiet().

static void libMesh::RBConstructionBase< LinearImplicitSystem >::get_global_max_error_pair	(	const Parallel::Communicator &	communicator,
		std::pair< numeric_index_type, Real > &	error_pair
	)

staticprotectedinherited

Static function to return the error pair (index,error) that is corresponds to the largest error on all processors.

const RBParameters& libMesh::RBConstruction::get_greedy_parameter ( unsigned int i )

inherited

Return the parameters chosen during the i^th step of the Greedy algorithm.

std::string libMesh::ReferenceCounter::get_info ( )

staticinherited

Gets a string containing the reference information.

Definition at line 47 of file reference_counter.C.

References libMesh::ReferenceCounter::_counts, and libMesh::Quality::name().

Referenced by libMesh::ReferenceCounter::print_info().

 {
 #if defined(LIBMESH_ENABLE_REFERENCE_COUNTING) && defined(DEBUG)
 
   std::ostringstream oss;
 
   oss << '\n'
       << " ---------------------------------------------------------------------------- \n"
       << "| Reference count information                                                |\n"
       << " ---------------------------------------------------------------------------- \n";
 
   for (Counts::iterator it = _counts.begin();
        it != _counts.end(); ++it)
     {
       const std::string name(it->first);
       const unsigned int creations    = it->second.first;
       const unsigned int destructions = it->second.second;
 
       oss << "| " << name << " reference count information:\n"
           << "|  Creations:    " << creations    << '\n'
           << "|  Destructions: " << destructions << '\n';
     }
 
   oss << " ---------------------------------------------------------------------------- \n";
 
   return oss.str();
 
 #else
 
   return "";
 
 #endif
 }

std::string libMesh::ReferenceCounter::get_info ( )

staticinherited

Gets a string containing the reference information.

Definition at line 47 of file reference_counter.C.

References libMesh::ReferenceCounter::_counts, and libMesh::Quality::name().

Referenced by libMesh::ReferenceCounter::print_info().

 {
 #if defined(LIBMESH_ENABLE_REFERENCE_COUNTING) && defined(DEBUG)
 
   std::ostringstream oss;
 
   oss << '\n'
       << " ---------------------------------------------------------------------------- \n"
       << "| Reference count information                                                |\n"
       << " ---------------------------------------------------------------------------- \n";
 
   for (Counts::iterator it = _counts.begin();
        it != _counts.end(); ++it)
     {
       const std::string name(it->first);
       const unsigned int creations    = it->second.first;
       const unsigned int destructions = it->second.second;
 
       oss << "| " << name << " reference count information:\n"
           << "|  Creations:    " << creations    << '\n'
           << "|  Destructions: " << destructions << '\n';
     }
 
   oss << " ---------------------------------------------------------------------------- \n";
 
   return oss.str();
 
 #else
 
   return "";
 
 #endif
 }

std::string libMesh::System::get_info ( ) const

inherited

Returns: A string containing information about the system.

Definition at line 1682 of file system.C.

Referenced by libMesh::System::read_parallel_data().

 {
   std::ostringstream oss;
 
 
   const std::string & sys_name = this->name();
 
   oss << "   System #"  << this->number() << ", \"" << sys_name << "\"\n"
       << "    Type \""  << this->system_type() << "\"\n"
       << "    Variables=";
 
   for (unsigned int vg=0; vg<this->n_variable_groups(); vg++)
     {
       const VariableGroup & vg_description (this->variable_group(vg));
 
       if (vg_description.n_variables() > 1) oss << "{ ";
       for (unsigned int vn=0; vn<vg_description.n_variables(); vn++)
         oss << "\"" << vg_description.name(vn) << "\" ";
       if (vg_description.n_variables() > 1) oss << "} ";
     }
 
   oss << '\n';
 
   oss << "    Finite Element Types=";
 #ifndef LIBMESH_ENABLE_INFINITE_ELEMENTS
   for (unsigned int vg=0; vg<this->n_variable_groups(); vg++)
     oss << "\""
         << Utility::enum_to_string<FEFamily>(this->get_dof_map().variable_group(vg).type().family)
         << "\" ";
 #else
   for (unsigned int vg=0; vg<this->n_variable_groups(); vg++)
     {
       oss << "\""
           << Utility::enum_to_string<FEFamily>(this->get_dof_map().variable_group(vg).type().family)
           << "\", \""
           << Utility::enum_to_string<FEFamily>(this->get_dof_map().variable_group(vg).type().radial_family)
           << "\" ";
     }
 
   oss << '\n' << "    Infinite Element Mapping=";
   for (unsigned int vg=0; vg<this->n_variable_groups(); vg++)
     oss << "\""
         << Utility::enum_to_string<InfMapType>(this->get_dof_map().variable_group(vg).type().inf_map)
         << "\" ";
 #endif
 
   oss << '\n';
 
   oss << "    Approximation Orders=";
   for (unsigned int vg=0; vg<this->n_variable_groups(); vg++)
     {
 #ifndef LIBMESH_ENABLE_INFINITE_ELEMENTS
       oss << "\""
           << Utility::enum_to_string<Order>(this->get_dof_map().variable_group(vg).type().order)
           << "\" ";
 #else
       oss << "\""
           << Utility::enum_to_string<Order>(this->get_dof_map().variable_group(vg).type().order)
           << "\", \""
           << Utility::enum_to_string<Order>(this->get_dof_map().variable_group(vg).type().radial_order)
           << "\" ";
 #endif
     }
 
   oss << '\n';
 
   oss << "    n_dofs()="             << this->n_dofs()             << '\n';
   oss << "    n_local_dofs()="       << this->n_local_dofs()       << '\n';
 #ifdef LIBMESH_ENABLE_CONSTRAINTS
   oss << "    n_constrained_dofs()=" << this->n_constrained_dofs() << '\n';
   oss << "    n_local_constrained_dofs()=" << this->n_local_constrained_dofs() << '\n';
 #endif
 
   oss << "    " << "n_vectors()="  << this->n_vectors()  << '\n';
   oss << "    " << "n_matrices()="  << this->n_matrices()  << '\n';
   //   oss << "    " << "n_additional_matrices()=" << this->n_additional_matrices() << '\n';
 
   oss << this->get_dof_map().get_info();
 
   return oss.str();
 }

ElemAssembly& libMesh::RBConstruction::get_inner_product_assembly ( )

inherited

Returns: A reference to the inner product assembly object

Referenced by libMesh::RBConstruction::get_delta_N().

SparseMatrix<Number>* libMesh::RBConstruction::get_inner_product_matrix ( )

inherited

Get a pointer to inner_product_matrix. Accessing via this function, rather than directly through the class member allows us to do error checking (e.g. inner_product_matrix is not defined in low-memory mode).

Referenced by libMesh::RBConstruction::is_quiet().

ElemAssembly& libMesh::TransientRBConstruction::get_L2_assembly ( )

Returns: A reference to the L2 assembly object

numeric_index_type libMesh::RBConstructionBase< LinearImplicitSystem >::get_last_local_training_index ( ) const

inherited

Get the last local index of the training parameters.

std::pair< unsigned int, Real > libMesh::ImplicitSystem::get_linear_solve_parameters ( ) const

virtualinherited

Returns: An integer corresponding to the upper iteration count limit and a Real corresponding to the convergence tolerance to be used in linear adjoint and/or sensitivity solves

Reimplemented in libMesh::NonlinearImplicitSystem, and libMesh::DifferentiableSystem.

Definition at line 1410 of file implicit_system.C.

References libMesh::Parameters::get(), libMesh::System::get_equation_systems(), libMesh::EquationSystems::parameters, and libMesh::Real.

Referenced by libMesh::ImplicitSystem::adjoint_solve(), libMesh::ImplicitSystem::sensitivity_solve(), libMesh::ImplicitSystem::system_type(), libMesh::ImplicitSystem::weighted_sensitivity_adjoint_solve(), and libMesh::ImplicitSystem::weighted_sensitivity_solve().

 {
   return std::make_pair(this->get_equation_systems().parameters.get<unsigned int>("linear solver maximum iterations"),
                         this->get_equation_systems().parameters.get<Real>("linear solver tolerance"));
 }

LinearSolver< Number > * libMesh::LinearImplicitSystem::get_linear_solver ( ) const

virtualinherited

Returns: A pointer to a linear solver appropriate for use in adjoint and/or sensitivity solves

Reimplemented from libMesh::ImplicitSystem.

Definition at line 353 of file linear_implicit_system.C.

References libMesh::LinearImplicitSystem::linear_solver.

Referenced by libMesh::LinearImplicitSystem::assemble().

 {
   return linear_solver.get();
 }

numeric_index_type libMesh::RBConstructionBase< LinearImplicitSystem >::get_local_n_training_samples ( ) const

inherited

Get the total number of training samples local to this processor.

SparseMatrix<Number>* libMesh::TransientRBConstruction::get_M_q ( unsigned int q )

Get a pointer to M_q.

const SparseMatrix< Number > & libMesh::ImplicitSystem::get_matrix ( const std::string & mat_name ) const

inherited

Returns: A const reference to this system's additional matrix named mat_name.

None of these matrices is involved in the solution process. Access is only granted when the matrix is already properly initialized.

Definition at line 255 of file implicit_system.C.

References libMesh::ImplicitSystem::_matrices.

Referenced by libMesh::NewmarkSystem::compute_matrix(), libMesh::EigenTimeSolver::solve(), and libMesh::NewmarkSystem::update_rhs().

 {
   // Make sure the matrix exists
   const_matrices_iterator pos = _matrices.find (mat_name);
 
   if (pos == _matrices.end())
     libmesh_error_msg("ERROR: matrix " << mat_name << " does not exist in this system!");
 
   return *(pos->second);
 }

SparseMatrix< Number > & libMesh::ImplicitSystem::get_matrix ( const std::string & mat_name )

inherited

Returns: A writable reference to this system's additional matrix named mat_name.

None of these matrices is involved in the solution process. Access is only granted when the matrix is already properly initialized.

Definition at line 268 of file implicit_system.C.

References libMesh::ImplicitSystem::_matrices.

 {
   // Make sure the matrix exists
   matrices_iterator pos = _matrices.find (mat_name);
 
   if (pos == _matrices.end())
     libmesh_error_msg("ERROR: matrix " << mat_name << " does not exist in this system!");
 
   return *(pos->second);
 }

virtual SparseMatrix<Number>& libMesh::TransientRBConstruction::get_matrix_for_output_dual_solves ( )

protectedvirtual

Override to return the L2 product matrix for output dual norm solves for transient state problems.

Reimplemented from libMesh::RBConstruction.

int libMesh::TransientRBConstruction::get_max_truth_solves ( ) const

inline

Get/set max_truth_solves, the maximum number of RB truth solves we are willing to compute in the transient case.

Note: In the steady state case, max_truth_solves is not needed since it is equivalent to Nmax.

Definition at line 208 of file transient_rb_construction.h.

References max_truth_solves.

208 { return max_truth_solves; }

libMesh::TransientRBConstruction::max_truth_solves

int max_truth_solves

Definition: transient_rb_construction.h:388

const MeshBase & libMesh::System::get_mesh ( ) const

inlineinherited

Returns: A constant reference to this systems's _mesh.

Definition at line 2012 of file system.h.

References libMesh::System::_mesh.

Referenced by libMesh::ExactSolution::_compute_error(), libMesh::HPCoarsenTest::add_projection(), libMesh::FEMSystem::assemble_qoi(), libMesh::FEMSystem::assemble_qoi_derivative(), libMesh::FEMSystem::assembly(), libMesh::System::calculate_norm(), DMCreateDomainDecomposition_libMesh(), DMCreateFieldDecomposition_libMesh(), DMlibMeshSetSystem_libMesh(), libMesh::WeightedPatchRecoveryErrorEstimator::estimate_error(), libMesh::JumpErrorEstimator::estimate_error(), libMesh::PatchRecoveryErrorEstimator::estimate_error(), libMesh::AdjointResidualErrorEstimator::estimate_error(), libMesh::ExactErrorEstimator::estimate_error(), libMesh::SystemSubsetBySubdomain::init(), libMesh::System::init_data(), libMesh::EigenSystem::init_matrices(), libMesh::ImplicitSystem::init_matrices(), libMesh::System::local_dof_indices(), libMesh::DofMap::max_constraint_error(), libMesh::FEMSystem::mesh_position_get(), libMesh::FEMSystem::mesh_position_set(), libMesh::BoundaryProjectSolution::operator()(), libMesh::petsc_auto_fieldsplit(), libMesh::System::point_gradient(), libMesh::System::point_hessian(), libMesh::System::point_value(), libMesh::FEMSystem::postprocess(), libMesh::ImplicitSystem::reinit(), libMesh::EigenSystem::reinit(), libMesh::HPSingularity::select_refinement(), libMesh::HPCoarsenTest::select_refinement(), libMesh::System::system_type(), and libMesh::System::zero_variable().

 {
   return _mesh;
 }

MeshBase & libMesh::System::get_mesh ( )

inlineinherited

Returns: A reference to this systems's _mesh.

Definition at line 2020 of file system.h.

References libMesh::System::_mesh.

 {
   return _mesh;
 }

unsigned int libMesh::RBParametrized::get_n_continuous_params ( ) const

inherited

Get the number of continuous parameters.

unsigned int libMesh::RBParametrized::get_n_discrete_params ( ) const

inherited

Get the number of discrete parameters.

unsigned int libMesh::RBParametrized::get_n_params ( ) const

inherited

Get the number of parameters.

unsigned int libMesh::RBTemporalDiscretization::get_n_time_steps ( ) const

inherited

Get/set the total number of time-steps.

numeric_index_type libMesh::RBConstructionBase< LinearImplicitSystem >::get_n_training_samples ( ) const

inherited

Get the total number of training samples.

unsigned int libMesh::RBConstruction::get_Nmax ( ) const

inlineinherited

Get/set Nmax, the maximum number of RB functions we are willing to compute.

Definition at line 199 of file rb_construction.h.

References libMesh::RBConstruction::Nmax, and libMesh::RBConstruction::set_Nmax().

199 { return Nmax; }

libMesh::RBConstruction::Nmax

unsigned int Nmax

Definition: rb_construction.h:692

SparseMatrix<Number>* libMesh::RBConstruction::get_non_dirichlet_Aq ( unsigned int q )

inherited

Get a pointer to non_dirichlet_Aq.

Referenced by libMesh::RBConstruction::is_quiet().

NumericVector<Number>* libMesh::RBConstruction::get_non_dirichlet_Fq ( unsigned int q )

inherited

Get a pointer to non-Dirichlet Fq.

Referenced by libMesh::RBConstruction::is_quiet().

SparseMatrix<Number>* libMesh::TransientRBConstruction::get_non_dirichlet_M_q ( unsigned int q )

Get a pointer to non_dirichlet_M_q.

NumericVector<Number>* libMesh::RBConstruction::get_non_dirichlet_output_vector	(	unsigned int	n,
		unsigned int	q_l
	)

inherited

Get a pointer to non-Dirichlet output vector.

Referenced by libMesh::RBConstruction::is_quiet().

bool libMesh::RBConstruction::get_normalize_rb_bound_in_greedy ( )

inlineinherited

Definition at line 193 of file rb_construction.h.

References libMesh::RBConstruction::normalize_rb_bound_in_greedy.

193 { return normalize_rb_bound_in_greedy; }

libMesh::RBConstruction::normalize_rb_bound_in_greedy

bool normalize_rb_bound_in_greedy

Definition: rb_construction.h:778

NumericVector<Number>* libMesh::RBConstruction::get_output_vector	(	unsigned int	n,
		unsigned int	q_l
	)

inherited

Get a pointer to the n^th output.

Referenced by libMesh::RBConstruction::is_quiet().

virtual void libMesh::RBConstruction::get_output_vectors ( std::map< std::string, NumericVector< Number > * > & all_vectors )

virtualinherited

Get a map that stores pointers to all of the vectors.

Referenced by libMesh::RBConstruction::is_quiet().

Real libMesh::RBParametrized::get_parameter_max ( const std::string & param_name ) const

inherited

Get maximum allowable value of parameter param_name.

Real libMesh::RBParametrized::get_parameter_min ( const std::string & param_name ) const

inherited

Get minimum allowable value of parameter param_name.

std::set<std::string> libMesh::RBParametrized::get_parameter_names ( ) const

inherited

Get a set that stores the parameter names.

const RBParameters& libMesh::RBParametrized::get_parameters ( ) const

inherited

Get the current parameters.

const RBParameters& libMesh::RBParametrized::get_parameters_max ( ) const

inherited

Get an RBParameters object that specifies the maximum allowable value for each parameter.

const RBParameters& libMesh::RBParametrized::get_parameters_min ( ) const

inherited

Get an RBParameters object that specifies the minimum allowable value for each parameter.

RBParameters libMesh::RBConstructionBase< LinearImplicitSystem >::get_params_from_training_set ( unsigned int index )

protectedinherited

Return the RBParameters in index index of training set.

Real libMesh::TransientRBConstruction::get_POD_tol ( ) const

inline

Get/set POD_tol

Definition at line 214 of file transient_rb_construction.h.

References POD_tol.

214 { return POD_tol; }

libMesh::TransientRBConstruction::POD_tol

Real POD_tol

Definition: transient_rb_construction.h:380

RBAssemblyExpansion& libMesh::RBConstruction::get_rb_assembly_expansion ( )

inherited

Returns: A reference to the rb_assembly_expansion object

virtual Real libMesh::RBConstruction::get_RB_error_bound ( )

protectedvirtualinherited

Returns: The RB error bound for the current parameters.

Used in the Greedy algorithm to select the next parameter.

Reimplemented in libMesh::RBEIMConstruction.

RBEvaluation& libMesh::RBConstruction::get_rb_evaluation ( )

inherited

Get a reference to the RBEvaluation object.

RBThetaExpansion& libMesh::RBConstruction::get_rb_theta_expansion ( )

inherited

Get a reference to the RBThetaExpansion object that that belongs to rb_eval.

Real libMesh::RBConstruction::get_rel_training_tolerance ( )

inlineinherited

Definition at line 179 of file rb_construction.h.

References libMesh::RBConstruction::rel_training_tolerance.

179 { return rel_training_tolerance; }

libMesh::RBConstruction::rel_training_tolerance

Real rel_training_tolerance

Definition: rb_construction.h:771

NumericVector< Number > & libMesh::System::get_sensitivity_rhs ( unsigned int i = 0 )

inherited

Returns: A reference to one of the system's sensitivity rhs vectors, by default the one corresponding to the first parameter. By default these vectors are built by the library, using finite differences, when assemble_residual_derivatives() is called.

When assembled, this vector should hold -(partial R / partial p_i)

Definition at line 1075 of file system.C.

References libMesh::System::get_vector().

Referenced by libMesh::ImplicitSystem::adjoint_qoi_parameter_sensitivity(), libMesh::System::project_solution_on_reinit(), and libMesh::ImplicitSystem::sensitivity_solve().

 {
   std::ostringstream sensitivity_rhs_name;
   sensitivity_rhs_name << "sensitivity_rhs" << i;
 
   return this->get_vector(sensitivity_rhs_name.str());
 }

const NumericVector< Number > & libMesh::System::get_sensitivity_rhs ( unsigned int i = 0 ) const

inherited

Returns: A reference to one of the system's sensitivity rhs vectors, by default the one corresponding to the first parameter.

Definition at line 1085 of file system.C.

References libMesh::System::get_vector().

 {
   std::ostringstream sensitivity_rhs_name;
   sensitivity_rhs_name << "sensitivity_rhs" << i;
 
   return this->get_vector(sensitivity_rhs_name.str());
 }

NumericVector< Number > & libMesh::System::get_sensitivity_solution ( unsigned int i = 0 )

inherited

Returns: A reference to one of the system's solution sensitivity vectors, by default the one corresponding to the first parameter.

Definition at line 930 of file system.C.

References libMesh::System::get_vector().

Referenced by libMesh::ImplicitSystem::forward_qoi_parameter_sensitivity(), libMesh::System::project_solution_on_reinit(), libMesh::ImplicitSystem::qoi_parameter_hessian(), and libMesh::ImplicitSystem::sensitivity_solve().

 {
   std::ostringstream sensitivity_name;
   sensitivity_name << "sensitivity_solution" << i;
 
   return this->get_vector(sensitivity_name.str());
 }

const NumericVector< Number > & libMesh::System::get_sensitivity_solution ( unsigned int i = 0 ) const

inherited

Returns: A reference to one of the system's solution sensitivity vectors, by default the one corresponding to the first parameter.

Definition at line 940 of file system.C.

References libMesh::System::get_vector().

 {
   std::ostringstream sensitivity_name;
   sensitivity_name << "sensitivity_solution" << i;
 
   return this->get_vector(sensitivity_name.str());
 }

ShellMatrix<Number>* libMesh::LinearImplicitSystem::get_shell_matrix ( )

inlineinherited

Returns: A pointer to the currently attached shell matrix, if any, or NULL else.

Definition at line 185 of file linear_implicit_system.h.

References libMesh::LinearImplicitSystem::_shell_matrix.

185 { return _shell_matrix; }

libMesh::LinearImplicitSystem::_shell_matrix

ShellMatrix< Number > * _shell_matrix

Definition: linear_implicit_system.h:203

unsigned int libMesh::RBTemporalDiscretization::get_time_step ( ) const

inherited

Get/set the current time-step.

const NumericVector< Number > & libMesh::System::get_vector ( const std::string & vec_name ) const

inherited

Returns: A const reference to this system's additional vector named vec_name. Access is only granted when the vector is already properly initialized.

Definition at line 788 of file system.C.

References libMesh::System::_vectors.

Referenced by libMesh::UniformRefinementEstimator::_estimate_error(), libMesh::UnsteadySolver::adjoint_advance_timestep(), libMesh::NewmarkSolver::advance_timestep(), libMesh::AdaptiveTimeSolver::advance_timestep(), libMesh::UnsteadySolver::advance_timestep(), libMesh::System::compare(), libMesh::NewmarkSolver::compute_initial_accel(), libMesh::UnsteadySolver::du(), libMesh::AdjointRefinementEstimator::estimate_error(), libMesh::System::get_adjoint_rhs(), libMesh::System::get_adjoint_solution(), libMesh::System::get_sensitivity_rhs(), libMesh::System::get_sensitivity_solution(), libMesh::System::get_weighted_sensitivity_adjoint_solution(), libMesh::System::get_weighted_sensitivity_solution(), libMesh::NewmarkSystem::initial_conditions(), libMesh::NewmarkSolver::project_initial_accel(), libMesh::SecondOrderUnsteadySolver::project_initial_rate(), libMesh::System::project_solution_on_reinit(), libMesh::SecondOrderUnsteadySolver::reinit(), libMesh::UnsteadySolver::reinit(), libMesh::MemorySolutionHistory::retrieve(), libMesh::UnsteadySolver::retrieve_timestep(), libMesh::TwostepTimeSolver::solve(), libMesh::TaoOptimizationSolver< T >::solve(), libMesh::NloptOptimizationSolver< T >::solve(), libMesh::FrequencySystem::solve(), libMesh::NewmarkSystem::update_rhs(), and libMesh::NewmarkSystem::update_u_v_a().

 {
   // Make sure the vector exists
   const_vectors_iterator pos = _vectors.find(vec_name);
 
   if (pos == _vectors.end())
     libmesh_error_msg("ERROR: vector " << vec_name << " does not exist in this system!");
 
   return *(pos->second);
 }

NumericVector< Number > & libMesh::System::get_vector ( const std::string & vec_name )

inherited

Returns: A writable reference to this system's additional vector named vec_name. Access is only granted when the vector is already properly initialized.

Definition at line 801 of file system.C.

References libMesh::System::_vectors.

 {
   // Make sure the vector exists
   vectors_iterator pos = _vectors.find(vec_name);
 
   if (pos == _vectors.end())
     libmesh_error_msg("ERROR: vector " << vec_name << " does not exist in this system!");
 
   return *(pos->second);
 }

const NumericVector< Number > & libMesh::System::get_vector ( const unsigned int vec_num ) const

inherited

Returns: A const reference to this system's additional vector number vec_num (where the vectors are counted starting with 0).

Definition at line 814 of file system.C.

References libMesh::libmesh_assert(), libMesh::System::vectors_begin(), and libMesh::System::vectors_end().

 {
   const_vectors_iterator v = vectors_begin();
   const_vectors_iterator v_end = vectors_end();
   unsigned int num = 0;
   while ((num<vec_num) && (v!=v_end))
     {
       num++;
       ++v;
     }
   libmesh_assert (v != v_end);
   return *(v->second);
 }

NumericVector< Number > & libMesh::System::get_vector ( const unsigned int vec_num )

inherited

Returns: A writable reference to this system's additional vector number vec_num (where the vectors are counted starting with 0).

Definition at line 830 of file system.C.

References libMesh::libmesh_assert(), libMesh::System::vectors_begin(), and libMesh::System::vectors_end().

 {
   vectors_iterator v = vectors_begin();
   vectors_iterator v_end = vectors_end();
   unsigned int num = 0;
   while ((num<vec_num) && (v!=v_end))
     {
       num++;
       ++v;
     }
   libmesh_assert (v != v_end);
   return *(v->second);
 }

NumericVector< Number > & libMesh::System::get_weighted_sensitivity_adjoint_solution ( unsigned int i = 0 )

inherited

Returns: A reference to one of the system's weighted sensitivity adjoint solution vectors, by default the one corresponding to the first qoi.

Definition at line 1015 of file system.C.

References libMesh::System::get_vector().

Referenced by libMesh::System::project_solution_on_reinit(), libMesh::ImplicitSystem::qoi_parameter_hessian_vector_product(), and libMesh::ImplicitSystem::weighted_sensitivity_adjoint_solve().

 {
   std::ostringstream adjoint_name;
   adjoint_name << "weighted_sensitivity_adjoint_solution" << i;
 
   return this->get_vector(adjoint_name.str());
 }

const NumericVector< Number > & libMesh::System::get_weighted_sensitivity_adjoint_solution ( unsigned int i = 0 ) const

inherited

Returns: A reference to one of the system's weighted sensitivity adjoint solution vectors, by default the one corresponding to the first qoi.

Definition at line 1025 of file system.C.

References libMesh::System::get_vector().

 {
   std::ostringstream adjoint_name;
   adjoint_name << "weighted_sensitivity_adjoint_solution" << i;
 
   return this->get_vector(adjoint_name.str());
 }

NumericVector< Number > & libMesh::System::get_weighted_sensitivity_solution ( )

inherited

Returns: A reference to the solution of the last weighted sensitivity solve

Definition at line 957 of file system.C.

References libMesh::System::get_vector().

Referenced by libMesh::System::project_solution_on_reinit(), libMesh::ImplicitSystem::qoi_parameter_hessian_vector_product(), and libMesh::ImplicitSystem::weighted_sensitivity_solve().

 {
   return this->get_vector("weighted_sensitivity_solution");
 }

const NumericVector< Number > & libMesh::System::get_weighted_sensitivity_solution ( ) const

inherited

Returns: A reference to the solution of the last weighted sensitivity solve

Definition at line 964 of file system.C.

References libMesh::System::get_vector().

 {
   return this->get_vector("weighted_sensitivity_solution");
 }

virtual bool libMesh::TransientRBConstruction::greedy_termination_test	(	Real	abs_greedy_error,
		Real	initial_greedy_error,
		int	count
	)

virtual

Function that indicates when to terminate the Greedy basis training.

Reimplemented from libMesh::RBConstruction.

bool libMesh::System::has_variable ( const std::string & var ) const

inherited

Returns: true if a variable named var exists in this System

Definition at line 1250 of file system.C.

References libMesh::System::_variable_numbers.

Referenced by libMesh::GMVIO::copy_nodal_solution(), and libMesh::System::project_solution_on_reinit().

 {
   return _variable_numbers.count(var);
 }

bool libMesh::ImplicitSystem::have_matrix ( const std::string & mat_name ) const

inlineinherited

Returns: true if this System has a matrix associated with the given name, false otherwise.

Definition at line 394 of file implicit_system.h.

Referenced by libMesh::ImplicitSystem::add_matrix(), and libMesh::EigenTimeSolver::init().

 {
   return (_matrices.count(mat_name));
 }

bool libMesh::System::have_vector ( const std::string & vec_name ) const

inlineinherited

Returns: true if this System has a vector associated with the given name, false otherwise.

Definition at line 2204 of file system.h.

References libMesh::System::_vectors.

Referenced by libMesh::System::add_vector(), libMesh::System::project_solution_on_reinit(), and libMesh::System::remove_vector().

 {
   return (_vectors.count(vec_name));
 }

bool& libMesh::System::hide_output ( )

inlineinherited

Returns: A writable reference to a boolean that determines if this system can be written to file or not. If set to true, then EquationSystems::write will ignore this system.

Definition at line 1655 of file system.h.

References libMesh::System::_hide_output, libMesh::System::discrete_var_norm(), end, libMesh::System::init_data(), libMesh::invalid_uint, libMesh::ReferenceCounter::n_objects(), libMesh::System::operator=(), libMesh::System::project_vector(), libMesh::System::read_SCALAR_dofs(), libMesh::System::read_serialized_blocked_dof_objects(), libMesh::System::read_serialized_vector(), libMesh::Real, and libMesh::System::System().

1655 { return _hide_output; }

libMesh::System::_hide_output

bool _hide_output

Definition: system.h:1988

bool libMesh::System::identify_variable_groups ( ) const

inlineinherited

Returns: true when VariableGroup structures should be automatically identified, false otherwise.

Definition at line 2180 of file system.h.

References libMesh::System::_identify_variable_groups.

Referenced by libMesh::System::add_variable(), and libMesh::System::project_solution_on_reinit().

 {
   return _identify_variable_groups;
 }

void libMesh::System::identify_variable_groups ( const bool ivg )

inlineinherited

Toggle automatic VariableGroup identification.

Definition at line 2188 of file system.h.

References libMesh::System::_identify_variable_groups.

 {
   _identify_variable_groups = ivg;
 }

void libMesh::ReferenceCounter::increment_constructor_count ( const std::string & name )

inlineprotectedinherited

Increments the construction counter. Should be called in the constructor of any derived class that will be reference counted.

Definition at line 185 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCounter::n_objects(), and libMesh::ReferenceCountedObject< RBParametrized >::ReferenceCountedObject().

 {
   Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
   std::pair<unsigned int, unsigned int> & p = _counts[name];
 
   p.first++;
 }

void libMesh::ReferenceCounter::increment_constructor_count ( const std::string & name )

inlineprotectedinherited

Increments the construction counter. Should be called in the constructor of any derived class that will be reference counted.

Definition at line 185 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCounter::n_objects(), and libMesh::ReferenceCountedObject< RBParametrized >::ReferenceCountedObject().

 {
   Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
   std::pair<unsigned int, unsigned int> & p = _counts[name];
 
   p.first++;
 }

void libMesh::ReferenceCounter::increment_destructor_count ( const std::string & name )

inlineprotectedinherited

Increments the destruction counter. Should be called in the destructor of any derived class that will be reference counted.

Definition at line 198 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCounter::n_objects(), and libMesh::ReferenceCountedObject< RBParametrized >::~ReferenceCountedObject().

 {
   Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
   std::pair<unsigned int, unsigned int> & p = _counts[name];
 
   p.second++;
 }

void libMesh::ReferenceCounter::increment_destructor_count ( const std::string & name )

inlineprotectedinherited

Increments the destruction counter. Should be called in the destructor of any derived class that will be reference counted.

Definition at line 198 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCounter::n_objects(), and libMesh::ReferenceCountedObject< RBParametrized >::~ReferenceCountedObject().

 {
   Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
   std::pair<unsigned int, unsigned int> & p = _counts[name];
 
   p.second++;
 }

void libMesh::System::init ( )

inherited

Initializes degrees of freedom on the current mesh. Sets the

Definition at line 230 of file system.C.

References libMesh::System::_basic_system_only, libMesh::System::init_data(), libMesh::System::is_initialized(), libMesh::libmesh_assert(), libMesh::System::n_vars(), and libMesh::System::user_initialization().

Referenced by libMesh::System::system().

 {
   // Calling init() twice on the same system currently works evil
   // magic, whether done directly or via EquationSystems::read()
   libmesh_assert(!this->is_initialized());
 
   // First initialize any required data:
   // either only the basic System data
   if (_basic_system_only)
     System::init_data();
   // or all the derived class' data too
   else
     this->init_data();
 
   // If no variables have been added to this system
   // don't do anything
   if (!this->n_vars())
     return;
 
   // Then call the user-provided intialization function
   this->user_initialization();
 }

virtual void libMesh::RBConstruction::init_context ( FEMContext & )

inlineprotectedvirtualinherited

Initialize the FEMContext prior to performing an element loop. Reimplement this in derived classes in order to call FE::get_*() as the particular physics requires.

Definition at line 672 of file rb_construction.h.

References libMesh::RBConstruction::check_convergence(), and libMesh::RBConstruction::get_convergence_assertion_flag().

672 {}

virtual void libMesh::TransientSystem< RBConstruction >::init_data ( )

protectedvirtualinherited

Initializes the member data fields associated with the system, so that, e.g., assemble() may be used.

Reimplemented from libMesh::RBConstructionBase< LinearImplicitSystem >.

void libMesh::ImplicitSystem::init_matrices ( )

protectedvirtualinherited

Initializes the matrices associated with this system.

Definition at line 106 of file implicit_system.C.

References libMesh::ImplicitSystem::_can_add_matrices, libMesh::ImplicitSystem::_matrices, libMesh::DofMap::attach_matrix(), libMesh::DofMap::compute_sparsity(), libMesh::System::get_dof_map(), libMesh::System::get_mesh(), libMesh::SparseMatrix< T >::initialized(), libMesh::DofMap::is_attached(), libMesh::libmesh_assert(), and libMesh::ImplicitSystem::matrix.

Referenced by libMesh::ImplicitSystem::init_data().

 {
   libmesh_assert(matrix);
 
   // Check for quick return in case the system matrix
   // (and by extension all the matrices) has already
   // been initialized
   if (matrix->initialized())
     return;
 
   // Get a reference to the DofMap
   DofMap & dof_map = this->get_dof_map();
 
   // no chance to add other matrices
   _can_add_matrices = false;
 
   // Tell the matrices about the dof map, and vice versa
   for (matrices_iterator pos = _matrices.begin();
        pos != _matrices.end(); ++pos)
     {
       SparseMatrix<Number> & m = *(pos->second);
       libmesh_assert (!m.initialized());
 
       // We want to allow repeated init() on systems, but we don't
       // want to attach the same matrix to the DofMap twice
       if (!dof_map.is_attached(m))
         dof_map.attach_matrix (m);
     }
 
   // Compute the sparsity pattern for the current
   // mesh and DOF distribution.  This also updates
   // additional matrices, \p DofMap now knows them
   dof_map.compute_sparsity (this->get_mesh());
 
   // Initialize matrices
   for (matrices_iterator pos = _matrices.begin();
        pos != _matrices.end(); ++pos)
     pos->second->init ();
 
   // Set the additional matrices to 0.
   for (matrices_iterator pos = _matrices.begin();
        pos != _matrices.end(); ++pos)
     pos->second->zero ();
 }

void libMesh::RBParametrized::initialize_parameters	(	const RBParameters &	mu_min_in,
		const RBParameters &	mu_max_in,
		const std::map< std::string, std::vector< Real > > &	discrete_parameter_values
	)

inherited

Initialize the parameter ranges and set current_parameters.

void libMesh::RBParametrized::initialize_parameters ( const RBParametrized & rb_parametrized )

inherited

Initialize the parameter ranges and set current_parameters.

virtual void libMesh::TransientRBConstruction::initialize_rb_construction	(	bool	skip_matrix_assembly = `false`,
		bool	skip_vector_assembly = `false`
	)

virtual

Allocate all the data structures necessary for the construction stage of the RB method. This function also performs matrix and vector assembly of the "truth" affine expansion.

Override to check that theta and assembly expansions are consistently sized.

Reimplemented from libMesh::RBConstruction.

virtual void libMesh::RBConstructionBase< LinearImplicitSystem >::initialize_training_parameters	(	const RBParameters &	mu_min,
		const RBParameters &	mu_max,
		unsigned int	n_training_parameters,
		std::map< std::string, bool >	log_param_scale,
		bool	deterministic = `true`
	)

virtualinherited

Initialize the parameter ranges and indicate whether deterministic or random training parameters should be used and whether or not we want the parameters to be scaled logarithmically.

virtual void libMesh::TransientRBConstruction::initialize_truth ( )

protectedvirtual

This function imposes a truth initial condition, defaults to zero initial condition if the flag nonzero_initialization is true.

bool libMesh::System::is_adjoint_already_solved ( ) const

inlineinherited

Accessor for the adjoint_already_solved boolean

Definition at line 372 of file system.h.

References libMesh::System::adjoint_already_solved.

Referenced by libMesh::ImplicitSystem::adjoint_qoi_parameter_sensitivity(), libMesh::AdjointRefinementEstimator::estimate_error(), libMesh::AdjointResidualErrorEstimator::estimate_error(), libMesh::ImplicitSystem::qoi_parameter_hessian(), and libMesh::ImplicitSystem::qoi_parameter_hessian_vector_product().

373 { return adjoint_already_solved;}

libMesh::System::adjoint_already_solved

bool adjoint_already_solved

Definition: system.h:1982

bool libMesh::RBParametrized::is_discrete_parameter ( const std::string & mu_name ) const

inherited

Is parameter mu_name discrete?

bool libMesh::System::is_initialized ( )

inlineinherited

Returns: true iff this system has been initialized.

Definition at line 2068 of file system.h.

References libMesh::System::_is_initialized.

Referenced by libMesh::System::add_variable(), libMesh::System::add_variables(), libMesh::System::init(), and libMesh::System::system().

 {
   return _is_initialized;
 }

bool libMesh::RBConstruction::is_quiet ( ) const

inlineinherited

Is the system in quiet mode?

Definition at line 213 of file rb_construction.h.

214 { return this->quiet_mode; }

libMesh::RBConstruction::quiet_mode

bool quiet_mode

Definition: rb_construction.h:704

bool libMesh::RBConstruction::is_rb_eval_initialized ( ) const

inherited

Returns: true if rb_eval is initialized. False, otherwise.

virtual void libMesh::RBConstruction::load_basis_function ( unsigned int i )

virtualinherited

Load the i^th RB function into the RBConstruction solution vector.

Reimplemented in libMesh::RBEIMConstruction.

Referenced by libMesh::RBConstruction::is_quiet().

virtual void libMesh::TransientRBConstruction::load_rb_solution ( )

virtual

Load the RB solution from the current time-level into the libMesh solution vector.

Reimplemented from libMesh::RBConstruction.

Referenced by set_delta_N().

virtual void libMesh::RBConstructionBase< LinearImplicitSystem >::load_training_set ( std::map< std::string, std::vector< Number > > & new_training_set )

virtualinherited

Overwrite the training parameters with new_training_set.

void libMesh::System::local_dof_indices	(	const unsigned int	var,
		std::set< dof_id_type > &	var_indices
	)		const

inherited

Fills the std::set with the degrees of freedom on the local processor corresponding the the variable number passed in.

Definition at line 1291 of file system.C.

References libMesh::MeshBase::active_local_elements_begin(), libMesh::MeshBase::active_local_elements_end(), libMesh::DofMap::dof_indices(), libMesh::DofMap::end_dof(), libMesh::DofMap::first_dof(), libMesh::System::get_dof_map(), and libMesh::System::get_mesh().

Referenced by libMesh::System::discrete_var_norm().

 {
   // Make sure the set is clear
   var_indices.clear();
 
   std::vector<dof_id_type> dof_indices;
 
   // Begin the loop over the elements
   MeshBase::const_element_iterator       el     =
     this->get_mesh().active_local_elements_begin();
   const MeshBase::const_element_iterator end_el =
     this->get_mesh().active_local_elements_end();
 
   const dof_id_type
     first_local = this->get_dof_map().first_dof(),
     end_local   = this->get_dof_map().end_dof();
 
   for ( ; el != end_el; ++el)
     {
       const Elem * elem = *el;
       this->get_dof_map().dof_indices (elem, dof_indices, var);
 
       for (std::size_t i=0; i<dof_indices.size(); i++)
         {
           dof_id_type dof = dof_indices[i];
 
           //If the dof is owned by the local processor
           if (first_local <= dof && dof < end_local)
             var_indices.insert(dof_indices[i]);
         }
     }
 }

void libMesh::TransientRBConstruction::mass_matrix_scaled_matvec	(	Number	scalar,
		NumericVector< Number > &	dest,
		NumericVector< Number > &	arg
	)

Perform a matrix-vector multiplication with the current mass matrix and store the result in dest.

dof_id_type libMesh::System::n_active_dofs ( ) const

inlineinherited

Returns: The number of active degrees of freedom for this System.

Definition at line 2196 of file system.h.

References libMesh::System::n_constrained_dofs(), and libMesh::System::n_dofs().

Referenced by libMesh::System::project_solution_on_reinit().

 {
   return this->n_dofs() - this->n_constrained_dofs();
 }

unsigned int libMesh::System::n_components ( ) const

inlineinherited

Returns: The total number of scalar components in the system's variables. This will equal n_vars() in the case of all scalar-valued variables.

Definition at line 2100 of file system.h.

References libMesh::System::_variables, libMesh::Variable::first_scalar_number(), and libMesh::Variable::n_components().

Referenced by libMesh::System::add_variables(), libMesh::WrappedFunction< Output >::operator()(), and libMesh::System::project_solution_on_reinit().

 {
   if (_variables.empty())
     return 0;
 
   const Variable & last = _variables.back();
   return last.first_scalar_number() + last.n_components();
 }

dof_id_type libMesh::System::n_constrained_dofs ( ) const

inherited

Returns: The total number of constrained degrees of freedom in the system.

Definition at line 150 of file system.C.

References libMesh::System::_dof_map.

Referenced by libMesh::System::get_info(), libMesh::System::n_active_dofs(), and libMesh::System::project_solution_on_reinit().

 {
 #ifdef LIBMESH_ENABLE_CONSTRAINTS
 
   return _dof_map->n_constrained_dofs();
 
 #else
 
   return 0;
 
 #endif
 }

dof_id_type libMesh::System::n_dofs ( ) const

inherited

Returns: The number of degrees of freedom in the system

Definition at line 143 of file system.C.

References libMesh::System::_dof_map.

Referenced by libMesh::System::add_vector(), libMesh::AdjointRefinementEstimator::estimate_error(), libMesh::System::get_info(), libMesh::SecondOrderUnsteadySolver::init_data(), libMesh::UnsteadySolver::init_data(), libMesh::System::init_data(), libMesh::OptimizationSystem::initialize_equality_constraints_storage(), libMesh::OptimizationSystem::initialize_inequality_constraints_storage(), libMesh::System::n_active_dofs(), libMesh::CondensedEigenSystem::n_global_non_condensed_dofs(), libMesh::System::project_solution_on_reinit(), libMesh::SecondOrderUnsteadySolver::reinit(), libMesh::UnsteadySolver::reinit(), libMesh::System::reinit(), and libMesh::System::restrict_vectors().

 {
   return _dof_map->n_dofs();
 }

unsigned int libMesh::LinearImplicitSystem::n_linear_iterations ( ) const

inlineinherited

Returns: The number of iterations taken for the most recent linear solve.

Definition at line 158 of file linear_implicit_system.h.

References libMesh::LinearImplicitSystem::_n_linear_iterations.

158 { return _n_linear_iterations; }

libMesh::LinearImplicitSystem::_n_linear_iterations

unsigned int _n_linear_iterations

Definition: linear_implicit_system.h:193

dof_id_type libMesh::System::n_local_constrained_dofs ( ) const

inherited

Returns: The number of constrained degrees of freedom on this processor.

Definition at line 165 of file system.C.

References libMesh::System::_dof_map.

Referenced by libMesh::System::get_info(), and libMesh::System::project_solution_on_reinit().

 {
 #ifdef LIBMESH_ENABLE_CONSTRAINTS
 
   return _dof_map->n_local_constrained_dofs();
 
 #else
 
   return 0;
 
 #endif
 }

dof_id_type libMesh::System::n_local_dofs ( ) const

inherited

Returns: The number of degrees of freedom local to this processor

Definition at line 180 of file system.C.

References libMesh::System::_dof_map, and libMesh::ParallelObject::processor_id().

Referenced by libMesh::System::add_vector(), libMesh::AdjointRefinementEstimator::estimate_error(), libMesh::System::get_info(), libMesh::SecondOrderUnsteadySolver::init_data(), libMesh::UnsteadySolver::init_data(), libMesh::System::init_data(), libMesh::OptimizationSystem::initialize_equality_constraints_storage(), libMesh::OptimizationSystem::initialize_inequality_constraints_storage(), libMesh::System::project_solution_on_reinit(), libMesh::SecondOrderUnsteadySolver::reinit(), libMesh::UnsteadySolver::reinit(), libMesh::System::reinit(), and libMesh::System::restrict_vectors().

 {
   return _dof_map->n_dofs_on_processor (this->processor_id());
 }

unsigned int libMesh::ImplicitSystem::n_matrices ( ) const

inlinevirtualinherited

Returns: The number of matrices handled by this system

Reimplemented from libMesh::System.

Definition at line 401 of file implicit_system.h.

References libMesh::ImplicitSystem::_matrices.

 {
   return cast_int<unsigned int>(_matrices.size());
 }

static unsigned int libMesh::ReferenceCounter::n_objects ( )

inlinestaticinherited

Prints the number of outstanding (created, but not yet destroyed) objects.

Definition at line 85 of file reference_counter.h.

References libMesh::ReferenceCounter::_n_objects, libMesh::ReferenceCounter::disable_print_counter_info(), libMesh::ReferenceCounter::enable_print_counter_info(), libMesh::ReferenceCounter::increment_constructor_count(), libMesh::ReferenceCounter::increment_destructor_count(), and libMesh::Quality::name().

Referenced by libMesh::System::hide_output(), and libMesh::LibMeshInit::LibMeshInit().

86 { return _n_objects; }

libMesh::ReferenceCounter::_n_objects

static Threads::atomic< unsigned int > _n_objects

Definition: reference_counter.h:132

static unsigned int libMesh::ReferenceCounter::n_objects ( )

inlinestaticinherited

Prints the number of outstanding (created, but not yet destroyed) objects.

Definition at line 85 of file reference_counter.h.

References libMesh::ReferenceCounter::_n_objects, libMesh::ReferenceCounter::disable_print_counter_info(), libMesh::ReferenceCounter::enable_print_counter_info(), libMesh::ReferenceCounter::increment_constructor_count(), libMesh::ReferenceCounter::increment_destructor_count(), and libMesh::Quality::name().

Referenced by libMesh::System::hide_output(), and libMesh::LibMeshInit::LibMeshInit().

86 { return _n_objects; }

libMesh::ReferenceCounter::_n_objects

static Threads::atomic< unsigned int > _n_objects

Definition: reference_counter.h:132

processor_id_type libMesh::ParallelObject::n_processors ( ) const

inlineinherited

Returns: The number of processors in the group.

Definition at line 93 of file parallel_object.h.

References libMesh::ParallelObject::_communicator, and libMesh::Parallel::Communicator::size().

94 { return cast_int<processor_id_type>(_communicator.size()); }

libMesh::Parallel::Communicator::size

unsigned int size() const

Definition: parallel.h:722

libMesh::ParallelObject::_communicator

const Parallel::Communicator & _communicator

Definition: parallel_object.h:105

unsigned int libMesh::System::n_variable_groups ( ) const

inlineinherited

Returns: The number of VariableGroup variable groups in the system

Definition at line 2092 of file system.h.

References libMesh::System::_variable_groups.

Referenced by libMesh::System::add_variable(), libMesh::FEMSystem::assembly(), libMesh::System::get_info(), libMesh::System::init_data(), and libMesh::System::project_solution_on_reinit().

 {
   return cast_int<unsigned int>(_variable_groups.size());
 }

unsigned int libMesh::System::n_vars ( ) const

inlineinherited

Returns: The number of variables in the system

Definition at line 2084 of file system.h.

References libMesh::System::_variables.

 {
   return cast_int<unsigned int>(_variables.size());
 }

unsigned int libMesh::System::n_vectors ( ) const

inlineinherited

Returns: The number of vectors (in addition to the solution) handled by this system This is the size of the _vectors map

Definition at line 2212 of file system.h.

References libMesh::System::_vectors.

Referenced by libMesh::ExplicitSystem::add_system_rhs(), libMesh::System::compare(), libMesh::System::get_info(), and libMesh::System::project_solution_on_reinit().

 {
   return cast_int<unsigned int>(_vectors.size());
 }

const std::string & libMesh::System::name ( ) const

inlineinherited

Returns: The system name.

Definition at line 1996 of file system.h.

References libMesh::System::_sys_name.

Referenced by libMesh::System::compare(), libMesh::WrappedFunction< Output >::component(), libMesh::ContinuationSystem::ContinuationSystem(), DMlibMeshSetUpName_Private(), libMesh::ExactErrorEstimator::estimate_error(), libMesh::ExactSolution::ExactSolution(), libMesh::ExactErrorEstimator::find_squared_element_error(), libMesh::System::get_info(), libMesh::ImplicitSystem::get_linear_solver(), libMesh::NewtonSolver::init(), libMesh::PetscDiffSolver::init(), libMesh::TimeSolver::init_data(), libMesh::WrappedFunction< Output >::operator()(), libMesh::petsc_auto_fieldsplit(), libMesh::TimeSolver::reinit(), libMesh::System::set_adjoint_already_solved(), libMesh::FrequencySystem::solve(), libMesh::LinearImplicitSystem::solve(), libMesh::NonlinearImplicitSystem::solve(), libMesh::FrequencySystem::system_type(), libMesh::System::user_assembly(), libMesh::System::user_constrain(), libMesh::System::user_initialization(), libMesh::System::user_QOI(), and libMesh::System::user_QOI_derivative().

 {
   return _sys_name;
 }

unsigned int libMesh::System::number ( ) const

inlineinherited

Returns: The system number.

Definition at line 2004 of file system.h.

References libMesh::System::_sys_number.

Referenced by libMesh::ExactSolution::_compute_error(), libMesh::System::add_variable(), libMesh::Sy

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Protected Types

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Static Protected Attributes

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation