Merged in feature/dirac_discr (pull request #17) Add multi-d dirac discretization with tests Approved-by: Vidar Stiernström <vidar.stiernstrom@it.uu.se> Approved-by: Martin Almquist <malmquist@stanford.edu> Approved-by: Jonatan Werpers <jonatan.werpers@it.uu.se>

function default_struct(s, val)
    if evalin('caller',sprintf('~exist(''%s'',''var'')',s))
        given = [];
    else
        given = evalin('caller', s);
    end

    final = copyWithDefault(given, val);
    assignin('caller', s, final);
end