Hello dear group, I have one question which is about how to build global stiffness matrix with cell stiffness matrix.
In our project,the density is a vector (actually we need a mu vector but in the code of deallii mu is a scalar)with different values,I would like to use the element of density vector to build cell matrix as written in the link below: https://www.dealii.org/current/doxygen/code-gallery/Quasi_static_Finite_strain_Compressible_Elasticity/cook_membrane.cc (starting from line 1782) And finally the cell stiffness matrix is used to build global stiffness matrix,but I don' know how the global stiffness matrix is bulit with cell stiffness matrix. https://www.dealii.org/current/doxygen/code-gallery/Quasi_static_Finite_strain_Compressible_Elasticity/cook_membrane.cc (starting from line 1599) //>>>>>>>>>>>>>>>>>>>>>>>>>> constraints.distribute_local_to_global( data.cell_matrix, data.cell_rhs, data.local_dof_indices, tangent_matrix, system_rhs); } //<<<<<<<<<<<<<<<<<<<<<<<<<<<< I have a idea that I use mu[0] to obtain the first cell_matrix and use mu[1] to get the second cell stiffness matrix ,step by step then I got final mu[m] to obtain the last cell stiffness matrix and in the end I use some function like distribute_local_to_global to construct the global stiffness matrix,and system rhs and even the solutions. I checked the code in AffineConstraints<number>::distribute_local_to_global <https://www.dealii.org/current/doxygen/deal.II/classAffineConstraints.html#a373fbdacd8c486e675b8d2bff8943192>. ( https://www.dealii.org/current/doxygen/deal.II/affine__constraints_8h_source.html ) But the hint information guides me to read the cm.templates.h file. ewcfp I was a little confused ,I did not find the file from dealii lib. could anyone provide any information or hint? Thanks in advance! Best regards Lance Could I -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/e3f16c37-7670-40ee-b0d2-821fdcb7b307n%40googlegroups.com.
