Hi, Wolfgang,

Thank you very much.

Best,
Xiang

On Friday, 8 May 2020 06:03:04 UTC-7, Wolfgang Bangerth wrote:
>
> On 5/7/20 10:25 PM, 孙翔 wrote: 
> > I have a .txt file which stores a permeability distribution in the 
> z-direction 
> > as follows 
> > 
> > z                 value 
> > 0                 0.1 
> > 1                 0.3 
> > 2                 0.5 
> > 3                 0.2 
> > 4                 0.8 
> > 5                 0.4 
> > 6                 0.7 
> > 7                 0.5 
> > 8                 0.6 
> > ...... 
> > 
> > Does Dealii have a function to read it into the code and apply the 
> > distribution in the model? Maybe, there is an interpolation function 
> that can 
> > help me get the value at quadrature points. if there is no such a 
> function, 
> > can I get the real coordinate at the quadrature point? I can wirte my 
> own 
> > interpolation function to get the coefficient at the quadrature point. 
>
> Xiang, 
> there is the InterpolatedTensorProductGridData that can do the 
> interpolation 
> for you. It's meant for arbitrary dimensions, but of course also works in 
> 1d. 
> But there are then two pieces missing that you need to provide yourself: 
>
> * You need to read the data into the arrays with which you create objects 
> of 
> the class mentioned above. 
>
> * You will create an object of that type with dim==1. In your code, you 
> will ask 
>    fe_values.quadrature_point(q) 
> to get the coordinates of the q'th quadrature point, of which you will 
> then 
> only use the last coordinate value to call the 1d function. 
>
> Best 
>   W. 
>
> -- 
> ------------------------------------------------------------------------ 
> Wolfgang Bangerth          email:                 bang...@colostate.edu 
> <javascript:> 
>                             www: http://www.math.colostate.edu/~bangerth/ 
>
>

-- 
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.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/9980a489-bf6f-4e2c-af87-5d2efba87930%40googlegroups.com.

Reply via email to