1 classdef LUPQ_Handle < handle
2 % handle
for passing l-u-p-q-factors to functions
12 function this = LUPQ_Handle(varargin)
14 [this.L, this.U, this.P, this.Q] = lu(varargin{1});
23 % matrix multiplication
24 function mat = mtimes(
this, mat)
25 mat = this.P' * (this.L * (this.U * (this.Q' * mat)));
28 % matrix left division
29 function mat = mldivide(this, mat)
30 mat = this.Q * (this.U \ (this.L \ (this.P * mat)));