- Timestamp:
- Sep 29, 2011, 6:29:37 PM (13 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- 172bbb269fb8ec8ad467f6d62ee5ea9c0151b3ba
- Parents:
- bebd746739b98e1580e2d7b75cf326c9350f2a7b
- git-author:
- Oleksandr Motsak <motsak@mathematik.uni-kl.de>2011-09-29 18:29:37+02:00
- git-committer:
- Mohamed Barakat <mohamed.barakat@rwth-aachen.de>2011-11-09 16:13:49+01:00
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/maps_ip.cc
rbebd746 r83a1714 32 32 33 33 /* debug output: Tok2Cmdname in maApplyFetch*/ 34 //#include "ipshell.h"34 #include "ipshell.h" 35 35 36 36 // define this if you want to use the fast_map routine for mapping ideals … … 55 55 int i; 56 56 int N = preimage_r->N; 57 //Print("N=%d what=%s ",N,Tok2Cmdname(what)); 58 //if (perm!=NULL) for(i=1;i<=N;i++) Print("%d -> %d ",i,perm[i]); 59 //PrintS("\n"); 60 //Print("P=%d ",P); 61 //if (par_perm!=NULL) for(i=0;i<P;i++) Print("%d -> %d ",i,par_perm[i]); 62 //PrintS("\n"); 57 #if 0 58 Print("N=%d what=%s ",N,Tok2Cmdname(what)); 59 if (perm!=NULL) for(i=1;i<=N;i++) Print("%d -> %d ",i,perm[i]); 60 PrintS("\n"); 61 Print("P=%d ",P); 62 if (par_perm!=NULL) for(i=0;i<P;i++) Print("%d -> %d ",i,par_perm[i]); 63 PrintS("\n"); 64 #endif 65 63 66 void *data=w->Data(); 64 67 res->rtyp = w->rtyp; … … 68 71 if (P!=0) 69 72 { 70 WerrorS("Sorry 'napPermNumber' was lost in the refactoring process (due to Frank): needs to be fixed"); 71 return TRUE; 72 #if 0 73 res->data=(void *)napPermNumber((number)data,par_perm,P, preimage_r); 73 // WerrorS("Sorry 'napPermNumber' was lost in the refactoring process (due to Frank): needs to be fixed"); 74 // return TRUE; 75 #if 1 76 // poly n_PermNumber(const number z, const int *par_perm, const int OldPar, const ring src, const ring dst); 77 res->data= (void *) n_PermNumber((number)data, par_perm, P, preimage_r, currRing); 74 78 #endif 75 79 res->rtyp=POLY_CMD; … … 80 84 else 81 85 { 82 res->data=(void *)nMap((number)data, preimage_r->cf, currRing->cf); 86 assume( nMap != NULL ); 87 88 number a = nMap((number)data, preimage_r->cf, currRing->cf); 89 90 83 91 if (nCoeff_is_Extension(currRing->cf)) 84 92 { 85 number a=(number)res->data; 86 number one=nInit(1); 87 number product = nMult(a,one ); 93 n_Normalize(a, currRing->cf); // ??? 94 /* 95 number a = (number)res->data; 96 number one = nInit(1); 97 number product = nMult(a, one ); 88 98 nDelete(&one); 89 99 nDelete(&a); 90 100 res->data=(void *)product; 101 */ 91 102 } 92 103 #ifdef LDEBUG 93 n Test((number) res->data);104 n_Test(a, currRing->cf); 94 105 #endif 106 res->data=(void *)a; 107 95 108 } 96 109 break; … … 101 114 else 102 115 if ((what==IMAP_CMD) || ((what==FETCH_CMD) /* && (nMap!=nCopy)*/)) 103 res->data=(void *)p_PermPoly((poly)data,perm,preimage_r,currRing, 104 nMap,par_perm,P); 116 res->data=(void *)p_PermPoly((poly)data,perm,preimage_r,currRing, nMap,par_perm,P); 105 117 else /*if (what==MAP_CMD)*/ 106 118 {
Note: See TracChangeset
for help on using the changeset viewer.