エキスパートの皆様に質問です。
Intel core i5にFedora23を入れて、SADを走らせようとしています。
コンパイルが終わり、サンプルプログラムも終了しているのですが、ちゃんと計算されていないようです。実行中につぎのようなエラーを吐いています。
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO
下記がサンプルプログラムの出力です。途中から主な数値がゼロをしめしています。どこかで数値の受け渡しに失敗しているようです。どこらからみたらいいでしょうか。
*** Welcome to SAD Ver.1.0.10.9.5k64 built at 2016-06-24 18:13:34 +0900 ***
*** Today: 18:18:31 Friday 06/24/2016 ***
OFF LOG ECHO;READ 77 ; 23
!
! ***** DEFINITION OF ELEMENTS IN MAIN LEVEL *****
!
;
DRIFT L1 = (L = 1)
;
BEND B = (L = 2)
;
QUAD QF = (L = 1 K1 = 0.1 )
QD = (L = 1 K1 = -0.1 )
QSF = (L = 1 K1 = 0.1 )
QSD = (L = 1 K1 = -0.1 )
QRF = (L = 1 K1 = 0.1 )
QRD = (L = 1 K1 = -0.1 )
;
SEXT SF =(L = 1 K2 = 0.1)
SD =(L = 1 K2 = -0.1)
;
MARK IP1 =(BETAX = 10 BETAY =10 EMIX = 4.0E-7
EMIY = 4.0E-7 DP = 0.01 )
;
CAVI CA1 = (L = 1 VOLT=1 MV HARM=100)
;
!
! ***** DEFINITION OF UNIT CELL IN MAIN LEVEL *****
!
! You need at least one LINE to start FFS, but
! other lines can be created in FFS.
!
;
LINE CELL = (IP1 QF L1 SF L1 B 3*L1 QD L1 SD L1 B 3*L1)
;
FFS USE=CELL;
*** SADScript Initialization: /home/mkuriki/SAD-master/Packages/init.n ***
*** Run time Environment: /home/mkuriki/SAD-master/Packages/init.local.n **
*
RFSW RADCOD RAD FLUC INTRA POL COD DAPER EMIOU CMPLO FOURI SMEAR
F F F T F F T F F F F T
Print[TimeUsed[]];! CPU Time used so far.
.030489999800920486
!
! Defining parameters:
!
nbends=24; ! number of bends per ring
nxcell=0.25; ! horizontal tune/cell
nycell=0.25; ! vertical tune/cell
Print[TimeUsed[]];! CPU Time used so far.
.030579999089241028
!
! ***** UNIT CELL MATCHING *****
!
CELL; ! peridic condition
B 2*Pi/nbends; ! settin bending angle to the BEND B
FIT; ! set fit point at end of line
NX nxcell; ! set fit condition NX
NY nycell; ! set fit condition NY
FREE Q*; ! set Q* (in this case QF and QD) as the matching
! variable
GO; ! start matching
Iterations Residual Method Reduction Variables
2 4.3472E-03 (NEWTON) 1.000 2
3 1.6466E-07 (NEWTON) 1.000 2
Matched. ( 1.6175E-15) DP = 0.01000 DP0 = 0.00000 ExponentOfResidual = 2.0 O
ffMomentumWeight = 1.000
$$$ f AX ####### # -1.496566 $$$ f BX ####### # 23.647704
$$$ f NX .25 1 .250000
$$$ f AY ####### # .467389 $$$ f BY ####### # 6.097759
$$$ f NY .25 1 .250000
$$$ f LENG ####### # 18.000000
! define dr as the drawing command (needs X-Window)
dr:=FFS["OUT 'a' DRAW BX BY & EX EY {BQ}*; TERM OUT; TDR 'a';"];
dr; ! draw optics of unit cell
SAVE; ! save the matching result (values of QF, QD) to
! keep them after swtch the beam line.
Print[TimeUsed[]];! CPU Time used so far.
.2547670006752014
!
! ***** DISPERSION SUPPRESSOR *****
!
unitcell=ExtractBeamLine[]; ! get the current BeamLine
! define a dispersion suppressor to insert rf
supp=BeamLine[QSF, 8*L1, QSD, 3*L1, B, 3*L1, QRF, 3*L1, CA1, 3*L1, QRD];
USE Join[unitcell, supp]; ! switch the beam line with suppressor
INS; ! now nonperiodic (a transport line)
QRD L 0.5; ! set the thickness of QD
FIT;
AX 0; ! set the end of line to be a symmetry point
AY 0;
EX 0; ! and dispersion-free
EPX 0;
Q*F MIN 0; ! set the lower limit of K1 of QRF, QSF
Q*D MAX 0; ! set the upper limit of K1 of QRD, QSD
FREE QS* QR*;
GO;
Iterations Residual Method Reduction Variables
3 258.4 (NEWTON) 6.2500E-02 4
4 118.5 (NEWTON) 0.2500 4
5 36.05 (NEWTON) 1.000 3
7 22.82 (NEWTON) 0.2044 3
8 4.008 (NEWTON) 0.8176 2
9 2.705 (NEWTON) 1.000 4
13 2.465 (NEWTON) 1.5625E-02 4
14 2.096 (NEWTON) 6.2500E-02 4
15 0.4112 (NEWTON) 0.2500 4
17 0.2073 (NEWTON) 0.2077 4
18 1.5378E-03 (NEWTON) 0.8306 4
19 1.4651E-06 (NEWTON) 1.000 4
Matched. ( 2.7823E-13) DP = 0.01000 DP0 = 0.00000 ExponentOfResidual = 2.0 O
ffMomentumWeight = 1.000
$$$ f AX 0 1 1.9485E-9 $$$ f BX ####### # 1.784152
$$$ f NX ####### # .716798
$$$ f AY 0 1 -5.273E-7 $$$ f BY ####### # 110.11073
$$$ f NY ####### # .792207
$$$ f EX 0 1 -2.011E-9 $$$ f EPX 0 1 -1.89E-10
$$$ f LENG ####### # 44.500000
dr; ! draw suppressor
SAVE;
Print[TimeUsed[]];! CPU Time used so far.
.2854500114917755
!
! ***** TOTAL RING *****
!
ncell=(nbends/2-2)/2; ! number of unitcells per half ring.
! Defining a half ring, removing the
! first QF and IP1 markers in unitcell
hring=BeamLine[IP1,-supp,Rest[ncell*Rest[unitcell]],supp];
USE Join[hring,-hring]; ! switch to full ring
CELL; ! now periodic condition again
CAL NX NY; ! set to display NX and NY after CAL
CAL ! calculate the optics of ring
Matched. ( 0.000 ) DP = 0.01000 DP0 = 0.00000 ExponentOfResidual = 2.0 O
ffMomentumWeight = 1.000
$$$ f AX ####### # -3.25E-15 $$$ f BX ####### # 1.784152
$$$ f NX ####### # 4.354010
$$$ f AY ####### # 3.553E-15 $$$ f BY ####### # 110.11071
$$$ f NY ####### # 4.615261
$$$ f LENG ####### # 284.00000
dr; ! draw ring
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Print[TimeUsed[]];! CPU Time used so far.
.33870500326156616
!
! ***** CHROMATICITY CORRECTION *****
!
! get the current tunes
{nx0, ny0}=Twiss[{"NX","NY"},"***"]/2/Pi;
FIT;
NX nx0 5; ! set the offmomentum tunes for 5 points
NY ny0 5; ! in -DP < dp/p0 < DP
DP=0.01; ! set the bandwidth |dp/p0|max = DP = 0.01
FREE S*; ! set S* (SF and SD) sextupoles as variable
GO;
Iterations Residual Method Reduction Variables
2 3.8142E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
19 3.8110E-04 (NEWTON) 1.000 2
*****qcod---> Overflow & closed orbit not found
***qmdiag---> Sum resonance: (TrA-TrB)^2 = NaN
Residual = 3.8110E-04 DP = 0.01000 DP0 = 0.00000 ExponentOfResidual = 2.0 O
ffMomentumWeight = 1.000
DP -.010000 -.005000 .000000 .005000 .010000
Res. 2.1016E-4 1.4048E-5 8.323E-29 8.4830E-6 1.4841E-4
$$$ f AX ####### # 1.786E-13 6.412E-15 -3.25E-15 3.691E-15 -4.88E-14
$$$ f BX ####### # 1.441761 1.605653 1.784152 1.987500 2.229826
$$$ f NX 4.35401 5 4.362037 4.356146 4.354010 4.355570 4.360806
$$$ f AY ####### # -1.94E-13 1.776E-15 3.553E-15 1.776E-15 3.020E-14
$$$ f BY ####### # 120.28086 114.82778 110.11071 105.30067 99.771916
$$$ f NY 4.61526 5 4.623219 4.617256 4.615261 4.616912 4.621898
$$$ f LENG ####### # 284.00000
FITP 1; ! set the off-momentum points to 1 (on-momentum only).
!
! Define Functions to display tunes as functions of dp/p0
nx[dp_]:=(DP0=dp;FFS["CALC"];Twiss["NX","***"]/2/Pi);
ny[dp_]:=(DP0=dp;FFS["CALC"];Twiss["NY","***"]/2/Pi);
!
! Plot tunes -3% < dp/p0 < 3%.
Plot[{nx[dp],ny[dp]}, {dp, -0.03, 0.03},
PlotLabel->"Tune Chromaticity",
FrameLabel->{{"`fD`np/p`d0`n","F X X"},{"`fn`n`dx`n, `fn`n`dy`n","GX X G
X X"}}];
! draw chromaticity
If[$DisplayFunction === CanvasDrawer, Update[]];
TYPE; ! print out all elements
;
DRIFT L1 =(L =1 )
;
BEND B =(L =2 ANGLE =.2617993877991494 )
;
QUAD QRD =(L =.5 K1 =-.10786211799536638 )
QRF =(L =1 K1 =.1945786529110958 )
QSD =(L =1 K1 =-.21812755206494053 )
QSF =(L =1 K1 =.15519425429347902 )
QD =(L =1 K1 =-.15606240584804387 )
QF =(L =1 K1 =.12696257666020425 )
;
SEXT SF =(L =1 K2 =.12077486557515338 )
SD =(L =1 K2 =-.2011189579437624 )
;
CAVI CA1 =(L =1 VOLT =1e+06 HARM =100 )
;
MARK IP1 =(AX =-8.574938606266326e-16 BX =2.1332114948917886 AY =-4.87
6713425966236e-16 BY =101.95488815812298
EX =-.14915073094352782 EPX =-6.386658675581085e-18 DX =-.00061712
54766648309 DPX =-3.394908253619185e-16
DDP =.008125 DP =.01 EMITX =4e-07 EMITY =4e-07 )
;
PrintBeamLine[ExtractBeamLine[]]; ! print out all beam line
BeamLine[
IP1, -QRD, -L1, -L1, -L1, -CA1,
-L1, -L1, -L1, -QRF, -L1, -L1,
-L1, -B, -L1, -L1, -L1, -QSD,
-L1, -L1, -L1, -L1, -L1, -L1,
-L1, -L1, -QSF, L1, SF, L1,
B, L1, L1, L1, QD, L1,
SD, L1, B, L1, L1, L1,
QF, L1, SF, L1, B, L1,
L1, L1, QD, L1, SD, L1,
B, L1, L1, L1, QF, L1,
SF, L1, B, L1, L1, L1,
QD, L1, SD, L1, B, L1,
L1, L1, QF, L1, SF, L1,
B, L1, L1, L1, QD, L1,
SD, L1, B, L1, L1, L1,
QF, L1, SF, L1, B, L1,
L1, L1, QD, L1, SD, L1,
B, L1, L1, L1, QSF, L1,
L1, L1, L1, L1, L1, L1,
L1, QSD, L1, L1, L1, B,
L1, L1, L1, QRF, L1, L1,
L1, CA1, L1, L1, L1, QRD,
-QRD, -L1, -L1, -L1, -CA1, -L1,
-L1, -L1, -QRF, -L1, -L1, -L1,
-B, -L1, -L1, -L1, -QSD, -L1,
-L1, -L1, -L1, -L1, -L1, -L1,
-L1, -QSF, -L1, -L1, -L1, -B,
-L1, -SD, -L1, -QD, -L1, -L1,
-L1, -B, -L1, -SF, -L1, -QF,
-L1, -L1, -L1, -B, -L1, -SD,
-L1, -QD, -L1, -L1, -L1, -B,
-L1, -SF, -L1, -QF, -L1, -L1,
-L1, -B, -L1, -SD, -L1, -QD,
-L1, -L1, -L1, -B, -L1, -SF,
-L1, -QF, -L1, -L1, -L1, -B,
-L1, -SD, -L1, -QD, -L1, -L1,
-L1, -B, -L1, -SF, -L1, -QF,
-L1, -L1, -L1, -B, -L1, -SD,
-L1, -QD, -L1, -L1, -L1, -B,
-L1, -SF, -L1, QSF, L1, L1,
L1, L1, L1, L1, L1, L1,
QSD, L1, L1, L1, B, L1,
L1, L1, QRF, L1, L1, L1,
CA1, L1, L1, L1, QRD, -IP1
];
Print[TimeUsed[]];! CPU Time used so far.
.45032501220703125
!
! ***** EMITTANCE CALCULATION *****
!
EMIOUT; ! turn on the output of matrices by EMIT.
EMIT; ! This (or Emittance[]) is necessary before tracking
Closed orbit:
x px/p0 y py/p0 z dp/p0
Entrance : .000000 .000000 .000000 .000000 .000000 .000000
Exit : .000000 .000000 .000000 .000000 .000000 .000000
Symplectic part of the transfer matrix:
x px/p0 y py/p0 z dp/p0
x : 1.000000 .000000 .000000 .000000 .000000 .000000
px/p0 : .000000 1.000000 .000000 .000000 .000000 .000000
y : .000000 .000000 1.000000 .000000 .000000 .000000
py/p0 : .000000 .000000 .000000 1.000000 .000000 .000000
z : .000000 .000000 .000000 .000000 1.000000 .000000
dp/p0 : .000000 .000000 .000000 .000000 .000000 1.000000
x px/p0 y py/p0 z dp/p0
x : 1.000000 .000000 .000000 .000000 .000000 .000000
px/p0 : .000000 1.000000 .000000 .000000 .000000 .000000
y : .000000 .000000 1.000000 .000000 .000000 .000000
py/p0 : .000000 .000000 .000000 1.000000 .000000 .000000
z : .000000 .000000 .000000 .000000 1.000000 .000000
dp/p0 : .000000 .000000 .000000 .000000 .000000 1.000000
X Px Y Py Z Pz
X : 1.000000 .000000 .000000 .000000 .000000 .000000
Px : .000000 1.000000 .000000 .000000 .000000 .000000
Y : .000000 .000000 1.000000 .000000 .000000 .000000
Py : .000000 .000000 .000000 1.000000 .000000 .000000
Z : .000000 .000000 .000000 .000000 1.000000 .000000
Pz : .000000 .000000 .000000 .000000 .000000 1.000000
Extended Twiss Parameters:
AX: .000000 BX: 1.000000 ZX: .000000 EX: .000000
PSIX: .000000 ZPX: .000000 EPX: .000000
R1: .000000 R2: .000000 AY: .000000 BY: 1.000000 ZY: .000000 EY: .000000
R3: .000000 R4: .000000 PSIY: .000000 ZPY: .000000 EPY: .000000
AZ: .000000 BZ: 1.000000
PSIZ: .000000
Units: B(X,Y,Z), E(X,Y), R2: m | PSI(X,Y,Z): radian | ZP(X,Y), R3: 1/m
Design momentum P0 = 1.0000000 GeV Revolution freq. f0 = 1055607.1 Hz
Energy loss per turn U0 = .0000000 MV Effective voltage Vc = NaN MV
Equilibrium position dz = NaN mm Momentum compact. alpha = .0000000
Orbit dilation dl = .0000000 mm Effective harmonic # h = NaN
Bucket height dV/P0 = NaN
Eigen values and eigen vectors:
Real: 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000
Imaginary: 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
Imag.tune: 0.0000000 0.0000000 0.0000000
Real tune: 0.0000000 0.0000000 0.0000000
X Px Y Py Z Pz
x : 1.000000 .000000 .000000 .000000 .000000 .000000
px/p0 : .000000 1.000000 .000000 .000000 .000000 .000000
y : .000000 .000000 1.000000 .000000 .000000 .000000
py/p0 : .000000 .000000 .000000 1.000000 .000000 .000000
z : .000000 .000000 .000000 .000000 1.000000 .000000
dp/p0 : .000000 .000000 .000000 .000000 .000000 1.000000
x px/p0 y py/p0 z dp/p0
X : 1.000000 .000000 .000000 .000000 .000000 .000000
Px : .000000 1.000000 .000000 .000000 .000000 .000000
Y : .000000 .000000 1.000000 .000000 .000000 .000000
Py : .000000 .000000 .000000 1.000000 .000000 .000000
Z : .000000 .000000 .000000 .000000 1.000000 .000000
Pz : .000000 .000000 .000000 .000000 .000000 1.000000
x px/p0 y py/p0 z dp/p0
x : 1.000000 .000000 .000000 .000000 .000000 .000000
px/p0 : .000000 1.000000 .000000 .000000 .000000 .000000
y : .000000 .000000 1.000000 .000000 .000000 .000000
py/p0 : .000000 .000000 .000000 1.000000 .000000 .000000
z : .000000 .000000 .000000 .000000 1.000000 .000000
dp/p0 : .000000 .000000 .000000 .000000 .000000 1.000000
Radiation part of the transfer matrix:
x px/p0 y py/p0 z dp/p0
x : NaN NaN NaN NaN NaN NaN
px/p0 : NaN NaN NaN NaN NaN NaN
y : NaN NaN NaN NaN NaN NaN
py/p0 : NaN NaN NaN NaN NaN NaN
z : NaN NaN NaN NaN NaN NaN
dp/p0 : NaN NaN NaN NaN NaN NaN
X Px Y Py Z Pz
X : NaN NaN NaN NaN NaN NaN
Px : NaN NaN NaN NaN NaN NaN
Y : NaN NaN NaN NaN NaN NaN
Py : NaN NaN NaN NaN NaN NaN
Z : NaN NaN NaN NaN NaN NaN
Pz : NaN NaN NaN NaN NaN NaN
Damping per one revolution:
X : NaN Y : NaN Z : NaN
Damping time (sec):
X : NaN Y : NaN Z : NaN
Tune shift due to radiation:
X : NaN Y : NaN Z : NaN
Damping partition number:
X : NaN Y : NaN Z : NaN
Beam matrix by radiation fluctuation:
x px/p0 y py/p0 z dp/p0
x NaN
px/p0 NaN NaN
y NaN NaN NaN
py/p0 NaN NaN NaN NaN
z NaN NaN NaN NaN NaN
dp/p0 NaN NaN NaN NaN NaN NaN
X Px Y Py Z Pz
X NaN
Px NaN NaN
Y NaN NaN NaN
Py NaN NaN NaN NaN
Z NaN NaN NaN NaN NaN
Pz NaN NaN NaN NaN NaN NaN
Equiliblium beam matrix:
X Px Y Py Z Pz
X .0000000
Px NaN .0000000
Y NaN NaN .0000000
Py NaN NaN NaN .0000000
Z NaN NaN NaN NaN .0000000
Pz NaN NaN NaN NaN NaN .0000000
x px/p0 y py/p0 z dp/p0
x NaN
px/p0 NaN NaN
y NaN NaN NaN
py/p0 NaN NaN NaN NaN
z NaN NaN NaN NaN NaN
dp/p0 NaN NaN NaN NaN NaN NaN
Emittance X = NaN m Emittance Y = NaN m
Emittance Z = NaN m Energy spread = NaN
Bunch Length = NaN mm Beam tilt = NaN rad
Beam size xi = NaN mm Beam size eta = NaN mm
Print[TimeUsed[]];! CPU Time used so far.
.45147401094436646
!
! ***** DYNAMIC APERTURE *****
!
! Dynamic aperture for 1000 turns,
DynamicApertureSurvey[{{0,100},{0,100},Range[-30,30,4]},
1000,Output->6];
Turns =1000 Maximum number of particles =224
Range Xmin: 0.000 Xmax: 100.000
(Ymin: 0.000 Ymax: 100.000)
Zmin: -30.000 Zmax: 30.000
Display: 100 turns/character
NZ 0----|----1----|----2----|----3----|----4----|----5
-30.00 0 . . . . .
-26.00 0 . . . . .
-22.00 0 . . . . .
-18.00 0 . . . . .
-14.00 0 . . . . .
-10.00 0 . . . . .
-6.00 0 . . . . .
-2.00 0 . . . . .
2.00 0 . . . . .
6.00 0 . . . . .
10.00 0 . . . . .
14.00 0 . . . . .
18.00 0 . . . . .
22.00 0 . . . . .
26.00 0 . . . . .
30.00 0 . . . . .
NZ 0----|----1----|----2----|----3----|----4----|----5
Score: 0
!
Print[TimeUsed[]];! CPU Time used so far.
.45451200008392334
ABORT ! stop SAD.
> エキスパートの皆様に質問です。
>
> Intel core i5にFedora23を入れて、SADを走らせようとしています。
> コンパイルが終わり、サンプルプログラムも終了しているのですが、ちゃんと計算されていないようです。実行中につぎのようなエラーを吐いています。
>
> Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO
>
このメッセージはGNU Fortran runtimeが出力しているもので実行途中の浮動小数点演算例外発生を知らせているだけで、SADによる出力ではありません
詳しくは、GNU Fortranのコンパイルオプションに関するマニュアルを参照のこと
https://gcc.gnu.org/gcc-4.9/changes.html
https://gcc.gnu.org/onlinedocs/gcc-4.9.0/gfortran/Debugging-Options.html
> エキスパートの皆様に質問です。
>
> Intel core i5にFedora23を入れて、SADを走らせようとしています。
> コンパイルが終わり、サンプルプログラムも終了しているのですが、ちゃんと計算されていないようです。実行中につぎのようなエラーを吐いています。
>
> Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO
>
> 下記がサンプルプログラムの出力です。途中から主な数値がゼロをしめしています。どこかで数値の受け渡しに失敗しているようです。どこらからみたらいいでしょうか。
>
付属のサンプルが正しく演算出来ない状況であれば、調査すべき主なポイントは以下の3つだと思われます
1. コンパイルオプションに依存した不具合の可能性
過度な最適化オプションを適用している場合は、コードの不具合(ill-formed code)やコンパイラの最適化バグの可能性が疑われます
2. 実行環境の不具合
ポインタ・オフセット演算(64bit環境の場合)や浮動小数点演算(libm)の環境依存の不具合の可能性が考えられます
3. ソースコードへのエラッタ混入
コードの改変に伴いエラッタが混入するケースは、割と多いです
実務的には、core dumpを伴わない不具合なので、正常に動作する組み合わせ(OS/コンパイラ/ソースリビジョン)と
不具合が発生する組み合わせの間で二分法を繰り返し原因部位を特定してゆく地道な作業が
一番確実です
情報どうもありがとうございます。
コンパイラーのオプション、バージョン等を変えてみて、
様子が変わるかどうか試してみます。
> > エキスパートの皆様に質問です。
> >
> > Intel core i5にFedora23を入れて、SADを走らせようとしています。
> > コンパイルが終わり、サンプルプログラムも終了しているのですが、ちゃんと計算されていないようです。実行中につぎのようなエラーを吐いています。
> >
> > Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO
> >
> > 下記がサンプルプログラムの出力です。途中から主な数値がゼロをしめしています。どこかで数値の受け渡しに失敗しているようです。どこらからみたらいいでしょうか。
> >
> 付属のサンプルが正しく演算出来ない状況であれば、調査すべき主なポイントは以下の3つだと思われます
>
> 1. コンパイルオプションに依存した不具合の可能性
> 過度な最適化オプションを適用している場合は、コードの不具合(ill-formed code)やコンパイラの最適化バグの可能性が疑われます
>
> 2. 実行環境の不具合
> ポインタ・オフセット演算(64bit環境の場合)や浮動小数点演算(libm)の環境依存の不具合の可能性が考えられます
>
> 3. ソースコードへのエラッタ混入
> コードの改変に伴いエラッタが混入するケースは、割と多いです
>
> 実務的には、core dumpを伴わない不具合なので、正常に動作する組み合わせ(OS/コンパイラ/ソースリビジョン)と
> 不具合が発生する組み合わせの間で二分法を繰り返し原因部位を特定してゆく地道な作業が
> 一番確実です
皆様、
sampleの実行結果をみて、計算に不具合があるように早合点してしまいましたが、htmlにあるsamplの実行結果と同じように見えます。
実際にDeckを用意してみて、動作を確認してみます。
> 情報どうもありがとうございます。
>
> コンパイラーのオプション、バージョン等を変えてみて、
> 様子が変わるかどうか試してみます。
>
> > > エキスパートの皆様に質問です。
> > >
> > > Intel core i5にFedora23を入れて、SADを走らせようとしています。
> > > コンパイルが終わり、サンプルプログラムも終了しているのですが、ちゃんと計算されていないようです。実行中につぎのようなエラーを吐いています。
> > >
> > > Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO
> > >
> > > 下記がサンプルプログラムの出力です。途中から主な数値がゼロをしめしています。どこかで数値の受け渡しに失敗しているようです。どこらからみたらいいでしょうか。
> > >
> > 付属のサンプルが正しく演算出来ない状況であれば、調査すべき主なポイントは以下の3つだと思われます
> >
> > 1. コンパイルオプションに依存した不具合の可能性
> > 過度な最適化オプションを適用している場合は、コードの不具合(ill-formed code)やコンパイラの最適化バグの可能性が疑われます
> >
> > 2. 実行環境の不具合
> > ポインタ・オフセット演算(64bit環境の場合)や浮動小数点演算(libm)の環境依存の不具合の可能性が考えられます
> >
> > 3. ソースコードへのエラッタ混入
> > コードの改変に伴いエラッタが混入するケースは、割と多いです
> >
> > 実務的には、core dumpを伴わない不具合なので、正常に動作する組み合わせ(OS/コンパイラ/ソースリビジョン)と
> > 不具合が発生する組み合わせの間で二分法を繰り返し原因部位を特定してゆく地道な作業が
> > 一番確実です