Re: [deal.II] Meshworker in dim=2, spacedim = 3

Many thanks. Andrea > On Apr 17, 2017, at 10:02 AM, Daniel Arndt > wrote: > > Andrea, > > This should be fixed upstream now. Have a look at PR #4265 [1]. > > Best, > Daniel > > [1] https://github.com/dealii/dealii/pull/4265 > > -- > The deal.II project is located at http://www.dealii.or

Re: [deal.II] Meshworker in dim=2, spacedim = 3

On 04/14/2017 09:47 AM, Bonito Andrea wrote: * * template voidintegration_loop(ITERATOR begin, typenameidentity::type end, DoFInfo &dof_info, IntegrationInfoBox &box, constLocalIntegrator &inte

[deal.II] Meshworker in dim=2, spacedim = 3

All: I am trying to have the meshworker going for a problem set on a hyper surface. As far as I can tell, it seems that this feature is already incorporated but cannot understand the implementation of *%* */Applications/deal.II.app/Contents/Resources/include/deal.II/meshworker/loop

[deal.II] MeshWorker clarifications

Dear all, I am trying to port the previous linear elastic code to MeshWorker, but I have some questions about this class. First of all, I have read this question , and I read that Me

[deal.II] Meshworker

we choose dim=2 template void Step12::integrate_face_term (DoFInfo& dinfo1, DoFInfo& dinfo2, CellInfo& info1, CellInfo& info2) { const FEValuesBase& fe_v = info1.fe_values(); const FEValuesBase& fe_v_neighbor = info2.fe_values();

[deal.II] meshworker getting gradient values

dim=2 template void Step12::integrate_cell_term (DoFInfo& dinfo, CellInfo& info) { const FEValuesBase& fe_v = info.fe_values(); const std::vector>& grad = info.gradients[0][0]; } When I call grad[0][0] it gives segmentation fault. Could any one tell how to access the gradie

Re: [deal.II] Meshworker

Sudarshan, there is a Triangulation::cell_iterator and a DoFHandler::cell_iterator and dinfo1.cell returns the first one while you are asking for the first one. Maybe have a look at [1] for the difference between the two. What exactly do you need the cell_iterator for? Best, Daniel [1] https:/

Re: [deal.II] Meshworker

I get the following error if I call cell iterator from dinfo1 ^ In file included from /Users/Sudarshan/Dropbox/ALWDG_code/2d_scalar_unsteady_legendre/dg.cc:26: In file included from /Applications/deal.II.app/Contents/Resources/include/deal.II/grid/tria.h:3631:

Re: [deal.II] Meshworker

Sudarshan, > In the following function I need to get the cell iterator of the > corresponding cell. > > void Step12::integrate_face_term (DoFInfo& dinfo1, DoFInfo& dinfo2, > { CellInfo& > info1, CellInfo& info2) > > typename DoFHandler::ce

[deal.II] Meshworker

In the following function I need to get the cell iterator of the corresponding cell. void Step12::integrate_face_term (DoFInfo& dinfo1, DoFInfo& dinfo2, { CellInfo& info1, CellInfo& info2) typename DoFHandler::cell_iterator cell=dinfo1