On Mon, 11 May 2009, Sebastiano Putoto wrote:
Dear Gretl users,
I know my question may seem (it actually is) very banal, but I can't figure out a way
to check if two matrices are equal.
The point is verifying idempotence and symmetry for an orthogonal projector. The
mathematical proof isn't required here: the point is showing that the actual
calculated
matrices are orthogonal projectors and therefore have the specified properties; I got to
show the numbers are the same. And apart from arbitrary extraction of submatrices to
prove they get out same, I can't think of anything 'scientific'.
Hope you will be able to help me with the matter.
Regards,
How about
check = sumc(sumr(abs(mat1-mat2))) < 1.0e-07
or something similar?
Riccardo (Jack) Lucchetti
Dipartimento di Economia
Università Politecnica delle Marche
r.lucchetti(a)univpm.it
http://www.econ.univpm.it/lucchetti