geometry_check.h

Go to the documentation of this file.

/***************************************************************************\
* Copyright (c) 2022, Claudio Pica, Sofie Martins                           *   
* All rights reserved.                                                      * 
\***************************************************************************/
 
#ifndef GEOMETRY_CHECK_H
#define GEOMETRY_CHECK_H
 
#include "geometry.h"
#include "error.h"
 
#ifdef CHECK_SPINOR_MATCHING
 
#define _TWO_SPINORS_MATCHING(s1, s2) error((s1)->type != (s2)->type, 1, __FILE__ ": ", "Spinor geometries don't agree!");
 
#define _ARRAY_SPINOR_MATCHING(s, n) \
    for (int _i = 0; _i < n; _i++)   \
        error((s)->type != ((s) + _i)->type, 1, __FILE__ ": ", "Spinors geometries don;t agree!");
 
#define _CHECK_GEOMETRY_MATCHING(s1, s2) error((s1)->type != (s2)->type, 1, __FILE__ ": ", "Field geometries don't agree!\n");
 
#define _CHECK_GEOMETRY_EO(s1, s2)                                                    \
    do {                                                                              \
        int pass1 = (s1)->type == &glattice && (s2)->type == &glattice;               \
        int pass2 = (s1)->type == &glat_odd && (s2)->type == &glat_even;              \
        int pass3 = (s1)->type == &glat_even && (s2)->type == &glat_odd;              \
        error(!(pass1 || pass2 || pass3), 1, __FILE__ ": ",                           \
              "Incorrect combination of geometries! "                                 \
              "Need even to odd, odd to even or both defined on the full lattice\n"); \
    } while (0)
 
#define _CHECK_GEOMETRY_FULL(s1)                                                               \
    do {                                                                                       \
        error(s1->type != &glattice, 1, __func__, "Field is not defined on all the lattice!"); \
    } while (0)
 
#define _CHECK_GEOMETRY_EVEN(s1)                                                \
    do {                                                                        \
        error(s1->type != &glat_even, 1, __func__, "Spinor needs to be even!"); \
    } while (0)
 
#define _CHECK_GEOMETRY_ODD(s1)                                               \
    do {                                                                      \
        error(s1->type != &glat_odd, 1, __func__, "Spinor needs to be odd!"); \
    } while (0)
 
#else
 
#define _TWO_SPINORS_MATCHING(s1, s2)
 
#define _ARRAY_SPINOR_MATCHING(s, n)
 
#define _CHECK_GEOMETRY_MATCHING(s1, s2)
 
#define _CHECK_GEOMETRY_EO(s1, s2)
 
#define _CHECK_GEOMETRY_FULL(s1)
 
#define _CHECK_GEOMETRY_EVEN(s1)
 
#define _CHECK_GEOMETRY_ODD(s1)
 
#endif
 
#endif