| 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]); |
| } |
} |
| } |
} |