search

JanJereczek / FastIsostasy.jl

Accelerated computation of glacial isostatic adjustment for laterally-variable solid-Earth
GNU General Public License v3.0
11 stars 0 forks source link

Compat with Float32; reduction of preallocated; remove jld2 + ncdatasets deps #125

Closed JanJereczek closed 3 months ago

JanJereczek commented 3 months ago

I repeatedly bump into memory issues when running high resolution (nx = ny > 600). This is partly addressed here by reducing the number of arrays defined in PreAllocated and by deleting the 3D arrays from LayeredEarth.

JanJereczek commented 3 months ago

Additionally, I removed the deps to JLD2 and to NCDatasets. The former is motivated by the fact that isostasy_data should not be julia-specific. The latter is motivated by the fact that we already need NetCDF.jl and hence we can easily remove NCDatasets.jl, which has the additional perc of being more legible.

JanJereczek commented 3 months ago

The code now includes a test with everything being Float32. This works and will avoid future issues about memory. In particular, this allows to run high resolution simulations on GPU!

JanJereczek commented 3 months ago

This addresses #119, #103, #120!