Metric given by means a positive definite hermitian matrix

> F:=matrix(4,4,[1,1+I,2,3-I,1-I,4,2-I,3*I,2,2+I,8,3-I,3+I,-3*I,3+I,50]):

> rank_signature(F,complejo);

> B := matrix([[1+I, 2-I], [4*I, 2], [5-I, 3+2*I],[0,2-3*I]]);

> S:=SCHMIDT(B,F); Returns an ortogonal (unitary metric defined by F) basis in Col(B)

Testing the output