public SparseLUSolver(int totalCases, double[,] matrix, double[] vector) { TotalCases = totalCases; CompressedColumnStorage <double> e = CompressedColumnStorage <double> .OfArray(matrix); SparseLUMatrix = SparseLU.Create(e, CSparse.ColumnOrdering.MinimumDegreeAtPlusA, 1.0); SparseLUVector = vector; SparseLUResultVector = Vector.Create(vector.Length, 1.0); SparseLUMethod(); }