bugGNU Octave - Bugs: bug #56541, [octave forge] (tisean) BIST...

 
 

bug #56541: [octave forge] (tisean) BIST failures in windows

Submitter:  John Donoghue <lostbard>
Submitted:  Mon 24 Jun 2019 04:06:36 PM UTC
   
 
Category:  Octave Package Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Incorrect Result
Status:  None Assigned to:  None
Originator Name:  lostbard Open/Closed:  * Open
Release:  * 5.1.0 Operating System:  * Microsoft Windows
Fixed Release:  None Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Sat 24 Oct 2020 09:05:56 PM UTC, comment #4: 

See also bug #54039.

 Markus Mützel <mmuetzel>
Group administrator
Mon 26 Aug 2019 04:50:21 PM UTC, comment #3: 

tisean for octave compiles with


gfortran -c -fPIC -g -O2 -std=legacy  -g    -freal-4-real-8 -O2

If I use that to compile TISEA_3.0.1 (instead of default -O0"


make FFLAGS="-O2 -std=legacy -freal-4-real-8" CFLAGS="-O2"

I get:


henon -l10 -A1.4 -B0.3 -X0.68587 -Y0.65876 -x10000

TISEAN 3.0.1 (C) R. Hegger, H. Kantz, T. Schreiber (1998-2007)

henon: Henon map
 writing to stdout
  0.69963836153225634      -0.21883607089751045
  0.24905780703219060       0.69963836153225634
   1.1230498007185195       0.24905780703219060
 -0.69101985474181138        1.1230498007185195
  0.66840312470920415      -0.69101985474181138
  0.16722621160801748       0.66840312470920415
   1.1613704892244840       0.16722621160801748
 -0.83812611505571910        1.1613704892244840
  0.36497360813359569      -0.83812611505571910
  0.56207413699560660       0.36497360813359569


which is still different from octave's results though...

Dmitri.
--

Dmitri A. Sergatskov <dasergatskov>
Mon 26 Aug 2019 02:38:58 PM UTC, comment #2: 

Perhaps related ?:

in linux:


octave>>  henon(10)
ans =

   0.67778   0.18938
   0.41367   0.67778
   0.96376   0.41367
  -0.17626   0.96376
   1.24563  -0.17626
  -1.22512   1.24563
  -0.72760  -1.22512
  -0.10869  -0.72760
   0.76518  -0.10869
   0.14769   0.76518


Using a compiled version of the Tisean_3.0.1 tools:


> ./henon -l10 -A1.4 -B0.3 -X0.68587 -Y0.65876 -x10000

TISEAN 3.0.1 (C) R. Hegger, H. Kantz, T. Schreiber (1998-2007)

./henon: Henon map
 writing to stdout
  0.653030694     -0.156294793
  0.356082797      0.653030694
   1.01839626      0.356082797
 -0.345158488       1.01839626
   1.13873076     -0.345158488
 -0.918938339       1.13873076
  0.159392521     -0.918938339
  0.688750148      0.159392521
  0.383690357      0.688750148
   1.00051939      0.383690357


This appears consistent in both Linux and windows

 John Donoghue <lostbard>
Group Member
Mon 24 Jun 2019 05:45:42 PM UTC, comment #1: 

I got similar errors (though not as gross) if I compile (on Fedora) with -mfpmath=387 flag:

+verbose+
test av_d2
*** test
 av_d2_res = [2.227339 0.02227155;2.077224 0.07391772;1.937226 0.2742541;1.806664 0.482474;1.684901 0.603453;1.571344 0.598633;1.465441 0.5230358;1.366675 0.4940035;1.274566 0.4025534;1.188664 0.4291942;1.108552 0.720235;1.033839 0.9861274;0.9641622 1.036355;0.8991808 0.8166605;0.8385791 0.6758107;0.7820616 0.6766519;0.7293533 0.7019168;0.6801973 0.7465489;0.6343543 0.7729089;0.5916009 0.7865262;0.551729 0.7587416;0.5145443 0.7644434;0.4798658 0.8011066;0.4475244 0.8391553;0.4173628 0.9049456;0.3892339 0.9205563;0.3630009 0.949196;0.3385358 0.9231388;0.3157196 0.9548996;0.2944412 0.9584774;0.2745968 0.9921522;0.2560899 0.9393911;0.2388303 0.9163786;0.2227339 0.8300805;0.2077224 0.8629162;0.1937226 0.865825;0.1806664 0.8647303;0.1684901 0.8767034;0.1571344 0.9034054;0.1465441 0.9768383];
 res = av_d2 (val);
 assert (cell2mat({res.d2}.')(1:length(av_d2_res),:), av_d2_res, -1e-6);
!!!!! test failed
ASSERT errors for:  assert (cell2mat ({res.d2}.') (1:length (av_d2_res), :),av_d2_res,-1e-6)

   Location  |  Observed  |  Expected  |  Reason
   (1,1)        2.2274       2.2273      Rel err 3.7778e-05 exceeds tol 1e-06 by 4e-05
   (2,1)        2.0773       2.0772      Rel err 3.7684e-05 exceeds tol 1e-06 by 4e-05
   (3,1)        1.9373       1.9372      Rel err 3.7713e-05 exceeds tol 1e-06 by 4e-05
   (4,1)        1.8067       1.8067      Rel err 3.7409e-05 exceeds tol 1e-06 by 4e-05
   (5,1)        1.685        1.6849      Rel err 3.7345e-05 exceeds tol 1e-06 by 4e-05
   (6,1)        1.5714       1.5713      Rel err 3.7544e-05 exceeds tol 1e-06 by 4e-05
   (7,1)        1.4655       1.4654      Rel err 3.7321e-05 exceeds tol 1e-06 by 4e-05
   (8,1)        1.3667       1.3667      Rel err 3.7454e-05 exceeds tol 1e-06 by 4e-05
   (9,1)        1.2746       1.2746      Rel err 3.7202e-05 exceeds tol 1e-06 by 4e-05
   (10,1)       1.1887       1.1887      Rel err 3.7637e-05 exceeds tol 1e-06 by 4e-05
   (11,1)       1.1086       1.1086      Rel err 3.7659e-05 exceeds tol 1e-06 by 4e-05
   (12,1)       1.0339       1.0338      Rel err 3.7955e-05 exceeds tol 1e-06 by 4e-05
   (13,1)       0.9642      0.96416      Rel err 3.7448e-05 exceeds tol 1e-06 by 4e-05
   (14,1)      0.89921      0.89918      Rel err 3.7553e-05 exceeds tol 1e-06 by 4e-05
   (15,1)      0.83861      0.83858      Rel err 3.7462e-05 exceeds tol 1e-06 by 4e-05
   (16,1)      0.78209      0.78206      Rel err 3.7559e-05 exceeds tol 1e-06 by 4e-05
   (17,1)      0.72938      0.72935      Rel err 3.7511e-05 exceeds tol 1e-06 by 4e-05
   (18,1)      0.68022       0.6802      Rel err 3.7554e-05 exceeds tol 1e-06 by 4e-05
   (19,1)      0.63438      0.63435      Rel err 3.752e-05 exceeds tol 1e-06 by 4e-05
   (20,1)      0.59162       0.5916      Rel err 3.7593e-05 exceeds tol 1e-06 by 4e-05
   (21,1)      0.55175      0.55173      Rel err 3.7555e-05 exceeds tol 1e-06 by 4e-05
   (22,1)      0.51456      0.51454      Rel err 3.7579e-05 exceeds tol 1e-06 by 4e-05
   (23,1)      0.47988      0.47987      Rel err 3.7446e-05 exceeds tol 1e-06 by 4e-05
   (24,1)      0.44754      0.44752      Rel err 3.756e-05 exceeds tol 1e-06 by 4e-05
   (25,1)      0.41738      0.41736      Rel err 3.7448e-05 exceeds tol 1e-06 by 4e-05
   (26,1)      0.38925      0.38923      Rel err 3.7562e-05 exceeds tol 1e-06 by 4e-05
   (27,1)      0.36301       0.363       Rel err 3.739e-05 exceeds tol 1e-06 by 4e-05
   (28,1)      0.33855      0.33854      Rel err 3.7551e-05 exceeds tol 1e-06 by 4e-05
   (29,1)      0.31573      0.31572      Rel err 3.7614e-05 exceeds tol 1e-06 by 4e-05
   (30,1)      0.29445      0.29444      Rel err 3.7464e-05 exceeds tol 1e-06 by 4e-05
   (31,1)      0.27461       0.2746      Rel err 3.7641e-05 exceeds tol 1e-06 by 4e-05
   (32,1)       0.2561      0.25609      Rel err 3.762e-05 exceeds tol 1e-06 by 4e-05
   (33,1)      0.23884      0.23883      Rel err 3.7611e-05 exceeds tol 1e-06 by 4e-05
   (34,1)      0.22274      0.22273      Rel err 3.7778e-05 exceeds tol 1e-06 by 4e-05
   (35,1)      0.20773      0.20772      Rel err 3.7684e-05 exceeds tol 1e-06 by 4e-05
   (36,1)      0.19373      0.19372      Rel err 3.7713e-05 exceeds tol 1e-06 by 4e-05
   (37,1)      0.18067      0.18067      Rel err 3.7409e-05 exceeds tol 1e-06 by 4e-05
   (38,1)       0.1685      0.16849      Rel err 3.7345e-05 exceeds tol 1e-06 by 4e-05
   (39,1)      0.15714      0.15713      Rel err 3.7544e-05 exceeds tol 1e-06 by 4e-05
   (40,1)      0.14655      0.14654      Rel err 3.7321e-05 exceeds tol 1e-06 by 4e-05
   (1,2)       0.20394      0.022272     Rel err 8.1569 exceeds tol 1e-06 by 8
   (2,2)       0.31823      0.073918     Rel err 3.3052 exceeds tol 1e-06 by 3
   (3,2)       0.38845      0.27425      Rel err 0.41637 exceeds tol 1e-06 by 0.4
   (4,2)       0.52847      0.48247      Rel err 0.095334 exceeds tol 1e-06 by 0.1
   (5,2)       0.61414      0.60345      Rel err 0.017707 exceeds tol 1e-06 by 0.02
   (6,2)       0.69605      0.59863      Rel err 0.16273 exceeds tol 1e-06 by 0.2
   (7,2)       0.66999      0.52304      Rel err 0.28097 exceeds tol 1e-06 by 0.3
   (8,2)       0.60489       0.494       Rel err 0.22447 exceeds tol 1e-06 by 0.2
   (9,2)       0.63678      0.40255      Rel err 0.58185 exceeds tol 1e-06 by 0.6
   (10,2)      0.56222      0.42919      Rel err 0.30995 exceeds tol 1e-06 by 0.3
   (11,2)      0.62336      0.72023      Rel err 0.13451 exceeds tol 1e-06 by 0.1
   (12,2)      0.59969      0.98613      Rel err 0.39188 exceeds tol 1e-06 by 0.4
   (13,2)      0.54431       1.0364      Rel err 0.47478 exceeds tol 1e-06 by 0.5
   (14,2)      0.55192      0.81666      Rel err 0.32417 exceeds tol 1e-06 by 0.3
   (15,2)      0.57677      0.67581      Rel err 0.14656 exceeds tol 1e-06 by 0.1
   (16,2)      0.67957      0.67665      Rel err 0.0043125 exceeds tol 1e-06 by 0.004
   (17,2)      0.78615      0.70192      Rel err 0.12 exceeds tol 1e-06 by 0.1
   (18,2)      0.77493      0.74655      Rel err 0.038013 exceeds tol 1e-06 by 0.04
   (19,2)      0.83007      0.77291      Rel err 0.073956 exceeds tol 1e-06 by 0.07
   (20,2)      0.74697      0.78653      Rel err 0.05029 exceeds tol 1e-06 by 0.05
   (21,2)      0.78251      0.75874      Rel err 0.031327 exceeds tol 1e-06 by 0.03
   (22,2)      0.80482      0.76444      Rel err 0.052815 exceeds tol 1e-06 by 0.05
   (23,2)      0.82803      0.80111      Rel err 0.033608 exceeds tol 1e-06 by 0.03
   (24,2)      0.84401      0.83916      Rel err 0.0057849 exceeds tol 1e-06 by 0.006
   (25,2)      0.80211      0.90495      Rel err 0.11364 exceeds tol 1e-06 by 0.1
   (26,2)      0.86179      0.92056      Rel err 0.063841 exceeds tol 1e-06 by 0.06
   (27,2)      0.86585       0.9492      Rel err 0.08781 exceeds tol 1e-06 by 0.09
   (28,2)      0.90356      0.92314      Rel err 0.021212 exceeds tol 1e-06 by 0.02
   (29,2)      0.85107       0.9549      Rel err 0.10873 exceeds tol 1e-06 by 0.1
   (30,2)      0.87128      0.95848      Rel err 0.090974 exceeds tol 1e-06 by 0.09
   (31,2)      0.87367      0.99215      Rel err 0.11942 exceeds tol 1e-06 by 0.1
   (32,2)       0.9279      0.93939      Rel err 0.012237 exceeds tol 1e-06 by 0.01
   (33,2)      0.97091      0.91638      Rel err 0.059508 exceeds tol 1e-06 by 0.06
   (34,2)      0.99549      0.83008      Rel err 0.19927 exceeds tol 1e-06 by 0.2
   (35,2)      0.92374      0.86292      Rel err 0.070485 exceeds tol 1e-06 by 0.07
   (36,2)      0.88883      0.86582      Rel err 0.026567 exceeds tol 1e-06 by 0.03
   (37,2)      0.83403      0.86473      Rel err 0.0355 exceeds tol 1e-06 by 0.04
   (38,2)      0.89994       0.8767      Rel err 0.026503 exceeds tol 1e-06 by 0.03
   (39,2)      0.93865      0.90341      Rel err 0.039013 exceeds tol 1e-06 by 0.04
   (40,2)       1.0153      0.97684      Rel err 0.039328 exceeds tol 1e-06 by 0.04
shared variables   scalar structure containing the fields:

    val =

      10x1 struct array containing the fields:

        dim
        c2
        d2
        h2

    c2d_res = [](0x0)
    c1_res = [](0x0)

-verbose-

+verbose+

test c2g.m
*** test
 c2g_res = [0.00255674845539033 0.00262371218 1.56393838;0.00274151703342795 0.00291572767 1.46379328;0.00293963821604848 0.00321952766 1.37964845;0.00315207871608436 0.00353588979 1.30904245;0.00337987160310149 0.00386583689 1.2498455;0.00362412468530238 0.00421059737 1.20023394;0.00388602935709059 0.00457157521 1.1586616;0.00416686292737722 0.00495032407 1.12383032;0.00446798978373408 0.00534852408 1.09465814;0.00479088025167584 0.00576799177 1.07024431;0.00513710314407945 0.00621065963 1.04983807;0.00550834601745009 0.00667859521 1.03281331;0.00590642029419541 0.00717399921 1.01864672;0.00633325986564159 0.00769920321 1.00690293;0.00679094903171062 0.00825670082 0.997214317;0.00728171039372683 0.00884913001 0.989262819;0.00780794164165854 0.00947930664 0.982763529;0.00837219785898924 0.0101501858 0.977451086;0.00897723622620106 0.0108648958 0.973080397;0.00962599366903305 0.011626686 0.969439983];
 val = d2 (henon (1000), 'd', 1, 't', 50, 'm', 5);
 res = c2g (val);
 assert (cell2mat ({res.g}.')(1:rows(c2g_res),:), c2g_res, -1e-5);
!!!!! test failed
ASSERT errors for:  assert (cell2mat ({res.g}.') (1:rows (c2g_res), :),c2g_res,-1e-5)

  Location  |  Observed  |  Expected  |  Reason
   (1,1)      0.0025568    0.0025567     Rel err 3.7466e-05 exceeds tol 1e-05 by 3e-05
   (2,1)      0.0027416    0.0027415     Rel err 3.7659e-05 exceeds tol 1e-05 by 3e-05
   (3,1)      0.0029397    0.0029396     Rel err 3.7876e-05 exceeds tol 1e-05 by 3e-05
   (4,1)      0.0031522    0.0031521     Rel err 3.7555e-05 exceeds tol 1e-05 by 3e-05
   (5,1)       0.00338     0.0033799     Rel err 3.7284e-05 exceeds tol 1e-05 by 3e-05
   (6,1)      0.0036243    0.0036241     Rel err 3.7519e-05 exceeds tol 1e-05 by 3e-05
   (7,1)      0.0038862     0.003886     Rel err 3.7712e-05 exceeds tol 1e-05 by 3e-05
   (8,1)       0.004167    0.0041669     Rel err 3.7466e-05 exceeds tol 1e-05 by 3e-05
   (9,1)      0.0044682     0.004468     Rel err 3.7634e-05 exceeds tol 1e-05 by 3e-05
   (10,1)     0.0047911    0.0047909     Rel err 3.7372e-05 exceeds tol 1e-05 by 3e-05
   (11,1)     0.0051373    0.0051371     Rel err 3.7511e-05 exceeds tol 1e-05 by 3e-05
   (12,1)     0.0055086    0.0055083     Rel err 3.7747e-05 exceeds tol 1e-05 by 3e-05
   (13,1)     0.0059066    0.0059064     Rel err 3.7508e-05 exceeds tol 1e-05 by 3e-05
   (14,1)     0.0063335    0.0063333     Rel err 3.7668e-05 exceeds tol 1e-05 by 3e-05
   (15,1)     0.0067912    0.0067909     Rel err 3.7363e-05 exceeds tol 1e-05 by 3e-05
   (16,1)      0.007282    0.0072817     Rel err 3.7591e-05 exceeds tol 1e-05 by 3e-05
   (17,1)     0.0078082    0.0078079     Rel err 3.7301e-05 exceeds tol 1e-05 by 3e-05
   (18,1)     0.0083725    0.0083722     Rel err 3.7544e-05 exceeds tol 1e-05 by 3e-05
   (19,1)     0.0089776    0.0089772     Rel err 3.7223e-05 exceeds tol 1e-05 by 3e-05
   (20,1)     0.0096264     0.009626     Rel err 3.748e-05 exceeds tol 1e-05 by 3e-05
   (1,2)      0.0027333    0.0026237     Rel err 0.041778 exceeds tol 1e-05 by 0.04
   (2,2)      0.0030363    0.0029157     Rel err 0.041358 exceeds tol 1e-05 by 0.04
   (3,2)      0.0033508    0.0032195     Rel err 0.040786 exceeds tol 1e-05 by 0.04
   (4,2)      0.0036775    0.0035359     Rel err 0.04005 exceeds tol 1e-05 by 0.04
   (5,2)      0.0040172    0.0038658     Rel err 0.039149 exceeds tol 1e-05 by 0.04
   (6,2)       0.004371    0.0042106     Rel err 0.038098 exceeds tol 1e-05 by 0.04
   (7,2)      0.0047404    0.0045716     Rel err 0.036919 exceeds tol 1e-05 by 0.04
   (8,2)      0.0051268    0.0049503     Rel err 0.035645 exceeds tol 1e-05 by 0.04
   (9,2)      0.0055321    0.0053485     Rel err 0.034314 exceeds tol 1e-05 by 0.03
   (10,2)     0.0059581     0.005768     Rel err 0.032961 exceeds tol 1e-05 by 0.03
   (11,2)      0.006407    0.0062107     Rel err 0.03162 exceeds tol 1e-05 by 0.03
   (12,2)     0.0068811    0.0066786     Rel err 0.030318 exceeds tol 1e-05 by 0.03
   (13,2)     0.0073826     0.007174     Rel err 0.029078 exceeds tol 1e-05 by 0.03
   (14,2)     0.0079141    0.0076992     Rel err 0.027918 exceeds tol 1e-05 by 0.03
   (15,2)     0.0084784    0.0082567     Rel err 0.026852 exceeds tol 1e-05 by 0.03
   (16,2)     0.0090782    0.0088491     Rel err 0.025892 exceeds tol 1e-05 by 0.03
   (17,2)     0.0097167    0.0094793     Rel err 0.025044 exceeds tol 1e-05 by 0.03
   (18,2)      0.010397     0.01015      Rel err 0.024318 exceeds tol 1e-05 by 0.02
   (19,2)      0.011123     0.010865     Rel err 0.023714 exceeds tol 1e-05 by 0.02
   (20,2)      0.011897     0.011627     Rel err 0.023234 exceeds tol 1e-05 by 0.02
   (1,3)        1.5591       1.5639      Rel err 0.0030755 exceeds tol 1e-05 by 0.003
   (2,3)        1.457        1.4638      Rel err 0.0046398 exceeds tol 1e-05 by 0.005
   (3,3)        1.3707       1.3796      Rel err 0.0065168 exceeds tol 1e-05 by 0.007
   (4,3)        1.2978       1.309       Rel err 0.0086211 exceeds tol 1e-05 by 0.009
   (5,3)        1.2363       1.2498      Rel err 0.010805 exceeds tol 1e-05 by 0.01
   (6,3)        1.1848       1.2002      Rel err 0.012889 exceeds tol 1e-05 by 0.01
   (7,3)        1.1416       1.1587      Rel err 0.014701 exceeds tol 1e-05 by 0.01
   (8,3)        1.1057       1.1238      Rel err 0.016114 exceeds tol 1e-05 by 0.02
   (9,3)        1.076        1.0947      Rel err 0.01706 exceeds tol 1e-05 by 0.02
   (10,3)       1.0515       1.0702      Rel err 0.017527 exceeds tol 1e-05 by 0.02
   (11,3)       1.0314       1.0498      Rel err 0.017546 exceeds tol 1e-05 by 0.02
   (12,3)       1.0151       1.0328      Rel err 0.017165 exceeds tol 1e-05 by 0.02
   (13,3)       1.0019       1.0186      Rel err 0.016443 exceeds tol 1e-05 by 0.02
   (14,3)      0.99136       1.0069      Rel err 0.015438 exceeds tol 1e-05 by 0.02
   (15,3)      0.98305      0.99721      Rel err 0.014201 exceeds tol 1e-05 by 0.01
   (16,3)      0.97662      0.98926      Rel err 0.012778 exceeds tol 1e-05 by 0.01
   (17,3)      0.97175      0.98276      Rel err 0.011208 exceeds tol 1e-05 by 0.01
   (18,3)      0.96814      0.97745      Rel err 0.0095291 exceeds tol 1e-05 by 0.01
   (19,3)       0.9655      0.97308      Rel err 0.0077881 exceeds tol 1e-05 by 0.008
   (20,3)      0.96357      0.96944      Rel err 0.0060551 exceeds tol 1e-05 by 0.006

-verbose-

Dmitri.
--

Dmitri A. Sergatskov <dasergatskov>
Mon 24 Jun 2019 04:06:36 PM UTC, original submission:  


When running the tests for tisean 0.2.3 in windows 64, using the mxe installer I get 10 failures for the tests:


>>>>> processing C:\Octave\OC328E~1.2\mingw64\share\octave\packages\tisean-0.2.3\av_d2.m
***** test
 av_d2_res = [2.227339 0.02227155;2.077224 0.07391772;1.937226 0.2742541;1.806664 0.482474;1.684901 0.603453;1.571344 0.598633;1.465441 0.5230358;1.366675 0.4940035;1.274566 0.4025534;1.188664 0.4291942;1.108552 0.720235;1.033839 0.9861274;0.9641622 1.036355;0.8991808 0.8166605;0.8385791 0.6758107;0.7820616 0.6766519;0.7293533 0.7019168;0.6801973 0.7465489;0.6343543 0.7729089;0.5916009 0.7865262;0.551729 0.7587416;0.5145443 0.7644434;0.4798658 0.8011066;0.4475244 0.8391553;0.4173628 0.9049456;0.3892339 0.9205563;0.3630009 0.949196;0.3385358 0.9231388;0.3157196 0.9548996;0.2944412 0.9584774;0.2745968 0.9921522;0.2560899 0.9393911;0.2388303 0.9163786;0.2227339 0.8300805;0.2077224 0.8629162;0.1937226 0.865825;0.1806664 0.8647303;0.1684901 0.8767034;0.1571344 0.9034054;0.1465441 0.9768383];
 res = av_d2 (val);
 assert (cell2mat({res.d2}.')(1:length(av_d2_res),:), av_d2_res, -1e-6);
!!!!! test failed
ASSERT errors for:  assert (cell2mat ({res.d2}.') (1:length (av_d2_res), :),av_d2_res,-1e-6)

  Location  |  Observed  |  Expected  |  Reason
   (1,2)       0.11731      0.022272     Rel err 4.2672 exceeds tol 1e-06 by 4
   (2,2)       0.12413      0.073918     Rel err 0.67927 exceeds tol 1e-06 by 0.7
   (3,2)       0.20589      0.27425      Rel err 0.24929 exceeds tol 1e-06 by 0.2
   (4,2)       0.34203      0.48247      Rel err 0.29109 exceeds tol 1e-06 by 0.3
   (5,2)        0.5021      0.60345      Rel err 0.16795 exceeds tol 1e-06 by 0.2
   (6,2)       0.57273      0.59863      Rel err 0.043275 exceeds tol 1e-06 by 0.04
   (7,2)       0.52025      0.52304      Rel err 0.0053256 exceeds tol 1e-06 by 0.005
   (8,2)       0.56711       0.494       Rel err 0.14799 exceeds tol 1e-06 by 0.1
   (9,2)       0.53054      0.40255      Rel err 0.31793 exceeds tol 1e-06 by 0.3
   (10,2)      0.62249      0.42919      Rel err 0.45036 exceeds tol 1e-06 by 0.5
   (11,2)      0.57009      0.72023      Rel err 0.20847 exceeds tol 1e-06 by 0.2
   (12,2)       0.6133      0.98613      Rel err 0.37807 exceeds tol 1e-06 by 0.4
   (13,2)      0.71159       1.0364      Rel err 0.31337 exceeds tol 1e-06 by 0.3
   (14,2)      0.72074      0.81666      Rel err 0.11745 exceeds tol 1e-06 by 0.1
   (15,2)      0.75801      0.67581      Rel err 0.12163 exceeds tol 1e-06 by 0.1
   (16,2)      0.68048      0.67665      Rel err 0.0056627 exceeds tol 1e-06 by 0.006
   (17,2)      0.67993      0.70192      Rel err 0.031319 exceeds tol 1e-06 by 0.03
   (18,2)      0.70543      0.74655      Rel err 0.055083 exceeds tol 1e-06 by 0.06
   (19,2)      0.85072      0.77291      Rel err 0.10068 exceeds tol 1e-06 by 0.1
   (20,2)      0.96862      0.78653      Rel err 0.23152 exceeds tol 1e-06 by 0.2
   (21,2)      0.92209      0.75874      Rel err 0.21529 exceeds tol 1e-06 by 0.2
   (22,2)      0.83271      0.76444      Rel err 0.089303 exceeds tol 1e-06 by 0.09
   (23,2)      0.82439      0.80111      Rel err 0.029066 exceeds tol 1e-06 by 0.03
   (24,2)      0.85831      0.83916      Rel err 0.022832 exceeds tol 1e-06 by 0.02
   (25,2)       0.8987      0.90495      Rel err 0.0068987 exceeds tol 1e-06 by 0.007
   (26,2)      0.94809      0.92056      Rel err 0.029905 exceeds tol 1e-06 by 0.03
   (27,2)       1.0606       0.9492      Rel err 0.11735 exceeds tol 1e-06 by 0.1
   (28,2)      0.93506      0.92314      Rel err 0.012912 exceeds tol 1e-06 by 0.01
   (29,2)      0.84373       0.9549      Rel err 0.11642 exceeds tol 1e-06 by 0.1
   (30,2)      0.79794      0.95848      Rel err 0.16749 exceeds tol 1e-06 by 0.2
   (31,2)       0.9104      0.99215      Rel err 0.082403 exceeds tol 1e-06 by 0.08
   (32,2)      0.96272      0.93939      Rel err 0.024829 exceeds tol 1e-06 by 0.02
   (33,2)       1.0368      0.91638      Rel err 0.1314 exceeds tol 1e-06 by 0.1
   (34,2)      0.99423      0.83008      Rel err 0.19776 exceeds tol 1e-06 by 0.2
   (35,2)       1.0314      0.86292      Rel err 0.19528 exceeds tol 1e-06 by 0.2
   (36,2)      0.98104      0.86582      Rel err 0.13307 exceeds tol 1e-06 by 0.1
   (37,2)      0.99293      0.86473      Rel err 0.14826 exceeds tol 1e-06 by 0.1
   (38,2)      0.94942       0.8767      Rel err 0.082946 exceeds tol 1e-06 by 0.08
   (39,2)      0.86051      0.90341      Rel err 0.04748 exceeds tol 1e-06 by 0.05
   (40,2)      0.85778      0.97684      Rel err 0.12188 exceeds tol 1e-06 by 0.1
shared variables   scalar structure containing the fields:

    val =

      10x1 struct array containing the fields:

        dim
        c2
        d2
        h2

    c2d_res = [](0x0)
    c1_res = [](0x0)
***** test
 av_d2_res_h2 = [1.965624 0.0536151;1.833148 0.07254237;1.7096 0.09431121;1.594379 0.1214435;1.486923 0.1495628;1.386709 0.1786416;1.29325 0.2090021;1.206089 0.2400019;1.124803 0.2694222;1.048995 0.2958831;0.978296 0.3231978;0.912362 0.3510871;0.8508719 0.3752123;0.7935261 0.3983067;0.740045 0.4201098;0.6901685 0.4411169;0.6436534 0.462419;0.6002733 0.4845437;0.5598169 0.5095501;0.5220871 0.5324418;0.4869002 0.55473;0.4540848 0.5774259;0.423481 0.6002299;0.3949398 0.6220529;0.3683222 0.6453185;0.3434985 0.6653961;0.3203478 0.6852458;0.2987574 0.706009;0.2786222 0.7280725;0.259844 0.7495944;0.2423314 0.7692112;0.225999 0.7882578;0.2107675 0.8074904;0.1965624 0.8228356;0.1833148 0.8381825];
 res = av_d2 (val, 'mindim', 2, 'maxdim', 7, 'a',4,'smoothall');
 assert (cell2mat({res.h2}.')(1:length(av_d2_res_h2),:), av_d2_res_h2,-1e-6);
 c2d_res = c2d (c1 (henon (1000)), 2);
!!!!! test failed
ASSERT errors for:  assert (cell2mat ({res.h2}.') (1:length (av_d2_res_h2), :),av_d2_res_h2,-1e-6)

  Location  |  Observed  |  Expected  |  Reason
   (1,2)       0.079127     0.053615     Rel err 0.47584 exceeds tol 1e-06 by 0.5
   (2,2)       0.10331      0.072542     Rel err 0.42411 exceeds tol 1e-06 by 0.4
   (3,2)       0.13024      0.094311     Rel err 0.381 exceeds tol 1e-06 by 0.4
   (4,2)         0.16       0.12144      Rel err 0.31749 exceeds tol 1e-06 by 0.3
   (5,2)       0.19073      0.14956      Rel err 0.27527 exceeds tol 1e-06 by 0.3
   (6,2)       0.22051      0.17864      Rel err 0.23435 exceeds tol 1e-06 by 0.2
   (7,2)        0.2519       0.209       Rel err 0.20525 exceeds tol 1e-06 by 0.2
   (8,2)       0.28101        0.24       Rel err 0.17085 exceeds tol 1e-06 by 0.2
   (9,2)       0.30876      0.26942      Rel err 0.146 exceeds tol 1e-06 by 0.1
   (10,2)      0.33706      0.29588      Rel err 0.13916 exceeds tol 1e-06 by 0.1
   (11,2)      0.36356       0.3232      Rel err 0.1249 exceeds tol 1e-06 by 0.1
   (12,2)       0.388       0.35109      Rel err 0.10515 exceeds tol 1e-06 by 0.1
   (13,2)       0.4113      0.37521      Rel err 0.096167 exceeds tol 1e-06 by 0.1
   (14,2)      0.43338      0.39831      Rel err 0.088053 exceeds tol 1e-06 by 0.09
   (15,2)      0.45463      0.42011      Rel err 0.08216 exceeds tol 1e-06 by 0.08
   (16,2)      0.47525      0.44112      Rel err 0.07738 exceeds tol 1e-06 by 0.08
   (17,2)      0.49657      0.46242      Rel err 0.073857 exceeds tol 1e-06 by 0.07
   (18,2)      0.51673      0.48454      Rel err 0.066428 exceeds tol 1e-06 by 0.07
   (19,2)      0.53564      0.50955      Rel err 0.051195 exceeds tol 1e-06 by 0.05
   (20,2)      0.55489      0.53244      Rel err 0.042153 exceeds tol 1e-06 by 0.04
   (21,2)      0.57633      0.55473      Rel err 0.038946 exceeds tol 1e-06 by 0.04
   (22,2)      0.59816      0.57743      Rel err 0.0359 exceeds tol 1e-06 by 0.04
   (23,2)      0.61961      0.60023      Rel err 0.032283 exceeds tol 1e-06 by 0.03
   (24,2)      0.64256      0.62205      Rel err 0.032966 exceeds tol 1e-06 by 0.03
   (25,2)      0.66539      0.64532      Rel err 0.031109 exceeds tol 1e-06 by 0.03
   (26,2)       0.6876       0.6654      Rel err 0.033371 exceeds tol 1e-06 by 0.03
   (27,2)      0.70788      0.68525      Rel err 0.03303 exceeds tol 1e-06 by 0.03
   (28,2)      0.72732      0.70601      Rel err 0.030182 exceeds tol 1e-06 by 0.03
   (29,2)      0.74665      0.72807      Rel err 0.025514 exceeds tol 1e-06 by 0.03
   (30,2)      0.76392      0.74959      Rel err 0.019107 exceeds tol 1e-06 by 0.02
   (31,2)      0.77778      0.76921      Rel err 0.011145 exceeds tol 1e-06 by 0.01
   (32,2)      0.79059      0.78826      Rel err 0.0029577 exceeds tol 1e-06 by 0.003
   (33,2)      0.80399      0.80749      Rel err 0.0043406 exceeds tol 1e-06 by 0.004
   (34,2)      0.81765      0.82284      Rel err 0.0062993 exceeds tol 1e-06 by 0.006
   (35,2)      0.83035      0.83818      Rel err 0.0093484 exceeds tol 1e-06 by 0.009
shared variables   scalar structure containing the fields:

    val =

      10x1 struct array containing the fields:

        dim
        c2
        d2
        h2

    c2d_res = [](0x0)
    c1_res = [](0x0)
***** test assert (!isequal(av_d2 (c2d_res),c2d_res(1)));
 c1_res = c1 (henon (1000));
!!!!! test failed
d2_c2d_c1_out must be the output of d2 or c2d
shared variables   scalar structure containing the fields:

    val =

      10x1 struct array containing the fields:

        dim
        c2
        d2
        h2

    c2d_res = [](0x0)
    c1_res = [](0x0)
***** test assert (!isequal(av_d2 (c1_res),c1_res(1)));
!!!!! test failed
d2_c2d_c1_out must be the output of d2 or c2d
shared variables   scalar structure containing the fields:

    val =

      10x1 struct array containing the fields:

        dim
        c2
        d2
        h2

    c2d_res = [](0x0)
    c1_res = [](0x0)
***** error <one field> warning("error", "Octave:tisean"); ...
                av_d2 (c2d_res,'smoothall');
!!!!! error failed.
Expected <one field>, but got <d2_c2d_c1_out must be the output of d2 or c2d>

***** error <one field> warning("error", "Octave:tisean"); ...
                av_d2 (c1_res,'smoothall');
!!!!! error failed.
Expected <one field>, but got <d2_c2d_c1_out must be the output of d2 or c2d>



>>>>> processing C:\Octave\OC328E~1.2\mingw64\share\octave\packages\tisean-0.2.3\c2g.m
***** test
 c2g_res = [0.00255674845539033 0.00262371218 1.56393838;0.00274151703342795 0.00291572767 1.46379328;0.00293963821604848 0.00321952766 1.37964845;0.00315207871608436 0.00353588979 1.30904245;0.00337987160310149 0.00386583689 1.2498455;0.00362412468530238 0.00421059737 1.20023394;0.00388602935709059 0.00457157521 1.1586616;0.00416686292737722 0.00495032407 1.12383032;0.00446798978373408 0.00534852408 1.09465814;0.00479088025167584 0.00576799177 1.07024431;0.00513710314407945 0.00621065963 1.04983807;0.00550834601745009 0.00667859521 1.03281331;0.00590642029419541 0.00717399921 1.01864672;0.00633325986564159 0.00769920321 1.00690293;0.00679094903171062 0.00825670082 0.997214317;0.00728171039372683 0.00884913001 0.989262819;0.00780794164165854 0.00947930664 0.982763529;0.00837219785898924 0.0101501858 0.977451086;0.00897723622620106 0.0108648958 0.973080397;0.00962599366903305 0.011626686 0.969439983];
 val = d2 (henon (1000), 'd', 1, 't', 50, 'm', 5);
 res = c2g (val);
 assert (cell2mat ({res.g}.')(1:rows(c2g_res),:), c2g_res, -1e-5);
!!!!! test failed
ASSERT errors for:  assert (cell2mat ({res.g}.') (1:rows (c2g_res), :),c2g_res,-1e-5)

  Location  |  Observed  |  Expected  |  Reason
   (19,3)      0.97309      0.97308      Rel err 1.0245e-05 exceeds tol 1e-05 by 2e-07
   (20,3)      0.96948      0.96944      Rel err 3.7241e-05 exceeds tol 1e-05 by 3e-05


>>>>> processing C:\Octave\OC328E~1.2\mingw64\share\octave\packages\tisean-0.2.3\ikeda.m
***** test
 ikd = [0.28571947432035805 0.64340009417266342;0.50841625820931591 -0.39972597757561129;0.87547944209544659 0.56858743316662752;0.70057389172681206 -0.89053147224601548;-1.9710645882970557E-002 -1.0429004993745194E-002;0.99222205960974108 -1.8501225510088638E-002;0.21620211688366087 -0.42823693940249830;1.3295528700280905 0.27892721893609196;1.0842998438472966 -1.2197367054961652;0.28153619668028618 -1.2810993040439602];
 res = ikeda (10);
 assert (res, ikd, 1e-16);
!!!!! test failed
ASSERT errors for:  assert (res,ikd,1e-16)

  Location  |  Observed  |  Expected  |  Reason
   (1,1)        1.0558      0.28572      Abs err 0.77008 exceeds tol 1e-16 by 0.8
   (2,1)       0.74872      0.50842      Abs err 0.2403 exceeds tol 1e-16 by 0.2
   (3,1)       -0.08221     0.87548      Abs err 0.95769 exceeds tol 1e-16 by 1
   (4,1)        1.2175      0.70057      Abs err 0.5169 exceeds tol 1e-16 by 0.5
   (5,1)       0.41677     -0.019711     Abs err 0.43648 exceeds tol 1e-16 by 0.4
   (6,1)       0.20145      0.99222      Abs err 0.79078 exceeds tol 1e-16 by 0.8
   (7,1)        0.5431       0.2162      Abs err 0.3269 exceeds tol 1e-16 by 0.3
   (8,1)        0.751        1.3296      Abs err 0.57855 exceeds tol 1e-16 by 0.6
   (9,1)       0.36374       1.0843      Abs err 0.72056 exceeds tol 1e-16 by 0.7
   (10,1)       1.1464      0.28154      Abs err 0.86491 exceeds tol 1e-16 by 0.9
   (1,2)       0.43243       0.6434      Abs err 0.21097 exceeds tol 1e-16 by 0.2
   (2,2)       -0.99561     -0.39973     Abs err 0.59589 exceeds tol 1e-16 by 0.6
   (3,2)       -0.2929      0.56859      Abs err 0.86149 exceeds tol 1e-16 by 0.9
   (4,2)       -0.16635     -0.89053     Abs err 0.72418 exceeds tol 1e-16 by 0.7
   (5,2)       -0.93961    -0.010429     Abs err 0.92918 exceeds tol 1e-16 by 0.9
   (6,2)       0.46704     -0.018501     Abs err 0.48554 exceeds tol 1e-16 by 0.5
   (7,2)       0.028187     -0.42824     Abs err 0.45642 exceeds tol 1e-16 by 0.5
   (8,2)       0.42137      0.27893      Abs err 0.14245 exceeds tol 1e-16 by 0.1
   (9,2)       -0.44253     -1.2197      Abs err 0.77721 exceeds tol 1e-16 by 0.8
   (10,2)      0.49432      -1.2811      Abs err 1.7754 exceeds tol 1e-16 by 2

>>>>> processing C:\Octave\OC328E~1.2\mingw64\share\octave\packages\tisean-0.2.3\lyap_spec.m
***** test
 hen = henon(1000)(:,1);
 lyap_spec_res = [97 1.475854e+00 4.569062e-01 -4.080476e-01 -9.133782e-01];
 res = lyap_spec (hen, 'm', 4, 'n', 100, 'k', 50, 'invert');
 assert (res, lyap_spec_res, -1e-6);
!!!!! test failed
ASSERT errors for:  assert (res,lyap_spec_res,-1e-6)

  Location  |  Observed  |  Expected  |  Reason
    (2)         1.4757       1.4759      Rel err 9.3901e-05 exceeds tol 1e-06 by 9e-05
    (3)        0.44349      0.45691      Rel err 0.029359 exceeds tol 1e-06 by 0.03
    (4)        -0.39322     -0.40805     Rel err 0.03634 exceeds tol 1e-06 by 0.04
    (5)        -0.91465     -0.91338     Rel err 0.0013961 exceeds tol 1e-06 by 0.001
***** test
 ik = ikeda(1000);
 [lyap, pars] = lyap_spec (ik, 'm',4, 'n',100, 'k', 50,'invert');
 res = {lyap, pars.rel_err, pars.abs_err, pars.nsize, pars.nno, pars.ky_dim};
 assert (res, tisean_res, 1e-6);
!!!!! test failed
ASSERT errors for:  assert (cond {i},expected {i},tol)

  Location  |  Observed  |  Expected  |  Reason
    (2)         2.0453       1.7247      Abs err 0.32055 exceeds tol 1e-06 by 0.3
    (3)         1.4319       1.3443      Abs err 0.087588 exceeds tol 1e-06 by 0.09
    (4)        0.50466      0.60553      Abs err 0.10087 exceeds tol 1e-06 by 0.1
    (5)        0.073399     0.11075      Abs err 0.037347 exceeds tol 1e-06 by 0.04
    (6)        -0.2113      -0.16793     Abs err 0.043373 exceeds tol 1e-06 by 0.04
    (7)        -0.39715     -0.48399     Abs err 0.086844 exceeds tol 1e-06 by 0.09
    (8)        -0.63888     -0.68313     Abs err 0.044249 exceeds tol 1e-06 by 0.04
    (9)        -1.3454      -1.2791      Abs err 0.066241 exceeds tol 1e-06 by 0.07

ASSERT errors for:  assert (cond {i},expected {i},tol)

  Location  |  Observed  |  Expected  |  Reason
    (1)        0.096786     0.098795     Abs err 0.0020091 exceeds tol 1e-06 by 0.002
    (2)        0.049795     0.056064     Abs err 0.0062684 exceeds tol 1e-06 by 0.006

ASSERT errors for:  assert (cond {i},expected {i},tol)

  Location  |  Observed  |  Expected  |  Reason
    (1)        0.048269     0.046455     Abs err 0.0018141 exceeds tol 1e-06 by 0.002
    (2)        0.03738      0.039155     Abs err 0.001775 exceeds tol 1e-06 by 0.002

ASSERT errors for:  assert (cond {i},expected {i},tol)

  Location  |  Observed  |  Expected  |  Reason
     ()        0.90642      0.70356      Abs err 0.20286 exceeds tol 1e-06 by 0.2
shared variables   scalar structure containing the fields:

    tisean_res =
    {
      [1,1] =

         97.00000    1.72473    1.34430    0.60553    0.11075   -0.16793   -0.48399   -0.68313   -1.27913

      [1,2] =

         0.098795   0.056064

      [1,3] =

         0.046455   0.039155

      [1,4] =  0.70356
      [1,5] =  50
      [1,6] =  8
    }


On my ubuntu machine I don't get the failures.

Some are obviously just tolerance differences, some seem larger than that.


John Donoghue <lostbard>
Group Member

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach Files:
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by mmuetzel (Posted a comment)
  • -email is unavailable- added by dasergatskov (Posted a comment)
  • -email is unavailable- added by lostbard (Submitted the item)

    •  

    There are 0 votes so far. Votes easily highlight which items people would like to see resolved in priority, independently of the priority of the item set by tracker managers.

    Only group members can vote.

     

    No changes have been made to this item

    Back to the top

    Copyright © 2024 Free Software Foundation, Inc.
    Verbatim copying and distribution of this entire article is permitted in any medium, provided this notice is preserved.
    The Levitating, Meditating, Flute-playing Gnu logo is a GNU GPL'ed image provided by the Nevrax Design Team.
    Page source code

    Powered by Savane 3.14-d582.
    Corresponding source code