diff options
Diffstat (limited to 'src/flow.c')
-rw-r--r-- | src/flow.c | 73 |
1 files changed, 65 insertions, 8 deletions
@@ -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); } |