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

2017-11-01 Thread Edith Sotelo 
Hi Denis,
Thank you for your answer but I still need to further understand some things :  
* I think I need to pass the point (xo) during assembly since it will different 
for different cells ( dof coordinate or center of the cell), and I do not see 
how I can pass the point when declaring the function (ex. Enrichment  
(param1,param2, etc...) ) to create the  FE_Enriched constructor, maybe I am 
missing something here.

*  The thing is that  during assembly the functions are already encapsulated  
in the “enrichments” vector and if I even  create a member function like 
Enrichment::get_point(Point & ) how do I apply it to the functions 
that are already in the “enrichments” vector ? Maybe I can make  the 
“enrichments” vector  public?...

Thank you for your help

Edith



> On Nov 1, 2017, at 5:29 PM, Denis Davydov  wrote:
> 
> Hi Edith,
> 
> f(x-x0) is just your enrichment function derived from dealii::Function. So 
> define a class which takes that extra parameter (or any number of them) and 
> use it with FE_Enriched.
> Assuming that your bilinear form needs gradients, then you would need to 
> implement gradients of your custom function so that FE_Enriched can work out 
> chain rule for you.
> 
> Cheers,
> Denis.
> 
> On Wednesday, November 1, 2017 at 5:27:10 AM UTC+1, Edith Sotelo wrote:
> 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 static point.I think it 
> has to be  during the assembly step, but FE_values call the enrich function 
> with a single argument (the quadrature point). 
> What sould be a way to do pass the addtional point?
> 
> Thank you for you help.
> 
> 
> -- 
> 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 you are subscribed to a topic in the Google 
> Groups "deal.II User Group" group.
> To unsubscribe from this topic, visit 
> https://groups.google.com/d/topic/dealii/fFbITWIMW5M/unsubscribe 
> .
> To unsubscribe from this group and all its topics, send an email to 
> dealii+unsubscr...@googlegroups.com 
> .
> For more options, visit https://groups.google.com/d/optout 
> .

Edith







-- 
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 you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Re: [deal.II] CellDataStorage with mesh refinement

2017-11-01 Thread Timo Heister 
> 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) from a parent cell to each child? Do you also have an example how
> this works?

You can transfer whatever you want and the "status" argument will tell
you if you have to convert your data to children or parents.

Examples? It is used in SolutionTransfer and ContinuousQuadratureDataTransfer:
https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_dealii_dealii_blob_60bfd3170a465002ac1404dc36a5ba9671254b71_source_distributed_solution-5Ftransfer.cc-23L84=DwIBaQ=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw=4DhDywRdwgpQcqLQhc9e3vXnNh9twLFXh6cAv-rYOtU=HDE5UAeI6pVJZbKAVdMR4mt-Hb-RJxrDcUsUjCX4myo=
 
https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_dealii_dealii_blob_a0073ac1c0e88a30702a52f6522cbfffe1834ed1_include_deal.II_base_quadrature-5Fpoint-5Fdata.h-23L735=DwIBaQ=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw=4DhDywRdwgpQcqLQhc9e3vXnNh9twLFXh6cAv-rYOtU=j5xv1X624c05TmUCkWaj5qBdYefoSr6JBbEZ1-PIhcU=
 

We also use it in ASPECT to transfer particles:
https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_geodynamics_aspect_blob_0eec247edd99cd782a8fc66b36b140ce7a4105ae_source_particle_particle-5Fhandler.cc-23L858=DwIBaQ=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw=4DhDywRdwgpQcqLQhc9e3vXnNh9twLFXh6cAv-rYOtU=OmRYoVuKEIL1dlaaCOgQLqNoCmHaHzJulO4dsMuAq2o=
 

-- 
Timo Heister
http://www.math.clemson.edu/~heister/

-- 
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 you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

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

2017-11-01 Thread 'Uwe Köcher' via deal . II User Group 
Hej Pawan,

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.
>

you can try to use the candi installer for deal.II
https://github.com/dealii/candi

With that, it was possible to compile deal.II with dependencies on a Mac. 
The installation even worked after the latest Xcode upgrade and the upgrade 
to HighSierra.

Kind regards
  Uwe

-- 
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 you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.