+
+ But if you supply a custom inner product, the result is
+ orthonormal with respect to that (and not the usual inner
+ product)::
+
+ sage: v1 = vector(QQ,(1,2,3))
+ sage: v2 = vector(QQ,(1,-1,6))
+ sage: v3 = vector(QQ,(2,1,-1))
+ sage: v = [v1,v2,v3]
+ sage: B = matrix(QQ, [ [6, 4, 2],
+ ....: [4, 5, 4],
+ ....: [2, 4, 9] ])
+ sage: ip = lambda x,y: (B*x).inner_product(y)
+ sage: norm = lambda x: ip(x,x)
+ sage: u = gram_schmidt(v,ip)
+ sage: all( norm(u_i) == 1 for u_i in u )
+ True
+ sage: ip(u[0],u[1]).is_zero()
+ True
+ sage: ip(u[0],u[2]).is_zero()
+ True
+ sage: ip(u[1],u[2]).is_zero()
+ True
+