PolyFEM
Loading...
Searching...
No Matches
polyfem::solver::LinearFilter Class Reference

Maps to average of neighboring. More...

#include <Parametrizations.hpp>

Inheritance diagram for polyfem::solver::LinearFilter:
[legend]
Collaboration diagram for polyfem::solver::LinearFilter:
[legend]

Public Member Functions

 LinearFilter (const mesh::Mesh &mesh, const double radius)
 
int size (const int x_size) const override
 Compute DOF of y given DOF of x.
 
int inverse_size (int y_size) const override
 Compute DOF of x given DOF of y.
 
Eigen::VectorXd eval (const Eigen::VectorXd &x) const override
 Eval y = f(x).
 
Eigen::VectorXd inverse_eval (const Eigen::VectorXd &y) const override
 Eval x = f^-1 (y).
 
Eigen::VectorXd apply_jacobian (const Eigen::VectorXd &grad, const Eigen::VectorXd &x) const override
 Apply jacobian for chain rule.
 
- Public Member Functions inherited from polyfem::solver::Parametrization
virtual ~Parametrization ()=default
 

Private Attributes

Eigen::SparseMatrix< double > tt_radius_adjacency
 
Eigen::VectorXd tt_radius_adjacency_row_sum
 

Detailed Description

Maps to average of neighboring.

Definition at line 159 of file Parametrizations.hpp.

Constructor & Destructor Documentation

◆ LinearFilter()

polyfem::solver::LinearFilter::LinearFilter ( const mesh::Mesh mesh,
const double  radius 
)

Definition at line 597 of file Parametrizations.cpp.

References polyfem::mesh::Mesh::bounding_box(), polyfem::mesh::Mesh::cell_barycenters(), polyfem::mesh::Mesh::face_barycenters(), polyfem::mesh::Mesh::is_volume(), tt_radius_adjacency, and tt_radius_adjacency_row_sum.

Here is the call graph for this function:

Member Function Documentation

◆ apply_jacobian()

Eigen::VectorXd polyfem::solver::LinearFilter::apply_jacobian ( const Eigen::VectorXd &  grad_full,
const Eigen::VectorXd &  x 
) const
overridevirtual

Apply jacobian for chain rule.

Let g(y) = g(f(x)). Given ∂g/∂y, compute ∂g/∂x = ∂g/∂y * ∂y/∂x.

Parameters
[in]grad_full∂g/∂y.
[in]xWhere ∂g/∂x is evaluated.
Returns
∂g/∂x.

Implements polyfem::solver::Parametrization.

Definition at line 661 of file Parametrizations.cpp.

References tt_radius_adjacency, tt_radius_adjacency_row_sum, and x.

◆ eval()

Eigen::VectorXd polyfem::solver::LinearFilter::eval ( const Eigen::VectorXd &  x) const
overridevirtual

Eval y = f(x).

Parameters
[in]xx.
Returns
y.

Implements polyfem::solver::Parametrization.

Definition at line 649 of file Parametrizations.cpp.

References tt_radius_adjacency, tt_radius_adjacency_row_sum, and x.

◆ inverse_eval()

Eigen::VectorXd polyfem::solver::LinearFilter::inverse_eval ( const Eigen::VectorXd &  y) const
overridevirtual

Eval x = f^-1 (y).

This is not a strict inverse in mathematical sense, one may choose "reasonable" x even if f is not one-to-one.

Parameters
[in]yy.
Returns
x.
Exceptions
std::runtime_errorThrow if inverse not implemented or is impossible.

Implements polyfem::solver::Parametrization.

Definition at line 655 of file Parametrizations.cpp.

References y.

◆ inverse_size()

int polyfem::solver::LinearFilter::inverse_size ( int  y_size) const
overridevirtual

Compute DOF of x given DOF of y.

Parameters
[in]y_sizeDOF of y.
Returns
DOF of x.
Exceptions
std::runtime_errorThrow if inverse not impossible or y_size is invalid.

Implements polyfem::solver::Parametrization.

Definition at line 634 of file Parametrizations.cpp.

References polyfem::log_and_throw_adjoint_error(), and tt_radius_adjacency.

Here is the call graph for this function:

◆ size()

int polyfem::solver::LinearFilter::size ( const int  x_size) const
overridevirtual

Compute DOF of y given DOF of x.

Parameters
[in]x_sizeThe DOF of x.
Returns
DOF of y.

Implements polyfem::solver::Parametrization.

Definition at line 644 of file Parametrizations.cpp.

Member Data Documentation

◆ tt_radius_adjacency

Eigen::SparseMatrix<double> polyfem::solver::LinearFilter::tt_radius_adjacency
private

Definition at line 171 of file Parametrizations.hpp.

Referenced by apply_jacobian(), eval(), inverse_size(), and LinearFilter().

◆ tt_radius_adjacency_row_sum

Eigen::VectorXd polyfem::solver::LinearFilter::tt_radius_adjacency_row_sum
private

Definition at line 172 of file Parametrizations.hpp.

Referenced by apply_jacobian(), eval(), and LinearFilter().


The documentation for this class was generated from the following files: