Dear all, I am having an issue that I would like to share and hopefullyI can get some feedback. In my code, I first do an assemble and I store the operators in the corresponding variables. Then, I have an* update()* function, where I use the pre-anssambled operators and change some values to compute different frequencies to build the *system_matrix*. Until here everything works.
I need to export the operators to perform some computations in a different environment. I used SparseMatrix::print_formatted() <https://www.dealii.org/current/doxygen/deal.II/classSparseMatrix.html#a6456ebe393c32b076eae0b5c9811e0c8> and verified that the exported values are the same as the ones in my dealii code. I did a sanity check, where I take the exported operators and solve the same system in in another environment (which is also validated). However, I get wrong results when using the exported operators. The single function in my dealii code that I am calling in the end of my *update()* function is *constraints.distribute_local_to_global(cell_matrix,* * cell_rhs,* * local_dof_indices,* * system_matrix,* * system_rhs);* *My question is, *if thee *distribute_local_to_global function * is performing something that I am missing, which can be changing the system_matrix so that I need to do something else with the exported operators when I create the *system_matrix *in the sanity check at the other environment. Thank you very much -- 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/29105e0c-ded0-4a65-82e6-99e15fdb8a74n%40googlegroups.com.
