Thank you guys for sharing all the info. >> Try in debug mode - I suspect you won't just see convergence failure, >> you'll see an assertion thrown when you try to access nonexistent data >> from a ghosted vector. > >Yes - you should. [...] > BTW: If you want to do this (access off-processor solution entries) _now_ > just use solution->localize() and pass in a std::vector...
This what I'm currently doing - localizing whatever information I need to a std::vector. Thus, I'm not getting an assertion failure in debug mode. > But that's only half of what Roman is asking for... adding to the send_list > will ensure that off-processor entries in the solution vector will be ghosted > onto the correct processors in current_local_solution. This is necessary for > assembly of residuals / jacobian entries for on-processor dofs... but doesn't > address the second part which is: new entries in the Jacobian. > > To achieve this you will need to modify the SparsityPattern objects. It's > been a while since I've been in that system... and don't really know where > you would start. You could think about adding some sort of a callback to get > extra dof pairs though... and see if you can get them put into the sparsity > pattern. Thanks for the pointer! If I have time, I'll see where I can go with this. Roman ------------------------------------------------------------------------------ Create and publish websites with WebMatrix Use the most popular FREE web apps or write code yourself; WebMatrix provides all the features you need to develop and publish your website. http://p.sf.net/sfu/ms-webmatrix-sf _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
