Ejemplos de FullMatrix::get_Value en C++ (Cpp)

Ejemplo n.º 1

Archivo: ConjugateGradient.cpp Proyecto: smoe1/ImplicitExplicit

void MatMult(const FullMatrix& A,
        const dTensor1& invec,
        dTensor1& outvec)
{


    for (int i=1; i<=A.get_NumRows(); i++)
    {
        int jmax = A.get_NZrow(i);
        double tmp = 0.0;

        for (int j=1; j<=jmax; j++)
        {
            int ind    = A.get_Index(i,j);
            double val = A.get_Value(i,j);

            tmp = tmp + val*invec.get(ind);
        }
        outvec.set(i, tmp );
    }

}