Good shooting! I introduced that bug in a recent change. I will fix it
for the next release.
Thank you for taking the time to run this down.
Henry Rich
On 7/14/2018 6:50 PM, Bill Heagy wrote:
It looks like a problem with "^":
2^0.5
1.41421
2^0.5j3
1.41421
For complex exponent with a real part of 0.5 (sqrt), it appears to
ignore the imaginary part.
On 2018-07-02 04:52 PM, Bill Heagy wrote:
I haven't checked it all out, but I suspect the problem is that the
Stirling formula is asymmptotic, and has a larger error near 1.
On 2018-07-02 12:20 PM, Bill Heagy wrote:
First four runs on 32bit linux:
1j7.26
1j8.53
1j_3.14
1j_2.67
64bit machine gives different imaginary parts, but real part=1.
Looks like 1j(whatever) may be causing the problem.
On 2018-07-02 07:45 AM, Henry Rich wrote:
When this fails, what is the value in z that gives the error?
BTW, Bill, that error with gcc that you reported a few weeks back
turned out to be a program error that showed up only on that
compiler. If you switch back to gcc now it should be gone.
Henry Rich
On 6/29/2018 9:59 PM, Bill Heagy wrote:
If I run test/g410.ijs often enough, I get:
100 RUN2 'g410'
.......
1 1 1 1 1 1 1 1
(g = gamma) x=: _4.1 _4.9 4.1 4.9 j./_0.4 _0.5 0.4 0.5 _4 _5 4 5
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
(g =!.1e_12 gamma) 1.5
1
(g =!.1e_12 gamma) 1r3 1r2 2r3
1 1 1
(g =!.1e_12 gamma) x=: 1e_2*_7e2+?2 20$2e3
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
(g =!.1e_12 gamma) x=: 3e_2*_7e2+?2 20$2e3
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
(g =!.1e_12 gamma) z=:j./1e_2*_7e2+?2 20$2e3
1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
|assertion failure
|[-103] /home/wheagy/git/jsource/test/g410.ijs
Which specific value is wrong changes from failure to failure, but
it's always this line, and it appears to happen about every 50 to
100 iterations. This happens on debian linux 32 and 64nonavx,
with the latest beta release and when compiled from the latest
version on github.
----------------------------------------------------------------------
For information about J forums see
http://www.jsoftware.com/forums.htm
---
This email has been checked for viruses by AVG.
https://www.avg.com
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
