Newer
Older
#include "cs.h"
///* solve Lx=b where x and b are dense. x=b on input, solution on output. */
//__global__
//void cs_lsolve_gpu (int n, const int *Lp, const int *Li, const double *Lx, double *x)
//{
// int p, j;
//
// // double *Lx ;
// // n = L->n ;
// // Lp = L->p ;
// // Li = L->i ;
// // Lx = L->x ;
//
// for (j = 0 ; j < n ; j++)
// {
// x [j] /= Lx [Lp [j]] ;
// for (p = Lp [j]+1 ; p < Lp [j+1] ; p++)
// {
// x [Li [p]] -= Lx [p] * x [j] ;
// }
// }
//
// //return (1);
//}