Metric given by means of a positive definite hermitian matrix

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

> rank_signature(F); In fact, F is positive definite

> S:=GRAM_SCHMIDT(B,F);

Testing the output