bugGNU Octave - Bugs: bug #53351, A = expint(x>32+yi) doesn't...

 
 

bug #53351: A = expint(x>32+yi) doesn't result any imaginary number

Submitter:  Mina Mourice <mourice>
Submitted:  Thu 15 Mar 2018 12:47:18 PM UTC
   
 
Category:  Octave Function Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Inaccurate Result
Status:  Fixed Assigned to:  None
Originator Name:  Open/Closed:  * Closed
Release:  * 4.2.1 Operating System:  * Any
Fixed Release:  None Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Fri 16 Mar 2018 09:11:18 AM UTC, comment #2: 

Dear Mina,
the new version of expint is already on the stable branch of Octave, and it solves this issue:

octave-cli-4.3.0+:1> version
ans = 4.3.0+
octave-cli-4.3.0+:2> expint(32.5 + 1i)
ans =  1.1811e-16 - 1.9663e-16i
octave-cli-4.3.0+:3> expint(44 + 1i)
ans =  9.0188e-22 - 1.4758e-21i

Michele Ginesi <m_ginesi>
Fri 16 Mar 2018 08:35:18 AM UTC, comment #1: 

Michele did a lot of work last summer on special functions. It is not yet on default branch, I think, but it should be in a few days. Michele, can you check if your job fixed this problem?

Marco Caliari <caliari>
Group Member
Thu 15 Mar 2018 12:47:18 PM UTC, original submission:  

In Exponential integral :
For complex values , the argument of the complex must be less than pi
This condition is satisfied in the condition of real number(x) > 33 and imaginary (y) = any number in a suitable range


Using MATLAB :

>> A = expint(32.5+1i)

A =

   1.1811e-16 - 1.9663e-16i

>> B = expint(44+1i)

B =

   9.0188e-22 - 1.4758e-21i

Using Octave GNU Release 4.2.1 :

>> A = expint(32.5+1i)

܀A =   1.1811e-016

>> B = expint(44+1i)

y㌒ऀB =   9.0188e-022

Mina Mourice <mourice>

 

(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 mtmiller (Updated the item)
  • -email is unavailable- added by m_ginesi (Posted a comment)
  • -email is unavailable- added by caliari (Posted a comment)
  • -email is unavailable- added by caliari
  • -email is unavailable- added by mourice (Submitted the item)
  • -email is unavailable- added by mourice (Mourice)
  •  

    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.

     

    Follow 4 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2018-03-16 mtmiller Open/ClosedOpen Closed
    2018-03-16 caliari StatusNone Fixed
    2018-03-16 caliari Carbon-Copy- Added -email is unavailable-
    2018-03-15 mourice Carbon-Copy- Added -email is unavailable-

    Back to the top

    Powered by Savane 3.13-02a9.
    Corresponding source code