Hello I'm solving a transient nonlinear implicit system in parallel. In each time step, a list of interacting element pairs is built, based on which the residual and jacobian are constructed. These element interactions spread over arbitrary processors, e.g. element 0 (processor 0) interacts with element 24 (processor 1), and in the next step the interaction happens between element 0 and 30 (processors 0 and 2), and so on. In other words, any node may interact with any other node on any processor at any time. The interaction pairs are recalculated after every nonlinear iteration step.
How can I implement this in libMesh? It works fine in serial mode, where I can just fill in ANY entries in the residual vector and jacobian matrix. The program fails to converge however in parallel mode, and I suppose this is because entries outside the (initially set) sparsity pattern are set, which then don't get communicated, because libMesh is not expecting arbitrary long range interactions. Will I need to deal with updated send_lists? How does this work? Thanks for your advice, Roman ------------------------------------------------------------------------------ Enable your software for Intel(R) Active Management Technology to meet the growing manageability and security demands of your customers. Businesses are taking advantage of Intel(R) vPro (TM) technology - will your software be a part of the solution? Download the Intel(R) Manageability Checker today! http://p.sf.net/sfu/intel-dev2devmar _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
