#! /usr/bin/env python

from openturns import *

TESTPREAMBLE()

matrix1 = IdentityMatrix(2)
print "matrix1 = ", repr(matrix1)

pt = NumericalPoint()
pt.add(5.)
pt.add(0.)
print "pt = ", repr(pt)

result = NumericalPoint()
result = matrix1.solveLinearSystem(pt)
print "result = ", repr(result)

determinant = matrix1.computeDeterminant()
print "determinant = %.6f" % determinant

ev = NumericalScalarCollection(2)
ev = matrix1.computeEigenValues()
print "ev = ", repr(ev)

if matrix1.isPositiveDefinite():
    isSPD = "true"
else:
    isSPD = "false"
print "isSPD = ", isSPD

matrix2 = matrix1.computeCholesky()
print "matrix2 = ", repr(matrix2)
