linear_spectral_fusion.h

Go to the documentation of this file.

 
#ifndef OAKFIELD_LINEAR_SPECTRAL_FUSION_H
#define OAKFIELD_LINEAR_SPECTRAL_FUSION_H
 
#include "oakfield/operator_split.h"
 
#ifdef __cplusplus
extern "C" {
#endif
 
struct SimContext;
 
typedef struct LinearSpectralFusionOperatorConfig {
    size_t field_index; 
    /* Dissipation: exp(dt * lambda(k)), lambda(k) = -viscosity * |k|^alpha. */
    double viscosity;           
    double alpha;               
    double dissipation_spacing; 
    /* Dispersion: exp(i * dt * omega(k)), omega(k) = coefficient * | |k| - k0 |^order.
     * Real fields use the projected real multiplier cos(dt * omega(k)).
     */
    double dispersion_coefficient; 
    double dispersion_order;       
    double dispersion_reference_k; 
    double dispersion_spacing;     
    /* Global phase: exp(i * dt * phase_rate). Real fields use cos(dt * phase_rate). */
    double phase_rate; 
} LinearSpectralFusionOperatorConfig;
 
SimResult sim_add_linear_spectral_fusion_operator(struct SimContext *context,
                                                  const LinearSpectralFusionOperatorConfig *config,
                                                  size_t *out_index);
 
SimResult sim_linear_spectral_fusion_config(struct SimContext *context, size_t operator_index,
                                            LinearSpectralFusionOperatorConfig *out_config);
 
SimResult sim_linear_spectral_fusion_update(struct SimContext *context, size_t operator_index,
                                            const LinearSpectralFusionOperatorConfig *config);
 
#ifdef __cplusplus
}
#endif
 
#endif /* OAKFIELD_LINEAR_SPECTRAL_FUSION_H */