[eigen] Complex symmetric Cholesky factorization |

[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]

*To*: eigen@xxxxxxxxxxxxxxxxxxx*Subject*: [eigen] Complex symmetric Cholesky factorization*From*: John Daily <John.Daily@xxxxxxxxxxxxxxxxx>*Date*: Thu, 13 Aug 2015 13:28:27 -0600

Greetings, certain symmetric complex (non-Hermitian) matrices admit an LL^T/LDL^T factorization[1] (rather than an LL^H LDL^H). I am working with sparse matrices and the code changes to support this are minimal and are all encompassed in Simplicial_Cholesky_impl.h, in the factorize method, as well as the updated solve calling transpose rather than selfadjointView. I have two questions:

1) Is there any desire to add this functionality to Eigen, or should I just keep these changes local? I would prefer to add this upstream is possible. The actual code changes are minimal, three lines or so.

2) If so what would be an acceptable way to do this, that would not cause surprise to users? It would seem that a separate template parameter or even function call would make this explicit, such as instead of .factorize() calling .factorizeComplexSymmetric() or something.

Thanks,

[1] Bereux, N., 2004, A cholesky algorithm for some complex symmetric systems: CMAP 515, Ecole Polytechnique.

(Available at http://www.cmap.polytechnique.fr/~natacha/CMAP TR 515.pdf).

**Follow-Ups**:**Re: [eigen] Complex symmetric Cholesky factorization***From:*Gael Guennebaud

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] eigen-3.2.5 and g++ 5.1.1** - Next by Date:
**[eigen] Complex symmetric Cholesky factorization** - Previous by thread:
**Re: [eigen] eigen-3.2.5 and g++ 5.1.1** - Next by thread:
**Re: [eigen] Complex symmetric Cholesky factorization**

Mail converted by MHonArc 2.6.19+ | http://listengine.tuxfamily.org/ |