--- nec2c.c.orig 2004-02-23 09:30:55 UTC +++ nec2c.c @@ -794,7 +794,7 @@ int main( int argc, char **argv ) continue; /* continue card input loop */ ntyp[idx]=3; - x11r[idx]=- tmp1; + x11r[idx]= - tmp1; continue; /* continue card input loop */ } @@ -1218,7 +1218,7 @@ int main( int argc, char **argv ) if( iperf != 1) { if( sig < 0.) - sig=- sig/(59.96*wlam); + sig= - sig/(59.96*wlam); epsc= cmplx( epsr, -sig*wlam*59.96); zrati=1./ csqrt( epsc); @@ -3197,7 +3197,7 @@ void conect( int ignd ) ix= jco[i]; if( ix <= 0) { - ix=- ix; + ix= - ix; ixx = ix-1; x[ixx]= xa; y[ixx]= ya; @@ -3931,7 +3931,7 @@ void efld( double xi, double yi, double zi, double ai, { if( ip == 1) ijx=1; - rfl=- rfl; + rfl= - rfl; salpr= salpj* rfl; zij= zi- rfl* zj; zp= xij* cabj+ yij* sabj+ zij* salpr; @@ -4032,7 +4032,7 @@ void efld( double xi, double yi, double zi, double ai, } else { - px=- yij/ xymag; + px= - yij/ xymag; py= xij/ xymag; cth= zij/ rmag; zrsin= csqrt(1.0 - zratx*zratx*(1.0 - cth*cth) ); @@ -4107,7 +4107,7 @@ void efld( double xi, double yi, double zi, double ai, /* displace observation point for thin wire approximation */ zij= zi+ zj; - salpr=- salpj; + salpr= - salpj; rhox= sabj* zij- salpr* yij; rhoy= salpr* xij- cabj* zij; rhoz= cabj* yij- sabj* xij; @@ -4123,7 +4123,7 @@ void efld( double xi, double yi, double zi, double ai, { rh= ai/ sqrt( rh); if( rhoz < 0.) - rh=- rh; + rh= - rh; xo= xi+ rh* rhox; yo= yi+ rh* rhoy; zo= zi+ rh* rhoz; @@ -4215,18 +4215,18 @@ void eksc( double s, double z, double rh, double xk, i gzp1= gp1* z1a; gzp2= gp2* z2a; *ezs= CONST1*(( gz2- gz1)* cs* xk-( gzp2+ gzp1)* ss); - *ezc=- CONST1*(( gz2+ gz1)* ss* xk+( gzp2- gzp1)* cs); + *ezc= - CONST1*(( gz2+ gz1)* ss* xk+( gzp2- gzp1)* cs); *erk= CONST1*( gp2- gp1)* rh; intx(- shk, shk, rhk, ij, &cint, &sint); - *ezk=- CONST1*( gzp2- gzp1+ xk* xk* cmplx( cint,- sint)); + *ezk= - CONST1*( gzp2- gzp1+ xk* xk* cmplx( cint,- sint)); gzp1= gzp1* z1a; gzp2= gzp2* z2a; if( rh >= 1.0e-10) { - *ers=- CONST1*(( gzp2+ gzp1+ gz2+ gz1)* + *ers= - CONST1*(( gzp2+ gzp1+ gz2+ gz1)* ss-( z2a* gz2- z1a* gz1)* cs*xk)/ rh; - *erc=- CONST1*(( gzp2- gzp1+ gz2- gz1)* + *erc= - CONST1*(( gzp2- gzp1+ gz2- gz1)* cs+( z2a* gz2+ z1a* gz1)* ss*xk)/ rh; return; } @@ -4304,16 +4304,16 @@ void ekscx( double bx, double s, double z, } *ezs= CONST1*(( gz2- gz1)* cs* xk-( gzp2+ gzp1)* ss); - *ezc=- CONST1*(( gz2+ gz1)* ss* xk+( gzp2- gzp1)* cs); - *ers=- CONST1*(( z2a* grp2+ z1a* grp1+ gr2+ gr1)*ss + *ezc= - CONST1*(( gz2+ gz1)* ss* xk+( gzp2- gzp1)* cs); + *ers= - CONST1*(( z2a* grp2+ z1a* grp1+ gr2+ gr1)*ss -( z2a* gr2- z1a* gr1)* cs* xk); - *erc=- CONST1*(( z2a* grp2- z1a* grp1+ gr2- gr1)*cs + *erc= - CONST1*(( z2a* grp2- z1a* grp1+ gr2- gr1)*cs +( z2a* gr2+ z1a* gr1)* ss* xk); *erk= CONST1*( grk2- grk1); intx(- shk, shk, rhk, ij, &cint, &sint); bk= b* xk; bk2= bk* bk*.25; - *ezk=- CONST1*( gzp2- gzp1+ xk* xk*(1.- bk2)* + *ezk= - CONST1*( gzp2- gzp1+ xk* xk*(1.- bk2)* cmplx( cint,- sint)-bk2*( gzz2- gzz1)); return; @@ -4373,10 +4373,10 @@ void etmns( double p1, double p2, double p3, double p4 set= sin( p3); pxl= cth* cph* cet- sph* set; pyl= cth* sph* cet+ cph* set; - pzl=- sth* cet; - wx=- sth* cph; - wy=- sth* sph; - wz=- cth; + pzl= - sth* cet; + wx= - sth* cph; + wy= - sth* sph; + wz= - cth; qx= wy* pzl- wz* pyl; qy= wz* pxl- wx* pzl; qz= wx* pyl- wy* pxl; @@ -4405,7 +4405,7 @@ void etmns( double p1, double p2, double p3, double p4 { for( i = 0; i < n; i++ ) { - arg=- TP*( wx* x[i]+ wy* y[i]+ wz* z[i]); + arg= - TP*( wx* x[i]+ wy* y[i]+ wz* z[i]); e[i]=-( pxl* cab[i]+ pyl* sab[i]+ pzl* salp[i])* cmplx( cos( arg), sin( arg)); } @@ -4415,11 +4415,11 @@ void etmns( double p1, double p2, double p3, double p4 tt1=( pyl* cph- pxl* sph)*( rrh- rrv); cx= rrv* pxl- tt1* sph; cy= rrv* pyl+ tt1* cph; - cz=- rrv* pzl; + cz= - rrv* pzl; for( i = 0; i < n; i++ ) { - arg=- TP*( wx* x[i]+ wy* y[i]- wz* z[i]); + arg= - TP*( wx* x[i]+ wy* y[i]- wz* z[i]); e[i]= e[i]-( cx* cab[i]+ cy* sab[i]+ cz* salp[i])* cmplx(cos( arg), sin( arg)); } @@ -4438,7 +4438,7 @@ void etmns( double p1, double p2, double p3, double p4 i++; i1 += 2; i2 = i1+1; - arg=- TP*( wx* px[i]+ wy* py[i]+ wz* pz[i]); + arg= - TP*( wx* px[i]+ wy* py[i]+ wz* pz[i]); tt1= cmplx( cos( arg), sin( arg))* psalp[i]* RETA; e[i2]=( qx* t1x[i]+ qy* t1y[i]+ qz* t1z[i])* tt1; e[i1]=( qx* t2x[i]+ qy* t2y[i]+ qz* t2z[i])* tt1; @@ -4459,7 +4459,7 @@ void etmns( double p1, double p2, double p3, double p4 i++; i1 += 2; i2 = i1+1; - arg=- TP*( wx* px[i]+ wy* py[i]- wz* pz[i]); + arg= - TP*( wx* px[i]+ wy* py[i]- wz* pz[i]); tt1= cmplx( cos( arg), sin( arg))* psalp[i]* RETA; e[i2]= e[i2]+( cx* t1x[i]+ cy* t1y[i]+ cz* t1z[i])* tt1; e[i1]= e[i1]+( cx* t2x[i]+ cy* t2y[i]+ cz* t2z[i])* tt1; @@ -4471,7 +4471,7 @@ void etmns( double p1, double p2, double p3, double p4 /* incident plane wave, elliptic polarization. */ tt1=-(CPLX_01)* p6; if( ipr == 3) - tt1=- tt1; + tt1= - tt1; if( n != 0) { @@ -4481,7 +4481,7 @@ void etmns( double p1, double p2, double p3, double p4 for( i = 0; i < n; i++ ) { - arg=- TP*( wx* x[i]+ wy* y[i]+ wz* z[i]); + arg= - TP*( wx* x[i]+ wy* y[i]+ wz* z[i]); e[i]=-( cx* cab[i]+ cy* sab[i]+ cz* salp[i])* cmplx( cos( arg), sin( arg)); } @@ -4491,11 +4491,11 @@ void etmns( double p1, double p2, double p3, double p4 tt2=( cy* cph- cx* sph)*( rrh- rrv); cx= rrv* cx- tt2* sph; cy= rrv* cy+ tt2* cph; - cz=- rrv* cz; + cz= - rrv* cz; for( i = 0; i < n; i++ ) { - arg=- TP*( wx* x[i]+ wy* y[i]- wz* z[i]); + arg= - TP*( wx* x[i]+ wy* y[i]- wz* z[i]); e[i]= e[i]-( cx* cab[i]+ cy* sab[i]+ cz* salp[i])* cmplx(cos( arg), sin( arg)); } @@ -4518,7 +4518,7 @@ void etmns( double p1, double p2, double p3, double p4 i++; i1 += 2; i2 = i1+1; - arg=- TP*( wx* px[i]+ wy* py[i]+ wz* pz[i]); + arg= - TP*( wx* px[i]+ wy* py[i]+ wz* pz[i]); tt2= cmplx( cos( arg), sin( arg))* psalp[i]* RETA; e[i2]=( cx* t1x[i]+ cy* t1y[i]+ cz* t1z[i])* tt2; e[i1]=( cx* t2x[i]+ cy* t2y[i]+ cz* t2z[i])* tt2; @@ -4539,7 +4539,7 @@ void etmns( double p1, double p2, double p3, double p4 i++; i1 += 2; i2 = i1+1; - arg=- TP*( wx* px[i]+ wy* py[i]- wz* pz[i]); + arg= - TP*( wx* px[i]+ wy* py[i]- wz* pz[i]); tt1= cmplx( cos( arg), sin( arg))* psalp[i]* RETA; e[i2]= e[i2]+( cx* t1x[i]+ cy* t1y[i]+ cz* t1z[i])* tt1; e[i1]= e[i1]+( cx* t2x[i]+ cy* t2y[i]+ cz* t2z[i])* tt1; @@ -4604,7 +4604,7 @@ void etmns( double p1, double p2, double p3, double p4 } /* if( arg >= 1.e-30) */ - arg=- TP* r; + arg= - TP* r; tt1= cmplx( cos( arg), sin( arg)); if( i < n ) @@ -4775,7 +4775,7 @@ complex double fbar( complex double p ) for( i = 1; i <= 100; i++ ) { - pow=- pow* zs/ (double)i; + pow= - pow* zs/ (double)i; term= pow/(2.* i+1.); sum= sum+ term; tms= creal( term* conj( term)); @@ -4793,7 +4793,7 @@ complex double fbar( complex double p ) if( creal( z) < 0.) { minus=1; - z=- z; + z= - z; } else minus=0; @@ -4804,13 +4804,13 @@ complex double fbar( complex double p ) for( i = 1; i <= 6; i++ ) { - term =- term*(2.*i -1.)* zs; + term = - term*(2.*i -1.)* zs; sum += term; } if( minus == 1) sum -= 2.* SP* z* cexp( z* z); - fbar=- sum; + fbar= - sum; return( fbar ); } @@ -4882,7 +4882,7 @@ void fblock( int nrow, int ncol, int imax, int ipsym ) { deter= ssx[i+j*nop]; ssx[i+(j+kk)*nop]= deter; - ssx[i+kk+(j+kk)*nop]=- deter; + ssx[i+kk+(j+kk)*nop]= - deter; ssx[i+kk+j*nop]= deter; } } @@ -4908,14 +4908,14 @@ void ffld( double thet, double phi, complex double zrsin, rrv, rrh, rrv1, rrh1, rrv2, rrh2; complex double tix, tiy, tiz, zscrn, ex, ey, ez, gx, gy, gz; - phx=- sin( phi); + phx= - sin( phi); phy= cos( phi); roz= cos( thet); rozs= roz; thx= roz* phy; - thy=- roz* phx; - thz=- sin( thet); - rox=- thz* phy; + thy= - roz* phx; + thz= - sin( thet); + rox= - thz* phy; roy= thz* phx; jump = FALSE; @@ -4954,11 +4954,11 @@ void ffld( double thet, double phi, zrsin= csqrt(1.- zrati2* zrati2* thz* thz); rrv2=-( roz- zrati2* zrsin)/( roz+ zrati2* zrsin); rrh2=( zrati2* roz- zrsin)/( zrati2* roz+ zrsin); - darg=- TP*2.* ch* roz; + darg= - TP*2.* ch* roz; } } /* if( ifar > 1) */ - roz=- roz; + roz= - roz; ccx= cix; ccy= ciy; ccz= ciz; @@ -5147,7 +5147,7 @@ void ffld( double thet, double phi, rfl=-1.; for( ip = 0; ip < ksymp; ip++ ) { - rfl=- rfl; + rfl= - rfl; rrz= roz* rfl; fflds( rox, roy, rrz, &cur[n], &gx, &gy, &gz); @@ -5161,9 +5161,9 @@ void ffld( double thet, double phi, if( iperf == 1) { - gx=- gx; - gy=- gy; - gz=- gz; + gx= - gx; + gy= - gy; + gz= - gz; } else { @@ -5285,7 +5285,7 @@ void gfld( double rho, double phi, double rz, thet= PI*.5; ffld( thet, phi, eth, epi); - arg=- TP* r; + arg= - TP* r; exa= cmplx( cos( arg), sin( arg))/ r; *eth= *eth* exa; *epi= *epi* exa; @@ -5296,10 +5296,10 @@ void gfld( double rho, double phi, double rz, /* computation of space and ground waves. */ u= ux; u2= u* u; - phx=- sin( phi); + phx= - sin( phi); phy= cos( phi); rx= rho* phy; - ry=- rho* phx; + ry= - rho* phx; cix=CPLX_00; ciy=CPLX_00; ciz=CPLX_00; @@ -5348,7 +5348,7 @@ void gfld( double rho, double phi, double rz, /* constant, sine, and cosine current distributions */ for( k = 0; k < 2; k++ ) { - rfl=- rfl; + rfl= - rfl; riz= rz- z[i]* rfl; rxyz= sqrt( rix* rix+ riy* riy+ riz* riz); rnx= rix/ rxyz; @@ -5409,7 +5409,7 @@ void gfld( double rho, double phi, double rz, } /* for( i = 0; i < n; i++ ) */ - arg=- TP* r; + arg= - TP* r; exa= cmplx( cos( arg), sin( arg)); cix= cix* exa; ciy= ciy* exa; @@ -5418,8 +5418,8 @@ void gfld( double rho, double phi, double rz, rny= ry/ r; rnz= rz/ r; thx= rnz* phy; - thy=- rnz* phx; - thz=- rho/ r; + thy= - rnz* phx; + thz= - rho/ r; *eth= cix* thx+ ciy* thy+ ciz* thz; *epi= cix* phx+ ciy* phy; *erd= cix* rnx+ ciy* rny+ ciz* rnz; @@ -5477,8 +5477,8 @@ void gwave( complex double *erv, complex double *ezv, cpp2=1.0e-20; cpp= sqrt( cpp2); - rk1=- TPJ* r1; - rk2=- TPJ* r2; + rk1= - TPJ* r1; + rk2= - TPJ* r2; t1=1. -u2* cpp2; t2= csqrt( t1); t3=(1. -1./ rk1)/ rk1; @@ -5542,7 +5542,7 @@ void gx( double zz, double rh, double xk, r= sqrt( r2); rkz= xk* r; *gz= cmplx( cos( rkz),- sin( rkz))/ r; - *gzp=- cmplx(1.0, rkz)* *gz/ r2; + *gzp= - cmplx(1.0, rkz)* *gz/ r2; return; } @@ -5580,7 +5580,7 @@ void gxx( double zz, double rh, double a, double a2, d if( ira != 1) { *g3=( *g3+ *gzp)* rh; - *gzp=- zz* c1* gz; + *gzp= - zz* c1* gz; if( rh <= 1.0e-10) { @@ -5596,11 +5596,11 @@ void gxx( double zz, double rh, double a, double a2, d } /* if( ira != 1) */ t2=.5* a; - *g2=- t2* c1* gz; + *g2= - t2* c1* gz; *g2p= t2* gz* c2/ r2; *g3= rh2* *g2p- a* gz* c1; *g2p= *g2p* zz; - *gzp=- zz* c1* gz; + *gzp= - zz* c1* gz; return; } @@ -5897,7 +5897,7 @@ void hintg( double xi, double yi, double zi ) for( ip = 1; ip <= ksymp; ip++ ) { - rfl=- rfl; + rfl= - rfl; rz= zi- zj* rfl; rsq= rx* rx+ ry* ry+ rz* rz; @@ -5925,12 +5925,12 @@ void hintg( double xi, double yi, double zi ) { if( iperf == 1) { - f1x=- f1x; - f1y=- f1y; - f1z=- f1z; - f2x=- f2x; - f2y=- f2y; - f2z=- f2z; + f1x= - f1x; + f1y= - f1y; + f1z= - f1z; + f2x= - f2x; + f2y= - f2y; + f2z= - f2z; } else { @@ -5944,7 +5944,7 @@ void hintg( double xi, double yi, double zi ) } else { - pxx=- ry/ xymag; + pxx= - ry/ xymag; pyy= rx/ xymag; cth= rz/ r; rrv= csqrt(1.- zrati* zrati*(1.- cth* cth)); @@ -5997,7 +5997,7 @@ void hsfld( double xi, double yi, double zi, double ai for( ip = 0; ip < ksymp; ip++ ) { - rfl=- rfl; + rfl= - rfl; salpr= salpj* rfl; zij= zi- rfl* zj; zp= xij* cabj+ yij* sabj+ zij* salpr; @@ -6061,7 +6061,7 @@ void hsfld( double xi, double yi, double zi, double ai } else { - px=- yij/ xymag; + px= - yij/ xymag; py= xij/ xymag; cth= zij/ rmag; rrv= csqrt(1.- zratx* zratx*(1.- cth* cth)); @@ -6137,7 +6137,7 @@ void hsflx( double s, double rh, double zpx, } else { - zp=- zpx; + zp= - zpx; hss=-1.; } @@ -6165,8 +6165,8 @@ void hsflx( double s, double rh, double zpx, t1= z1* ekr1/ r1; t2= z2a* ekr2/ r2; *hps=( cdk*( ekr2- ekr1)- CPLX_01* sdk*( t2+ t1))* hss; - *hpc=- sdk*( ekr2+ ekr1)- CPLX_01* cdk*( t2- t1); - cons=- CPLX_01/(2.* TP* rh); + *hpc= - sdk*( ekr2+ ekr1)- CPLX_01* cdk*( t2- t1); + cons= - CPLX_01/(2.* TP* rh); *hps= cons* *hps; *hpc= cons* *hpc; return; @@ -6807,7 +6807,7 @@ void move( double rox, double roy, double roz, double yx= sph* cth; yy= sph* sth* sps+ cph* cps; yz= sph* sth* cps- cph* sps; - zx=- sth; + zx= - sth; zy= cth* sps; zz= cth* cps; @@ -7338,15 +7338,15 @@ void netwk( complex double *cm, complex double *cmb, y12r=0.; y12i=1./( x11r[j]* sin( y22r)); y11r= x12r[j]; - y11i=- y12i* cos( y22r); + y11i= - y12i* cos( y22r); y22r= x22r[j]; y22i= y11i+ x22i[j]; y11i= y11i+ x12i[j]; if( ntyp[j] != 2) { - y12r=- y12r; - y12i=- y12i; + y12r= - y12r; + y12i= - y12i; } } /* if( ntyp[j] <= 1) */ @@ -8013,7 +8013,7 @@ void patch( int nx, int ny, if( xa >= 1.0e-6) { - t1x[mi]=- ynv/ xa; + t1x[mi]= - ynv/ xa; t1y[mi]= xnv/ xa; t1z[mi]=0.; } @@ -8256,10 +8256,10 @@ void subph( int nx, int ny ) psalp[mia]= saln; if( ix == 2) - yt=- yt; + yt= - yt; if( (ix == 1) || (ix == 3) ) - xt=- xt; + xt= - xt; mia++; } @@ -8470,7 +8470,7 @@ void qdsrc( int is, complex double v, complex double * if( ipr < 0 ) { - ipr=- ipr; + ipr= - ipr; ipr--; if( -icon1[ipr-1] != jp1 ) ind1=2; @@ -8828,7 +8828,7 @@ void rdpat( void ) stilta= sin( tilta); tstor1= tstor1* stilta* stilta; tstor2= tstor2* stilta* cos( tilta); - emajr2=- tstor1+ tstor2+ ethm2; + emajr2= - tstor1+ tstor2+ ethm2; eminr2= tstor1- tstor2+ ephm2; if( eminr2 < 0.) eminr2=0.; @@ -9148,7 +9148,7 @@ void readgm( char *gm, int *i1, int *i2, double *x1, d (line_buf[ line_idx] > '9')) && (line_buf[ line_idx] != '+') && (line_buf[ line_idx] != '-') ) - if( (line_buf[line_idx] == '\0') ) + if( line_buf[line_idx] == '\0') { *i1= iarr[0]; *i2= iarr[1]; @@ -9212,7 +9212,7 @@ void readgm( char *gm, int *i1, int *i2, double *x1, d (line_buf[ line_idx] != '+') && (line_buf[ line_idx] != '-') && (line_buf[ line_idx] != '.') ) - if( (line_buf[line_idx] == '\0') ) + if( line_buf[line_idx] == '\0') { *i1= iarr[0]; *i2= iarr[1]; @@ -9341,7 +9341,7 @@ void readmn( char *gm, int *i1, int *i2, int *i3, int (line_buf[ line_idx] > '9')) && (line_buf[ line_idx] != '+') && (line_buf[ line_idx] != '-') ) - if( (line_buf[line_idx] == '\0') ) + if( line_buf[line_idx] == '\0') { *i1= iarr[0]; *i2= iarr[1]; @@ -9407,7 +9407,7 @@ void readmn( char *gm, int *i1, int *i2, int *i3, int (line_buf[ line_idx] != '+') && (line_buf[ line_idx] != '-') && (line_buf[ line_idx] != '.') ) - if( (line_buf[line_idx] == '\0') ) + if( line_buf[line_idx] == '\0') { *i1= iarr[0]; *i2= iarr[1]; @@ -9538,10 +9538,10 @@ void reflc( int ix, int iy, int iz, int itx, int nop ) x[nx]= x[i]; y[nx]= y[i]; - z[nx]=- e1; + z[nx]= - e1; x2[nx]= x2[i]; y2[nx]= y2[i]; - z2[nx]=- e2; + z2[nx]= - e2; itagi= itag[i]; if( itagi == 0) @@ -9587,14 +9587,14 @@ void reflc( int ix, int iy, int iz, int itx, int nop ) px[nx]= px[i]; py[nx]= py[i]; - pz[nx]=- pz[i]; + pz[nx]= - pz[i]; t1x[nx]= t1x[i]; t1y[nx]= t1y[i]; - t1z[nx]=- t1z[i]; + t1z[nx]= - t1z[i]; t2x[nx]= t2x[i]; t2y[nx]= t2y[i]; - t2z[nx]=- t2z[i]; - psalp[nx]=- psalp[i]; + t2z[nx]= - t2z[i]; + psalp[nx]= - psalp[i]; pbi[nx]= pbi[i]; } @@ -9637,10 +9637,10 @@ void reflc( int ix, int iy, int iz, int itx, int nop ) } x[nx]= x[i]; - y[nx]=- e1; + y[nx]= - e1; z[nx]= z[i]; x2[nx]= x2[i]; - y2[nx]=- e2; + y2[nx]= - e2; z2[nx]= z2[i]; itagi= itag[i]; @@ -9686,15 +9686,15 @@ void reflc( int ix, int iy, int iz, int itx, int nop ) } px[nx]= px[i]; - py[nx]=- py[i]; + py[nx]= - py[i]; pz[nx]= pz[i]; t1x[nx]= t1x[i]; - t1y[nx]=- t1y[i]; + t1y[nx]= - t1y[i]; t1z[nx]= t1z[i]; t2x[nx]= t2x[i]; - t2y[nx]=- t2y[i]; + t2y[nx]= - t2y[i]; t2z[nx]= t2z[i]; - psalp[nx]=- psalp[i]; + psalp[nx]= - psalp[i]; pbi[nx]= pbi[i]; } /* for( i = m2; i <= m; i++ ) */ @@ -9738,10 +9738,10 @@ void reflc( int ix, int iy, int iz, int itx, int nop ) stop(-1); } - x[nx]=- e1; + x[nx]= - e1; y[nx]= y[i]; z[nx]= z[i]; - x2[nx]=- e2; + x2[nx]= - e2; y2[nx]= y2[i]; z2[nx]= z2[i]; itagi= itag[i]; @@ -9786,16 +9786,16 @@ void reflc( int ix, int iy, int iz, int itx, int nop ) stop(-1); } - px[nx]=- px[i]; + px[nx]= - px[i]; py[nx]= py[i]; pz[nx]= pz[i]; - t1x[nx]=- t1x[i]; + t1x[nx]= - t1x[i]; t1y[nx]= t1y[i]; t1z[nx]= t1z[i]; - t2x[nx]=- t2x[i]; + t2x[nx]= - t2x[i]; t2y[nx]= t2y[i]; t2z[nx]= t2z[i]; - psalp[nx]=- psalp[i]; + psalp[nx]= - psalp[i]; pbi[nx]= pbi[i]; } @@ -10106,11 +10106,11 @@ void sbf( int i, int is, double *aa, double *bb, doubl if( jcox != 0 ) { if( jcox < 0 ) - jcox=- jcox; + jcox= - jcox; else { - sig=- sig; - jend=- jend; + sig= - sig; + jend= - jend; } jcoxx = jcox-1; @@ -10135,7 +10135,7 @@ void sbf( int i, int is, double *aa, double *bb, doubl { *aa= aj/ sd* sig; *bb= aj/(2.* cdh); - *cc=- aj/(2.* sdh)* sig; + *cc= - aj/(2.* sdh)* sig; june= iend; } @@ -10161,14 +10161,14 @@ void sbf( int i, int is, double *aa, double *bb, doubl } /* if( jcox != i ) */ else if( jcox == is) - *bb=- *bb; + *bb= - *bb; if( iend == 1) break; } /* if( jcox != 0 ) */ - pm=- pp; + pm= - pp; pp=0.; njun1= jsno; @@ -10220,9 +10220,9 @@ void sbf( int i, int is, double *aa, double *bb, doubl if( june == 1) { - *aa=- *aa* qp; + *aa= - *aa* qp; *bb= *bb* qp; - *cc=- *cc* qp; + *cc= - *cc* qp; if( i != is) return; } @@ -10273,9 +10273,9 @@ void sbf( int i, int is, double *aa, double *bb, doubl } else { - *aa=- *aa* qp; + *aa= - *aa* qp; *bb= *bb* qp; - *cc=- *cc* qp; + *cc= - *cc* qp; } if( i != is) @@ -10319,7 +10319,7 @@ void sflds( double t, complex double *e ) { rhx= rhx/ rho; rhy= rhy/ rho; - phx=- rhy; + phx= - rhy; phy= rhx; } @@ -10668,11 +10668,11 @@ void tbf( int i, int icap ) if( jcox != 0 ) { if( jcox < 0 ) - jcox=- jcox; + jcox= - jcox; else { - sig=- sig; - jend=- jend; + sig= - sig; + jend= - jend; } jcoxx = jcox-1; @@ -10696,7 +10696,7 @@ void tbf( int i, int icap ) pp= pp- omc/ sd* aj; ax[jsnox]= aj/ sd* sig; bx[jsnox]= aj/(2.* cdh); - cx[jsnox]=- aj/(2.* sdh)* sig; + cx[jsnox]= - aj/(2.* sdh)* sig; if( jcox != i) { @@ -10719,14 +10719,14 @@ void tbf( int i, int icap ) } /* if( jcox != i) */ else - bx[jsnox] =- bx[jsnox]; + bx[jsnox] = - bx[jsnox]; if( iend == 1) break; } /* if( jcox != 0 ) */ - pm=- pp; + pm= - pp; pp=0.; njun1= jsno; @@ -10799,9 +10799,9 @@ void tbf( int i, int icap ) for( iend = 0; iend < njun2; iend++ ) { - ax[iend]=- ax[iend]* qp; + ax[iend]= - ax[iend]* qp; bx[iend]= bx[iend]* qp; - cx[iend]=- cx[iend]* qp; + cx[iend]= - cx[iend]* qp; } jsno= jsnop+1; @@ -10855,9 +10855,9 @@ void tbf( int i, int icap ) jend= njun1; for( iend = jend; iend < jsno; iend++ ) { - ax[iend]=- ax[iend]* qp; + ax[iend]= - ax[iend]* qp; bx[iend]= bx[iend]* qp; - cx[iend]=- cx[iend]* qp; + cx[iend]= - cx[iend]* qp; } jsno= jsnop+1; @@ -10947,9 +10947,9 @@ void trio( int j ) do { if( jcox < 0 ) - jcox=- jcox; + jcox= - jcox; else - jend=- jend; + jend= - jend; jcoxx = jcox-1; if( jcox != j) @@ -11036,9 +11036,9 @@ void unere( double xob, double yob, double zob ) if( ipgnd == 2) { - zr=- zr; - t1zr=- t1zr; - t2zr=- t2zr; + zr= - zr; + t1zr= - t1zr; + t2zr= - t2zr; } rx= xob- xj; @@ -11058,7 +11058,7 @@ void unere( double xob, double yob, double zob ) } r= sqrt( r2); - tt1=- TP* r; + tt1= - TP* r; tt2= tt1* tt1; rt= r2* r; er= cmplx( sin( tt1),- cos( tt1))*( CONST2* s); @@ -11078,12 +11078,12 @@ void unere( double xob, double yob, double zob ) if( iperf == 1) { - exk=- exk; - eyk=- eyk; - ezk=- ezk; - exs=- exs; - eys=- eys; - ezs=- ezs; + exk= - exk; + eyk= - eyk; + ezk= - ezk; + exs= - exs; + eys= - eys; + ezs= - ezs; return; } @@ -11097,7 +11097,7 @@ void unere( double xob, double yob, double zob ) } else { - px=- ry/ xymag; + px= - ry/ xymag; py= rx/ xymag; cth= rz/ sqrt( xymag* xymag+ rz* rz); rrv= csqrt(1.- zrati* zrati*(1.- cth* cth));