Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0 |

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

*To*: eigen@xxxxxxxxxxxxxxxxxxx*Subject*: Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0*From*: Benoit Jacob <jacob.benoit.1@xxxxxxxxx>*Date*: Sat, 17 Jul 2010 11:39:00 -0400*Dkim-signature*: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:received:in-reply-to :references:date:message-id:subject:from:to:content-type :content-transfer-encoding; bh=kLMIzjd0JJujMxCRzIGIwUdAv6BY9E66m5lfTqoqm/0=; b=JUIWeyRIwulekx8uNBxQoZK6bZgN2bvlWTLW0ek+YUhqFs1JgPQZWbXYAKucYhvXqY ZgfJoyrAimWJrNxOLPXp7Lo/Z3f7ZJb6dYtWjfL3ExPV+9wMIWMYR2WH9wWiiQvvalI/ ZSL/I8Vw4p3MDXG6+G5V1Wioq1cPyMzgbdWYQ=*Domainkey-signature*: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type:content-transfer-encoding; b=kOJFkDwspD0YyDNiTqAF3maKPyMO+Ywl9nbwVoheaT2+lh9dLhfymdvZqYteSd7ekX GYaa7XjR6jjwfskXOYq8edARxB2k1Y3RosPUBHBe4o8UZiKQOOrDO+fdJaDMXiVNZTjM S76j+NNd3tGJj5f3aMAb/f7x6rZE0t5TaBl0M=

2010/7/17 Aron Ahmadia <aja2111@xxxxxxxxxxxx>: > Hi List, > > Just responding to this discussion on condition number functions and > their availability in similar packages to Eigen.. > > LAPACK's expert driver routines do calculate the condition number of > the matrix (http://www.netlib.org/lapack/single/). Sorry I shouldn't have written that "LAPACK offers any condition number", instead I should have written that LAPACK doesn't standardize on a single notion of condition number. LAPACK does exposes the condition number of certain decompositions in certain expert solvers. However, these condition numbers are mutually inequivalent in different expert solvers, and are not the notion of condition number that you describe below. For example, the first link on that page mentioning condition numbers: http://www.netlib.org/lapack/single/sgesvx.f is actually doing "LU condition numbers". Namely, it is calling SGECON, http://www.netlib.org/lapack/single/sgecon.f which is doing this: * An estimate is obtained for norm(inv(A)), and the reciprocal of the * condition number is computed as * RCOND = 1 / ( norm(A) * norm(inv(A)) ). * * Arguments * ========= * * NORM (input) CHARACTER*1 * Specifies whether the 1-norm condition number or the * infinity-norm condition number is required: * = '1' or 'O': 1-norm; * = 'I': Infinity-norm. This is not the same as the condition number computed from the quotient of singular values (that would be the above formula with the spectral norm instead of L1 or Linf norm, which can be very different) > In Numerical > Linear Algebra, I would consider the condition number's definition > under the 2-norm as the largest singular value divided by the smallest > to be well-accepted (see Trefethen and Bau, "Numerical Linear > Algebra", Strang, "Introduction to Linear Algebra")

**Follow-Ups**:**Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0***From:*Benoit Jacob

**Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0***From:*Benoit Jacob

**References**:**[eigen] Problem inverting a Matrix4f with Eigen 2.0.0***From:*Robert Lupton the Good

**Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0***From:*Christoph Hertzberg

**Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0***From:*Gael Guennebaud

**Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0***From:*Robert Lupton the Good

**Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0***From:*Benoit Jacob

**Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0***From:*Aron Ahmadia

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0** - Next by Date:
**Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0** - Previous by thread:
**Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0** - Next by thread:
**Re: [eigen] Problem inverting a Matrix4f with Eigen 2.0.0**

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