File complexscalarkernels.h

File List > code_source > cosmolattice > include > CosmoInterface > evolvers > kernels > complexscalarkernels.h

Go to the documentation of this file

#ifndef COSMOINTERFACE_EVOLVERS_KERNELS_COMPLEXSCALARKERNELS_H
#define COSMOINTERFACE_EVOLVERS_KERNELS_COMPLEXSCALARKERNELS_H

/* This file is part of CosmoLattice, available at www.cosmolattice.net .
   Copyright Daniel G. Figueroa, Adrien Florio, Francisco Torrenti and Wessel Valkenburg.
   Released under the MIT license, see LICENSE.md. */

// File info: Main contributor(s): Daniel G. Figueroa, Adrien Florio, Francisco Torrenti,  Year: 2020

#include "CosmoInterface/definitions/potential.h"
#include "CosmoInterface/definitions/gaugederivatives.h"
#include "CosmoInterface/evolvers/kernels/kernelstypes.h"

namespace TempLat
{
  // @label:complexscalarkernels_class
  class ComplexScalarKernels
  {
  public:
    // Put public methods here. These should change very little over time.
    ComplexScalarKernels() = delete;

    // Equations of motion:
    template <class Model, int N> static auto get(Model &model, Tag<N> n, KernelsTypes::EoM<Model> eom)
    {
      // Returns kernel for complex scalars (formed by the covariant laplacian and potential derivative terms):

      // return ZeroType();
      return pow(model.aI, 1 + model.alpha) * GaugeDerivatives::covLaplacianCS(model, n) -
             pow(model.aI, 3 + model.alpha) / 2 * Potential::derivCS(model, n);
    }

    template <class Model, int N> static auto get_momentum(Model &model, Tag<N> n, KernelsTypes::EoM<Model> eom)
    {
      // Returns momentum for complex scalars:
      return pow(model.aI, model.alpha - 3) * model.piCS(n);
    }

    template <class Model, int N> static auto get(Model &model, Tag<N> n, KernelsTypes::Diffusion<Model> diffusion)
    {
      // Returns momentum for complex scalars:
      return Complexify(ZeroType(), ZeroType());
    }

    template <class Model, int N> static auto get_momentum(Model &model, Tag<N> n, KernelsTypes::Diffusion<Model> diffusion)
    {
      // Returns momentum for complex scalars:
      return ComplexScalarKernels::get(model, n, KernelsTypes::EoM<Model>());
    }

    // Default function returns EoM kernels, for backward compatibility.
    template <class Model, int N> static auto get(Model &model, Tag<N> n)
    {
      return ComplexScalarKernels::get(model, n, KernelsTypes::EoM<Model>());
    }
  };
  // @endlabel
} // namespace TempLat

#endif