linint.h

Go to the documentation of this file.

#ifndef LININT_H
#define LININT_H
#include "../ascot5.h"
#include "../offload.h"
#include "../spline/interp.h"

typedef struct {
    int n_x;     
    int bc_x;    
    real x_min;  
    real x_max;  
    real x_grid; 
    real* c;     
} linint1D_data;

typedef struct {
    int n_x;     
    int n_y;     
    int bc_x;    
    int bc_y;    
    real x_min;  
    real x_max;  
    real x_grid; 
    real y_min;  
    real y_max;  
    real y_grid; 
    real* c;     
} linint2D_data;

typedef struct {
    int n_x;     
    int n_y;     
    int n_z;     
    int bc_x;    
    int bc_y;    
    int bc_z;    
    real x_min;  
    real x_max;  
    real x_grid; 
    real y_min;  
    real y_max;  
    real y_grid; 
    real z_min;  
    real z_max;  
    real z_grid; 
    real* c;     
} linint3D_data;
 
void linint1D_init(linint1D_data* str, real* c,
                   int n_x, int bc_x,
                   real x_min, real x_max);
 
void linint2D_init(linint2D_data* str, real* c,
                   int n_x, int n_y, int bc_x, int bc_y,
                   real x_min, real x_max,
                   real y_min, real y_max);
 
void linint3D_init(linint3D_data* str, real* c,
                   int n_x, int n_y, int n_z,
                   int bc_x, int bc_y, int bc_z,
                   real x_min, real x_max,
                   real y_min, real y_max,
                   real z_min, real z_max);
 
GPU_DECLARE_TARGET_SIMD_UNIFORM(str)
int linint1D_eval_f(real* f, linint1D_data* str, real x);
DECLARE_TARGET_END
GPU_DECLARE_TARGET_SIMD_UNIFORM(str)
int linint2D_eval_f(real* f, linint2D_data* str, real x, real y);
DECLARE_TARGET_END
GPU_DECLARE_TARGET_SIMD_UNIFORM(str)
int linint3D_eval_f(real* f, linint3D_data* str,
                    real x, real y, real z);
DECLARE_TARGET_END
#endif