> Here 0 are the zero matrix, and * presents the matrix we want to > assemble. My question is how to give the finite elment system? I use > > FESystem<dim> fe; > > fe (FE_Q<dim>(degree), 1, > FE_Q<dim>(degree), 1, > FE_Q<dim>(degree), 1, > FE_Q<dim>(degree), dim) > > It is wrong because it is not supported by deal ii.
Indeed, it seems as if nobody ever needed the ability to concatenate four groups of elements. The easiest solution would probably be to add a constructor to FESystem that takes 4 groups of elements. If you take a look at fe_system.cc, you'll see that all you really need to do is extend a few helper functions from 3 to 4 arguments -- if you'd like to go this direction, we'd be happy to help you with it (including adding whatever code you write to the next release!). Best Wolfgang -- ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
