view reshapeRowMaj.m @ 724:a9e8c9d71307 feature/poroelastic

Modify penalty for Neumann in Head2d so that data is for normal derivative and not u_x or u_y.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 14 Mar 2018 13:04:25 -0700
parents 5dad3b754c18
children
line wrap: on
line source

% Reshapes a matrix as if it was stored in row major order.
function B = reshapeRowMaj(A, m)
    D = length(m);

    if D == 1
        m = [m 1];
        D = 2;
    end

    % Reshape and reverse order of indecies
    B = permute(reshape(permute(A, ndims(A):-1:1), rot90(m,2)), D:-1:1);
end