Opened 7 years ago

Closed 7 years ago

#1428 closed enhancement (fixed)

[with patch, with positive review] add SVD method to matrix_complex_double_dense

Reported by: mhansen Owned by: mhansen
Priority: major Milestone: sage-2.10
Component: linear algebra Keywords:
Cc: dfdeshom@… Merged in:
Authors: Reviewers:
Report Upstream: Work issues:
Branch: Commit:
Dependencies: Stopgaps:

Description

Hallo!

I tried to compute the SVD of a complex matrix (module
matrix.matrix_complex_double_dense), but I didn't found a function to
do so. However, real matrices (module matrix.matrix_real_double_dense)
support it. Is there really no way to compute a complex SVD? If I
remember correctly, at least the underlying library GSL supports
complex SVDs... What would I have to do to integrate those functions
into Sage?

Sander

Basically all that needs to be done is to translate the real code over to the complex case.

Attachments (1)

1428.patch (9.1 KB) - added by mhansen 7 years ago.

Download all attachments as: .zip

Change History (5)

Changed 7 years ago by mhansen

comment:1 Changed 7 years ago by mhansen

  • Owner changed from was to mhansen
  • Status changed from new to assigned
  • Summary changed from add SVD method to matrix_complex_double_dense to [with patch] add SVD method to matrix_complex_double_dense

comment:2 Changed 7 years ago by rlm

  • Summary changed from [with patch] add SVD method to matrix_complex_double_dense to [with patch, needs review] add SVD method to matrix_complex_double_dense

comment:3 Changed 7 years ago by dfdeshom

  • Cc dfdeshom@… added
  • Summary changed from [with patch, needs review] add SVD method to matrix_complex_double_dense to [with patch, with positive review] add SVD method to matrix_complex_double_dense

The patch looks great. I would suggest making an option that would just return S, instead of the tuple (U,S,V') since people that use this method tend to care more about S than anything else.

comment:4 Changed 7 years ago by mabshoff

  • Resolution set to fixed
  • Status changed from assigned to closed

Merged in Sage 2.10.alpha3

Note: See TracTickets for help on using tickets.