Lines Matching refs:mkj
5375 giac_double * mkj=&m[k][i],*col0j=col0+i,*col1j=col1+i,*col2j=col2+i,*col3j=col3+i,*col0end=col0+k; in do_thread_double_linv() local
5380 mkj += shift; in do_thread_double_linv()
5387 for (;col0j<col0end;++mkj,++col3j,++col2j,++col1j,++col0j){ in do_thread_double_linv()
5388 giac_double tmp=(*mkj); in do_thread_double_linv()
5408 giac_double * mkj=&m[k][i],*colj=col+i,*colend=col+k; in do_thread_double_linv() local
5409 for (;colj<colend;++mkj,++colj) in do_thread_double_linv()
5410 res -= (*mkj)*(*colj); in do_thread_double_linv()
5455 …giac_double * mkj=&m[k][lastnon0pos],*col0j=col0+lastnon0pos,*colend=col0+k,*col1j=col1+lastnon0po… in do_thread_double_lu2inv() local
5456 for (;col0j>colend;--mkj,--col5j,--col4j,--col3j,--col2j,--col1j,--col0j){ in do_thread_double_lu2inv()
5457 giac_double tmp=*mkj; in do_thread_double_lu2inv()
5466 giac_double tmp=*mkj; in do_thread_double_lu2inv()
5479 giac_double * mkj=&m[k][n-1],*colj=col+n-1,*colend=col+k; in do_thread_double_lu2inv() local
5480 for (;colj>colend;--mkj,--colj){ in do_thread_double_lu2inv()
5481 res -= (*mkj)*(*colj); in do_thread_double_lu2inv()
5483 *colj=res/(*mkj); in do_thread_double_lu2inv()
5502 …giac_double * mkj=&m[k][lastnon0pos],*col0j=col0+lastnon0pos,*colend=col0+k,*col1j=col1+lastnon0po… in do_thread_double_lu2inv() local
5503 for (;col0j>colend;--mkj,--col5j,--col4j,--col3j,--col2j,--col1j,--col0j){ in do_thread_double_lu2inv()
5504 giac_double tmp=*mkj; in do_thread_double_lu2inv()
5513 *col0j=res0/(*mkj); in do_thread_double_lu2inv()
5514 *col1j=res1/(*mkj); in do_thread_double_lu2inv()
5515 *col2j=res2/(*mkj); in do_thread_double_lu2inv()
5516 *col3j=res3/(*mkj); in do_thread_double_lu2inv()
5517 *col4j=res4/(*mkj); in do_thread_double_lu2inv()
5518 *col5j=res5/(*mkj); in do_thread_double_lu2inv()
5525 giac_double * mkj=&m[k][cstart],*colj=col+cstart,*colend=col+k; in do_thread_double_lu2inv() local
5526 for (;colj>colend;--mkj,--colj){ in do_thread_double_lu2inv()
5527 res -= (*mkj)*(*colj); in do_thread_double_lu2inv()
5529 *colj=res/(*mkj); in do_thread_double_lu2inv()
5543 …giac_double * mkj=&m[k][cstart],*col0j=col0+cstart,*colend=col0+cend,*col1j=col1+cstart,*col2j=col… in do_thread_double_lu2inv() local
5544 for (;col0j>=colend;--mkj,--col5j,--col4j,--col3j,--col2j,--col1j,--col0j){ in do_thread_double_lu2inv()
5545 giac_double tmp=*mkj; in do_thread_double_lu2inv()
5566 giac_double * mkj=&m[k][cstart],*colj=col+cstart,*colend=col+cend; in do_thread_double_lu2inv() local
5567 for (;colj>=colend;--mkj,--colj){ in do_thread_double_lu2inv()
5568 res -= (*mkj)*(*colj); in do_thread_double_lu2inv()
5691 const int * mkj=&m[k+l][c]; in int_linsolve_l() local
5694 for (;aj<ak;++mkj,++aj) in int_linsolve_l()
5695 res -= longlong(*mkj)*(*aj); in int_linsolve_l()
5716 const int * mkj=&m[k+l][c]; in int_linsolve_l4() local
5725 for (;a0j<a0k;++mkj,++a0j,++a1j,++a2j,++a3j){ in int_linsolve_l4()
5726 longlong tmp=*mkj; in int_linsolve_l4()
5748 const double * mkj=&m[k+l][c]; in double_linsolve_l() local
5751 for (;aj<ak;++mkj,++aj) in double_linsolve_l()
5752 res -= (*mkj)*(*aj); in double_linsolve_l()
5773 const double * mkj=&m[k+l][c]; in double_linsolve_l4() local
5782 for (;a0j<a0k;++mkj,++a0j,++a1j,++a2j,++a3j){ in double_linsolve_l4()
5783 double tmp=*mkj; in double_linsolve_l4()
6007 int * mkj=&m[k][i],*col0j=col0+i,*col1j=col1+i,*col2j=col2+i,*col3j=col3+i,*col0end=col0+k; in int_lu2inv() local
6008 for (;col0j<col0end;++mkj,++col3j,++col2j,++col1j,++col0j){ in int_lu2inv()
6009 longlong tmp=(*mkj); in int_lu2inv()
6029 int * mkj=&m[k][i],*colj=col+i,*colend=col+k; in int_lu2inv() local
6030 for (;colj<colend;++mkj,++colj) in int_lu2inv()
6031 res -= longlong(*mkj)*(*colj); in int_lu2inv()
6052 …int * mkj=&m[k][n-1],*col0j=col0+n-1,*colend=col0+k,*col1j=col1+n-1,*col2j=col2+n-1,*col3j=col3+n-… in int_lu2inv() local
6053 for (;col0j>colend;--mkj,--col3j,--col2j,--col1j,--col0j){ in int_lu2inv()
6054 longlong tmp=*mkj; in int_lu2inv()
6061 int tmp=invmod(*mkj,p); in int_lu2inv()
6072 int * mkj=&m[k][n-1],*colj=col+n-1,*colend=col+k; in int_lu2inv() local
6073 for (;colj>colend;--mkj,--colj){ in int_lu2inv()
6074 res -= longlong(*mkj)*(*colj); in int_lu2inv()
6076 *colj=((res % p)*invmod(*mkj,p))%p; in int_lu2inv()