Re: [deal.II] How to use ManifoldMapping

I know that there could be cancellation errors and other magic stuff… :) … I was joking Btw, I checked your code carefully and I do not see errors in the implementation. I re-checked it carefully together with Luca and we do not see errors in the implementation. Now, I would like to isolate t

Re: [deal.II] How to use ManifoldMapping

Dear Praveen, I noticed that the divergence theorem is safe. :) Mapping Manifold: \int_{\partial K} v \cdot n = -8.51746e-21 \int_{K} \div v = -3.74547e-21 Mapping Q: \int_{\partial K} v \cdot n = -5.70692e-22 \int_{K} \div v = 3.09366e-21 BTW, as y

Re: [deal.II] Re: Segmentation fault (core dumped) error

Dear Bruno, I did what you have recommended and the problem is solved. Thanks. Ehsan On Friday, May 20, 2016 at 6:01:14 PM UTC+2, Bruno Turcksin wrote: > > 2016-05-20 11:55 GMT-04:00 Bruno Turcksin >: > > looking more at the code I don't think that this is your problem. Is > > your code paralle

[deal.II] Re: Post-processing of solution (computing stress distribution)

David, if you look at step-33 On Monday, May 23, 2016 at 6:59:27 AM UTC-4, David wrote: > > Thanks for your help. But I have already looked at step 29 -33, the > question is I don't quite understand it. Just look at this small piece of > code: > > Tensor<2,dim> >

[deal.II] Re: Post-processing of solution (computing stress distribution)

Hi, Daniel: Thanks for your help. But I have already looked at step 29 -33, the question is I don't quite understand it. Just look at this small piece of code: Tensor<2,dim> grad_u; for (unsigned int d=0; dhttp://www

[deal.II] Re: Post-processing of solution (computing stress distribution)

Hey David, Yes, this won't work since FEValues need a cell to reinit to get mapping > data. I looded at the class DataPostprocessor as you mentioned, but don't > quite understand the parameters passed to compute_derived_quantities_vector > method. I don't understand what exactly is the paramete

Re: [deal.II] Re: Hanging node constraints generate new nonzero locations in parallel.

Dear Martin, It is clear now. Thank you very much. Best regards, Ce Qin 2016-05-23 16:14 GMT+08:00 Martin Kronbichler : > Dear Ce, > > It must have been coincidence that it worked. I guess it depends on the > data structures used by the PETSc matrices. If you were using deal.II > matrices, I g

Re: [deal.II] Re: Hanging node constraints generate new nonzero locations in parallel.

Dear Ce, It must have been coincidence that it worked. I guess it depends on the data structures used by the PETSc matrices. If you were using deal.II matrices, I guess that the serial case should have failed as soon as constraints were presented. Best, Martin On 23.05.2016 10:12, Ce Qin wr

Re: [deal.II] Re: Hanging node constraints generate new nonzero locations in parallel.

Dear Martin, Yes, it works now. Could you please explain why it works in serial or with a small number of processes? Best regards, Ce Qin 2016-05-23 15:27 GMT+08:00 Martin Kronbichler : > Dear Ce, > > I faced a strange problem. With relatively few cells and non-conforming >> triangulation, the

[deal.II] Re: Hanging node constraints generate new nonzero locations in parallel.

Dear Ce, I faced a strange problem. With relatively few cells and non-conforming > triangulation, the assembling process generates new nonzero locations in > parallel with relatively large amount of processes . It works well in > serial or with a small number of processes or with conforming >

[deal.II] Re: "Make test" failure following Trilinos installation with Mesquite

FYI, I created a FAQ page for using deal.II with Spack: https://github.com/dealii/dealii/wiki/deal.II-in-Spack However given your compiler issue, i think in this case it will be easier to use Luca's prebuilt deal.II suite. Regards, Denis. -- The deal.II project is located at http://www.dealii