// find strides for efficient lexsort
auto init_strides() const {
return compute_strides(shape);
}
vertical_advection_reference(repository &repo) : m_repo(repo), m_domain(repo.domain()), m_halo(repo.halo()) {
compute_strides(m_domain, m_strides);
}
repository(IJKSize domain, IJKSize halo)
: m_domain(domain), m_halo(halo), m_field_size(domain.m_i * domain.m_j * domain.m_k) {
compute_strides(m_domain, m_strides);
}