wilsonflow.h

/***************************************************************************\
* Copyright (c) 2011, Agostino Patella                                      *   
* All rights reserved.                                                      * 
\***************************************************************************/
 
/*******************************************************************************
*
* File wilsonflow.h
* 
*******************************************************************************/
 
#ifndef WILSONFLOW_H
#define WILSONFLOW_H
 
#include "suN_types.h"
#include "hr_complex.h"
#include "spinor_field.h"
#include "suN.h"
#include "data_storage.h"
 
#ifdef __cplusplus
extern "C" {
#endif
 
typedef enum { EUL = 0, RK3 = 1, RK3_ADAPTIVE = 2 } WF_integrator_type;
 
void WF_initialize();
void WF_free();
 
void WF_set_bare_anisotropy(double *chi);
 
double max_distance(suNg_field *V, suNg_field *Vprime);
 
void WilsonFlow1(suNg_field *V, const double epsilon);
void WilsonFlow3(suNg_field *V, const double epsilon);
int WilsonFlow3_adaptative(suNg_field *V, double *epsilon, double *epsilon_new, double *delta);
 
data_storage_array *WF_update_and_measure(WF_integrator_type wft, suNg_field *V, double *tmax, double *eps, double *delta,
                                          int nmeas, storage_switch swc);
 
#ifdef __cplusplus
}
#endif
#endif /* WILSONFLOW_H */