/// <summary>
///
/// </summary>
/// <param name="vector">The vector is always a column entity!!</param>
/// <returns></returns>
public _nVector MultiplyByVector(_nVector vector)
{
double[,] a = new double[m, n];
for (int j = 1; j <= n; j++)
{
for (int i = 1; i <= m; i++)
{
a[i - 1, j - 1] = GetValueAt(j, i) * vector.V(i);
}
}
var vec = new List <double>();
for (int i = 0; i < n; i++)
{
double temp = 0;
for (int j = 0; j < m; j++)
{
temp += a[j, i];
}
vec.Add(temp);
}
return(new _nVector(vec));
}
public double Norm(_nVector vec)
{
double sum = 0;
foreach (var component in vec.vector)
{
sum += Math.Pow(component, 2);
}
return(Math.Sqrt(sum));
}
public double DotProductWithAnotherVectorTranspose(_nVector vec)
{
double magnitude = 0;
int ctr = 0;
foreach (var component in vector)
{
magnitude += component * vec.vector[ctr];
ctr++;
}
return(magnitude);
}
public _nVector Subtract(_nVector nVector)
{
var vec = new List <double>();
int cntr = 0;
foreach (var value in vector)
{
vec.Add(value - nVector.vector[cntr]);
cntr++;
}
return(new _nVector(vec));
}
public void AddToUHouseholder(_nVector U, _nMatrix matrix, int i, int dim) // Householder Matrices are correct and functioning
{
var uMatrix = U.MatrixTypeMultiplyWithMe();
var resultafterUMatrixTimes2 = uMatrix.MultiplyByScalar(2);
var identity = new nIdentityMatrix();
var finalHouseholderMatrix = identity.ReturnNIdentityMatrix(i)
.SubtractAnotherMatrix(resultafterUMatrixTimes2);
var matrixToAdd = identity.ReturnNIdentityMatrix(dim);
matrixToAdd.InsertSubMatrix(finalHouseholderMatrix, i);
UHouseholderMatrices.Add(matrixToAdd);
//var result = matrixToAdd.MultiplyByAnotherMatrix(matrix);
}
public void AddToVHouseholder(_nVector V, _mnMatrix mnMatrix, int i, int dim)
{
var vMatrix = V.MatrixTypeMultiplyWithMe();
var resultafterUMatrixTimes2 = vMatrix.MultiplyByScalar(2);
var identity = new nIdentityMatrix();
var finvlHouseholderMatrix = identity.ReturnNIdentityMatrix(i - 1)
.SubtractAnotherMatrix(resultafterUMatrixTimes2);
var matrixToAdd = identity.ReturnNIdentityMatrix(dim - 1);
matrixToAdd.InsertSubMatrix(finvlHouseholderMatrix, i - 1);
VHouseholderMatrices.Add(matrixToAdd);
//var maybe = matrixToAdd.MultiplyByAnotherMatrix(matrixToAdd.ReturnTranspose(matrixToAdd));
//var result = mnMatrix.MultiplyMeByNMatrix(matrixToAdd);
}
public _nVector GetSubRow(int m, int k)
{
var mRows = new List <_nVector>();
int skipRows = n - m, skipCols = n - k;
var nVector = new List <double>();
for (int j = skipCols; j < n; j++)
{
nVector.Add(rows[skipRows].vector[j]);
}
var tempNVector = new _nVector(nVector);
return(tempNVector);
}
public _mnMatrix GetSubMatrix(int mm, int k)
{
var mRows = new List <_nVector>();
int skipRows = m - mm, skipCols = m - k;
for (int i = skipRows; i < m; i++)
{
var nVector = new List <double>();
for (int j = skipCols; j < m; j++)
{
nVector.Add(rows[i].vector[j]);
}
var tempNVector = new _nVector(nVector);
mRows.Add(tempNVector);
}
return(new _mnMatrix(mRows, mm, k));
}
public _nMatrix GetSubMatrix(int m)
{
var mRows = new List <_nVector>();
int skipRows = n - m, skipCols = n - m;
for (int i = skipRows; i < n; i++)
{
var nVector = new List <double>();
for (int j = skipRows; j < n; j++)
{
nVector.Add(rows[i].vector[j]);
}
var tempNVector = new _nVector(nVector);
mRows.Add(tempNVector);
}
return(new _nMatrix(mRows, m));
}
public _mnMatrix mnMatrixTypeMultiplyWithAnotherVector(_nVector vec)
{
var componentList = new List <double>();
var dictionary = new List <_nVector>();
int i = 0, j = 0;
foreach (var component in vector)
{
var vctor = new List <double>();
foreach (var component2 in vec.vector)
{
vctor.Add(component * component2);
}
dictionary.Add(new _nVector(vctor));
i++;
}
return(new _mnMatrix(dictionary, i, j));
}
public _nMatrix MatrixTypeMultiplyWithAnotherVector(_nVector vec)
{
var componentList = new List <double>();
int n = this.vector.Count();
double[,] matrixComponents = new double[n, n];
var dictionary = new List <_nVector>();
int i = 0, j = 0;
foreach (var component in vector)
{
var vctor = new List <double>();
foreach (var component2 in vec.vector)
{
vctor.Add(component * component2);
}
dictionary.Add(new _nVector(vctor));
i++;
}
return(new _nMatrix(dictionary, n));
}
public _nMatrix MultiplyByAnotherMatrix(_nMatrix matrix)
{
_nVector[] _row = new _nVector[n];
double[] a = new double[n];
for (int k = 1; k <= n; k++)
{
for (int j = 1; j <= n; j++)
{
for (int i = 1; i <= n; i++)
{
var tmp1 = GetValueAt(k, i); var tmp2 = matrix.GetValueAt(i, j);
a[j - 1] += GetValueAt(k, i) * matrix.GetValueAt(i, j);
}
var tmp3 = a[j - 1];
}
_row[k - 1] = new _nVector(a.ToList());
a = new double[n];
}
return(new _nMatrix(_row.ToList(), n));
}
public _mnMatrix MultiplyMeByNMatrix(_nMatrix matrix)
{
_nVector[] _row = new _nVector[n];
var mm = matrix.rows.Count();
double[] a = new double[n];
for (int k = 1; k <= n; k++)
{
for (int j = 1; j <= mm; j++)
{
for (int i = 1; i <= mm; i++)
{
//var tmp1 = GetValueAt(k, i); var tmp2 = matrix.GetValueAt(i, j);
a[j - 1] += GetValueAt(k, i) * matrix.GetValueAt(i, j);
}
// var tmp3 = a[j - 1];
}
_row[k - 1] = new _nVector(a.ToList());
a = new double[n];
}
return(new _mnMatrix(_row.ToList(), n, mm));
}
public _nMatrix ReturnNIdentityMatrix(int n)
{
var matrixRows = new List <_nVector>();
for (int i = 0; i < n; i++)
{
var row = new List <double>();
for (int j = 0; j < n; j++)
{
if (j == i)
{
row.Add(1);
}
else
{
row.Add(0);
}
}
var vectorRow = new _nVector(row);
matrixRows.Add(vectorRow);
}
return(new _nMatrix(matrixRows, n));
}
