Changeset 6cb81b in git
- Timestamp:
- Dec 6, 2018, 4:51:01 PM (5 years ago)
- Branches:
- (u'spielwiese', '4a9821a93ffdc22a6696668bd4f6b8c9de3e6c5f')
- Children:
- 138ecdd4aa5eadc67fa825adb892b78fd942cb32
- Parents:
- ab4286e2b8df65d60f51280bf2f1bf56bf2c88dc
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
libpolys/polys/matpol.cc
rab4286 r6cb81b 227 227 matrix c = mpNew(m,q); 228 228 229 for (i= 1; i<=m; i++)230 { 231 for (k= 1; k<=p; k++)229 for (i=0; i<m; i++) 230 { 231 for (k=0; k<p; k++) 232 232 { 233 233 poly aik; 234 if ((aik=MATELEM (a,i,k))!=NULL)235 { 236 for (j= 1; j<=q; j++)234 if ((aik=MATELEM0(a,i,k))!=NULL) 235 { 236 for (j=0; j<q; j++) 237 237 { 238 238 poly bkj; 239 if ((bkj=MATELEM (b,k,j))!=NULL)239 if ((bkj=MATELEM0(b,k,j))!=NULL) 240 240 { 241 poly *cij=&(MATELEM (c,i,j));242 poly s = pp_Mult_qq(aik /*MATELEM (a,i,k)*/, bkj/*MATELEM(b,k,j)*/, R);243 (*cij)/*MATELEM (c,i,j)*/ = p_Add_q((*cij) /*MATELEM(c,i,j)*/ ,s, R);241 poly *cij=&(MATELEM0(c,i,j)); 242 poly s = pp_Mult_qq(aik /*MATELEM0(a,i,k)*/, bkj/*MATELEM0(b,k,j)*/, R); 243 (*cij)/*MATELEM0(c,i,j)*/ = p_Add_q((*cij) /*MATELEM0(c,i,j)*/ ,s, R); 244 244 } 245 245 } 246 246 } 247 // pNormalize(t);248 // MATELEM(c,i,j) = t;249 247 } 250 248 }
Note: See TracChangeset
for help on using the changeset viewer.