version 1.6, 2002/09/08 10:49:50 |
version 1.8, 2004/06/12 07:29:46 |
|
|
/* $OpenXM: OpenXM/src/kan96xx/Kan/poly3.c,v 1.5 2002/02/09 06:21:02 takayama Exp $ */ |
/* $OpenXM: OpenXM/src/kan96xx/Kan/poly3.c,v 1.7 2004/02/23 09:03:42 takayama Exp $ */ |
#include <stdio.h> |
#include <stdio.h> |
#include "datatype.h" |
#include "datatype.h" |
#include "extern2.h" |
#include "extern2.h" |
Line 31 void initT(void) { |
|
Line 31 void initT(void) { |
|
DList = (int *)sGC_malloc(sizeof(int)*Lsize*N0); |
DList = (int *)sGC_malloc(sizeof(int)*Lsize*N0); |
MList = (int *)sGC_malloc(sizeof(int)*Lsize*N0); |
MList = (int *)sGC_malloc(sizeof(int)*Lsize*N0); |
*/ |
*/ |
DList = (int *)malloc(sizeof(int)*Lsize*N0); |
DList = (int *)sGC_malloc(sizeof(int)*Lsize*N0); |
MList = (int *)malloc(sizeof(int)*Lsize*N0); |
MList = (int *)sGC_malloc(sizeof(int)*Lsize*N0); |
|
|
RList = (POLY *)sGC_malloc(sizeof(POLY)*Lsize); |
RList = (POLY *)sGC_malloc(sizeof(POLY)*Lsize); |
RListRoot = (POLY *)sGC_malloc(sizeof(POLY)*Lsize); |
RListRoot = (POLY *)sGC_malloc(sizeof(POLY)*Lsize); |
Line 117 void monomialMult_diff(e,f) |
|
Line 117 void monomialMult_diff(e,f) |
|
tmp.e[i] = f->m->e[i]; |
tmp.e[i] = f->m->e[i]; |
} |
} |
if ((!weightedHomogenization) && Homogenize) { |
if ((!weightedHomogenization) && Homogenize) { |
tmp.e[0].D += EList[k]; /* homogenization. |
if (Homogenize == 3) { |
e[0].D will be added later. */ |
tmp.e[0].D += EList[k]/2; /* Double homogenization. */ |
|
tmp.e[0].x += EList[k]/2; /* Dx x = x Dx + h H */ |
|
}else{ |
|
tmp.e[0].D += EList[k]; /* homogenization. |
|
e[0].D will be added later. */ |
|
} |
}else if (weightedHomogenization && Homogenize) { |
}else if (weightedHomogenization && Homogenize) { |
tmp.e[0].D += EList[k]/2 ; /* homogenization. Weight is (1,0) (special). |
tmp.e[0].D += EList[k]/2 ; /* homogenization. Weight is (1,0) (special). |
*/ |
*/ |