IndexedBareField.hpp

Go to the documentation of this file.

// -*- C++ -*-
/***************************************************************************
 *
 * The IPPL Framework
 * 
 * This program was prepared by PSI. 
 * All rights in the program are reserved by PSI.
 * Neither PSI nor the author(s)
 * makes any warranty, express or implied, or assumes any liability or
 * responsibility for the use of this software
 *
 * Visit www.amas.web.psi for more details
 *
 ***************************************************************************/
 
// -*- C++ -*-
/***************************************************************************
 *
 * The IPPL Framework
 * 
 *
 * Visit http://people.web.psi.ch/adelmann/ for more details
 *
 ***************************************************************************/
 
// include files
#include "Field/IndexedBareField.h"
#include "Utility/IpplInfo.h"
#include "Utility/PAssert.h"
 
#include "FieldLayout/FieldLayout.h"
 
#include <vector>
#include <iostream>
 
//----------------------------------------------------------------------
 
//MWERKS: moved this member template into class definition, in 
//MWERKS: IndexedBareField.h
// template<class T, unsigned D1>
// template<unsigned int D2>
// bool
// IndexedBareFieldIterator<T,D1>::plugBase(const NDIndex<D2>& i)
 
 
 
template< class T, unsigned Dim, unsigned Brackets >
void 
IndexedBareField<T,Dim,Brackets>::write(std::ostream& out)
{
  
  
  // make sure we have the right number of brackets
  PInsist(Dim == Brackets,
          "Field not fully indexed in IndexedBareField::write!!");
  NDIndex<Dim> testIndex;
  for (unsigned d=0; d<Dim; d++)
    testIndex[d] = Indexes[d];
 
  // make a BareField which will store the subset
  FieldLayout<Dim> subfl(testIndex);
  BareField<T,Dim> subset(subfl);
 
  // assign values to this subfield
  //  my_indexed_assign(subset, A, testIndex);
  assign(subset[testIndex], A[testIndex]);
 
  // finally, print out the subfield
  out << subset;
}
 
 
 
template<class T, unsigned Dim, unsigned Brackets>
void
IndexedBareField<T,Dim,Brackets>::getsingle(T& r)
{
  
  
  A.getsingle(Indexes, r);
}
 
// Return a single value.
template<class T, unsigned Dim, unsigned Brackets>
T
IndexedBareField<T,Dim,Brackets>::get()
{
  T r;
  
  
  getsingle(r);
  return r;
}
 
/***************************************************************************
 * $RCSfile: IndexedBareField.cpp,v $   $Author: adelmann $
 * $Revision: 1.1.1.1 $   $Date: 2003/01/23 07:40:26 $
 * IPPL_VERSION_ID: $Id: IndexedBareField.cpp,v 1.1.1.1 2003/01/23 07:40:26 adelmann Exp $ 
 ***************************************************************************/