// 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); }