| Re: [eigen] eigen blas/lapack + suitesparse |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] eigen blas/lapack + suitesparse
- From: Rhys Ulerich <rhys.ulerich@xxxxxxxxx>
- Date: Fri, 17 Jun 2011 14:10:48 -0500
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:in-reply-to:references:from:date :message-id:subject:to:content-type; bh=Un9a35k0ovVlGLNJsVIv44bhYtMCdov4vjlr7CFNPm0=; b=mtNRHDVUoyg2V7iAGwzV5Nz8nEZgJk9kXoxw2lo0RBVjkvLDdFeFUzd5uOIuoLbISb f76km7WvlqKR9vg0VIe6ibT2XeHfdCk7M1yRwTpoo050wCeTuJyRM66otf6NH2++60Wr tBeAZ13zLNNea/1uX7Rodz56IVKBE4XrXIL+M=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :content-type; b=OXHqFBqmTQ2Ub4jbf4V4T+Vm8gc3+xuyCP3WvhEk+ERiQQQnzZ+QBnW2MBnJgG2ZKN XAfFNQ8MMMkSeed0AkEi4cIE04MjLXdgM3bpinocbXriZAw6En2VxpLMdEYpRKWltKIh 6Cwm1Vhnb5CT55jVVPUnvFHWMR9D2Fpa8w4rY=
On Fri, Jun 17, 2011 at 12:52 PM, Gael Guennebaud
<gael.guennebaud@xxxxxxxxx> wrote:
> yes, it is complete enough for all direct solver libraries that I've
> tried.
Was the Eigen BLAS intended to handle negative strides? I've noticed
code like the following where incx < 0 looks like it behaves
incorrectly:
int EIGEN_BLAS_FUNC(scal)(int *n, RealScalar *palpha, RealScalar *px, int *incx)
{
if(*n<=0) return 0;
Scalar* x = reinterpret_cast<Scalar*>(px);
Scalar alpha = *reinterpret_cast<Scalar*>(palpha);
if(*incx==1) vector(x,*n) *= alpha;
else vector(x,*n,std::abs(*incx)) *= alpha;
return 0;
}
Has anyone run into such a problem in practice?
- Rhys