Ian Jauslin
summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
Diffstat (limited to 'src/flow.c')
-rw-r--r--src/flow.c73
1 files changed, 65 insertions, 8 deletions
diff --git a/src/flow.c b/src/flow.c
index b294bb0..10c5ae0 100644
--- a/src/flow.c
+++ b/src/flow.c
@@ -1,5 +1,5 @@
/*
-Copyright 2015 Ian Jauslin
+Copyright 2015-2022 Ian Jauslin
Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
@@ -25,14 +25,20 @@ limitations under the License.
#include "array.h"
#include "coefficient.h"
#include "rcc.h"
+#include "grouped_polynomial.h"
// compute flow numerically, no exponentials
-int numerical_flow(Grouped_Polynomial flow_equation, RCC init, Labels labels, int niter, int display_mode){
+int numerical_flow(Grouped_Polynomial flow_equation, RCC init, Grouped_Polynomial postprocess_flow_equation, Labels labels, int niter, int display_mode){
// running coupling contants
RCC rccs=init;
+ // for printing
+ RCC rcc_print;
int i,j;
+ // init printing rcc
+ init_RCC(&rcc_print, rccs.length);
+
if(display_mode==DISPLAY_NUMERICAL){
// print labels
printf("%5s ","n");
@@ -42,9 +48,25 @@ int numerical_flow(Grouped_Polynomial flow_equation, RCC init, Labels labels, in
printf("\n\n");
// print initial values
+ RCC_cpy_noinit(rccs,&rcc_print);
+ if(postprocess_flow_equation.length>0){
+ // ignore constants
+ for(j=0;j<rcc_print.length;j++){
+ if(rcc_print.indices[j]<0){
+ rcc_print.values[j]=1.;
+ }
+ }
+ evaleq(rcc_print, rcc_print, postprocess_flow_equation);
+ }
printf("%5d ",0);
- for(j=0;j<rccs.length;j++){
- printf("% 14.7Le ",rccs.values[j]);
+ for(j=0;j<rcc_print.length;j++){
+ // use constants from rcc
+ if(rcc_print.indices[j]<0){
+ printf("% 14.7Le ",rccs.values[j]);
+ }
+ else{
+ printf("% 14.7Le ",rcc_print.values[j]);
+ }
}
printf("\n");
}
@@ -52,12 +74,32 @@ int numerical_flow(Grouped_Polynomial flow_equation, RCC init, Labels labels, in
for(i=0;i<niter;i++){
// compute a single step
step_flow(&rccs, flow_equation);
- // convert ls to alphas
+
+ // print
+ if(postprocess_flow_equation.length>0){
+ RCC_cpy_noinit(rccs,&rcc_print);
+ // ignore constants
+ for(j=0;j<rcc_print.length;j++){
+ if(rcc_print.indices[j]<0){
+ rcc_print.values[j]=1.;
+ }
+ }
+ evaleq(rcc_print, rcc_print, postprocess_flow_equation);
+ }
+ else{
+ RCC_cpy_noinit(rccs,&rcc_print);
+ }
if(display_mode==DISPLAY_NUMERICAL){
// print the result
printf("%5d ",i+1);
- for(j=0;j<rccs.length;j++){
- printf("% 14.7Le ",rccs.values[j]);
+ for(j=0;j<rcc_print.length;j++){
+ // use constants from rcc
+ if(rcc_print.indices[j]<0){
+ printf("% 14.7Le ",rccs.values[j]);
+ }
+ else{
+ printf("% 14.7Le ",rcc_print.values[j]);
+ }
}
printf("\n");
}
@@ -74,9 +116,24 @@ int numerical_flow(Grouped_Polynomial flow_equation, RCC init, Labels labels, in
}
if(display_mode==DISPLAY_FINAL){
- RCC_print(rccs);
+ if(postprocess_flow_equation.length>0){
+ RCC_cpy_noinit(rccs,&rcc_print);
+ // ignore constants
+ for(j=0;j<rcc_print.length;j++){
+ if(rcc_print.indices[j]<0){
+ rcc_print.values[j]=1.;
+ }
+ }
+ evaleq(rcc_print, rcc_print, postprocess_flow_equation);
+ }
+ else{
+ RCC_cpy_noinit(rccs,&rcc_print);
+ }
+ RCC_print(rcc_print);
}
+ free_RCC(rcc_print);
+
return(0);
}