--- DarthVader/TrackerLevel2/src/F77/tricircle.f 2009/08/04 14:01:40 1.4 +++ DarthVader/TrackerLevel2/src/F77/tricircle.f 2014/01/16 15:30:00 1.5 @@ -31,12 +31,16 @@ integer npoints !fit number of points real dep(npoints),indep(npoints) !dependent and independent variables - real angle(npoints) !angle between the tangent line in the input points and +c real angle(npoints) !angle between the tangent line in the input points and +c ! the independent variable axis +c real residual(npoints) !residuals +c real chi !sum of squared residuals +c real xc,zc,radius !circle parameters + double precision angle(npoints) !angle between the tangent line in the input points and ! the independent variable axis - - real residual(npoints) !residuals - real chi !sum of squared residuals - real xc,zc,radius !circle parameters + double precision residual(npoints) !residuals EM GCC4.7 + double precision chi !sum of squared residuals EM GCC4.7 + double precision xc,zc,radius !circle parameters EM GCC4.7 integer eflag !error flag =1 if the procedure fails @@ -51,7 +55,9 @@ integer ifail !=-1 if singular matrix error, =0 if not singular integer jfail !=0 if determinant can be evaluated, =-1 if determinat is probably too small, =+1 if too large - parameter (ibig=1.e4) !just a number greater than C(npoints,3) + integer ibig ! EM GCC4.7 +c parameter (ibig=1.e4) !just a number greater than C(npoints,3) + parameter (ibig=10000) !just a number greater than C(npoints,3) EM GCC 4.7 double precision xxc(ibig),zzc(ibig),rrr(ibig) !centres and radii to be averaged double precision tmp1,tmp2,tmp(npoints) !temp variables