Dear Wolfgang, Jean-Paul and all I followed your suggestions and the code run. Thank for your helps.
But I still a concern that the function class (the Coefficient class in step 6) will be called in many time in the main code. Do we have another approach, for example store the heterogeneous permeability belonged to each cell, to avoid wasted computational cost? Best regard Lam DANG University of Orleans, France On Wednesday, January 25, 2017 at 2:47:18 PM UTC+1, Wolfgang Bangerth wrote: > > On 01/25/2017 05:37 AM, Jean-Paul Pelteret wrote: > > Dear Lam, > > > > There are a number of ways that you can do this. One approach would be > to > > define a class that derives from the function class > > <https://www.dealii.org/developer/doxygen/deal.II/classFunction.html>; > this is > > alluded to in step-8 > > < > https://www.dealii.org/developer/doxygen/deal.II/step_8.html#ElasticProblemassemble_system>, > > > > but here the parameters are considered constant. For this approach you > would > > need to describe the material parameters by position, rather than > material_id. > > Indeed, this is what step-6 does. (For the Laplace equation, but if you > see it > happen there, you will know how it works for the elasticity equation as > well.) > > Best > W. > > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@colostate.edu > <javascript:> > 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.
