public AMatrix <double[, ]> Multiply(AMatrix <double[, ]> element_0, double other)
{
return(new MatrixKozzion <double>(this, Multiply(element_0.Data, other)));
}
public AMatrix <double[, ]> Multiply(IReadOnlyList <double> operant_0, AMatrix <double[, ]> operant_1)
{
throw new NotImplementedException();
}
public override AMatrix <MatrixType> Integrate(IFunction <AMatrix <MatrixType>, AMatrix <MatrixType> > differntial, AMatrix <MatrixType> inital_value, double step_size)
{
AMatrix <MatrixType> k1 = differntial.Compute(inital_value) * (step_size / 6.0);
AMatrix <MatrixType> k2 = differntial.Compute(inital_value + (k1 * 0.5)) * (step_size / 3.0);
AMatrix <MatrixType> k3 = differntial.Compute(inital_value + (k2 * 0.5)) * (step_size / 3.0);
AMatrix <MatrixType> k4 = differntial.Compute(inital_value + k3) * (step_size / 6.0);
AMatrix <MatrixType> increment = k1 + k2 + k3 + k4;
return(inital_value + increment);
}
public AMatrix <double[, ]> Multiply(AMatrix <double[, ]> element_0, AMatrix <double[, ]> element_1)
{
return(new MatrixKozzion <double>(this, Multiply(element_0.Data, element_1.Data)));
}
public AMatrix <double[, ]> SelectRows(AMatrix <double[, ]> operant_0, IReadOnlyList <int> row_indexes)
{
throw new NotImplementedException();
}
public AMatrix <double[, ]> SelectColumns(AMatrix <double[, ]> operant_0, IList <int> column_indexes)
{
throw new NotImplementedException();
}
public AMatrix <MatrixType>[] Integrate(IFunction <AMatrix <MatrixType>, AMatrix <MatrixType> > differential_equation, AMatrix <MatrixType> initial_value, double step_size, int step_count)
{
AMatrix <MatrixType>[] steps = new AMatrix <MatrixType> [step_count + 1]; // Add initial step
steps[0] = initial_value;
for (int time_step = 0; time_step < step_count; time_step++)
{
steps[time_step + 1] = Integrate(differential_equation, steps[time_step], step_size);
}
return(steps);
}
public AMatrix <double[, ]> ComputeRoot(AMatrix <double[, ]> operant_0)
{
throw new NotImplementedException();
}
public SVD(AMatrix <MatrixType> u, AMatrix <MatrixType> s, AMatrix <MatrixType> vt)
{
this.U = u;
this.S = s;
this.VT = vt;
}
public abstract AMatrix <MatrixType> Integrate(IFunction <AMatrix <MatrixType>, AMatrix <MatrixType> > differential_equation, AMatrix <MatrixType> value, double step_size);
public DimensionReductionPCA(IDataContext data_context, AMatrix <MatrixType> projection)
{
this.DataContext = data_context;
this.algebra = projection.Algebra;
this.projection = projection;
}
public AMatrix <Matrix <double> > Solve(AMatrix <Matrix <double> > A, AMatrix <Matrix <double> > b, AMatrix <Matrix <double> > initial, int max_iter)
{
return(new MatrixMathNet(A.Data.Solve(b.Data)));
}
public AMatrix <Matrix <double> > Solve(AMatrix <Matrix <double> > A, AMatrix <Matrix <double> > b)
{
return(new MatrixMathNet(A.Data.Solve(b.Data)));
}
public AMatrix <float[, ]> ComputeRoot(AMatrix <float[, ]> covariance_matrix)
{
throw new NotImplementedException();
}
public AMatrix <double[, ]> Add(AMatrix <double[, ]> element_0, AMatrix <double[, ]> element_1)
{
return(new MatrixKozzion <double>(this, Add(element_0.Data, element_1.Data)));
}
public double[,] ToArray2D(AMatrix <double[, ]> operant_0)
{
throw new NotImplementedException();
}
public AMatrix <double[, ]> Add(AMatrix <double[, ]> element_0, double other)
{
return(new MatrixKozzion <double>(this, Add(element_0.Data, other)));
}
public SVD <MatrixType> ComputeSVD <MatrixType>(AMatrix <MatrixType> data)
{
throw new NotImplementedException();
}
public AMatrix <double[, ]> Subtract(AMatrix <double[, ]> element_0, AMatrix <double[, ]> element_1)
{
return(new MatrixKozzion <double>(this, Subtract(element_0.Data, element_1.Data)));
}
public SVD <double[, ]> ComputeSVD(AMatrix <double[, ]> data)
{
throw new NotImplementedException();
}
public AMatrix <double[, ]> Subtract(AMatrix <double[, ]> element_0, double other)
{
return(new MatrixKozzion <double>(this, Subtract(element_0.Data, other)));
}
public AMatrix <double[, ]> Add(AMatrix <double[, ]> operant_0, IList <double> operant_1)
{
throw new NotImplementedException();
}
public IFunction <double, double> GetIRF(double[] signal)
{
AMatrix <MatrixType> weights = solver.Solve(forward_matrix, algebra.Create(signal, true));
return(new FunctionWeigthed <double, double>(new AlgebraRealFloat64(), basis_function_list, weights.ToArray1DFloat64()));
}
