Is there a way to use files.dat.ref?
As far as I know they are the reference results that you must compare against the results of the simulation with your version of CCX.
1 Like
If they are different, as in my case, then what to do?
You should investigate the reason of the difference. There is an script that make the comparation, but again it dosen´t tells you what is the error, only count the number of mismatching results. The script is called “compare” (without extension), don´t know if it works in Windows.
Of course I should and I did with diff that spitted out a lot of mismatch lines. I have no idea where to go further.
Which test had you try it? Can you copy and past the results (frd / dat)?
cubenewt test
cubenewt.dat
displacements (vx,vy,vz) for set NALL and time 0.9125000E-03
1 5.758572E-32 5.435292E-33 5.414972E-16
2 2.249961E-31 -1.227343E-31 5.414972E-16
3 -1.198752E-31 2.125827E-31 5.414972E-16
4 -1.373028E-31 -1.536206E-31 5.414972E-16
5 -6.699346E-32 1.283106E-31 5.414972E-16
6 -2.164502E-31 -2.293819E-31 5.414972E-16
7 2.138184E-31 2.496039E-31 5.414972E-16
8 4.819658E-32 -5.039250E-33 5.414972E-16
9 1.410040E-31 1.727160E-32 5.414972E-16
10 -1.396722E-32 -3.165642E-32 5.414972E-16
11 -1.283061E-31 -2.414720E-32 5.414972E-16
12 2.982393E-33 2.383543E-32 5.414972E-16
13 -1.543961E-32 6.608370E-32 5.414972E-16
14 1.533440E-32 -1.768496E-31 5.414972E-16
15 -3.364718E-33 2.319976E-31 5.414972E-16
16 8.644271E-33 -7.842873E-32 5.414972E-16
17 -1.420066E-31 4.718563E-32 5.414972E-16
18 2.762394E-32 2.564612E-32 5.414972E-16
19 1.312882E-31 -3.488285E-32 5.414972E-16
20 -2.081149E-32 -2.187946E-33 5.414972E-16
21 -1.204639E-31 -1.852048E-31 -5.414972E-16
22 -1.186037E-31 2.039159E-31 -5.414972E-16
23 5.477700E-32 -2.651656E-32 -5.414972E-16
24 6.569672E-33 3.334881E-33 -5.414972E-16
25 1.474045E-31 8.499714E-32 -5.414972E-16
26 2.572726E-31 5.036445E-32 -5.414972E-16
27 -2.015051E-31 -3.349787E-31 -5.414972E-16
28 -2.738900E-32 1.813812E-31 -5.414972E-16
29 -1.194154E-31 -9.124975E-33 -5.414972E-16
30 9.876246E-33 -3.823637E-32 -5.414972E-16
31 3.054211E-32 8.198990E-33 -5.414972E-16
32 2.247752E-32 3.805802E-32 -5.414972E-16
33 -9.611187E-33 -4.948125E-32 -5.414972E-16
34 -3.109875E-32 1.277655E-31 -5.414972E-16
35 2.594461E-32 -1.813868E-31 -5.414972E-16
36 6.924623E-33 9.172302E-32 -5.414972E-16
37 2.024544E-31 -2.815267E-32 -5.414972E-16
38 -6.722632E-33 7.695846E-32 -5.414972E-16
39 -1.145785E-31 2.496694E-32 -5.414972E-16
40 -1.824648E-32 -7.831867E-32 -5.414972E-16
displacements (vx,vy,vz) for set NALL and time 0.7588672E-02
1 9.569320E-30 2.261410E-29 3.745088E-14
2 4.047669E-29 -3.263358E-29 3.745088E-14
3 2.114798E-30 3.518758E-29 3.745088E-14
4 -6.076564E-29 -4.077744E-29 3.745088E-14
5 -1.008912E-30 1.085543E-29 3.745088E-14
6 -3.696503E-29 5.714924E-30 3.745088E-14
7 7.473368E-30 5.270769E-29 3.745088E-14
8 6.336176E-29 -4.110549E-29 3.745088E-14
9 2.762196E-29 -6.380963E-30 3.745088E-14
10 -1.396900E-29 1.416379E-30 3.745088E-14
11 -3.088982E-29 -2.973413E-30 3.745088E-14
12 1.040598E-29 4.178258E-30 3.745088E-14
13 3.713273E-30 1.801588E-29 3.745088E-14
14 1.504706E-30 -1.217824E-29 3.745088E-14
15 4.966234E-30 4.252395E-29 3.745088E-14
16 -9.833971E-32 -4.236518E-29 3.745088E-14
17 -1.638800E-29 7.229767E-30 3.745088E-14
18 9.738463E-30 -1.020959E-29 3.745088E-14
19 3.385314E-29 7.191163E-30 3.745088E-14
20 -1.226633E-29 2.974390E-30 3.745088E-14
21 -3.895896E-29 -9.880591E-30 -3.745088E-14
22 -1.267292E-29 4.300529E-29 -3.745088E-14
23 1.449084E-29 -2.145340E-30 -3.745088E-14
24 2.123560E-30 -3.244199E-31 -3.745088E-14
25 7.073837E-29 4.665816E-29 -3.745088E-14
26 1.040346E-29 -5.454777E-29 -3.745088E-14
27 -2.500628E-29 -1.227283E-28 -3.745088E-14
28 1.562745E-29 7.893146E-29 -3.745088E-14
29 -2.060586E-29 1.018007E-29 -3.745088E-14
30 -3.831429E-30 -8.564149E-30 -3.745088E-14
31 5.512950E-30 -1.440151E-30 -3.745088E-14
32 9.400589E-30 8.597588E-30 -3.745088E-14
33 1.281297E-29 1.624346E-29 -3.745088E-14
34 5.771796E-30 -7.916543E-30 -3.745088E-14
35 -9.748903E-30 -6.140114E-29 -3.745088E-14
36 3.977937E-30 4.033922E-29 -3.745088E-14
37 4.578097E-29 -3.438268E-31 -3.745088E-14
38 2.196488E-30 4.871939E-29 -3.745088E-14
39 -7.483647E-30 -2.169731E-29 -3.745088E-14
40 -1.622472E-29 -3.849025E-29 -3.745088E-14
displacements (vx,vy,vz) for set NALL and time 0.1000000E-01
1 1.219985E-29 -6.315731E-30 6.503250E-14
2 -6.053538E-30 -7.367495E-30 6.503250E-14
3 -7.765785E-30 2.977484E-29 6.503250E-14
4 -3.229449E-29 -1.685612E-29 6.503250E-14
5 5.257581E-30 -2.037498E-29 6.503250E-14
6 -3.622925E-31 3.174031E-29 6.503250E-14
7 5.739304E-30 4.616445E-29 6.503250E-14
8 3.056938E-29 -2.729276E-29 6.503250E-14
9 6.438768E-30 -4.409615E-30 6.503250E-14
10 -1.628106E-30 -3.086243E-30 6.503250E-14
11 -2.216738E-29 -4.262317E-31 6.503250E-14
12 4.805331E-30 7.436297E-30 6.503250E-14
13 8.352997E-30 -7.440745E-30 6.503250E-14
14 -2.737490E-30 1.809100E-29 6.503250E-14
15 3.372792E-30 3.235973E-29 6.503250E-14
16 -9.602432E-31 -2.768438E-29 6.503250E-14
17 5.813281E-30 8.960808E-30 6.503250E-14
18 -2.346139E-30 -1.861642E-29 6.503250E-14
19 1.601707E-29 7.033963E-30 6.503250E-14
20 -9.493372E-30 9.886231E-30 6.503250E-14
21 -3.407485E-29 -3.069507E-30 -6.503250E-14
22 3.738162E-30 4.227570E-29 -6.503250E-14
23 -2.295752E-29 -1.173253E-29 -6.503250E-14
24 -4.689795E-29 -2.500658E-29 -6.503250E-14
25 6.831822E-29 1.258067E-28 -6.503250E-14
26 2.935210E-29 -9.762193E-29 -6.503250E-14
27 -9.320664E-30 -1.535444E-28 -6.503250E-14
28 2.234005E-29 9.262142E-29 -6.503250E-14
29 -1.085318E-29 2.348920E-29 -6.503250E-14
30 7.429297E-31 -9.243263E-30 -6.503250E-14
31 -3.538091E-29 -2.071254E-29 -6.503250E-14
32 1.411527E-29 7.874307E-30 -6.503250E-14
33 2.050826E-29 5.788122E-29 -6.503250E-14
34 1.073638E-29 -3.116048E-29 -6.503250E-14
35 -9.389085E-30 -7.994200E-29 -6.503250E-14
36 -1.482867E-29 3.650389E-29 -6.503250E-14
37 5.315029E-29 8.783182E-30 -6.503250E-14
38 -5.749191E-30 6.903344E-29 -6.503250E-14
39 6.056541E-30 -2.350474E-29 -6.503250E-14
40 -2.123549E-29 -6.170550E-29 -6.503250E-14
cubenewt.frd
1C
1UUSER
1UDATE 26.september.2023
1UTIME 01:39:27
1UHOST
1UPGM CalculiX
1UVERSION Version 2.17
1UCOMPILETIME Thu Jul 23 21:43:50 CEST 2020
1UDIR
1UDBN
1UMAT 1EL
2C 40 1
-1 1 1.00000E+00 1.00000E+00 0.00000E+00
-1 2 1.00000E+00 0.00000E+00 0.00000E+00
-1 3 1.00000E+00 0.00000E+00 1.00000E+00
-1 4 1.00000E+00 1.00000E+00 1.00000E+00
-1 5 0.00000E+00 1.00000E+00 0.00000E+00
-1 6 0.00000E+00 0.00000E+00 0.00000E+00
-1 7 0.00000E+00 0.00000E+00 1.00000E+00
-1 8 0.00000E+00 1.00000E+00 1.00000E+00
-1 9 1.00000E+00 5.00000E-01 0.00000E+00
-1 10 1.00000E+00 0.00000E+00 5.00000E-01
-1 11 1.00000E+00 5.00000E-01 1.00000E+00
-1 12 1.00000E+00 1.00000E+00 5.00000E-01
-1 13 5.00000E-01 1.00000E+00 0.00000E+00
-1 14 5.00000E-01 0.00000E+00 0.00000E+00
-1 15 5.00000E-01 0.00000E+00 1.00000E+00
-1 16 5.00000E-01 1.00000E+00 1.00000E+00
-1 17 0.00000E+00 5.00000E-01 0.00000E+00
-1 18 0.00000E+00 0.00000E+00 5.00000E-01
-1 19 0.00000E+00 5.00000E-01 1.00000E+00
-1 20 0.00000E+00 1.00000E+00 5.00000E-01
-1 21 1.00000E+00 1.00000E+00 2.00000E+00
-1 22 1.00000E+00 0.00000E+00 2.00000E+00
-1 23 1.00000E+00 0.00000E+00 3.00000E+00
-1 24 1.00000E+00 1.00000E+00 3.00000E+00
-1 25 0.00000E+00 1.00000E+00 2.00000E+00
-1 26 0.00000E+00 0.00000E+00 2.00000E+00
-1 27 0.00000E+00 0.00000E+00 3.00000E+00
-1 28 0.00000E+00 1.00000E+00 3.00000E+00
-1 29 1.00000E+00 5.00000E-01 2.00000E+00
-1 30 1.00000E+00 0.00000E+00 2.50000E+00
-1 31 1.00000E+00 5.00000E-01 3.00000E+00
-1 32 1.00000E+00 1.00000E+00 2.50000E+00
-1 33 5.00000E-01 1.00000E+00 2.00000E+00
-1 34 5.00000E-01 0.00000E+00 2.00000E+00
-1 35 5.00000E-01 0.00000E+00 3.00000E+00
-1 36 5.00000E-01 1.00000E+00 3.00000E+00
-1 37 0.00000E+00 5.00000E-01 2.00000E+00
-1 38 0.00000E+00 0.00000E+00 2.50000E+00
-1 39 0.00000E+00 5.00000E-01 3.00000E+00
-1 40 0.00000E+00 1.00000E+00 2.50000E+00
-3
3C 2 1
-1 1 4 0 1
-2 1 2 3 4 5 6 7 8 9 10
-2 11 12 13 14 15 16 17 18 19 20
-1 2 4 0 1
-2 21 22 23 24 25 26 27 28 29 30
-2 31 32 33 34 35 36 37 38 39 40
-3
9999
Using 2.21 with the SPOOLES solver, I get exactly the same output as the reference.
Those displacements are all effectively zero and there’s always numerical error in zeros which would depend on the solver or compiler or optimization settings or some other detail. It probably doesn’t indicate any problem.
It’s not an easy problem to compare floating point values when there’s any non-determinism since whether they are “the same” or not depends on the problem. But I don’t think anything in CCX works on the scale of 1e-14 so it should be pretty safe.
1 Like