14-numcxx-lapack-sharedptr.cxx File Reference

Include dependency graph for 14-numcxx-lapack-sharedptr.cxx:

Go to the source code of this file.

Functions
int	main (void)

Function Documentation

int main

(

void

)

Examples:

14-numcxx-lapack-sharedptr.cxx.

Definition at line 12 of file 14-numcxx-lapack-sharedptr.cxx.

{
    int n=5000;
    auto pA=numcxx::DMatrix::create(n,n);
    auto pF=numcxx::DArray1::create(n);
    auto pU=numcxx::DArray1::create(n);
    
    // Here, we need to derefernce the smart pointer to references
    // in order to use the () operator in a convenient way.
    auto &A=*pA;
    auto &F=*pF;
    auto &U=*pU;
    
    F=1.0;
    for (int i=0;i<n;i++)
    {
        A(i,i)=3.0;
        if (i>0) A(i,i-1)=-1;
        if (i<n-1) A(i,i+1)=-1;
    }
    
    auto pLapack=numcxx::DSolverLapackLU::create(pA);
    pLapack->solve(U,F);
    double residual=normi(A*U-F);
    
    std::cout << "residual:" << residual << std::endl;
}

Here is the call graph for this function: