lpsolution.inc

Go to the documentation of this file.

#ifndef ABA_LPSOLUTION_INC
#define ABA_LPSOLUTION_INC

#include "abacus/master.h"
#include "abacus/lpsolution.h"
#include "abacus/active.h"
#include "abacus/sub.h"

  template <class BaseType, class CoType> 
  ABA_LPSOLUTION<BaseType, CoType>::ABA_LPSOLUTION(ABA_SUB *sub, 
        bool primalVariables, ABA_ACTIVE<BaseType, CoType> *active )
  :  
    master_(sub->master_), 
    nVarCon_(primalVariables?sub->nVar():sub->nCon()), 
    idSub_(sub->id_), 
    idLp_(sub->nIter_), 
    #ifdef PARALLLEL
    idProc_(0),         
    #endif
    zVal_(master_,nVarCon_) 
  { 


  if(active==0)
    active_=0;
  else   
      active_=new ABA_ACTIVE<BaseType, CoType>(*active);

  double *p;
  if(primalVariables==true)
    p=sub->xVal_;
  else
    p=sub->yVal_;

  for(int i=0;i<nVarCon_;i++)
    zVal_[i]=p[i];

}

  template <class BaseType, class CoType> 
  ABA_LPSOLUTION<BaseType, CoType>::ABA_LPSOLUTION(ABA_MASTER *master) 
  :  
    master_(master), 
    nVarCon_(0), 
    idSub_(0), 
    idLp_(0), 
    #ifdef PARALLLEL
    idProc_(0),         
    #endif
    zVal_(master_,0) 
  { 
    active_=0;
    
  }

  template <class BaseType, class CoType>
  ABA_LPSOLUTION<BaseType, CoType>::ABA_LPSOLUTION(const ABA_LPSOLUTION<BaseType, CoType> &rhs)
  :  
    master_(rhs.master_), 
    nVarCon_(rhs.nVarCon_), 
    idSub_(rhs.idSub_), 
    idLp_(rhs.idLp_), 
    #ifdef PARALLLEL
    idProc_(rhs.idProc_),       
    #endif
    zVal_(rhs.zVal_) 
  {
   if(rhs.active_)
     active_=new ABA_ACTIVE<BaseType,CoType>(*rhs.active_);
  }

  template <class BaseType, class CoType>
  ABA_LPSOLUTION<BaseType, CoType>::~ABA_LPSOLUTION()
  {
    delete active_;
  }

  template <class BaseType, class CoType>
  ostream &operator<<(ostream &os, const ABA_LPSOLUTION<BaseType, CoType> &rhs)
  {
    double machEps=rhs.master_->machineEps();
    os << rhs.idSub_ << " ";
    os << rhs.idLp_ << " ";
    os << rhs.nVarCon_ << " ";
    for(int i=0;i<rhs.nVarCon_;i++) 
      if(fabs(rhs.zVal_[i])<machEps)
        os << "0.0 ";
      else
        os << rhs.zVal_[i] << " " ;
    os << endl;
    return os;
  }

  template <class BaseType, class CoType>
  int ABA_LPSOLUTION<BaseType, CoType>::nVarCon() const
  {
    return nVarCon_;
  }

  template <class BaseType, class CoType>
  double* ABA_LPSOLUTION<BaseType, CoType>::zVal() 
  {
    return &(zVal_[0]);
  }

  template <class BaseType, class CoType>
  ABA_ACTIVE<BaseType, CoType> * ABA_LPSOLUTION<BaseType, CoType>::active() 
  {
    return active_;
  }

  template <class BaseType, class CoType>
  int ABA_LPSOLUTION<BaseType, CoType>::idSub() const
  {
    return idSub_;
  }

  template <class BaseType, class CoType>
  int ABA_LPSOLUTION<BaseType, CoType>::idLp() const
  {
    return idLp_;
  }

  #ifdef PARALLLEL
  template <class BaseType, class CoType>
  int ABA_LPSOLUTION<BaseType, CoType>::idProc() const
  {
    return idProc_;
  }
  #endif


#endif   // ABA_LPSOLUTION_INC

---