[Trilinos-Users] FW: teuchos

Day, David dmday at sandia.gov
Tue Dec 11 11:22:13 MST 2012



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/b2ca709c/attachment-0001.html 


More information about the Trilinos-Users mailing list