[deal.II] How to add addtional point to enrichment functions (FE_Enriched)

2017-10-31 Thread Edith Sotelo
Hello, I would like to implement enrichment functions of this type: f(x-xo) where xo is a user defined point. It could be the center of the cell or the dof coordinate. I would like a suggestion about how to implement this, because I do not see how to pass it in the constructor since it is not a

Re: [deal.II] Re: VectorTools::point_value() throw exception point not locally owned

2017-10-31 Thread Wolfgang Bangerth
On 10/31/2017 12:12 PM, Bruno Turcksin wrote: The above code seems to output the "ERROR" statement in this code, which means the GridTools::find_active_cell_around_point() will return a different cell. It is understandable that to return the cell that a vertex belongs to is a

[deal.II] Re: VectorTools::point_value() throw exception point not locally owned

2017-10-31 Thread Yiyang Zhang
Hello Bruno, I have opened a github issue, https://github.com/dealii/dealii/issues/5367. Also attached a test code in that webpage. Best, Yiyang -- 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 ---

[deal.II] Re: VectorTools::point_value() throw exception point not locally owned

2017-10-31 Thread Bruno Turcksin
Yiyang, On Tuesday, October 31, 2017 at 10:56:00 AM UTC-4, Yiyang Zhang wrote: > The above code seems to output the "ERROR" statement in this code, which > means the GridTools::find_active_cell_around_point() will return a > different cell. > It is understandable that to return the cell that a

Re: [deal.II] deal.II on Mac OS Sierra

2017-10-31 Thread Pawan Takhar
Thanks for compiling the newer version of .dmg file, Timo. I tried compiling 8.5.1 .tar.gz file, but make failed on that as well with newer versions of MacOS. I look forward to your .dmg file. Regards, Pawan -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum opt

[deal.II] VectorTools::point_value() throw exception point not locally owned

2017-10-31 Thread Yiyang Zhang
template void SharedTria::test_local_vertex() { for(auto cell = this->tria.begin_active(); cell != this->tria.end(); ++ cell){ if(cell->is_locally_owned()){ for(auto i = 0; i::vertices_per_cell; ++i ){ auto located_cell = GridTools:: find_ac

Re: [deal.II] Re: Read in tethex grid

2017-10-31 Thread Bruno Turcksin
Felix, 2017-10-31 7:07 GMT-04:00 Felix : > > However, what I would try is to use gmsh to directly produce a mesh of >> quadrangles. >> > > Do you mean that I should avoid this "iso.msh" (consist of triangles) and > directly construct the isosurface consisting of quadrangle and save it as > .msh?

Re: [deal.II] CellDataStorage with mesh refinement

2017-10-31 Thread 'Frederik S.' via deal.II User Group
Hi Timo, Thank you for your answer, this sounds really good! I tried to use the class, but I still have some Problems. In the example you are just transferring an integer variable. Is it also possible to transfer a whole CellDataStorage element (or at least some packed vector with ist elements)

Re: [deal.II] CellDataStorage with mesh refinement

2017-10-31 Thread 'Frederik S.' via deal.II User Group
Hi Jean-Paul, Thank you for your answer! I have a truly discontinuous case, the case you mentioned with internal yield variables. Do you think it will be easier to use the method Timo mentioned "from scratch" or will it be easier to copy the source files of ContinuousQuadratureDataTransfer clas

[deal.II] Re: Read in tethex grid

2017-10-31 Thread Felix
> However, what I would try is to use gmsh to directly produce a mesh of > quadrangles. > Do you mean that I should avoid this "iso.msh" (consist of triangles) and directly construct the isosurface consisting of quadrangle and save it as .msh? -- The deal.II project is located at http://ww