Fixup1.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2002 Ralf Deiterding
// Brandenburgische Universitaet Cottbus

#ifndef AMROC_FIXUP1_H
#define AMROC_FIXUP1_H

#include "FixupBase.h"

template <class VectorType, class FixupType>
class FixupOps<VectorType,FixupType,1> {
public:
  FixupOps() {}

  void AlignBBox(BBox &bb, const int s) {
    int DAGH_Base = DAGH_X;
#ifdef DEBUG_PRINT
    assert (s>=0 && s<2*1);
#endif
    int dim_1 = 1;
    if (s % 2 != 0) 
      DAGH_Base = DAGH_Base | DAGH_Dim1ToTop;
    gdbAlignBBox(dim_1, bb, DAGH_Base - s/2);
    bb.rank = dim_1;
    bb.lower().rank = dim_1;
    bb.upper().rank = dim_1;
    bb.stepsize().rank = dim_1;
  } 

  void copy_to(GridData<FixupType,1> &target, const BBox &to,
               const GridData<VectorType,1> &source, const BBox &fromwhere, const int s) {   
    BBox tobb(to * target.bbox()); 
    BBox frombb(fromwhere * source.bbox()); AlignBBox(frombb, s);
    equals_to(target(tobb.lower(0)),source(frombb.lower(0)));
  }
  void copy_to(GridData<FixupType,1> &target,  
               const GridData<VectorType,1> &source, const BBox &fromwhere, const int s) {   
    BBox tobb(target.bbox()); 
    BBox frombb(fromwhere * source.bbox()); AlignBBox(frombb, s);
    equals_to(target(tobb.lower(0)),source(frombb.lower(0)));
  }
  
  void copy_from(GridData<VectorType,1> &target, const BBox &towhere, 
                 const GridData<FixupType,1> &source, const BBox &from, const int s) {    
    BBox tobb(towhere * target.bbox()); AlignBBox(tobb, s);
    BBox frombb(from * source.bbox()); 
    equals_from(target(tobb.lower(0)),source(frombb.lower(0)));
  }
  void copy_from(GridData<VectorType,1> &target, const BBox &towhere, 
                 const GridData<FixupType,1> &source, const int s) {    
    BBox tobb(towhere * target.bbox()); AlignBBox(tobb, s);
    BBox frombb(source.bbox()); 
    equals_from(target(tobb.lower(0)),source(frombb.lower(0)));
  }
  
  void add_to(GridData<FixupType,1> &target, const BBox &to, 
              const GridData<VectorType,1> &source, const BBox &fromwhere, const int s) {    
    BBox tobb(to * target.bbox()); 
    BBox frombb(fromwhere * source.bbox()); AlignBBox(frombb, s);
    plus_to(target(tobb.lower(0)),source(frombb.lower(0)));
  }
  void add_to(GridData<FixupType,1> &target,
              const GridData<VectorType,1> &source, const BBox &fromwhere, const int s) {    
    BBox tobb(target.bbox()); 
    BBox frombb(fromwhere * source.bbox()); AlignBBox(frombb, s);
    plus_to(target(tobb.lower(0)),source(frombb.lower(0)));
  }
  
  void add_from(GridData<VectorType,1> &target, const BBox &towhere, 
                const GridData<FixupType,1> &source, const BBox &from, const int s) {    
    BBox tobb(towhere * target.bbox()); AlignBBox(tobb, s);
    BBox frombb(from * source.bbox()); 
    plus_from(target(tobb.lower(0)),source(frombb.lower(0)));
  }
  void add_from(GridData<VectorType,1> &target, const BBox &towhere, 
                const GridData<FixupType,1> &source, const int s) {    
    BBox tobb(towhere * target.bbox()); AlignBBox(tobb, s);
    BBox frombb(source.bbox()); 
    plus_from(target(tobb.lower(0)),source(frombb.lower(0)));
  }
  
  void add_to(GridData<FixupType,1> &target, const BBox &to, 
              const GridData<VectorType,1> &source, const BBox &from) {   
    BBox tobb(to * target.bbox()); 
    BBox frombb(from * source.bbox()); 
    plus_to(target(tobb.lower(0)),source(frombb.lower(0)));
  }
  void add_to(GridData<FixupType,1> &target, 
              const GridData<VectorType,1> &source, const BBox &where) {   
    BBox tobb(target.bbox()); 
    BBox frombb(where * source.bbox());
    plus_to(target(tobb.lower(0)),source(frombb.lower(0)));
  }
};

template <class VectorType, class FixupType>
class Fixup<VectorType,FixupType,1> : public FixupBase<VectorType,FixupType,1> {
  typedef GridData<VectorType,1>  vec_grid_data_type;
  typedef GridData<VectorType,minus_1<1>::dim> ld_vec_grid_data_type;
  typedef GridData<FixupType,minus_1<1>::dim>  ld_fixup_grid_data_type;

public:
  Fixup() : FixupBase<VectorType,FixupType,1>() {}  

protected:
  void debug_print(vec_grid_data_type &gd, const BBox &where) {
#ifdef DEBUG_PRINT
    int lower = std::max(where.lower(0),gd.lower(0));
    int upper = std::min(where.upper(0),gd.upper(0));
    for (register int i=lower; i<=upper; i+=where.stepsize(0)) 
      ( comm_service::log() << "[" << i << "]=" 
        << gd(i)(0) << " " ).flush();
#endif
  }    

  void debug_print_ldv(ld_vec_grid_data_type &gd, const BBox &where) {
#ifdef DEBUG_PRINT
    ( comm_service::log() << "[" << gd.bbox().lower(0) << "]=" 
      << gd(gd.bbox().lower(0))(0) << " " ).flush();
    ( comm_service::log() << "\n" ).flush();
#endif
  }    
  void debug_print_ldv(ld_vec_grid_data_type &gd) {
    debug_print_ldv(gd, gd.bbox());
  }

  void debug_print_ld(ld_fixup_grid_data_type &gd, const BBox &where) {
#ifdef DEBUG_PRINT
    ( comm_service::log() << "[" << gd.bbox().lower(0) << "]=" 
      << gd(gd.bbox().lower(0))(0) << " " ).flush();
    ( comm_service::log() << "\n" ).flush();
#endif
  }    
  void debug_print_ld(ld_fixup_grid_data_type &gd) {
    debug_print_ld(gd, gd.bbox());
  }
};


#endif

---