home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / SMOOTH.C < prev    next >
Encoding:
C/C++ Source or Header  |  1988-02-02  |  7.9 KB  |  320 lines
  1. # include <stdio.h>
  2. # include <math.h>
  3.  
  4.  
  5.  
  6. void DataSmoothSg(float *dataset,int numdat,int smoothnum,
  7.                  int derivnum, float *smoothdata)
  8. {
  9.   int normcoef[3][5];
  10.   int smoothcoef[3][5][13];
  11.   int i;
  12.   int j;
  13.   int k;
  14.   int smoothdeg;
  15.   int startconv;
  16.   int stopconv;
  17.   int numcoef;
  18.  
  19.   normcoef[0][0] =35;
  20.   normcoef[0][1]=21;
  21.   normcoef[0][2]=231; 
  22.   normcoef[0][3]=429; 
  23.   normcoef[0][4]=143; 
  24.   normcoef[1][0]=10;
  25.   normcoef[1][1]=28;
  26.   normcoef[1][2]=60; 
  27.   normcoef[1][3]=110; 
  28.   normcoef[1][4]=182; 
  29.   normcoef[2][0]=7;
  30.   normcoef[2][1]=42; 
  31.   normcoef[2][2]=462; 
  32.   normcoef[2][3]=429; 
  33.   normcoef[2][4]=1001;
  34.   
  35.  
  36.   smoothcoef[0][0][0] = -3;
  37.   smoothcoef[0][0][1]=12;
  38.   smoothcoef[0][0][2]=17;
  39.   smoothcoef[0][0][3]=12; 
  40.   smoothcoef[0][0][4]=-3;
  41.   smoothcoef[0][0][5]=0;
  42.   smoothcoef[0][0][6]=0; 
  43.   smoothcoef[0][0][7]=0; 
  44.   smoothcoef[0][0][8]=0; 
  45.   smoothcoef[0][0][9]=0; 
  46.   smoothcoef[0][0][10]=0; 
  47.   smoothcoef[0][0][11]=0; 
  48.   smoothcoef[0][0][12]=0;
  49.   
  50.  
  51.   smoothcoef[0][1][0]=-2; 
  52.   smoothcoef[0][1][1]=3; 
  53.   smoothcoef[0][1][2]=6; 
  54.   smoothcoef[0][1][3]=7; 
  55.   smoothcoef[0][1][4]=6; 
  56.   smoothcoef[0][1][5]=3; 
  57.   smoothcoef[0][1][6]=-2; 
  58.   smoothcoef[0][1][7]=0; 
  59.   smoothcoef[0][1][8]=0; 
  60.   smoothcoef[0][1][9]=0;
  61.   smoothcoef[0][1][11]=0; 
  62.   smoothcoef[0][1][11]=0; 
  63.   smoothcoef[0][1][12]=0;
  64.  
  65.   smoothcoef[0][2][0]=-21;
  66.   smoothcoef[0][2][1]=14;
  67.   smoothcoef[0][2][2]=39; 
  68.   smoothcoef[0][2][3]=54; 
  69.   smoothcoef[0][2][4]=59; 
  70.   smoothcoef[0][2][5]=54; 
  71.   smoothcoef[0][2][6]=39; 
  72.   smoothcoef[0][2][7]=14; 
  73.   smoothcoef[0][2][8]=-21; 
  74.   smoothcoef[0][2][9]=0; 
  75.   smoothcoef[0][2][10]=0; 
  76.   smoothcoef[0][2][11]=0; 
  77.   smoothcoef[0][2][12]=0;
  78.  
  79.  
  80.  
  81.   smoothcoef[0][3][0]=-36; 
  82.   smoothcoef[0][3][1]=9; 
  83.   smoothcoef[0][3][2]=44;
  84.   smoothcoef[0][3][3]=69; 
  85.   smoothcoef[0][3][4]=84;
  86.   smoothcoef[0][3][5]=89; 
  87.   smoothcoef[0][3][6]=84;
  88.   smoothcoef[0][3][7]=69; 
  89.   smoothcoef[0][3][8]=44; 
  90.   smoothcoef[0][3][9]=9; 
  91.   smoothcoef[0][3][10]=-36; 
  92.   smoothcoef[0][3][11]=0; 
  93.   smoothcoef[0][3][12]=0; 
  94.  
  95.  
  96.   smoothcoef[0][4][0]=-11; 
  97.   smoothcoef[0][4][1]=0; 
  98.   smoothcoef[0][4][2]=9; 
  99.   smoothcoef[0][4][3]=16; 
  100.   smoothcoef[0][4][4]=21; 
  101.   smoothcoef[0][4][5]=24; 
  102.   smoothcoef[0][4][6]=25; 
  103.   smoothcoef[0][4][7]=24; 
  104.   smoothcoef[0][4][8]=21; 
  105.   smoothcoef[0][4][9]=16;
  106.   smoothcoef[0][4][10]=9;
  107.   smoothcoef[0][4][11]=0; 
  108.   smoothcoef[0][4][12]=-11;
  109.  
  110.  
  111.   smoothcoef[1][0][0]=-2; 
  112.   smoothcoef[1][0][1]=-1; 
  113.   smoothcoef[1][0][2]=0; 
  114.   smoothcoef[1][0][3]=1; 
  115.   smoothcoef[1][0][4]=2; 
  116.   smoothcoef[1][0][5]=0; 
  117.   smoothcoef[1][0][6]=0; 
  118.   smoothcoef[1][0][7]=0; 
  119.   smoothcoef[1][0][8]=0; 
  120.   smoothcoef[1][0][9]=0; 
  121.   smoothcoef[1][0][10]=0;
  122.   smoothcoef[1][0][11]=0; 
  123.   smoothcoef[1][0][12]=0;
  124.  
  125.   smoothcoef[1][1][0]=-3;
  126.   smoothcoef[1][1][1]=-2; 
  127.   smoothcoef[1][1][2]=-1; 
  128.   smoothcoef[1][1][3]=0; 
  129.   smoothcoef[1][1][4]=1; 
  130.   smoothcoef[1][1][5]=2; 
  131.   smoothcoef[1][1][6]=3; 
  132.   smoothcoef[1][1][7]=0; 
  133.   smoothcoef[1][1][8]=0;
  134.   smoothcoef[1][1][9]=0; 
  135.   smoothcoef[1][1][10]=0; 
  136.   smoothcoef[1][1][11]=0; 
  137.   smoothcoef[1][1][12]=0;
  138.  
  139.   smoothcoef[1][2][0]=-4; 
  140.   smoothcoef[1][2][1]=-3; 
  141.   smoothcoef[1][2][2]=-2; 
  142.   smoothcoef[1][2][3]=-1; 
  143.   smoothcoef[1][2][4]=0; 
  144.   smoothcoef[1][2][5]=1; 
  145.   smoothcoef[1][2][6]=2; 
  146.   smoothcoef[1][2][7]=3; 
  147.   smoothcoef[1][2][8]=4; 
  148.   smoothcoef[1][2][0]=0; 
  149.   smoothcoef[1][2][10]=0; 
  150.   smoothcoef[1][2][11]=0; 
  151.   smoothcoef[1][2][12]=0;
  152.   
  153.    
  154.   smoothcoef[1][3][0]=-5; 
  155.   smoothcoef[1][3][1]=-4; 
  156.   smoothcoef[1][3][2]=-3;
  157.   smoothcoef[1][3][3]=-2; 
  158.   smoothcoef[1][3][4]=-1; 
  159.   smoothcoef[1][3][5]=0; 
  160.   smoothcoef[1][3][6]=1; 
  161.   smoothcoef[1][3][7]=2; 
  162.   smoothcoef[1][3][8]=3; 
  163.   smoothcoef[1][3][9]=4; 
  164.   smoothcoef[1][3][10]=5; 
  165.   smoothcoef[1][3][11]=0; 
  166.   smoothcoef[1][3][12]=0;
  167.   
  168.    
  169.   smoothcoef[1][4][0]=-6; 
  170.   smoothcoef[1][4][1]=-5; 
  171.   smoothcoef[1][4][2]=-4; 
  172.   smoothcoef[1][4][3]=-3; 
  173.   smoothcoef[1][4][4]=-2; 
  174.   smoothcoef[1][4][5]=-1; 
  175.   smoothcoef[1][4][6]=0; 
  176.   smoothcoef[1][4][7]=1; 
  177.   smoothcoef[1][4][8]=2; 
  178.   smoothcoef[1][4][9]=3; 
  179.   smoothcoef[1][4][10]=4;
  180.   smoothcoef[1][4][11]=5; 
  181.   smoothcoef[1][4][12]=6; 
  182.   
  183.   
  184.   smoothcoef[2][0][0]=2; 
  185.   smoothcoef[2][0][1]=-1; 
  186.   smoothcoef[2][0][2]=-2; 
  187.   smoothcoef[2][0][3]=-1; 
  188.   smoothcoef[2][0][4]=2; 
  189.   smoothcoef[2][0][5]=0; 
  190.   smoothcoef[2][0][6]=0; 
  191.   smoothcoef[2][0][7]=0; 
  192.   smoothcoef[2][0][8]=0; 
  193.   smoothcoef[2][0][9]=0; 
  194.   smoothcoef[2][0][10]=0; 
  195.   smoothcoef[2][0][11]=0; 
  196.   smoothcoef[2][0][12]=0;
  197.   
  198.    
  199.   smoothcoef[2][1][0]=5; 
  200.   smoothcoef[2][1][1]=0; 
  201.   smoothcoef[2][1][2]=-3; 
  202.   smoothcoef[2][1][3]=-4;
  203.   smoothcoef[2][1][4]=-3; 
  204.   smoothcoef[2][1][5]=0; 
  205.   smoothcoef[2][1][6]=5; 
  206.   smoothcoef[2][1][7]=0; 
  207.   smoothcoef[2][1][8]=0; 
  208.   smoothcoef[2][1][9]=0; 
  209.   smoothcoef[2][1][10]=0; 
  210.   smoothcoef[2][1][11]=0; 
  211.   smoothcoef[2][1][12]=0;
  212.   
  213.    
  214.   smoothcoef[2][2][0]=28; 
  215.   smoothcoef[2][2][1]=7; 
  216.   smoothcoef[2][2][2]=-8; 
  217.   smoothcoef[2][2][3]=-17; 
  218.   smoothcoef[2][2][4]=-20; 
  219.   smoothcoef[2][2][5]=-17; 
  220.   smoothcoef[2][2][6]=-8; 
  221.   smoothcoef[2][2][7]=7; 
  222.   smoothcoef[2][2][8]=28; 
  223.   smoothcoef[2][2][9]=0; 
  224.   smoothcoef[2][2][10]=0; 
  225.   smoothcoef[2][2][11]=0;
  226.   smoothcoef[2][2][12]=0;
  227.   
  228.    
  229.   smoothcoef[2][3][0]=15; 
  230.   smoothcoef[2][3][1]=6; 
  231.   smoothcoef[2][3][2]=-1; 
  232.   smoothcoef[2][3][3]=-6; 
  233.   smoothcoef[2][3][4]=-9; 
  234.   smoothcoef[2][3][5]=-10; 
  235.   smoothcoef[2][3][6]=-9; 
  236.   smoothcoef[2][3][7]=-6; 
  237.   smoothcoef[2][3][8]=-1; 
  238.   smoothcoef[2][3][9]=6; 
  239.   smoothcoef[2][3][10]=15; 
  240.   smoothcoef[2][3][11]=0; 
  241.   smoothcoef[2][3][12]=0; 
  242.   
  243.   
  244.   smoothcoef[2][4][0]=22; 
  245.   smoothcoef[2][4][1]=11; 
  246.   smoothcoef[2][4][2]=2; 
  247.   smoothcoef[2][4][3]=-5; 
  248.   smoothcoef[2][4][4]=-10;
  249.   smoothcoef[2][4][5]=-13;
  250.   smoothcoef[2][4][6]=-14;
  251.   smoothcoef[2][4][7]=-13;
  252.   smoothcoef[2][4][8]=-10;
  253.   smoothcoef[2][4][9]=-5;
  254.   smoothcoef[2][4][10]=2;
  255.   smoothcoef[2][4][11]=11;
  256.   smoothcoef[2][4][12]=22;
  257.  
  258.      if ( (smoothnum >= 1) && (smoothnum <= 5) && (derivnum >= 0) && (derivnum <= 2) ) {
  259.       numcoef = 2 * smoothnum + 3;
  260.       smoothdeg = (numcoef - 1) / 2;
  261.       startconv = smoothdeg;
  262.       stopconv = numdat - smoothdeg - 1;
  263.       for ( i = 0; i <= numdat-1; ++i ) {
  264.          smoothdata[i] = 0.0;
  265.       }
  266.       for ( i = startconv; i <= stopconv; ++i ) {
  267.          for ( j = 0; j <= numcoef - 1; ++j ) {
  268.             smoothdata[i] = smoothdata[i] +
  269.      dataset[i - smoothdeg + j] * smoothcoef[derivnum][smoothnum-1][j];
  270.    
  271.          }
  272.          smoothdata[i] = smoothdata[i] / normcoef[derivnum][smoothnum-1];
  273.       }
  274.       j = smoothdeg-1;
  275.       k = numdat - smoothdeg;
  276.       for ( i = 0; i <= smoothdeg-1; ++i ) {
  277.          smoothdata[j] = smoothdata[j + 1];
  278.          smoothdata[k] = smoothdata[k - 1];
  279.          j = j - 1;
  280.          k = k + 1;
  281.       }
  282.    }
  283.  }
  284.  
  285. void DataSmoothWeights(float *dataset,int numdat,int smoothnum,
  286.                        float *weights,float wdivisor,float *smoothdata)
  287. {
  288.    int i;
  289.    int j;
  290.    int k;
  291.    int smoothdeg;
  292.    int startconv;
  293.    int stopconv;
  294.    int numcoef;
  295.    if ( (smoothnum >= 1) && (fmod(smoothnum,2) != 0) ) {
  296.       smoothdeg = smoothnum / 2;
  297.       startconv = smoothdeg;
  298.       stopconv = numdat - smoothdeg -1;
  299.       for ( i = 0; i <= numdat-1; ++i ) {
  300.          smoothdata[i] = 0;
  301.       }
  302.       for ( i = startconv; i <= stopconv; ++i ) {
  303.          for ( j = 0; j <= smoothnum - 1; ++j ) {
  304.             smoothdata[i] = smoothdata[i] +
  305.             dataset[i - smoothdeg + j] * weights[j];
  306.          }
  307.          smoothdata[i] = smoothdata[i] / wdivisor;
  308.       }
  309.       j = smoothdeg-1;
  310.       k = numdat - smoothdeg;
  311.       for ( i = 1; i <= smoothdeg; ++i ) {
  312.          smoothdata[j] = smoothdata[j + 1];
  313.          smoothdata[k] = smoothdata[k - 1];
  314.          j = j - 1;
  315.          k = k + 1;
  316.       }
  317.    }
  318. }
  319.  
  320.