[eigen] How to speed up svd? |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: [eigen] How to speed up svd?
- From: Hugh Perkins <hughperkins@xxxxxxxxx>
- Date: Fri, 5 Oct 2012 20:25:59 +0800
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=mime-version:sender:date:x-google-sender-auth:message-id:subject :from:to:content-type; bh=oE0P7F0tms6xqk7uhwf1EWBoP78pJSrosX01nhwhvNU=; b=pUrPbCBnmPV/1XcoNsLLkGSJWzxVP7DyO21uMka2bQ2Iqm9oqHsKAhXOA3896OuJ6t rhrHri0V37VvgY/R8s20wrSrTFEAgguOrWUkQTO6VBCrhS3K1+OupN93wEoT7RS/sJRG r6ZVNDItHsLD8gt9SpDt4IHVcAbqRBNZ/OKBPQypwXAHEyG+cNnLChppBOUUAKYKG/ta VuSAxk46yFQ1Go3tvDXpOAl9wD/dYZqGmGy1yqg8EyLO+5CHyIS6qFO835j2Nqp51UzS Prl+TRyMvvYWETzgS0cmbyJ7yoak1NTqNXMWsUCTjuITRGYv2aRk2ux6/+yCeXMN1hQ7 K/7A==
Hi,
I run SVD on a 500*1000 dense double matrix of random number, with
thin U and V, and it takes about 26 seconds. The code looks like:
JacobiSVD<MatrixXd,HouseholderQRPreconditioner> svd(In,
ComputeThinU | ComputeThinV);
On matlab it takes about 1.1 seconds, though admittedly the solution
doesn't seem terribly accurate.
What are possible reasons for this disparity, and is there anything I
can do to improve the situation?
Hugh