Changeset 19


Ignore:
Timestamp:
30.10.2008 21:04:31 (11 years ago)
Author:
marcus
Message:

r4 | svn | 2008-01-15 14:21:11 +0100 (Di, 15 Jan 2008) | 1 line

two SQR macros in one command line will not work in gcc 4 any more

File:
1 edited

Legend:

Unmodified
Added
Removed
  • SH_SHM/trunk/source/numres.c

    r16 r19  
    33 *      ======== 
    44 * 
    5  * version 15, 22-May-2006 
     5 * version 16, 15-Jan-2008 
    66 * 
    77 * routines from numerical recipes 
     
    482482                for (j=1;j<=(n-k);j++) { 
    483483                        num += wk1[j]*wk2[j]; 
    484                         denom += SQR(wk1[j])+SQR(wk2[j]); 
     484                        /*denom += SQR(wk1[j])+SQR(wk2[j]); K.S.:dont trust two SQRs*/ 
     485                        denom += wk1[j]*wk1[j] + wk2[j]*wk2[j]; 
    485486                } 
    486487                cof[k]=2.0*num/denom; 
     
    576577                } 
    577578                nr_four1(w1,mm,1); 
    578                 p[1] += (SQR(w1[1])+SQR(w1[2])); 
     579                /*p[1] += (SQR(w1[1])+SQR(w1[2]));*/ /* K.S. doesn't work with gcc 4 */ 
     580                p[1] += ( w1[1]*w1[1] + w1[2]*w1[2] ); 
    579581                for (j=2;j<=m;j++) { 
    580582                        j2=j+j; 
    581                         p[j] += (SQR(w1[j2])+SQR(w1[j2-1]) 
    582                                 +SQR(w1[m44-j2])+SQR(w1[m43-j2])); 
     583                        p[j] += ( w1[j2]*w1[j2] + w1[j2-1]*w1[j2-1] 
     584                                + w1[m44-j2]*w1[m44-j2] + w1[m43-j2]*w1[m43-j2] ); 
    583585                } 
    584586                den += sumw; 
Note: See TracChangeset for help on using the changeset viewer.