Re: [eigen] Complex symmetric Cholesky factorization |

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

*To*: eigen <eigen@xxxxxxxxxxxxxxxxxxx>*Subject*: Re: [eigen] Complex symmetric Cholesky factorization*From*: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>*Date*: Tue, 1 Sep 2015 09:13:21 +0200*Dkim-signature*: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :content-type; bh=FLuAkZGRg1XSExo3CY4bCTpa/uTnvBkdH6niOORns2k=; b=GIcWjhBf+BjmVvBOYOYBSgJ2B/bt8XTsqn13aP6CeGrVaYbZ7DHmyHa4ZUzMmzrgQq 2ZOGHcQKboLV5xtO3pu/kyi23/muEx+Hq2k+EePdExSBS4/btG4b4+itIOMtnqTb0o4y arBNGvYizHHZgMBANuLoYqRgmrr91xxU1Nwx8fJ2umf0XZgXA5zTA5G+bc/1qEX2B5/Q wjUJE5GQqHlo1I1lig7EDUNEV81uKrIWRU5zLBLB0efbUzsbIrzgsjOugYXT+zz0O6oU t1AY4rKbgPdIj+6PpE5JhIlv4S9ekw3T1gSBErGbIde6UuYtl920iHqo/W7Ym+IR/qUL 6q7A==

Yes, it is for MKL's Pardiso.

gael

On Mon, Aug 17, 2015 at 10:10 PM, John Daily <John.Daily@xxxxxxxxxxxxxxxxx> wrote:

Thanks.Does Eigen's Pardiso wrapper only work with the MKL version of Pardiso?On 17 August 2015 at 13:46, Gael Guennebaud <gael.guennebaud@xxxxxxxxx> wrote:On Mon, Aug 17, 2015 at 6:33 PM, John Daily <John.Daily@xxxxxxxxxxxxxxxxx> wrote:We already have a SelfAdjoint and Symmetric enums to distinguish between the two cases, so perhaps an idea would be to allow to pass Symmetric to the second template parameter, like:[Simplicial]L[D]LT<MatrixType, {{Lower,Upper,Lower|Upper}[|Symmetric] , Symmetric}>I don't quite follow your notation here for the Symmetric enum, did you mean one of those to be SelfAdjoint as written?The trick is that the default is Lower which is also equal to 0. Therefore, one can simply pass "Symmetric" as a synonym of "Lower|Symmetric". This is why I put Symmetric twice. So basically, the test for symmetric instead of self-adjoint should be "(Options&Symmetric)==Symmetric" (after renaming the UpLo parameter to Options).Actually, this solution is already the one we adopted for PardisoLLT<> and PardisoLDLT<>.I'll take a look at that. I'll try and make a branch achieving this and will let you know when I have something that can be reviewed.Great!In a second time, for consistency, we should also add a symmetricView<>() method to [Sparse]MatrixBase<> returning a SelfAdjointView<...,{Upper,Lower}[|Symmetric]>. Likewise, we could also think about adding templated variants of adjoint() / selfadjointView() to ease writing generic code with true compile-time branches...I don't have enough familiarity with the sparse internals to know how difficult this part will be. I'll start with the above.Sure, this part might be a bit tricky.gael

**Messages sorted by:**[ date | thread ]- Next by Date:
**Re: [eigen] Performance of SparseMatrix** - Next by thread:
**Re: [eigen] Performance of SparseMatrix**

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