diff options
Diffstat (limited to 'src/flow_mpfr.c')
| -rw-r--r-- | src/flow_mpfr.c | 63 |
1 files changed, 55 insertions, 8 deletions
diff --git a/src/flow_mpfr.c b/src/flow_mpfr.c index e03c130..1701b15 100644 --- a/src/flow_mpfr.c +++ b/src/flow_mpfr.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. @@ -32,14 +32,21 @@ limitations under the License. #include "coefficient.h" #include "flow.h" #include "rcc_mpfr.h" +#include "grouped_polynomial.h" // compute flow numerically -int numerical_flow_mpfr(Grouped_Polynomial flow_equation, RCC_mpfr init, Labels labels, int niter, int display_mode){ +int numerical_flow_mpfr(Grouped_Polynomial flow_equation, RCC_mpfr init, Grouped_Polynomial postprocess_flow_equation, Labels labels, int niter, int display_mode){ // running coupling contants RCC_mpfr rccs=init; int i,j; + // for printing + RCC_mpfr rcc_print; + + + // init printing rcc + init_RCC_mpfr(&rcc_print, rccs.length); if(display_mode==DISPLAY_NUMERICAL){ // print labels @@ -50,9 +57,25 @@ int numerical_flow_mpfr(Grouped_Polynomial flow_equation, RCC_mpfr init, Labels printf("\n\n"); // print initial values + RCC_mpfr_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){ + mpfr_set_ui(rcc_print.values[j], 1, MPFR_RNDN); + } + } + evaleq_mpfr(rcc_print, rccs, postprocess_flow_equation); + } printf("%5d ",0); - for(j=0;j<rccs.length;j++){ - mpfr_printf("% 14.7Re ",rccs.values[j]); + for(j=0;j<rcc_print.length;j++){ + // use constants from rcc + if(rcc_print.indices[j]<0){ + mpfr_printf("% 14.7Re ",rccs.values[j]); + } + else{ + mpfr_printf("% 14.7Re ",rcc_print.values[j]); + } } printf("\n"); } @@ -60,12 +83,29 @@ int numerical_flow_mpfr(Grouped_Polynomial flow_equation, RCC_mpfr init, Labels for(i=0;i<niter;i++){ // compute a single step step_flow_mpfr(&rccs, flow_equation); - // convert ls to alphas + + // print + RCC_mpfr_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){ + mpfr_set_ui(rcc_print.values[j], 1, MPFR_RNDN); + } + } + evaleq_mpfr(rcc_print, rccs, postprocess_flow_equation); + } if(display_mode==DISPLAY_NUMERICAL){ // print the result printf("%5d ",i+1); - for(j=0;j<rccs.length;j++){ - mpfr_printf("% 14.7Re ",rccs.values[j]); + for(j=0;j<rcc_print.length;j++){ + // use constants from rcc + if(rcc_print.indices[j]<0){ + mpfr_printf("% 14.7Re ",rccs.values[j]); + } + else{ + mpfr_printf("% 14.7Re ",rcc_print.values[j]); + } } printf("\n"); } @@ -82,9 +122,16 @@ int numerical_flow_mpfr(Grouped_Polynomial flow_equation, RCC_mpfr init, Labels } if(display_mode==DISPLAY_FINAL){ - RCC_mpfr_print(rccs); + if(postprocess_flow_equation.length>0){ + evaleq_mpfr(rcc_print, rccs, postprocess_flow_equation); + } + else{ + rcc_print=rccs; + } + RCC_mpfr_print(rcc_print); } + free_RCC_mpfr(rcc_print); return(0); } |