[Trilinos-Users] teuchos

[Thornquist, Heidi K]

Hi David,

When I upgraded Anasazi's Thyra interface, our call to Thyra::scale went from:

Thyra::scale(alpha,&mv);

to

Thyra::scale (alpha, Teuchos::inOutArg (mv));

It looks like you are just grabbing the first column of the multivector and scaling that vector.  That could be dangerous if there is more than one column in Y_r.

Does that make sense?

Thanks,
Heidi

________________________________
From: trilinos-users-bounces at software.sandia.gov [trilinos-users-bounces at software.sandia.gov] on behalf of Day, David [dmday at sandia.gov]
Sent: Tuesday, December 11, 2012 11:22 AM
To: trilinos-users at software.sandia.gov
Subject: [Trilinos-Users] FW: teuchos



From: <Heroux>, Michael A <maherou at sandia.gov<mailto:maherou at sandia.gov>>
Date: Tuesday, December 11, 2012 11:20 AM
To: David Day <dmday at sandia.gov<mailto:dmday at sandia.gov>>
Subject: Re: teuchos

David,

I think you might have sent this to trilinos-users-bounce instead of trilinos-users.  Please try again.

Mike

From: David Day <dmday at sandia.gov<mailto:dmday at sandia.gov>>
Date: Tuesday, December 11, 2012 11:52 AM
To: "trilinos-users-bounces at software.sandia.gov<mailto:trilinos-users-bounces at software.sandia.gov>" <trilinos-users-bounces at software.sandia.gov<mailto:trilinos-users-bounces at software.sandia.gov>>
Subject: teuchos

I am trying to upgrade from an old version of Trilinos to a new version of Trilinos.


void aXDpbY(
    Real alpha,
    const Thyra::MultiVectorBase<Real> &X_r,
    const Thyra::MultiVectorBase<Real> &X_i,
    std::vector<Real> &D_r,
    std::vector<Real> &D_i,
    Real beta,
    Thyra::MultiVectorBase<Real> &Y_r,
    Thyra::MultiVectorBase<Real> &Y_i)
{
  // (x_r + i*x_i)*(d_r + i*d_i)
  // = (x_r*d_r - x_i*d_i)
  // + i*(x_r*d_i + x_i*d_r)

  // Y_r = beta*Y_r + alpha*X_r*D_r - alpha*X_i*D_i
  //Thyra::scale(beta,&Y_r);   // old Trilinos   ,  this is what used to work
  Thyra::scale(beta,Y_r.col(0).ptr());  //   This is what compiles now.
…


Is this the right way to pass a multi – vector to Thyra::scale ?

--David Day

-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20121211/ddcd93fd/attachment.html