16 T
matrix_inner_product(
const Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic> &A,
const Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic> &B)
18 return (A.array() * B.array()).sum();
16 T
matrix_inner_product(
const Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic> &A,
const Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic> &B) {
…}
21 template <
typename T,
int rows,
int cols,
int option,
int maxRow,
int maxCol>
22 T
determinant(
const Eigen::Matrix<T, rows, cols, option, maxRow, maxCol> &mat)
24 assert(mat.rows() == mat.cols());
28 else if (mat.rows() == 2)
29 return mat(0, 0) * mat(1, 1) - mat(0, 1) * mat(1, 0);
30 else if (mat.rows() == 3)
31 return mat(0, 0) * (mat(1, 1) * mat(2, 2) - mat(1, 2) * mat(2, 1)) - mat(0, 1) * (mat(1, 0) * mat(2, 2) - mat(1, 2) * mat(2, 0)) + mat(0, 2) * (mat(1, 0) * mat(2, 1) - mat(1, 1) * mat(2, 0));
22 T
determinant(
const Eigen::Matrix<T, rows, cols, option, maxRow, maxCol> &mat) {
…}
38 Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic, 0, 3, 3>
inverse(
const Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic, 0, 3, 3> &mat)
40 assert(mat.rows() == mat.cols());
44 Eigen::Matrix<T, 1, 1> inv;
45 inv(0, 0) = T(1.) / mat(0);
49 else if (mat.rows() == 2)
51 Eigen::Matrix<T, 2, 2> inv;
53 inv(0, 0) = mat(1, 1) / det;
54 inv(0, 1) = -mat(0, 1) / det;
55 inv(1, 0) = -mat(1, 0) / det;
56 inv(1, 1) = mat(0, 0) / det;
60 else if (mat.rows() == 3)
62 Eigen::Matrix<T, 3, 3> inv;
64 inv(0, 0) = (-mat(1, 2) * mat(2, 1) + mat(1, 1) * mat(2, 2)) / det;
65 inv(0, 1) = (mat(0, 2) * mat(2, 1) - mat(0, 1) * mat(2, 2)) / det;
66 inv(0, 2) = (-mat(0, 2) * mat(1, 1) + mat(0, 1) * mat(1, 2)) / det;
67 inv(1, 0) = (mat(1, 2) * mat(2, 0) - mat(1, 0) * mat(2, 2)) / det;
68 inv(1, 1) = (-mat(0, 2) * mat(2, 0) + mat(0, 0) * mat(2, 2)) / det;
69 inv(1, 2) = (mat(0, 2) * mat(1, 0) - mat(0, 0) * mat(1, 2)) / det;
70 inv(2, 0) = (-mat(1, 1) * mat(2, 0) + mat(1, 0) * mat(2, 1)) / det;
71 inv(2, 1) = (mat(0, 1) * mat(2, 0) - mat(0, 0) * mat(2, 1)) / det;
72 inv(2, 2) = (-mat(0, 1) * mat(1, 0) + mat(0, 0) * mat(1, 1)) / det;
78 Eigen::Matrix<T, 1, 1> inv;
38 Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic, 0, 3, 3>
inverse(
const Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic, 0, 3, 3> &mat) {
…}
85 Eigen::SparseMatrix<double> I(rows, cols);
91 Eigen::VectorXd
flatten(
const Eigen::MatrixXd &X);
94 Eigen::MatrixXd
unflatten(
const Eigen::VectorXd &
x,
int dim);
101 Eigen::SparseMatrix<double>
lump_matrix(
const Eigen::SparseMatrix<double> &M);
111 const int reduced_size,
112 const std::vector<int> &removed_vars,
123 const Eigen::MatrixXd &in,
124 const Eigen::VectorXi &in_to_out,
126 const int block_size = 1);
135 const Eigen::MatrixXd &out,
136 const Eigen::VectorXi &in_to_out,
138 const int block_size = 1);
145 const Eigen::MatrixXi &in,
146 const Eigen::VectorXi &index_mapping);
148 template <
typename DstMat,
typename SrcMat>
154 dst.resize(dst.rows(), src.cols());
155 assert(dst.cols() == src.cols());
156 dst.conservativeResize(dst.rows() + src.rows(), dst.cols());
157 dst.bottomRows(src.rows()) = src;
160 template <
typename DstMat>
163 assert(dst.cols() > 0);
164 if (n_zero_rows == 0)
166 dst.conservativeResize(dst.rows() + n_zero_rows, dst.cols());
167 dst.bottomRows(n_zero_rows).setZero();
172 const Eigen::MatrixXd &A,
173 const Eigen::MatrixXd &b,
174 const std::vector<int> &local_to_global,
176 Eigen::MatrixXd &bout);
180 const Eigen::MatrixXd &A,
181 const Eigen::MatrixXd &b,
182 const std::vector<int> &local_to_global,
184 Eigen::MatrixXd &bout);
188 const std::vector<long> &shape,
189 const std::vector<int> &rows,
190 const std::vector<int> &cols,
191 const std::vector<double> &
vals,
192 const Eigen::MatrixXd &b,
193 const std::vector<int> &local_to_global,
195 Eigen::MatrixXd &bout);
199 const std::vector<long> &shape,
200 const std::vector<int> &rows,
201 const std::vector<int> &cols,
202 const std::vector<double> &
vals,
203 const Eigen::MatrixXd &b,
204 const std::vector<int> &local_to_global,
206 Eigen::MatrixXd &bout);
ElementAssemblyValues vals
T matrix_inner_product(const Eigen::Matrix< T, Eigen::Dynamic, Eigen::Dynamic > &A, const Eigen::Matrix< T, Eigen::Dynamic, Eigen::Dynamic > &B)
Eigen::SparseMatrix< double > lump_matrix(const Eigen::SparseMatrix< double > &M)
Lump each row of a matrix into the diagonal.
void show_matrix_stats(const Eigen::MatrixXd &M)
Eigen::MatrixXd reorder_matrix(const Eigen::MatrixXd &in, const Eigen::VectorXi &in_to_out, int out_blocks=-1, const int block_size=1)
Reorder row blocks in a matrix.
Eigen::SparseMatrix< double > sparse_identity(int rows, int cols)
void vector2matrix(const Eigen::VectorXd &vec, Eigen::MatrixXd &mat)
void scatter_matrix(const int n_dofs, const int dim, const Eigen::MatrixXd &A, const Eigen::MatrixXd &b, const std::vector< int > &local_to_global, StiffnessMatrix &Aout, Eigen::MatrixXd &bout)
Eigen::Matrix< T, Eigen::Dynamic, Eigen::Dynamic, 0, 3, 3 > inverse(const Eigen::Matrix< T, Eigen::Dynamic, Eigen::Dynamic, 0, 3, 3 > &mat)
Eigen::MatrixXd unreorder_matrix(const Eigen::MatrixXd &out, const Eigen::VectorXi &in_to_out, int in_blocks=-1, const int block_size=1)
Undo the reordering of row blocks in a matrix.
Eigen::MatrixXd unflatten(const Eigen::VectorXd &x, int dim)
Unflatten rowwises, so every dim elements in x become a row.
void append_rows_of_zeros(DstMat &dst, const size_t n_zero_rows)
void append_rows(DstMat &dst, const SrcMat &src)
void scatter_matrix_col(const int n_dofs, const int dim, const Eigen::MatrixXd &A, const Eigen::MatrixXd &b, const std::vector< int > &local_to_global, StiffnessMatrix &Aout, Eigen::MatrixXd &bout)
Eigen::VectorXd flatten(const Eigen::MatrixXd &X)
Flatten rowwises.
Eigen::MatrixXi map_index_matrix(const Eigen::MatrixXi &in, const Eigen::VectorXi &index_mapping)
Map the entrys of an index matrix to new indices.
void full_to_reduced_matrix(const int full_size, const int reduced_size, const std::vector< int > &removed_vars, const StiffnessMatrix &full, StiffnessMatrix &reduced)
Map a full size matrix to a reduced one by dropping rows and columns.
T determinant(const Eigen::Matrix< T, rows, cols, option, maxRow, maxCol > &mat)
Eigen::SparseMatrix< double, Eigen::ColMajor > StiffnessMatrix