[BACK]Return to poly3.c CVS log [TXT][DIR] Up to [local] / OpenXM / src / kan96xx / Kan

Diff for /OpenXM/src/kan96xx/Kan/poly3.c between version 1.6 and 1.9

version 1.6, 2002/09/08 10:49:50 version 1.9, 2020/10/06 11:33:46
Line 1 
Line 1 
 /* $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.8 2004/06/12 07:29:46 takayama Exp $ */
 #include <stdio.h>  #include <stdio.h>
 #include "datatype.h"  #include "datatype.h"
 #include "extern2.h"  #include "extern2.h"
   
   void outputTable(void);
   
 int Homogenize = 1;  int Homogenize = 1;
   
 #define I(i,j) (i*N0+j)  #define I(i,j) (i*N0+j)
Line 31  void initT(void) {
Line 33  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 119  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).
                                                                   */                                                                    */
Line 257  POLY mpMult_difference_org(POLY f,POLY g)
Line 264  POLY mpMult_difference_org(POLY f,POLY g)
   return(r);    return(r);
 }  }
   
 outputTable() {  void outputTable() {
   int i,j;    int i,j;
   printf("Maxv = %d Plist=%d\n",Maxv,Plist);    printf("Maxv = %d Plist=%d\n",Maxv,Plist);
   for (i=0; i<Maxv; i++) printf("%5d",V[i]);    for (i=0; i<Maxv; i++) printf("%5d",V[i]);
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.9
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>