67 int const idx_bhe_unknowns,
68 Eigen::MatrixBase<SingleUnknownMatrixType>
const& matBHE_loc_R,
69 Eigen::MatrixBase<RMatrixType>& R_matrix,
70 Eigen::MatrixBase<RPiSMatrixType>& R_pi_s_matrix,
71 Eigen::MatrixBase<RSMatrixType>& R_s_matrix)
73 switch (idx_bhe_unknowns)
76 R_matrix.block(0, 2 * NPoints, NPoints, NPoints) +=
78 R_matrix.block(2 * NPoints, 0, NPoints, NPoints) +=
81 R_matrix.block(0, 0, NPoints, NPoints) +=
83 R_matrix.block(2 * NPoints, 2 * NPoints, NPoints, NPoints) +=
87 R_matrix.block(NPoints, 3 * NPoints, NPoints, NPoints) +=
89 R_matrix.block(3 * NPoints, NPoints, NPoints, NPoints) +=
92 R_matrix.block(NPoints, NPoints, NPoints, NPoints) +=
94 R_matrix.block(3 * NPoints, 3 * NPoints, NPoints, NPoints) +=
98 R_matrix.block(2 * NPoints, 3 * NPoints, NPoints, NPoints) +=
100 R_matrix.block(3 * NPoints, 2 * NPoints, NPoints, NPoints) +=
103 R_matrix.block(2 * NPoints, 2 * NPoints, NPoints, NPoints) +=
106 R_matrix.block(3 * NPoints, 3 * NPoints, NPoints, NPoints) +=
111 R_s_matrix.template block<NPoints, NPoints>(0, 0).noalias() +=
114 R_pi_s_matrix.block(2 * NPoints, 0, NPoints, NPoints) +=
116 R_pi_s_matrix.block(3 * NPoints, 0, NPoints, NPoints) +=
118 R_matrix.block(2 * NPoints, 2 * NPoints, NPoints, NPoints) +=
120 R_matrix.block(3 * NPoints, 3 * NPoints, NPoints, NPoints) +=
125 "Error!!! In the function BHE_1U::assembleRMatrices, "
126 "the index of bhe unknowns is out of range! ");