@@ -247,7 +247,6 @@ typedef struct Namelist Namelist;
247247#define s_copy (A ,B ,C ,D ) { int __i,__m; for (__i=0, __m=f2cmin((C),(D)); __i<__m && (B)[__i] != 0; ++__i) (A)[__i] = (B)[__i]; }
248248#define sig_die (s , kill ) { exit(1); }
249249#define s_stop (s , n ) {exit(0);}
250- static char junk [] = "\n@(#)LIBF77 VERSION 19990503\n" ;
251250#define z_abs (z ) (cabs(Cd(z)))
252251#define z_exp (R , Z ) {pCd(R) = cexp(Cd(Z));}
253252#define z_sqrt (R , Z ) {pCd(R) = csqrt(Cd(Z));}
@@ -256,7 +255,7 @@ static char junk[] = "\n@(#)LIBF77 VERSION 19990503\n";
256255#define myceiling (w ) {ceil(w)}
257256#define myhuge (w ) {HUGE_VAL}
258257//#define mymaxloc_(w,s,e,n) {if (sizeof(*(w)) == sizeof(double)) dmaxloc_((w),*(s),*(e),n); else dmaxloc_((w),*(s),*(e),n);}
259- #define mymaxloc (w ,s ,e ,n ) {dmaxloc_(w,*(s),*(e),n)}
258+ // #define mymaxloc(w,s,e,n) {dmaxloc_(w,*(s),*(e),n)}
260259
261260/* procedure parameter types for -A and -C++ */
262261
@@ -267,6 +266,7 @@ typedef logical (*L_fp)(...);
267266typedef logical (* L_fp )();
268267#endif
269268
269+ #if 0
270270static float spow_ui (float x , integer n ) {
271271 float pow = 1.0 ; unsigned long int u ;
272272 if (n != 0 ) {
@@ -378,6 +378,7 @@ static integer smaxloc_(float *w, integer s, integer e, integer *n)
378378 if (w [i - 1 ]> m ) mi = i ,m = w [i - 1 ];
379379 return mi - s + 1 ;
380380}
381+ #endif
381382
382383static inline void cdotc_ (complex * z , integer * n_ , complex * x , integer * incx_ , complex * y , integer * incy_ ) {
383384 integer n = * n_ , incx = * incx_ , incy = * incy_ , i ;
@@ -411,7 +412,7 @@ static inline void cdotc_(complex *z, integer *n_, complex *x, integer *incx_, c
411412 zdotc += conjf (Cf (& x [i * incx ])) * Cf (& y [i * incy ]);
412413 }
413414 }
414- pCf ( z ) = zdotc ;
415+ z -> r = creal ( zdotc ); z -> i = cimag ( zdotc ) ;
415416}
416417#endif
417418static inline void zdotc_ (doublecomplex * z , integer * n_ , doublecomplex * x , integer * incx_ , doublecomplex * y , integer * incy_ ) {
0 commit comments