Re: [deal.II] Re: where and how should I use the MeshWorker class? how should I do to implement DG without MeshWorker?

2016-08-29 Thread yoczhang
Dear Daniel, Thank you, I have looked at it, it may take some time to learn. It's a good news for me. 在 2016年8月29日星期一 UTC+8下午5:23:27,Daniel Arndt写道: > > Yongchao, > > If you are still interested in using MeshWorker, there are some more > examples (and in particular Stokes) at > https://bitbucket

Re: [deal.II] Re: How to apply a point source condition (point boundary condition) in Discontinuous Galerkin methods

2016-08-29 Thread Jiaqi ZHANG
Thanks Bruno, I should have searched a little bit more. Best, Jiaqi 2016-08-29 11:00 GMT-04:00 Bruno Turcksin : > Jiaqi, > > to create a point source, you can use VectorTools::create_point_source_vector > (http://dealii.org/developer/doxygen/deal.II/namespaceVectorTools.html# > ad03b858b1a3b5900

[deal.II] Re: How to apply a point source condition (point boundary condition) in Discontinuous Galerkin methods

2016-08-29 Thread Bruno Turcksin
Jiaqi, to create a point source, you can use VectorTools::create_point_source_vector (http://dealii.org/developer/doxygen/deal.II/namespaceVectorTools.html#ad03b858b1a3b59003a76f6224e67efc7) Best, Bruno On Monday, August 29, 2016 at 10:39:27 AM UTC-4, Jiaqi Zhang wrote: > > Dear All, > > I wa

Re: [deal.II] Re: assembling individual block step 43

2016-08-29 Thread Daniel Arndt
> > > There actually *is* a function that takes two DoFHandler objects, but it > doesn't take a constraint matrix or the last boolean argument. > [...] > You can of course see what variants of the make_sparsity_pattern function > there are by looking here: >https://www.dealii.org/8.4.1/do

[deal.II] How to apply a point source condition (point boundary condition) in Discontinuous Galerkin methods

2016-08-29 Thread Jiaqi Zhang
Dear All, I was trying to use DG to solve a Hamilton-Jacobi equation, and the problem has a point source condition for example, the center of a unit square pointCenter = 0. First I thought of using ZeroSourcePoint.add_line(cell->vertex_dof_index(v,0)); but later found that DG elements aren't d

Re: [deal.II] Re: assembling individual block step 43

2016-08-29 Thread Wolfgang Bangerth
DoFTools::make_sparsity_pattern (saturation_dof_handler, darcy_dof_handler, dsp, saturation_constraints, false); ssemble_matrix_H.reinit(dsp); dealii is complaining and is saying /home/franckm/apps/candi/deal.II-toolchain/deal.II-v

Re: [deal.II] Re: where and how should I use the MeshWorker class? how should I do to implement DG without MeshWorker?

2016-08-29 Thread Daniel Arndt
Yongchao, If you are still interested in using MeshWorker, there are some more examples (and in particular Stokes) at https://bitbucket.org/guidokanschat/amandus Best, Daniel Am Montag, 29. August 2016 11:15:22 UTC+2 schrieb yocz...@gmail.com: > > Dear Wolfgang , > Thanks for your advice, I ha

Re: [deal.II] Re: where and how should I use the MeshWorker class? how should I do to implement DG without MeshWorker?

2016-08-29 Thread yoczhang
Dear Wolfgang , Thanks for your advice, I have read the 6.2.1 version of step-12, and I know how to deal DG without using MeshWorker. What's more there is the 6.2.1 version of step-12 without using MeshWorker, if someone needs. https://www.dealii.org/6.2.1/doxygen/deal.II/step_12.html 在 2016年8