Quick Start

Minimal example: simple shear

using LinearAlgebra
using DRex

# Create an olivine mineral with 3500 grains
mineral = Mineral(
    phase   = olivine,
    fabric  = olivine_A,
    regime  = matrix_dislocation,
    n_grains = 3500,
    seed    = 8816,
)

# Analytical simple-shear velocity gradient (strain rate 1e-15 s⁻¹)
_, get_velocity_gradient = simple_shear_2d("X", "Z", 1e-15)

# Integrate CPO forward along a (stationary) pathline
params = default_params()
timestamps = range(0, 1e15, length=26)   # 25 steps

deformation_gradient = Matrix{Float64}(I, 3, 3)
for t in eachindex(timestamps)[2:end]
    global deformation_gradient = update_orientations!(
        mineral, params, deformation_gradient,
        get_velocity_gradient,
        (timestamps[t-1], timestamps[t], _ -> zeros(3)),
    )
end

# Post-process
orientations_resampled, _ = resample_orientations(mineral.orientations, mineral.fractions)
m_index = misorientation_index(orientations_resampled[end], orthorhombic)
println("Final M-index: ", round(m_index; sigdigits=3))

Key types

Key functions

See the API Reference for the full list.