On 3/19/19 10:39 PM, Yidong ZHAO wrote: > > While using a new kind of shape function, degrees of freedom on each cell > will > couple to the degrees of freedom on other cells connected to the current one > by a common face or common vertex. > > From the description of make_flux_sparsity_pattern(), it can solve the > problem > when two elements are connected by a common face, as depicted in > elements_connected_by_common_line.png. > > But if it's the case as elements_connected_by_common_vertex.png, is there any > way to make_vertex_sparsity_pattern() ?
Yidong, there is currently no function in deal.II that does this. But it should not be too hard for you to write one (in which case I would suggest you contribute it to deal.II! :) ). In essence, the functions that build a sparsity pattern _simulate_ what will happen during assembly. They do not compute individual matrix elements, but investigate what entries will be written to, and those are the ones for which elements of the matrix are allocated in the sparsity pattern. You will need to do something similar for your case. It would probably be useful for you to take a look at how the existing functions are implemented. Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- 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. For more options, visit https://groups.google.com/d/optout.