Add matfaust.lazylinop.LazyLinearOperator ctor and update unit tests consequently.

misc/test/src/Matlab/LazyLinearOpFFTTest.m

@@ -20,7 +20,8 @@ classdef LazyLinearOpFFTTest < LazyLinearOpTest
 		function instantiateTestFaust(this)
 			import matfaust.rand
             import matfaust.lazylinop.*
-			this.lop = LazyLinearOp.create_from_funcs(@(x) fft(x), @(x) 8 * ifft(x), [8, 8]);
+			%this.lop = LazyLinearOp.create_from_funcs(@(x) fft(x), @(x) 8 * ifft(x), [8, 8]);
+            this.lop = LazyLinearOperator([8, 8], 'matmat', @(x) fft(x), 'rmatvec', @(x) 8 * ifft(x))
             this.from_func = true;
             this.lopA = full(this.lop);
             F2 = rand(8, 8);

wrapper/matlab/+matfaust/+lazylinop/LazyLinearOperator.m

+function L = LazyLinearOperator(shape, varargin)
+    import matfaust.lazylinop.LazyLinearOp
+    p = inputParser;
+    validDtype = @(dtype) any(strcmp(dtype, {'complex', 'double', 'single', 'undefined'}));
+    addOptional(p, 'dtype', 'undefined', validDtype)
+    addOptional(p, 'matvec', 'none', @(l) isa(l, 'function_handle'))
+    addOptional(p, 'rmatvec', 'none', @(l) isa(l, 'function_handle'))
+    addOptional(p, 'matmat', 'none', @(l) isa(l, 'function_handle'))
+    addOptional(p, 'rmatmat', 'none', @(l) isa(l, 'function_handle'))
+
+    parse(p, varargin{:})
+    dtype = p.Results.dtype;
+    matvec = p.Results.matvec;
+    rmatvec = p.Results.rmatvec;
+    matmat = p.Results.matmat;
+    rmatmat = p.Results.rmatmat;
+
+    if strcmp(matvec, 'none') && strcmp(matmat, 'none')
+        error(['At least a matvec or a matmat function must be passed in'
+        ' kwargs'])
+    end
+
+    % define a matmat function from matvec in case matmat is not provided
+    function out = matmat_(M, matvec_)
+        out = zeros(shape(1), size(M, 2));
+        for i=1:size(M, 2)
+            out(:, i) = matvec_(M(:, i));
+        end
+    end
+
+    if strcmp(matmat, 'none')
+        matmat = @(M) matmat_(M, matvec);
+    end
+
+    if strcmp(rmatmat, 'none') && ~ strcmp(rmatvec, 'none')
+        rmatmat = @(M) matmat_(M, rmatvec);
+    end
+
+    L = LazyLinearOp.create_from_funcs(matmat, rmatmat, shape);
+end