TSparseMtx.Diag Method

Copies the kk-th diagonal of TSparseMtx to the Vec object.

Pascal

procedure Diag(const Vec: TVec; kk: integer = 0); overload;

Copies the kk-th diagonal of TSparseMtx to the Vec object. If k = 0 then the main diagonal matrix is obtained, if k < 0 then the k-th subdiagonal matrix is obtained and if k > 0 the then the k-th super diagonal is copied to Vec. The Vec parameters size is set automatically.

Example

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; }

TSparseMtx Class, TSparseMtx Members, Sparse Namespace, public, Example

What do you think about this topic? Send feedback!