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

2017-04-17 Thread Daniel Arndt
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.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You

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

2017-04-15 Thread Bonito Andrea
done, thanks. On Friday, April 14, 2017 at 4:39:45 PM UTC-5, Wolfgang Bangerth wrote: > > On 04/14/2017 12:55 PM, Bonito Andrea wrote: > > Unfortunately, I am using the dmg version and do not have access to the > > source files... > > > > I might have to convince Luca to try it. > >

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

2017-04-14 Thread Wolfgang Bangerth
On 04/14/2017 12:55 PM, Bonito Andrea wrote: Unfortunately, I am using the dmg version and do not have access to the source files... I might have to convince Luca to try it. Alternatively, can you create a small program that shows the compiler error? It doesn't have to do anything useful,

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

2017-04-14 Thread Bonito Andrea
Unfortunately, I am using the dmg version and do not have access to the source files... I might have to convince Luca to try it. Andrea On Friday, April 14, 2017 at 10:47:22 AM UTC-5, Bonito Andrea wrote: > > All: > > I am trying to have the meshworker going for a problem set on a hyper >

[deal.II] Re: MeshWorker clarifications

2017-03-10 Thread Franco Milicchio
Just one more question about MeshWorker. I am trying to figure how I could assemble simultaneously more than one form. I've found that step 39 uses two different assemblers. My guess is that I could use a MatrixSimple assembler for bilinear forms, and ResidualSimple for linear forms. Is this

[deal.II] Re: MeshWorker clarifications

2017-03-10 Thread Franco Milicchio
It was all my mistake in copy/pasting the weak form. The two matrices now are the same. Sorry for the noise. -- 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

[deal.II] Re: MeshWorker clarifications

2017-03-09 Thread Daniel Arndt
Franco, on the first glance your code looks good. Can you reduce the size of your problem to a single cell? Is there still a difference then? Best, Daniel Am Mittwoch, 8. März 2017 13:04:01 UTC+1 schrieb Franco Milicchio: > > I have made available the snippet for this example here >

[deal.II] Re: MeshWorker clarifications

2017-03-07 Thread Franco Milicchio
Thanks for the answer, Daniel. On Monday, March 6, 2017 at 7:49:53 PM UTC+1, Daniel Arndt wrote: So my first question is, should I avoid using this class and implement >> parallel loops by hand (via TBB or other means)? >> > > "amandus"[1] is in fact based on MeshWorker. If you are trying to

[deal.II] Re: MeshWorker clarifications

2017-03-06 Thread Daniel Arndt
Franco, First of all, I have read this question > , > > and I read that MeshWorker is sort of abandoned: > > the issue with MeshWorker is that when it was implemented, we didn't

Re: [deal.II] Re: Meshworker

2016-11-22 Thread sudarshan kumar
Yes the FEValues defined for each dell dinfo1.cell and dinfo2.cell provides all the interior information of the cells( however I did not succeed in the accuracy test so far, but it works ). In fact, apart from the interface values I need solution values at the interior of each cell,

Re: [deal.II] Re: Meshworker

2016-11-22 Thread Wolfgang Bangerth
On 11/22/2016 06:13 AM, sudarshan kumar wrote: Thanks a lot, yes it works. In fact I am working on time dependent problem, and so in my method I need an adjacent interior value when I am at the face. Thanks a lot. But don't the FEValues objects you get for the face provide these values? I

[deal.II] Re: Meshworker

2016-11-20 Thread Daniel Arndt
Sudarshan, When I am at the face I need the solution at the interior quadrature > points as well. However from dinfo1.indices we could get the the dof > indices for the cell1. But is it possible to get the shape value at the > interior quadrature points ? > You can of course work

[deal.II] Re: Meshworker

2016-11-18 Thread Sudarshan Kumar
Dear Daniel, When I am at the face I need the solution at the interior quadrature points as well. However from dinfo1.indices we could get the the dof indices for the cell1. But is it possible to get the shape value at the interior quadrature points ? Thanks a lot On Friday,