uses MtxExpr, Sparse, MtxVecEdit; procedure Example; var x,b: Vector; SparseA: TMtxSparse; begin SparseA := TMtxSparse.Create; try // load data SparseA.LoadFromMatrixFile('system.mtx'); b.LoadFromFile('coefficients.Vec'); // set solution size x.Size(b); // define iterative method - in this case Jacobi iterative method SparseA.RefinementIters := 4; SparseA.IterativeMethod := itmJAC; // solve SparseA.SolveIterative(b,x); // view solution ViewValues(x); finally SparseA.Free; end; end;

#include "MtxExpr.hpp" #include "Sparse.hpp" #include "MtxVecEdit.hpp" void __fastcall Example(); { TMtxSparse *SparseA = new TMtxSparse(); sVector x,b; try { // load data SparseA->LoadFromMatrixFile("system.mtx"); b.LoadFromFile("coefficients.Vec"); // set solution size x.Size(b); // define iterative method - in this case Jacobi iterative method SparseA->RefinementIters = 4; SparseA->IterativeMethod = itmJAC; // solve SparseA->SolveIterative(b,x); // view solution ViewValues(x); } __finally { delete SparseA; } }