24 #include "kernels_trace.h"
26 #ifndef DOXYGEN_SHOULD_SKIP_THIS
27 #define MAXSIZEOFBLOCKS 64
28 static double done = 1.0;
29 static double mdone = -1.0;
71 pastix_int_t *nbpivots,
80 if ( fabs(*Akk) < criterion ) {
81 (*Akk) = (double)criterion;
88 pastix_print_error(
"Negative diagonal term\n" );
95 cblas_dscal(n-k-1, ( alpha ), Amk, 1 );
100 cblas_dsyr(CblasColMajor, CblasLower,
147 pastix_int_t *nbpivots,
150 pastix_int_t k, blocknbr, blocksize, matrixsize;
151 double *tmp,*tmp1,*tmp2;
154 blocknbr = pastix_iceil( n, MAXSIZEOFBLOCKS );
156 for (k=0; k<blocknbr; k++) {
158 blocksize = pastix_imin(MAXSIZEOFBLOCKS, n-k*MAXSIZEOFBLOCKS);
159 tmp = A+(k*MAXSIZEOFBLOCKS)*(lda+1);
164 if ((k*MAXSIZEOFBLOCKS+blocksize) < n) {
166 tmp1 = tmp + blocksize;
167 tmp2 = tmp1 + blocksize * lda;
169 matrixsize = n-(k*MAXSIZEOFBLOCKS+blocksize);
175 cblas_dtrsm(CblasColMajor,
176 CblasRight, CblasLower,
177 CblasTrans, CblasNonUnit,
178 matrixsize, blocksize,
183 cblas_dsyrk(CblasColMajor, CblasLower, CblasNoTrans,
184 matrixsize, blocksize,
185 (
double)mdone, tmp1, lda,
186 (
double)done, tmp2, lda);
224 pastix_int_t ncols, stride;
225 pastix_int_t nbpivots = 0;
226 pastix_fixdbl_t time, flops;
229 double criterion = solvmtx->diagthreshold;
231 time = kernel_trace_start( PastixKernelPOTRF );
234 stride = (cblk->
cblktype & CBLK_LAYOUT_2D) ? ncols : cblk->
stride;
240 if ( cblk->
cblktype & CBLK_COMPRESSED ) {
246 assert( lrL->
rk == -1 );
247 assert( stride == lrL->
rkmax );
253 flops = FLOPS_DPOTRF( ncols );
254 kernel_trace_start_lvl2( PastixKernelLvl2POTRF );
256 kernel_trace_stop_lvl2( flops );
258 kernel_trace_stop( cblk->
fblokptr->
inlast, PastixKernelPOTRF, ncols, 0, 0, flops, time );
261 pastix_atomic_add_32b( &(solvmtx->nbpivots), nbpivots );
296 pastix_int_t nbpivots;
301 cblk, L, L, &(solvmtx->lowrank) );
342 pastix_int_t nbpivots;
350 for( ; blok < lblk; blok++ )
352 fcblk = solvmtx->cblktab + blok->
fcblknm;
354 if ( fcblk->
cblktype & CBLK_FANIN ) {
361 work, lwork, &(solvmtx->lowrank) );