Orthonormalization of a matrix

General Discussion
1

Hello dear FreeFem++ community,

I want to orthonormalize a rectangular real matrix A of size mxn in such a way to satisfy the property of a column-orthogonal matrix :

(A)^T A = Id_n

Is there any way to do it using FreeFem++?

Thank you in advance.