version 1.6, 2001/10/09 01:36:14 |
version 1.10, 2015/08/08 14:19:41 |
|
|
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* |
* |
* $OpenXM: OpenXM_contrib2/asir2000/engine/up_lm.c,v 1.5 2001/09/03 07:01:06 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/up_lm.c,v 1.9 2015/08/06 10:01:52 fujimoto Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include <math.h> |
#include <math.h> |
Line 251 void crup_lm(ModNum **f,int d,int *mod,int index,N m,N |
|
Line 251 void crup_lm(ModNum **f,int d,int *mod,int index,N m,N |
|
k[j] += c2*f[i][j]; |
k[j] += c2*f[i][j]; |
} |
} |
uiarraytoup(sum,len,d,&s); |
uiarraytoup(sum,len,d,&s); |
GC_free(sum_b); |
GCFREE(sum_b); |
|
|
u = UPALLOC(d); |
u = UPALLOC(d); |
for ( j = 0; j <= d; j++ ) { |
for ( j = 0; j <= d; j++ ) { |
#if 1 |
#if 1 |
a = (UL)floor(k[j]); |
a = (UL)floor(k[j]); |
#if defined(i386) || defined(__alpha) || defined(VISUAL) |
#if defined(i386) || defined(__alpha) || defined(VISUAL) || defined(__MINGW32__) || defined(__MINGW64__) || defined(__x86_64) |
au = ((unsigned int *)&a)[1]; |
au = ((unsigned int *)&a)[1]; |
al = ((unsigned int *)&a)[0]; |
al = ((unsigned int *)&a)[0]; |
#else |
#else |
Line 500 void hybrid_powertabup(UP f,UP xp,UP *tab) |
|
Line 500 void hybrid_powertabup(UP f,UP xp,UP *tab) |
|
invmodup(t,f->d,&invf); |
invmodup(t,f->d,&invf); |
|
|
for ( i = 2; i < d; i++ ) { |
for ( i = 2; i < d; i++ ) { |
if ( debug_up ) |
if ( debug_up ){ |
fprintf(stderr,"."); |
fprintf(stderr,"."); |
hybrid_mulup(FF_GFP,tab[i-1],xp,&t); |
#if defined(__MINGW32__) || defined(__MINGW64__) |
hybrid_rembymulup_special(FF_GFP,t,f,invf,&tab[i]); |
fflush(stderr); |
|
#endif |
|
} |
|
hybrid_mulup(FF_GFP,tab[i-1],xp,&t); |
|
hybrid_rembymulup_special(FF_GFP,t,f,invf,&tab[i]); |
} |
} |
} |
} |
|
|
Line 523 void powertabup(UP f,UP xp,UP *tab) |
|
Line 527 void powertabup(UP f,UP xp,UP *tab) |
|
invmodup(t,f->d,&invf); |
invmodup(t,f->d,&invf); |
|
|
for ( i = 2; i < d; i++ ) { |
for ( i = 2; i < d; i++ ) { |
if ( debug_up ) |
if ( debug_up ){ |
fprintf(stderr,"."); |
fprintf(stderr,"."); |
kmulup(tab[i-1],xp,&t); |
#if defined(__MINGW32__) || defined(__MINGW64__) |
rembymulup_special(t,f,invf,&tab[i]); |
fflush(stderr); |
|
#endif |
|
} |
|
kmulup(tab[i-1],xp,&t); |
|
rembymulup_special(t,f,invf,&tab[i]); |
} |
} |
} |
} |