``L.m'' and ``K.m''

Both L() and K() are diagonal operators based on the values of G2. We thus need only check whether the input has a size corresponding to G2 or G2c and to then use the values from the corresponding array. This can be accomplished by adding the statement ``global gbl_G2c'' to the top of each of the above functions and employing constructions ``if size(in,1)==length(gbl_G2c) [code block using G2c] else [code block using G2] endif'' which execute, as appropriate, either your current code or your current code with G2c replacing G2.



Tomas Arias 2006-04-26