uses MtxExpr, Sparse; procedure Example; var asp, bsp: TSparseMtx; csp, dsp: TSparseMtx; d: Vector; begin // ... // #1 : split asp into lower(bsp), upper(csp) and diagonal(d) sections bsp.LowerTriangle(asp); csp.UpperTriangle(asp); asp.Diag(d); // now combine lower(bsp), upper, and diagonal(d) back dsp.AddSplit(bsp,d,csp); // should be equal if not asp.Equal(dsp) then Eraise('Not Equal!'); end;

#include "MtxExpr.hpp" #include "Sparse.hpp" void __fastcall Example() { sVector d; TSparseMtx *asp = new TSparseMtx(); TSparseMtx *bsp = new TSparseMtx(); TSparseMtx *csp = new TSparseMtx(); TSparseMtx *dsp = new TSparseMtx(); // ... // #1 : split asp into lower(bsp), upper(csp) and diagonal(d) sections bsp->LowerTriangle(asp); csp->UpperTriangle(asp); asp->Diag(d); // now combine lower(bsp), upper, and diagonal(d) back dsp->AddSplit(bsp,d,csp); // should be equal if (!asp->Equal(dsp)) ERaise("Not Equal!"); delete asp; delete bsp; delete csp; delete dsp; }