[aspect-devel] log file

[Thorsten Becker]

As Timo said, the code is presumably loading a spherical harmonic expansion via a library package to evaluate a single point, each time for each integration point. Besides efficiency questions, there's no point in doing that, of course, if the initial models are not parametrized as spherical harmonics. I think it would be preferable to load a 3d grid once and then have a simple (eg bi linear) interpolation routine spit out values at desired locations. Even with a clunky importation (eg via slices of Netcdf grids, read by each processor, as we do for citcom), this shouldn't take more than a few tens of seconds after the I/O  for a ~1million elements. 

Thorsten W Becker
geodynamics.usc.edu

> On Jul 10, 2015, at 23:36, Wolfgang Bangerth <bangerth at tamu.edu> wrote:
> 
>> On 07/10/2015 02:40 PM, Timo Heister wrote:
>> 
>> I assume this is from your initial condition: S40RTS, right? I haven't
>> looked at it closely, but it looks like the construction and spline
>> interpolation is done for every single quadrature point. One should
>> really precompute as much as possible here.
> 
> Indeed, and likely the same with the Savani initial temperature model. I recall Shangxin telling me that it took very long to run through his code. I don't know what he is doing in there -- I don't think that the slowdown is in the ASPECT code itself, but in the initial temperature code.
> 
> Best
> W.
> 
> 
> -- 
> ------------------------------------------------------------------------
> Wolfgang Bangerth               email:            bangerth at math.tamu.edu
>                                www: http://www.math.tamu.edu/~bangerth/
> 
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel