125 ring dring =
IDRING( dringHdl );
126 ring sring =
IDRING( sringHdl );
130 WerrorS(
"rings must have same characteristic" );
133 if ( (sring->OrdSgn != 1) || (dring->OrdSgn != 1) )
135 WerrorS(
"only works for global orderings" );
138 if ( sring->N != dring->N )
140 WerrorS(
"rings must have same number of variables" );
145 WerrorS(
"rings must have same number of parameters" );
148 if ( state !=
FglmOk )
return state;
152 int npar =
rPar(sring);
155 pperm= (
int *)
omAlloc0( (npar+1)*
sizeof( int ) );
159 dring->names, nvar,
rParameter(dring), npar, vperm, pperm,
161 for (
k= nvar; (
k > 0) && (state ==
FglmOk);
k-- )
164 WerrorS(
"variable names do not agree" );
167 for (
k= npar-1; (
k >= 0) && (state ==
FglmOk);
k-- )
170 WerrorS(
"parameter names do not agree" );
175 if ( state !=
FglmOk )
return state;
177 if ( sring->qideal !=
NULL )
179 if ( dring->qideal ==
NULL )
181 Werror(
"%s is a qring, current ring not", sringHdl->
id );
189 for (
k=
IDELEMS( sring->qideal )-1;
k >= 0;
k-- )
190 (sqind->m)[
k]=
p_PermPoly( (sring->qideal->m)[
k], vperm, sring,
192 ideal sqindred =
kNF( dring->qideal,
NULL, sqind );
193 if ( !
idIs0( sqindred ) )
195 WerrorS(
"the quotients do not agree" );
201 if ( state !=
FglmOk )
return state;
203 int * dsvperm = (
int *)
omAlloc0( (nvar+1)*
sizeof( int ) );
205 dsvperm,
NULL, sring->cf->type);
208 for (
k=
IDELEMS( dring->qideal )-1;
k >= 0;
k-- )
209 (dqins->m)[
k]=
p_PermPoly( (dring->qideal->m)[
k], dsvperm, sring,
211 ideal dqinsred =
kNF( sring->qideal,
NULL, dqins );
212 if ( !
idIs0( dqinsred ) )
214 WerrorS(
"the quotients do not agree" );
220 if ( state !=
FglmOk )
return state;
224 if ( dring->qideal !=
NULL )
226 Werror(
"current ring is a qring, %s not", sringHdl->
id );
static FORCE_INLINE nMapFunc n_SetMap(const coeffs src, const coeffs dst)
set the mapping function pointers for translating numbers from src to dst
number(* nMapFunc)(number a, const coeffs src, const coeffs dst)
maps "a", which lives in src, into dst
void WerrorS(const char *s)
#define idDelete(H)
delete an ideal
BOOLEAN idIs0(ideal h)
returns true if h is the zero ideal
poly kNF(ideal F, ideal Q, poly p, int syzComp, int lazyReduce)
void maFindPerm(char const *const *const preim_names, int preim_n, char const *const *const preim_par, int preim_p, char const *const *const names, int n, char const *const *const par, int nop, int *perm, int *par_perm, n_coeffType ch)
#define omFreeSize(addr, size)
poly p_PermPoly(poly p, const int *perm, const ring oldRing, const ring dst, nMapFunc nMap, const int *par_perm, int OldPar, BOOLEAN use_mult)
VAR ring currRing
Widely used global variable which specifies the current polynomial ring for Singular interpreter and ...
void Werror(const char *fmt,...)
static char const ** rParameter(const ring r)
(r->cf->parameter)
static int rPar(const ring r)
(r->cf->P)
ideal idInit(int idsize, int rank)
initialise an ideal / module