function diffOpsStruct = constructDiffOps(method, varargin)
	switch method
	case {'bdf','backwards-difference-formula'}
		diffOpsStruct = rv.diffops.constructDiffOpsBdf(varargin{:});
	case {'ms','multi-stage'}
		diffOpsStruct = rv.diffops.constructDiffOpsMultiStage(varargin{:});
	case {'mg','multi-grid'}
		diffOpsStruct = rv.diffops.constructDiffOpsMultiGrid(varargin{:});
end