BackupPatchOverlap.py

Go to the documentation of this file.

# This file is part of the Peano project. For conditions of distribution and
# use, please see the copyright notice at www.peano-framework.org
from peano4.solversteps.ActionSet import ActionSet
 
 
 
class BackupPatchOverlap(ActionSet):
  """
  
  This class plugs into the first access to a face and copies the data over
  
  patch_overlap_in: peano4.datamodel.Patch
     Patch data associated with cell.
        
  patch_overlap_out: peano4.datamodel.Patch
     Patch data associated with faces.
  
  invoke_in_touch_first: Boolean
    If you set this one to True, then I do the copying in touchFaceFirstTime. 
    Otherwise, I do it in touchFaceLastTime
    
  guard: Boolean 
    Predicate to determine when to invoke this thing. Pass in
    true if you wanna be it called always
  
  """
  
  
  def __init__(self,patch_overlap_in,patch_overlap_out,invoke_in_touch_first,guard="",additional_includes=""):
    super(BackupPatchOverlap,self).__init__(descend_invocation_order=1,parallel=False)
    self.d = {}
    self.d[ "UNKNOWNS" ]           = str(patch_overlap_in.no_of_unknowns)
    self.d[ "DOFS_PER_AXIS" ]      = str(patch_overlap_in.dim[1])
    self.d[ "OVERLAP" ]            = str(int(patch_overlap_in.dim[0]/2))
    self.d[ "FACES_ACCESSOR_IN" ]  = "fineGridFace"  + patch_overlap_in.name
    self.d[ "FACES_ACCESSOR_OUT" ] = "fineGridFace"  + patch_overlap_out.name
    self.d[ "GUARD_PREDICATE" ]    = guard
    
    self.invoke_in_touch_first = invoke_in_touch_first
    self._additional_includes  = additional_includes
 
 
  def set_guard_predicate(self,guard):
    self.d[ "GUARD_PREDICATE" ]    = guard
    
 
  def get_constructor_body(self):
    return ""
 
  
  def get_destructor_body(self):
    return ""
 
 
  def get_body_of_getGridControlEvents(self):
    return "  return std::vector< peano4::grid::GridControlEvent >();\n" 
 
 
  def get_action_set_name(self):
    return __name__.replace(".py", "").replace(".", "_")
 
 
  def user_should_modify_template(self):
    return False
 
 
  __Template = """
    if ({GUARD_PREDICATE}) {{
      logTraceInWith3Arguments( "--backup patch overlap from {FACES_ACCESSOR_IN} into {FACES_ACCESSOR_OUT}--(...)", marker.toString(), {FACES_ACCESSOR_OUT}.value, {FACES_ACCESSOR_IN}.value );
      int counter = 0;
      dfore(k,{DOFS_PER_AXIS},0,0) {{
        for (int i=0; i<{OVERLAP}*2; i++) {{
          for (int j=0; j<{UNKNOWNS}; j++) {{
            // assertion4( {FACES_ACCESSOR_IN}.value[counter]=={FACES_ACCESSOR_IN}.value[counter], k,i,j, marker.toString() );
            {FACES_ACCESSOR_OUT}.value[counter] = {FACES_ACCESSOR_IN}.value[counter];
            counter++;
          }}
        }}
      }}
      logTraceOutWith1Argument( "--backup patch overlap from {FACES_ACCESSOR_IN} into {FACES_ACCESSOR_OUT}--(...)", {FACES_ACCESSOR_IN}.value[0] );
    }}
"""
 
 
  def get_body_of_operation(self,operation_name):
    result = "\n"
    if self.invoke_in_touch_first and operation_name==ActionSet.OPERATION_TOUCH_FACE_FIRST_TIME:
      result = self.__Template.format(**self.d)
      pass 
    if not self.invoke_in_touch_first and operation_name==ActionSet.OPERATION_TOUCH_FACE_LAST_TIME:
      result = self.__Template.format(**self.d)
      pass 
    return result
 
 
  def get_attributes(self):
    return ""
 
 
  def get_includes(self):
    return """
#include "peano4/utils/Loop.h"
""" + self._additional_includes