Merge with feature/poroelastic. Use only the changes made to multiblock.DiffOp and scheme.Elastic2dVariable. DiffOp.getBoundaryOperator/Quadrature now use scheme methods instead of propeties.

function benchmark(func, N)
    default_arg('N',100);

    tic
    profile on

    for i = 1:N
        func();
    end

    profile viewer
end