From 3f0510629e422e979b57d3f93791937912a4183a Mon Sep 17 00:00:00 2001
From: Ian Jauslin <ian@jauslin.org>
Date: Tue, 14 Jun 2022 09:26:07 +0200
Subject: Update to v1.5.

  The update to version 1.5 is rather substantial, and introduces some minor
  backward-incompatibilities:
    * The header "#!symbols" has been replaced by "#!virtual_fields"
    * Multiplying polynomials using the '*' symbol is no longer supported (or,
      rather, the symbolic capabilities of meankondo were enhanced, and the
      syntax has been changed).
    * 'meantools exp' has been removed (its functionality is now handled by
      other means)
    * 'meantoolds derive' has been replaced by 'meantools differentiate'

  * The symbolic capabilities were enhanced: polynomials can now be
    multiplied, added, exponentiated, and their logarithms can be taken
    directly in the configuration file.

  * The flow equation can now be processed after being computed using the
    various "#!postprocess_*" entries.

  * Deprecated kondo_preprocess.

  * Compute the mean using an LU decomposition if possible.

  * More detailed checks for syntax errors in configuration file.

  * Check that different '#!group' entries are indeed uncorrelated.

  * New flags in meankondo: '-p' and '-A'.

  * New tool: meantools expand.

  * Improve conversion to LaTeX using meantools-convert

  * Assign terms randomly to different threads.

  * Created vim files to implement syntax highlighting for configuration
    files.

  * Multiple bug fixes
---
 src/flow_mpfr.c | 63 +++++++++++++++++++++++++++++++++++++++++++++++++--------
 1 file changed, 55 insertions(+), 8 deletions(-)

(limited to 'src/flow_mpfr.c')

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);
 }
 
-- 
cgit v1.2.3-70-g09d2