Bug#796780: paw-common: Erroneous minuit in paw in Debian 7
Mikhail Ukhanov
Mikhail.Ukhanov at ihep.ru
Mon Aug 24 12:57:23 UTC 2015
Package: paw-common
Version: 1:2.14.04.dfsg.2-9
Severity: important
Tags: newcomer
Dear Maintainer,
The following kumac gives erroneous result in Debian 7 (wheezy)
--------------------------------------------------------------------------------------------------------------------------
$ uname -a
Linux xxx 3.16.0-4-586 #1 Debian 3.16.7-ckt11-1 (2015-05-24) i686 GNU/Linux
$ cat test2.kumac
*** compatibility test
SET *
OPT *
OPT NBOX
TRACE ON
OPT LOGY
VECT/CRE DX(14) R 14*0.0
VECT/CRE X(14) R 14*0
VECT/CRE Y(14) R 14*0
VECT/CRE W(14) R 14*0
VECT/CRE DY(14) R 14*0
APPLI COMIS QUIT
subroutine filvec
VECTOR X,Y,DY,W
DO J=1,14
X(J)=1.+0.4*(J-1)
Y(J)=EXP(8.-3.0*X(J))
W(J)=EXP(8.-3.0*X(J))
DY(J)=0.1*Y(J)
PRINT *,J,X(J),Y(J),W(J),DY(J)
ENDDO
W(1)=0.9*W(1)
END
CALL FILVEC
END
QUIT
NULL 0 7 1.e-07 1.E+03
HPLOT/ERROR X Y DX DY 14 21
VECT/FIT X Y DY E S
WAIT
NULL 0 7 1.e-07 1.E+03
HPLOT/ERROR X W DX DY 14 21
VECT/FIT X W DY E S
--------------------------------------------------------------------------------------------------------------------------
The correct behavior of the minuit is observed on a similar linux box:
$ uname -a
Linux TestBench 3.2.0-4-686-pae #1 SMP Debian 3.2.68-1+deb7u1 i686 GNU/Linux
but the paw version in this case was not updated since jessie (Debian 6)
The examples of the paw output are given in the attacments.
Best regards,
Mikhail Ukhanov.
-- System Information:
Debian Release: 8.1
APT prefers stable
APT policy: (500, 'stable')
Architecture: amd64 (x86_64)
Kernel: Linux 3.16.0-4-amd64 (SMP w/4 CPU cores)
Locale: LANG=ru_RU.UTF-8, LC_CTYPE=ru_RU.UTF-8 (charmap=UTF-8)
Shell: /bin/sh linked to /bin/dash
Init: systemd (via /run/systemd/system)
Versions of packages paw-common depends on:
ii cernlib-base 20061220+dfsg3-4.1
ii netbase 5.3
Versions of packages paw-common recommends:
ii libc6-dev [libc-dev] 2.19-18
ii paw [paw-binary] 1:2.14.04.dfsg.2-9+b1
ii paw++ [paw-binary] 1:2.14.04.dfsg.2-9+b1
paw-common suggests no packages.
-- no debconf information
-------------- next part --------------
$ paw
******************************************************
* *
* W E L C O M E to P A W *
* *
* Version 2.14/04 12 January 2004 *
* *
******************************************************
Workstation type (?=HELP) <CR>=1 :
Version 1.29/04 of HIGZ started
*** No default PAWLOGON file "/home/ukhanov/.pawlogon.kumac" found
PAW > exec test2.kumac
1. 8: OPT LOGY
1. 9:
1. 10:
1. 11: VECT/CRE DX(14) R 14*0.0
1. 12: VECT/CRE X(14) R 14*0
1. 13: VECT/CRE Y(14) R 14*0
1. 14: VECT/CRE W(14) R 14*0
1. 15: VECT/CRE DY(14) R 14*0
1. 16:
1. 17: APPLI COMIS QUIT
1. 18: subroutine filvec
1. 19: VECTOR X,Y,DY,W
1. 20: DO J=1,14
1. 21: X(J)=1.+0.4*(J-1)
1. 22: Y(J)=EXP(8.-3.0*X(J))
1. 23: W(J)=EXP(8.-3.0*X(J))
1. 24: DY(J)=0.1*Y(J)
1. 25: PRINT *,J,X(J),Y(J),W(J),DY(J)
1. 26: ENDDO
1. 27: W(1)=0.9*W(1)
1. 28: END
1. 29: CALL FILVEC
1. 30: END
1. 31: QUIT
1 1.000000 148.4132 148.4132 14.84132
2 1.400000 44.70119 44.70119 4.470119
3 1.800000 13.46374 13.46374 1.346374
4 2.200000 4.055198 4.055198 0.4055198
5 2.600000 1.221403 1.221403 0.1221403
6 3.000000 0.3678795 0.3678795 0.3678795E-01
7 3.400000 0.1108031 0.1108031 0.1108031E-01
8 3.800000 0.3337328E-01 0.3337328E-01 0.3337328E-02
9 4.200000 0.1005184E-01 0.1005184E-01 0.1005184E-02
10 4.600000 0.3027551E-02 0.3027551E-02 0.3027551E-03
11 5.000000 0.9118820E-03 0.9118820E-03 0.9118820E-04
12 5.400000 0.2746534E-03 0.2746534E-03 0.2746534E-04
13 5.800000 0.8272394E-04 0.8272394E-04 0.8272395E-05
14 6.200000 0.2491600E-04 0.2491600E-04 0.2491600E-05
1. 32:
1. 33: NULL 0 7 1.e-07 1.E+03
1. 34: HPLOT/ERROR X Y DX DY 14 21
1. 35: VECT/FIT X Y DY E S
MINUIT RELEASE 96.03 INITIALIZED. DIMENSIONS 100/ 50 EPSMAC= 0.89E-15
**********
** 1 **SET EPS 0.1000E-06
**********
FLOATING-POINT NUMBERS ASSUMED ACCURATE TO 0.100E-06
**********
** 2 **SET ERR 1.000
**********
1
**********************************************
* *
* Function minimization by SUBROUTINE HFITV *
* Variable-metric method *
* ID = 0 CHOPT = S *
* *
**********************************************
Convergence when estimated distance to minimum (EDM) .LT. 0.10E+01
FCN= 0.1684534E-09 FROM MIGRAD STATUS=CONVERGED 24 CALLS 25 TOTAL
EDM= 0.42E-03 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 P1 8.0000 0.17001 0.0000 0.14353
2 P2 -3.0000 0.43085E-01 -0.0000 1.5079
CHISQUARE = 0.1404E-10 NPFIT = 14
1. 36: WAIT
Type <CR> to continue or Q to quit
1. 37:
1. 38: NULL 0 7 1.e-07 1.E+03
1. 39: HPLOT/ERROR X W DX DY 14 21
1. 40: VECT/FIT X W DY E S
1
**********************************************
* *
* Function minimization by SUBROUTINE HFITV *
* Variable-metric method *
* ID = 0 CHOPT = S *
* *
**********************************************
Convergence when estimated distance to minimum (EDM) .LT. 0.10E+01
FCN= 0.7456951 FROM MIGRAD STATUS=FAILED 24 CALLS 25 TOTAL
EDM= 0.11E-01 STRATEGY= 1 ERR MATRIX APPROXIMATE
EXT PARAMETER APPROXIMATE STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 P1 7.9654 0.17893 0.0000 -0.75687
2 P2 -2.9925 0.44889E-01 -0.0000 1.7350
CHISQUARE = 0.6214E-01 NPFIT = 14
PAW > q
Exiting from PAW.
-------------- next part --------------
$ paw
******************************************************
* *
* W E L C O M E to P A W *
* *
* Version 2.14/04 12 January 2004 *
* *
******************************************************
Workstation type (?=HELP) <CR>=1 :
Version 1.29/04 of HIGZ started
*** No default PAWLOGON file "/home/ukhanov/.pawlogon.kumac" found
PAW > exec test2.kumac
1. 8: OPT LOGY
1. 9:
1. 10:
1. 11: VECT/CRE DX(14) R 14*0.0
1. 12: VECT/CRE X(14) R 14*0
1. 13: VECT/CRE Y(14) R 14*0
1. 14: VECT/CRE W(14) R 14*0
1. 15: VECT/CRE DY(14) R 14*0
1. 16:
1. 17: APPLI COMIS QUIT
1. 18: subroutine filvec
1. 19: VECTOR X,Y,DY,W
1. 20: DO J=1,14
1. 21: X(J)=1.+0.4*(J-1)
1. 22: Y(J)=EXP(8.-3.0*X(J))
1. 23: W(J)=EXP(8.-3.0*X(J))
1. 24: DY(J)=0.1*Y(J)
1. 25: PRINT *,J,X(J),Y(J),W(J),DY(J)
1. 26: ENDDO
1. 27: W(1)=0.9*W(1)
1. 28: END
1. 29: CALL FILVEC
1. 30: END
1. 31: QUIT
1 1.000000 148.4132 148.4132 14.84132
2 1.400000 44.70119 44.70119 4.470119
3 1.800000 13.46374 13.46374 1.346374
4 2.200000 4.055198 4.055198 0.4055198
5 2.600000 1.221403 1.221403 0.1221403
6 3.000000 0.3678795 0.3678795 0.3678795E-01
7 3.400000 0.1108031 0.1108031 0.1108031E-01
8 3.800000 0.3337328E-01 0.3337328E-01 0.3337328E-02
9 4.200000 0.1005184E-01 0.1005184E-01 0.1005184E-02
10 4.600000 0.3027551E-02 0.3027551E-02 0.3027551E-03
11 5.000000 0.9118820E-03 0.9118820E-03 0.9118820E-04
12 5.400000 0.2746534E-03 0.2746534E-03 0.2746534E-04
13 5.800000 0.8272394E-04 0.8272394E-04 0.8272395E-05
14 6.200000 0.2491600E-04 0.2491600E-04 0.2491600E-05
1. 32:
1. 33: NULL 0 7 1.e-07 1.E+03
1. 34: HPLOT/ERROR X Y DX DY 14 21
1. 35: VECT/FIT X Y DY E S
MINUIT RELEASE 96.03 INITIALIZED. DIMENSIONS 100/ 50 EPSMAC= 0.89E-15
**********
** 1 **SET EPS 0.1000E-06
**********
FLOATING-POINT NUMBERS ASSUMED ACCURATE TO 0.100E-06
**********
** 2 **SET ERR 1.000
**********
1
**********************************************
* *
* Function minimization by SUBROUTINE HFITV *
* Variable-metric method *
* ID = 0 CHOPT = S *
* *
**********************************************
Convergence when estimated distance to minimum (EDM) .LT. 0.10E+01
FCN= 0.2913009E-11 FROM MIGRAD STATUS=CONVERGED 34 CALLS 35 TOTAL
EDM= 0.16E-03 STRATEGY= 1 ERR MATRIX NOT POS-DEF
EXT PARAMETER APPROXIMATE STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 P1 8.0000 1.1162 0.10119E-01 0.20447E-01
2 P2 -3.0000 0.91753 0.37947E-02 0.53462E-02
CHISQUARE = 0.9999E+04 NPFIT = 2
1. 36: WAIT
Type <CR> to continue or Q to quit
1. 37:
1. 38: NULL 0 7 1.e-07 1.E+03
1. 39: HPLOT/ERROR X W DX DY 14 21
1. 40: VECT/FIT X W DY E S
1
**********************************************
* *
* Function minimization by SUBROUTINE HFITV *
* Variable-metric method *
* ID = 0 CHOPT = S *
* *
**********************************************
Convergence when estimated distance to minimum (EDM) .LT. 0.10E+01
FCN= 0.3664465E-04 FROM MIGRAD STATUS=CONVERGED 44 CALLS 45 TOTAL
EDM= 0.42E-03 STRATEGY= 1 ERR MATRIX NOT POS-DEF
EXT PARAMETER APPROXIMATE STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 P1 7.6286 1.1671 0.19299E-01 0.29995E-01
2 P2 -2.7344 0.94331 0.69175E-02 0.67855E-01
CHISQUARE = 0.9999E+04 NPFIT = 2
PAW > q
Exiting from PAW.
More information about the debian-science-maintainers
mailing list