Changeset e0cb6c in git
- Timestamp:
- Oct 4, 2010, 10:47:47 AM (14 years ago)
- Branches:
- (u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '38dfc5131670d387a89455159ed1e071997eec94')
- Children:
- d322f7f690342c6a4c286ae7a938ccccacacb7f3
- Parents:
- a990f52adef334d4f73400116e91e4d0ec764bcd
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/gring.cc
ra990f5 re0cb6c 989 989 int k,m; 990 990 int rN=r->N; 991 matrix cMT=r->GetNC()->MT[UPMATELEM(j,i,rN)]; /* cMT=current MT */ 991 const int cMTindex = UPMATELEM(j,i,rN); 992 matrix cMT=r->GetNC()->MT[cMTindex]; /* cMT=current MT */ 992 993 993 994 poly x=pOne();p_SetExp(x,j,1,r);p_Setm(x,r); … … 1011 1012 // t=p_Copy(MATELEM(cMT,k-1,1),r); 1012 1013 t = gnc_mm_Mult_p(y,t,r); 1014 cMT=r->GetNC()->MT[cMTindex]; // since multiplication can change the MT table... 1015 assume( t != NULL ); 1016 #ifdef PDEBUG 1017 p_Test(t,r); 1018 #endif 1013 1019 MATELEM(cMT,k,1) = nc_p_CopyPut(t,r); 1014 1020 // omCheckAddr(cMT->m); … … 1025 1031 { 1026 1032 t = nc_p_CopyGet(MATELEM(cMT,a,m-1),r); 1033 assume( t != NULL ); 1027 1034 // t=p_Copy(MATELEM(cMT,a,m-1),r); 1028 1035 t = gnc_p_Mult_mm(t,x,r); 1036 cMT=r->GetNC()->MT[cMTindex]; // since multiplication can change the MT table... 1037 #ifdef PDEBUG 1038 p_Test(t,r); 1039 #endif 1029 1040 MATELEM(cMT,a,m) = nc_p_CopyPut(t,r); 1030 1041 // MATELEM(cMT,a,m) = t; … … 1036 1047 p_Delete(&x,r); 1037 1048 p_Delete(&y,r); 1038 // t=MATELEM(cMT,a,b); 1039 t= nc_p_CopyGet(MATELEM(cMT,a,b),r); 1049 t=MATELEM(cMT,a,b); 1050 assume( t != NULL ); 1051 1052 t= nc_p_CopyGet(t,r); 1053 #ifdef PDEBUG 1054 p_Test(t,r); 1055 #endif 1040 1056 // return(p_Copy(t,r)); 1041 1057 /* since the last computed element was cMT[a,b] */
Note: See TracChangeset
for help on using the changeset viewer.