Re: [eigen] Sparse Matrix Support |

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

*To*: eigen@xxxxxxxxxxxxxxxxxxx*Subject*: Re: [eigen] Sparse Matrix Support*From*: Trevor Irons <trevorirons@xxxxxxxxx>*Date*: Fri, 18 Jun 2010 09:06:34 -0600*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; bh=WuGeqDHaclCslktQPGTQA57oDnxX5V/trin+9L1L/cU=; b=P9BvhUaEhOdPp9YwBpixdQoFv8LN0gjoy1oBOnlQ9DVd7Hr1NmifAu+IAHZWWj/tJi tEE6rKKEV92t+HIRDIZuvPyJHVhWlvBJq6fYfGbylUqppKtqfULHAhPJAloKyPdNnNMa sld2sJaouoPTJEbiX4tLY3fI0Dz7Pw/QbJEoM=*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; b=kjn+YtSoFALwNvxF2EIdBpBHpyhjn+1rSPYFykXyVSFzeAf8NTt1/ywZHtNZcv9v52 JFTlSLpVUYVqgZbo0Y2H7V9ZEa+qe0g65kDWBnVEstkZ8v3om7YQrGV2k1/QbNfPFYU0 UFNOu5SYjAjqC0sjCsDiHHVyom1KNSW+wtAZY=

> And last but not the least, have you given any thought to supportingThis is very tricky and I'm still not sure how to implement it the

> block sparse matrices?

most efficient and simple way. Indeed, the idea is of course to reuse

the SparseMatrix class and let it stores small fixed size dense Matrix

for the scalar type. For instance, let's assume a 300x300 block matrix

represented as a 100x100 standard sparse matrix of Matrix3d. The pb is

that depending on the context, sometimes we want to see it as a

300x300 matrix of double, and sometimes as a 100x100 matrix of

Matrix3d blocks. So this probably requires a lot of work to get it

right.

It would also be really useful to support block sparse where the blocks are themselves sparse too. You run into this a lot in solutions of pde's.

**References**:**[eigen] Sparse Matrix Support***From:*Sameer Agarwal

**Re: [eigen] Sparse Matrix Support***From:*Gael Guennebaud

**Re: [eigen] Sparse Matrix Support***From:*Sameer Agarwal

**Re: [eigen] Sparse Matrix Support***From:*Gael Guennebaud

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] documentation: the long tutorial** - Next by Date:
**Re: [eigen] Sparse Matrix Support** - Previous by thread:
**Re: [eigen] Sparse Matrix Support** - Next by thread:
**Re: [eigen] Sparse Matrix Support**

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