Matrix-Vector Operations in IronPython QuickStart Sample

Illustrates how to perform operations that involve both matrices and vectors in IronPython.

View this sample in: C# Visual Basic F#

```Python
import numerics

from Extreme.Mathematics import *

# The Vector and DenseMatrix classes reside in the 
# Extreme.Mathematics.LinearAlgebra namespace.
from Extreme.Mathematics.LinearAlgebra import *

# Illustrates operations on DenseMatrix objects and combined
# operations on Vector and DenseMatrix objects from the
# Extreme.Mathematics.LinearAlgebra namespace of Numerics.NET.

# For details on the basic workings of Vector 
# objects, including constructing, copying and
# cloning vectors, see the BasicVectors QuickStart
# Sample.
#
# For details on the basic workings of DenseMatrix
# objects, including constructing, copying and
# cloning vectors, see the BasicVectors QuickStart
# Sample.
#
# Let's create some vectors to work with.
v1 = Vector([1, 2, 3, 4, 5])
v2 = Vector([1, -2, 3, -4, 5])
print "v1 = {0:.4f}".format(v1)
print "v2 = {0:.4f}".format(v2)

# Also, here are a couple of matrices.
# We start out with a 5x5 identity matrix:
m1 = DenseMatrix.GetIdentity(5)
# Now we use the GetDiagonal method and combine it
# with the SetValue method of the Vector class to
# set some of the off-diagonal elements:
m1.GetDiagonal(1).SetValue(2)
m1.GetDiagonal(2).SetValue(3)
m1.GetDiagonal(-1).SetValue(4)
print "m1 = {0:.4f}".format(m1)
# We define our second matrix by hand:
m2 = Matrix([[1,1,1,1,1], [2,3,4,8,-1], [3,5,9,27,1], [4,7,16,64,-1], [5,9,25,125,1]])
print "m2 = {0:.4f}".format(m2)
print 
			
#
# Matrix arithmetic
#

# The Matrix class defines operator overloads for
# addition, subtraction, and multiplication of
# matrices.

# Addition:
print "Matrix arithmetic:"
m = m1 + m2
print "m1 + m2 = {0:.4f}".format(m)
# Subtraction:
m = m1 - m2
print "m1 - m2 = {0:.4f}".format(m)
# Multiplication is the True matrix product:
m = m1 * m2
print "m1 * m2 = {0:.4f}".format(m)
print 

#
# Matrix-Vector products
#

# The DenseMatrix class defines overloaded addition, # subtraction, and multiplication operators 
# for vectors and matrices:
print "Matrix-vector products:"
v = m1 * v1
print "m1 v1 = {0:.4f}".format(v)
# You can also multiply a vector by a matrix on the right.
# This is equivalent to multiplying on the left by the 
# transpose of the matrix:
v = v1 * m1
print "v1 m1 = {0:.4f}".format(v)
			
# Now for some methods of the Vector class that
# involve matrices:
# Add a product of a matrix and a vector:
v.Add(m1, v1)
print "v + m1 v1 = {0:.4f}".format(v)
# Or add a scaled product:
v.Add(-2, m1, v2)
print "v - 2 m1 v2 = {0:.4f}".format(v)
# You can also use static Subtract methods:
v.Subtract(m1, v1)
print "v - m1 v1 = {0:.4f}".format(v)
print 

#
# Matrix norms
#
print "Matrix norms"
# Matrix norms are not as easily defined as
# vector norms. Three matrix norms are available.
# 1. The one-norm through the OneNorm property:
a = m2.OneNorm()
print "OneNorm of m2 = {0:.4f}".format(a)
# 2. The infinity norm through the 
#    InfinityNorm property:
a = m2.InfinityNorm()
print "InfinityNorm of m2 = {0:.4f}".format(a)
# 3. The Frobenius norm is often used because it
#    is easy to calculate.
a = m2.FrobeniusNorm()
print "FrobeniusNorm of m2 = {0:.4f}".format(a)
print 

# The trace of a matrix is the sum of its diagonal
# elements. It is returned by the Trace property:
a = m2.Trace()
print "Trace(m2) = {0:.4f}".format(a)

# The Transpose method returns the transpose of a 
# matrix. This transposed matrix shares element storage
# with the original matrix. Use the CloneData method
# to give the transpose its own data storage.
m = m2.Transpose()
print "Transpose(m2) = {0:.4f}".format(m)
```