imat5.cpp

Go to the documentation of this file.

/*
 * $Id$
 *
 * Author: David Fournier
 * Copyright (c) 2008-2012 Regents of the University of California
 */
#include "fvar.hpp"

imatrix& imatrix::operator=(const imatrix& m1)
 {
   if (allocated(*this))
   {
     if (rowmin() != m1.rowmin() || rowmax() != m1.rowmax() ||
       colmin() != m1.colmin() || colmax() != m1.colmax() )
     {
       cerr << rowmin() << endl;
       cerr << rowmax() << endl;
       cerr << m1.rowmin() << endl;
       cerr << m1.rowmax() << endl;
       cerr << colmin() << endl;
       cerr << colmax() << endl;
       cerr << m1.colmin() << endl;
       cerr << m1.colmax() << endl;
       cerr << " Incompatible array bounds in "
       "imatrix& operator=(const imatrix&)\n";
       ad_exit(21);
     }

     if (m != m1.m)            // check for condition that both matrices
     {                         // point to the same object
       for (int i=rowmin();i<=rowmax();i++)
       {
         *(m+i) = m1[i];
       }
     }
   }
   else
   {
     index_min=m1.index_min;
     index_max=m1.index_max;
     shape=m1.shape;
     if (shape)
     {
       (shape->ncopies)++;
     }
     m = m1.m;
   }
   return(*this);
 }

 imatrix& imatrix::operator= (const int n)
 {
   int mmin=indexmin();
   int mmax=indexmax();
   for (int i=mmin;i<=mmax;i++)
   {
     (*this)(i) = n;
   }
   return(*this);
 }