ad_structure.hpp

Go to the documentation of this file.

#pragma once

#include "../include/datatype_structure.hpp"

namespace AD{
  void StartRecording();

  void StopRecording();

  void RegisterInput(su2double &data);

  void RegisterOutput(su2double &data);

  void ClearAdjoints();

  void ComputeAdjoint();

  void Reset();

  void ResetInput(su2double &data);

  void SetPreaccIn(const su2double &data);

  void SetPreaccIn(const su2double* data, const int size);

  void SetPreaccIn(const su2double* const *data, const int size_x, const int size_y);

  void StartPreacc();

  void SetPreaccOut(su2double &data);

  void SetPreaccOut(su2double* data, const int size);

  void SetPreaccOut(su2double** data, const int size_x, const int size_y);

  void EndPreacc();

}

/*--- Macro to begin and end sections with a passive tape ---*/

#ifdef CODI_REVERSE_TYPE
#define AD_BEGIN_PASSIVE         \
  if(AD::globalTape.isActive()) {\
     AD::globalTape.setPassive();\
     AD::Status = true;          \
  }
#define AD_END_PASSIVE           \
  if(AD::Status) {               \
     AD::globalTape.setActive(); \
     AD::Status = false;         \
  }
#else
#define AD_BEGIN_PASSIVE
#define AD_END_PASSIVE
#endif

/*--- If we compile under OSX we have to overload some of the operators for
 *   complex numbers to avoid the use of the standard operators
 *  (they use a lot of functions that are only defined for doubles) ---*/

#ifdef __APPLE__

namespace std{
  template<>
  su2double abs(const complex<su2double>& x);

  template<>
  complex<su2double> operator/(const complex<su2double>& x, const complex<su2double>& y);

  template<>
  complex<su2double> operator*(const complex<su2double>& x, const complex<su2double>& y);
}
#endif
#include "ad_structure.inl"