ygrenier / SwissEphNet

Swiss Ephemeris for .Net. A Swiss Ephemeris portage for .Net.
Other
77 stars 35 forks source link

Problems using files jpl: DE406 and DE431 #7

Closed blucobalto closed 10 years ago

blucobalto commented 10 years ago

Hi, I think there is a problem using the files jpl: DE406 and DE431. The results are not correct.

I have tested the results jpl using this: 1) SweWin.exe (swissEph); 2) http://www.astro.com/swisseph/swetest.htm and comparing them with SweWin.exe (SwissEphNet).

Thanks a lot, blu

ygrenier commented 10 years ago

Hi,

Could you indicate some values for testing it's more easy for me.

Thanks,

Yanos

blucobalto commented 10 years ago

Hi, for the test: 1) Run Swewin.exe (SwissEphNet). 2) Leave the default data; 3) Select "Swiss Eph"; 4) Run. 5) All ok; 6) Select "DE406"; 7) Run.

The values ​​do not change. I expect, however, that the values ​​of the last decimal places are changed. I tested the problem with SweWin.exe (SwissEph). In this case the last decimal places change.

Thanks a lot, Antonio

ygrenier commented 10 years ago

Hi,

I can't reproduce the issue ! I have the same behavior between the two applications.

Perhaps it's a problem with the ephemeris path that I just correct (a few minutes).

Now the application search the ephemeris files in "c:\sweph\ephe" path like original application.

Try this.

Thanks.

Yanos

blucobalto commented 10 years ago

Hi, I downloaded the latest version from Github. I use the default path: "c:\sweph\ephe".

I have two problems with SweWin.exe (SwissEphNet): 1) If I choose "Moshier" or "Swiss Eph" or "DE406" I always get the same results. I expect that the last decimal places are different.

2) I tried to debug with the parameter "DE406". In the file: SweJpl.cs - function: swi_open_jpl_file the line "retc = read_const_jpl (ss, ref serr);" returns -2.

I hope I have not created any confusion.

Thanks. Antonio

ygrenier commented 10 years ago

How I'm a stupid guy !! I correct the default search path in SweWin, look at the last code version.

Normally SwissEph and JPL have same behavior and Moshier have differents results now.

Sorry, I think I was so tired yesterday when I define the default path :(

Yanos

blucobalto commented 10 years ago

No problem. Stuff happens. You're doing so much already.

The problem between Moshier and "Swiss Eph / DE406" is solved.

The problem between the "Swiss Eph" and "DE406" remains. Test If I enter the following parameters: Date: 07/05/2014 Time: 10:31:49 Type: JPL Ephemeris DE406 I get:

SweWin (original) Saturn 20 sc 16' 5.9163 2°29' 5.2078 8.901383263 -0° 4'30.0268 4.4696

SweWin (Net) Saturn 20 sc 16' 6.0006 2°29' 5.2541 8,901385122 -0° 4'30.0267 4,4696

Thanks a lot, Antonio

ygrenier commented 10 years ago

When I do the test I have exactely the same values :

Original : Saturn 20 sc 16' 6.0006 2°29' 5.2541 8.901385122 -0° 4'30.0267 4.4696 .Net : Saturn 20 sc 16' 6.0006 2°29' 5.2541 8,901385122 -0° 4'30.0267 4,4696

For all planets. The only différences i have between Original et .Net version is for Moshier method on Mean/True node and Mean/Osc. apogee.

Perhaps we don't have the same files that explain the difference calculations between our tests ! Personnly I use the last files downloaded from "ftp://ftp.astro.ch/pub/swisseph/ephe/" but only the files, the Subfolders are missing on my computer and I added directly three asteroid files. Here my list :

03/05/2014  20:37               235 .message
03/05/2014  20:34           119 336 fixstars.cat
03/05/2014  20:36           107 400 fixstars.cat.1
03/05/2014  20:40           101 457 fixstars.txt
03/05/2014  20:39             2 568 LICENSE.TXT
03/05/2014  20:41             2 538 md5sum.txt
03/05/2014  20:37             6 249 README.TXT
04/05/2014  19:20           105 656 se00433s.se1
04/05/2014  18:19            41 407 se03045s.se1
04/05/2014  19:21            23 009 se07066s.se1
03/05/2014  20:35           225 432 seasm06.se1
03/05/2014  20:37           222 642 seasm102.se1
03/05/2014  20:35           222 363 seasm108.se1
03/05/2014  20:33           222 388 seasm114.se1
03/05/2014  20:36           226 922 seasm12.se1
03/05/2014  20:33           222 577 seasm120.se1
03/05/2014  20:40           222 724 seasm126.se1
03/05/2014  20:35           223 235 seasm132.se1
03/05/2014  20:41           228 057 seasm18.se1
03/05/2014  20:34           228 975 seasm24.se1
03/05/2014  20:34           229 539 seasm30.se1
03/05/2014  20:40           229 518 seasm36.se1
03/05/2014  20:35           228 933 seasm42.se1
03/05/2014  20:34           227 954 seasm48.se1
03/05/2014  20:38           226 670 seasm54.se1
03/05/2014  20:40           225 519 seasm60.se1
03/05/2014  20:40           224 589 seasm66.se1
03/05/2014  20:40           223 760 seasm72.se1
03/05/2014  20:36           223 297 seasm78.se1
03/05/2014  20:36           222 929 seasm84.se1
03/05/2014  20:39           222 879 seasm90.se1
03/05/2014  20:36           222 694 seasm96.se1
03/05/2014  20:39         7 528 290 seasnam.txt
03/05/2014  20:38         1 682 841 seasnam2.txt
03/05/2014  20:33           224 354 seas_00.se1
03/05/2014  20:35           223 809 seas_06.se1
03/05/2014  20:37           229 041 seas_102.se1
03/05/2014  20:33           229 615 seas_108.se1
03/05/2014  20:36           229 816 seas_114.se1
03/05/2014  20:41           223 396 seas_12.se1
03/05/2014  20:38           229 960 seas_120.se1
03/05/2014  20:35           230 536 seas_126.se1
03/05/2014  20:37           229 995 seas_132.se1
03/05/2014  20:34           229 328 seas_138.se1
03/05/2014  20:37           228 501 seas_144.se1
03/05/2014  20:36           227 639 seas_150.se1
03/05/2014  20:37           226 805 seas_156.se1
03/05/2014  20:34           225 935 seas_162.se1
03/05/2014  20:34           223 002 seas_18.se1
03/05/2014  20:35           222 673 seas_24.se1
03/05/2014  20:34           222 208 seas_30.se1
03/05/2014  20:40           221 567 seas_36.se1
03/05/2014  20:34           221 190 seas_42.se1
03/05/2014  20:41           221 138 seas_48.se1
03/05/2014  20:40           221 392 seas_54.se1
03/05/2014  20:40           221 950 seas_60.se1
03/05/2014  20:37           222 906 seas_66.se1
03/05/2014  20:39           223 798 seas_72.se1
03/05/2014  20:35           224 916 seas_78.se1
03/05/2014  20:41           226 330 seas_84.se1
03/05/2014  20:40           227 170 seas_90.se1
03/05/2014  20:36           228 105 seas_96.se1
03/05/2014  20:36           107 546 sefstars.txt
03/05/2014  20:34         1 332 945 semom06.se1
03/05/2014  20:33         1 457 499 semom102.se1
03/05/2014  20:34         1 466 886 semom108.se1
03/05/2014  20:40         1 476 216 semom114.se1
03/05/2014  20:39         1 347 197 semom12.se1
03/05/2014  20:37         1 490 638 semom120.se1
03/05/2014  20:40         1 499 802 semom126.se1
03/05/2014  20:36         1 514 331 semom132.se1
03/05/2014  20:36         1 356 379 semom18.se1
03/05/2014  20:35         1 365 357 semom24.se1
03/05/2014  20:33         1 371 740 semom30.se1
03/05/2014  20:40         1 380 820 semom36.se1
03/05/2014  20:41         1 390 210 semom42.se1
03/05/2014  20:35         1 398 880 semom48.se1
03/05/2014  20:37         1 403 463 semom54.se1
03/05/2014  20:36         1 409 454 semom60.se1
03/05/2014  20:36         1 414 798 semom66.se1
03/05/2014  20:34         1 423 009 semom72.se1
03/05/2014  20:40         1 427 325 semom78.se1
03/05/2014  20:36         1 433 946 semom84.se1
03/05/2014  20:41         1 443 026 semom90.se1
03/05/2014  20:35         1 453 178 semom96.se1
03/05/2014  20:41         1 322 594 semo_00.se1
03/05/2014  20:39         1 316 045 semo_06.se1
03/05/2014  20:39         1 411 422 semo_102.se1
03/05/2014  20:41         1 416 577 semo_108.se1
03/05/2014  20:39         1 420 469 semo_114.se1
03/05/2014  20:36         1 305 931 semo_12.se1
03/05/2014  20:40         1 424 162 semo_120.se1
03/05/2014  20:37         1 432 625 semo_126.se1
03/05/2014  20:34         1 442 619 semo_132.se1
03/05/2014  20:39         1 450 003 semo_138.se1
03/05/2014  20:38         1 459 433 semo_144.se1
03/05/2014  20:38         1 474 021 semo_150.se1
03/05/2014  20:37         1 495 817 semo_156.se1
03/05/2014  20:38         1 505 928 semo_162.se1
03/05/2014  20:39         1 304 771 semo_18.se1
03/05/2014  20:33         1 308 123 semo_24.se1
03/05/2014  20:41         1 312 208 semo_30.se1
03/05/2014  20:35         1 326 257 semo_36.se1
03/05/2014  20:35         1 339 659 semo_42.se1
03/05/2014  20:37         1 349 937 semo_48.se1
03/05/2014  20:37         1 358 320 semo_54.se1
03/05/2014  20:35         1 366 885 semo_60.se1
03/05/2014  20:40         1 375 610 semo_66.se1
03/05/2014  20:34         1 383 983 semo_72.se1
03/05/2014  20:41         1 391 143 semo_78.se1
03/05/2014  20:37         1 397 403 semo_84.se1
03/05/2014  20:36         1 402 879 semo_90.se1
03/05/2014  20:36         1 407 105 semo_96.se1
03/05/2014  20:40             4 832 seorbel.txt
03/05/2014  20:37               389 sepl.memo
03/05/2014  20:41           484 266 seplm06.se1
03/05/2014  20:39           630 636 seplm102.se1
03/05/2014  20:34           641 587 seplm108.se1
03/05/2014  20:36           650 870 seplm114.se1
03/05/2014  20:40           483 990 seplm12.se1
03/05/2014  20:35           656 541 seplm120.se1
03/05/2014  20:34           662 357 seplm126.se1
03/05/2014  20:38           665 907 seplm132.se1
03/05/2014  20:35           488 142 seplm18.se1
03/05/2014  20:35           489 846 seplm24.se1
03/05/2014  20:34           493 416 seplm30.se1
03/05/2014  20:35           500 226 seplm36.se1
03/05/2014  20:41           513 210 seplm42.se1
03/05/2014  20:34           521 081 seplm48.se1
03/05/2014  20:35           531 230 seplm54.se1
03/05/2014  20:34           541 501 seplm60.se1
03/05/2014  20:38           553 441 seplm66.se1
03/05/2014  20:39           565 716 seplm72.se1
03/05/2014  20:37           582 598 seplm78.se1
03/05/2014  20:34           609 773 seplm84.se1
03/05/2014  20:41           617 973 seplm90.se1
03/05/2014  20:38           625 668 seplm96.se1
03/05/2014  20:33           483 871 sepl_00.se1
03/05/2014  20:34           482 953 sepl_06.se1
03/05/2014  20:34           566 462 sepl_102.se1
03/05/2014  20:37           574 871 sepl_108.se1
03/05/2014  20:40           589 817 sepl_114.se1
03/05/2014  20:38           483 757 sepl_12.se1
03/05/2014  20:40           598 274 sepl_120.se1
03/05/2014  20:39           605 675 sepl_126.se1
03/05/2014  20:38           619 007 sepl_132.se1
03/05/2014  20:37           623 354 sepl_138.se1
03/05/2014  20:35           628 335 sepl_144.se1
03/05/2014  20:36           644 186 sepl_150.se1
03/05/2014  20:33           652 910 sepl_156.se1
03/05/2014  20:35           658 153 sepl_162.se1
03/05/2014  20:34           484 055 sepl_18.se1
03/05/2014  20:40           483 402 sepl_24.se1
03/05/2014  20:39           484 373 sepl_30.se1
03/05/2014  20:40           484 686 sepl_36.se1
03/05/2014  20:40           483 628 sepl_42.se1
03/05/2014  20:34           485 517 sepl_48.se1
03/05/2014  20:34           488 956 sepl_54.se1
03/05/2014  20:33           493 649 sepl_60.se1
03/05/2014  20:36           501 136 sepl_66.se1
03/05/2014  20:37           508 815 sepl_72.se1
03/05/2014  20:34           516 226 sepl_78.se1
03/05/2014  20:38           530 714 sepl_84.se1
03/05/2014  20:34           543 316 sepl_90.se1
03/05/2014  20:39           554 703 sepl_96.se1
             166 fichier(s)      119 028 228 octets

Yanos

blucobalto commented 10 years ago

Maybe I understand the problem. The software SwissEph reads three types of files.

Swiss Ephemeris The core part of Swiss Ephemeris is a compression of the JPL-DE431 ephemeris, All Swiss Ephemeris files have the file suffix. Se1.

Moshier This is a semi-analytical approximation of the JPL planetary and lunar ephemerides DE404, developed by Steve Moshier.

JPL This is the full precision state-of-the-art ephemeris. It provides the highest precision and is the basis of the Astronomical Almanac. Filename is: de406.eph or de431.eph.

Problem When I select the type in SweWin.exe "JPL DE406" the software to read the file data jpl de406.eph. If not found try to use the Swiss Eph file (*.se1). I noticed that in your folder does not exist the file de406.eph (194MB) or de431.eph (2700MB). That's why you see the same results.

Antonio

ygrenier commented 10 years ago

OK it's that I supposed.

I looking for that tomorrow.

Thanks.

Yanos

ygrenier commented 10 years ago

OK I just fix the DE406 file reading.

I waiting to download the DE431 file for more tests.

Yanos

ygrenier commented 10 years ago

Now I just fix an error for reading the DE431 file that very large and crash the applications.

My tests are now same than the test page. So I close this issue.