interp_method_config_internal.h

Go to the documentation of this file.

// Copyright (c) 2026 The YAC Authors
//
// SPDX-License-Identifier: BSD-3-Clause
 
#ifndef INTERP_METHOD_CONFIG_INTERNAL_H
#define INTERP_METHOD_CONFIG_INTERNAL_H
 
#include <mpi.h>
#include "interpolation/interp_stack_config.h"
 
struct yac_interp_method_config;
struct interp_method;
 
struct yac_interp_method_config_vtable {
  void (*delete)(struct yac_interp_method_config *config);
 
  struct yac_interp_method_config *(*copy)
    (struct yac_interp_method_config const *config);
 
  int (*compare)(void const *a, void const *b);
 
  size_t (*get_pack_size)(
    struct yac_interp_method_config const *config, MPI_Comm comm);
 
  void (*pack)(
    struct yac_interp_method_config const *config,
    void *buffer, int buffer_size, int *position, MPI_Comm comm);
 
  enum yac_interpolation_list (*get_type)(void);
 
  struct yac_param *(*get_param)
    (struct yac_interp_method_config const *config);
 
  struct interp_method *(*generate)(struct yac_interp_method_config const *config);
};
 
struct yac_interp_method_config {
  struct yac_interp_method_config_vtable const * vtable;
};
 
// Comparison helper macros
#define CHECK_BASIC(FIELD) \
  if ((a->FIELD) != (b->FIELD)) \
    return ((a->FIELD) > (b->FIELD)) - ((a->FIELD) < (b->FIELD))
 
#define CHECK_INT CHECK_BASIC
#define CHECK_SIZE_T CHECK_BASIC
#define CHECK_DOUBLE CHECK_BASIC
#define CHECK_ENUM CHECK_BASIC
 
#define CHECK_ANGLE(FIELD) \
  if (fabs((a->FIELD) - (b->FIELD)) > yac_angle_tol) \
    return ((a->FIELD) > (b->FIELD)) - ((a->FIELD) < (b->FIELD))
 
#define RETURN_ON_DIFFERENT(expr) \
  do { int ret; if ((ret = (expr))) return ret; } while(0)
 
// Use memcmp for correct handling of NaNs in fixed values
#define CHECK_DOUBLE_NAN(FIELD) \
  RETURN_ON_DIFFERENT(memcmp(&(a->FIELD), &(b->FIELD), sizeof(a->FIELD)))
 
#define CHECK_STRUCT(PTR_FIELD, CMP_ROUTINE) \
  RETURN_ON_DIFFERENT(CMP_ROUTINE(a->PTR_FIELD, b->PTR_FIELD))
 
#define CHECK_NULL_POINTER(FIELD) \
  RETURN_ON_DIFFERENT((a->FIELD == NULL) - (b->FIELD == NULL))
 
#define CHECK_STRING(FIELD) \
  do { \
    CHECK_NULL_POINTER(FIELD); \
    if (a->FIELD != NULL) RETURN_ON_DIFFERENT(strcmp(a->FIELD, b->FIELD)); \
  } while(0)
 
#define DEF_INTERP_METHOD_CONFIG_TYPE(TYPE) \
  const struct yac_interp_method_config_ ## TYPE *config_ ## TYPE = \
    (const struct yac_interp_method_config_ ## TYPE *)config;
 
#define DEF_INTERP_METHOD_CONFIG_COMPARE_TYPE(TYPE) \
  const struct yac_interp_method_config_ ## TYPE *a = \
    (const struct yac_interp_method_config_ ## TYPE *)a_; \
  const struct yac_interp_method_config_ ## TYPE *b = \
    (const struct yac_interp_method_config_ ## TYPE *)b_;
 
#endif // INTERP_METHOD_CONFIG_INTERNAL_H