 Timestamp:
 Oct 4, 2010, 10:47:47 AM (14 years ago)
 Branches:
 (u'spielwiese', '4a9821a93ffdc22a6696668bd4f6b8c9de3e6c5f')
 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,k1,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,m1),r); 1033 assume( t != NULL ); 1027 1034 // t=p_Copy(MATELEM(cMT,a,m1),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.