[vsipl++] [patch] matvec: dot, trans, kron

[Don McCoy]

Stefan Seefeld wrote:

> Jules Bergmann wrote:
>
>> Uh, the return type for herm is a bit more complex...  Maybe Stefan 
>> can suggest a type to use.  If not, go ahead and keep the current 
>> function.
>
>
> what about
>
> template <typename T, typename Block>
> typename Unary_func_view<conj_functor,
>   typename const_Matrix<complex<T>,
>                         Block>::transpose_type>::result_type
> herm(const_Matrix<complex<T>, Block> m) VSIP_NOTHROW
> {
>   typedef typename const_Matrix<complex<T>, Block>::transpose_type 
> transpose_type;
>   typedef Unary_func_view<conj_functor, transpose_type> functor_type;
>   return functor_type::apply(m.transpose());
> }
>
> This assumes the conj_functor is already defined (through the macro
> machinery in fns_elementwise.hpp that defines the conj function).
>
> Regards,
>         Stefan


This worked.  Thank you Stefan.  Function trans() is updated as per your 
suggestion Jules.  Thank you also.

Retested with icc 8.0 and gcc 3.4.0.  Checked in.
_____

Changelog:

    Added support for dot, trans and kron functions in [math.matvec]
    * src/vsip/math.hpp: included impl/matvec.hpp
    * src/vsip/impl/matvec.hpp: new file
    * tests/matvec.cpp: new file

-- 

Don McCoy
CodeSourcery, LLC



-------------- next part --------------
A non-text attachment was scrubbed...
Name: mv.diff
Type: text/x-patch
Size: 7432 bytes
Desc: not available
URL: <http://sourcerytools.com/pipermail/vsipl++/attachments/20050919/8f03df59/attachment.bin>