
Boost : 
From: walter_at_[hidden]
Date: 20011123 13:13:11
 In boost_at_y <mailto:boost_at_y>..., Peter Schmitteckert (boost) <
boost_at_s <mailto:boost_at_s>...> wrote:
> Hello,
>
> On Wednesday 21 November 2001 11:58, walter_at_g <mailto:walter_at_g>...
wrote:
>
> > You've lost me. Could you please explain or give a reference?
> Please see below.
Thanks.
> 
> Outline of my most cpu intensive part
> 
>
> In my application I have to iterativly diagonalize
> large sparse matrices which involve the MatrixVector
> product of a matrix C with a vector x.
> To improve performance one can use a representation
> of C where the vector space of C is a tensor product of two
> vector spaces V and W.
> C is now a sum of operators A_i \otimes B_i, where A_i (B_i) acts on
> V (W) only. The total dimension of C is equal the product of the
> dimensions of V and W. A basis of the vector space of the product
space of
> V and W can be represented by a dyadic product of basis states of V
and W,
> i.e. if v * w^T.
I guess, this is a rather domain specific vector space decomposition
(if dim (C) is prime for example, the decomposition seems to be
useless).
> This representation has the advantage that the (sparse) matrix
vector
> multiplication can be represented by BLAS3 operations instead of
BLAS2
> operations using dense matrices.
Ok.
> In case you're not lost again, The spaces V and W themselves can
> be represented by by direct sums of subspaces V_l (W_k).
> The matrices A can now be represente by blocks A_lm, where
> A_lm is a mapping from V_l to V_m.
I'll believe that ;)
> If you're still interested in details you may look at my thesis,
> <http://www.Physik.UniAugsburg.DE/~peters/thesis/index.html>
> chapter 5.5 .
I guess, your list of matrices conceptually can be viewed as tensor
or as (domain specific?) special sparse matrix format. But may be,
that's a matter of taste.
Regards
Joerg
Boost list run by bdawes at acm.org, gregod at cs.rpi.edu, cpdaniel at pacbell.net, john at johnmaddock.co.uk