bugGNU Octave - Bugs: bug #60577, Detect (sparse) transposed lower...

 
 

bug #60577: Detect (sparse) transposed lower triangular linear systems

Submitter:  Marco Caliari <caliari>
Submitted:  Tue 11 May 2021 06:55:54 AM UTC
   
 
Category:  Interpreter Severity:  2 - Minor
Priority:  5 - Normal Item Group:  Performance
Status:  None Assigned to:  None
Originator Name:  Open/Closed:  * Open
Release:  * 6.1.0 Operating System:  * GNU/Linux
Fixed Release:  None Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Tue 11 May 2021 06:55:54 AM UTC, original submission:  

If I run the following code with M sparse or full


m = 500;
M = toeplitz(sparse([1,1],[1,2],[4,1],1,m));
R = chol(M); % try it with chol(full(M));
b = randn(m,1);
tic
for i = 1:10000
  R\b;
end
toc
tic
for i = 1:10000
  R'\b;
end
toc


I observe the transposed system 1.5x slower than the untransposed one in the sparse case, while there is no difference in the full case. I have two possible explanations:
1) in the full case the interpreter detects that R' is lower triangular and avoids transposition
2) in the full case the computational cost hides the cost related to transposition

Would it be possible to detect in any case that R' is lower triangular and avoid transposition?

Marco Caliari <caliari>
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 caliari (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

    Powered by Savane 3.13-d3ae.
    Corresponding source code