LinearEquation.h

/********************************************************************/
/* Copyright (c) 2015 DG Technologies Inc. and Yuzi Mizuno          */
/* All rights reserved.                                             */
/********************************************************************/
#ifndef _MGLINEAREQUATION_HH_
#define _MGLINEAREQUATION_HH_


class MGMatrix;
class MGBPointSeq;

void factorizeLU(
    MGMatrix& W,
    int* id     
);



void solveLU(
    const MGMatrix& W,
    const int* id,  
    const MGBPointSeq& A,
    MGBPointSeq& X  
);

void solveSymetricTridiagonal(
    MGBPointSeq& W,
    const MGBPointSeq& A,
    MGBPointSeq& X  
);



int factorizeBandLU(
    MGBPointSeq& W,
    int nlower  
);



void solveBandLU(
    const MGBPointSeq& W,
    int nlower, 
    const MGBPointSeq& A,
    MGBPointSeq& X  
);


void factorizeCholeLU(
    MGBPointSeq& W 
);



void solveCholeLU(
    const MGBPointSeq& W,
    const MGBPointSeq& A,
    MGBPointSeq& X  
);
 // end of ALGORITHM group
#endif