interp_operator_fixed.c

Go to the documentation of this file.

// Copyright (c) 2024 The YAC Authors
//
// SPDX-License-Identifier: BSD-3-Clause
 
#include <string.h>
 
#include "interpolation/operators/interp_operator_fixed.h"
#include "interp_operator.h"
#include "utils_core.h"
#include "interpolation/interpolation_utils.h"
#include "collection_selection_internal.h"
#include "interpolation/operators/interp_operator_internal.h"
 
static int yac_interp_operator_fixed_is_source(
  struct yac_interp_operator * interp);
static int yac_interp_operator_fixed_is_target(
  struct yac_interp_operator * interp);
static void yac_interp_operator_fixed_execute(
  struct yac_interp_operator * interp,
  double *** src_fields, double *** src_frac_masks, double ** tgt_field,
  double frac_mask_fallback_value, double scale_factor, double scale_summand);
static void yac_interp_operator_fixed_execute_put(
  struct yac_interp_operator * interp,
  double *** src_fields, double *** src_frac_masks,
  int is_target, double frac_mask_fallback_value,
  double scale_factor, double scale_summand);
static void yac_interp_operator_fixed_execute_get(
  struct yac_interp_operator * interp, double ** tgt_field,
  double frac_mask_fallback_value, double scale_factor, double scale_summand);
static void yac_interp_operator_fixed_execute_get_async(
  struct yac_interp_operator * interp, double ** tgt_field,
  double frac_mask_fallback_value, double scale_factor, double scale_summand);
static enum YAC_INTERP_TEST_STATUS yac_interp_operator_fixed_execute_test(
  struct yac_interp_operator * interp);
static void yac_interp_operator_fixed_execute_wait(
  struct yac_interp_operator * interp);
static struct yac_interp_operator * yac_interp_operator_fixed_copy(
  struct yac_interp_operator * interp);
static void yac_interp_operator_fixed_delete(
  struct yac_interp_operator * interp);
 
static struct yac_interp_operator_vtable const interpolation_fixed_vtable = {
  .is_source         = yac_interp_operator_fixed_is_source,
  .is_target         = yac_interp_operator_fixed_is_target,
  .execute           = yac_interp_operator_fixed_execute,
  .execute_put       = yac_interp_operator_fixed_execute_put,
  .execute_get       = yac_interp_operator_fixed_execute_get,
  .execute_get_async = yac_interp_operator_fixed_execute_get_async,
  .execute_put_test  = yac_interp_operator_fixed_execute_test,
  .execute_get_test  = yac_interp_operator_fixed_execute_test,
  .execute_wait      = yac_interp_operator_fixed_execute_wait,
  .copy              = yac_interp_operator_fixed_copy,
  .delete            = yac_interp_operator_fixed_delete,
};
 
struct yac_interp_operator_fixed {
 
  struct yac_interp_operator_vtable const * vtable;
  struct yac_collection_selection * collection_selection;
 
  double value;
  size_t * pos;
  size_t count;
 
  int ref_count;
};
 
struct yac_interp_operator * yac_interp_operator_fixed_new(
  struct yac_collection_selection const * collection_selection,
  double value, size_t count, size_t const * pos) {
 
  struct yac_interp_operator_fixed * fixed = xmalloc(1 * sizeof(*fixed));
 
  fixed->vtable = &interpolation_fixed_vtable;
  fixed->collection_selection =
    yac_collection_selection_copy(collection_selection);
  fixed->value = value;
  fixed->count = count;
  fixed->pos = COPY_DATA(pos, count);
  fixed->ref_count = 1;
  
  return (struct yac_interp_operator *)fixed;
}
 
static int yac_interp_operator_fixed_is_source(
  struct yac_interp_operator * interp) {
  UNUSED(interp);
 
  return 0;
}
 
static int yac_interp_operator_fixed_is_target(
  struct yac_interp_operator * interp) {
  return ((struct yac_interp_operator_fixed*)interp)->count > 0;
}
 
static void yac_interp_operator_fixed_execute_get_(
  struct yac_interp_operator * interp, double ** tgt_field) {
  struct yac_interp_operator_fixed * fixed =
    (struct yac_interp_operator_fixed*)interp;
 
  double const value = fixed->value;
  size_t * restrict pos = fixed->pos;
  size_t const count = fixed->count;
  size_t const collection_size =
    yac_collection_selection_get_collection_size(
      fixed->collection_selection);
  size_t const * collection_indices =
    yac_collection_selection_get_indices(
      fixed->collection_selection);
 
  for (size_t l = 0; l < collection_size; ++l) {
    size_t l_ = (collection_indices != NULL)?collection_indices[l]:l;
    for (size_t j = 0; j < count; ++j) {
      tgt_field[l_][pos[j]] = value;
    }
  }
}
 
static void yac_interp_operator_fixed_execute_get(
  struct yac_interp_operator * interp, double ** tgt_field,
  double frac_mask_fallback_value, double scale_factor, double scale_summand) {
 
  UNUSED(frac_mask_fallback_value);
  UNUSED(scale_factor);
  UNUSED(scale_summand);
 
  yac_interp_operator_fixed_execute_get_(interp, tgt_field);
}
 
static void yac_interp_operator_fixed_execute_get_async(
  struct yac_interp_operator * interp, double ** tgt_field,
  double frac_mask_fallback_value, double scale_factor, double scale_summand) {
 
  yac_interp_operator_fixed_execute_get(
    interp, tgt_field, frac_mask_fallback_value, scale_factor, scale_summand);
}
 
static void yac_interp_operator_fixed_execute(
  struct yac_interp_operator * interp,
  double *** src_fields, double *** src_frac_masks, double ** tgt_field,
  double frac_mask_fallback_value, double scale_factor, double scale_summand) {
 
  UNUSED(src_fields);
  UNUSED(src_frac_masks);
  UNUSED(frac_mask_fallback_value);
  UNUSED(scale_factor);
  UNUSED(scale_summand);
 
  yac_interp_operator_fixed_execute_get_(interp, tgt_field);
}
 
static void yac_interp_operator_fixed_execute_put(
  struct yac_interp_operator * interp,
  double *** src_fields, double *** src_frac_masks,
  int is_target, double frac_mask_fallback_value,
  double scale_factor, double scale_summand) {
 
  UNUSED(interp);
  UNUSED(src_fields);
  UNUSED(src_frac_masks);
  UNUSED(is_target);
  UNUSED(frac_mask_fallback_value);
  UNUSED(scale_factor);
  UNUSED(scale_summand);
  return;
}
 
static enum YAC_INTERP_TEST_STATUS yac_interp_operator_fixed_execute_test(
  struct yac_interp_operator * interp) {
  UNUSED(interp);
  return YAC_INTERP_COMPLETE;
}
 
static void yac_interp_operator_fixed_execute_wait(
  struct yac_interp_operator * interp) {
  UNUSED(interp);
  return;
}
 
static struct yac_interp_operator * yac_interp_operator_fixed_copy(
  struct yac_interp_operator * interp) {
  struct yac_interp_operator_fixed * fixed =
    (struct yac_interp_operator_fixed*)interp;
 
  fixed->ref_count++;
 
  return interp;
}
 
static void yac_interp_operator_fixed_delete(
  struct yac_interp_operator * interp) {
  if (interp == NULL) return;
 
  struct yac_interp_operator_fixed * fixed =
    (struct yac_interp_operator_fixed*)interp;
 
  if(--(fixed->ref_count)) return;
 
  yac_collection_selection_delete(fixed->collection_selection);
  free(fixed->pos);
  free(fixed);
}