commit 2653c1f473a47f6af6e5e5547a2ab7a3e34392fb
parent 7020270c63400828533363a2ac2824c8e1237b3e
Author: Christophe Coustet <christophe.coustet@meso-star.com>
Date: Thu, 30 Apr 2020 13:20:49 +0200
Remove obsolete post-process
Diffstat:
| D | pp/green2xslx.pl | | | 580 | ------------------------------------------------------------------------------- |
1 file changed, 0 insertions(+), 580 deletions(-)
diff --git a/pp/green2xslx.pl b/pp/green2xslx.pl
@@ -1,580 +0,0 @@
-#!/usr/bin/perl
-# Copyright (C) 2018-2020 |Meso|Star>
-#
-# This program is free software: you can redistribute it and/or modify
-# it under the terms of the GNU General Public License as published by
-# the Free Software Foundation, either version 3 of the License, or
-# (at your option) any later version.
-#
-# This program is distributed in the hope that it will be useful,
-# but WITHOUT ANY WARRANTY; without even the implied warranty of
-# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
-# GNU General Public License for more details.
-#
-# You should have received a copy of the GNU General Public License
-# along with this program. If not, see <http://www.gnu.org/licenses/>.
-
-use strict;
-use warnings;
-
-use Excel::Writer::XLSX;
-use Data::Dumper;
-
-#
-# Read and parse STDIN
-#
-
-# Drop warnings and text until leading '---BEGIN GREEN---'
-my $found_start = 0;
-my $line;
-while ($line = <STDIN>) {
- chomp $line;
- if ($line eq "---BEGIN GREEN---") { $found_start = 1; last;}
-}
-die 'No green found!' unless $found_start;
-
-# Read counters
-my $solids_rk=0;
-my $fluids_rk=1;
-my $tbounds_rk=2;
-my $hbounds_rk=3;
-my $fbounds_rk=4;
-my $ok_rk=5;
-my $failed_rk=6;
-$line = <STDIN>;
-chomp $line;
-die "Unexpected content found ($line)!" unless $line eq "# #solids #fluids #t_boundaries #h_boundaries #f_boundaries #ok #failures";
-$line = <STDIN>;
-chomp $line;
-my @counts = $line =~ /(\d+)/g;
-
-# Need at least 1 successful sample to proceed!
-die "No successful samples in this green ($counts[$failed_rk] ) failed samples\n" unless $counts[$ok_rk];
-
-my $last_id = -1;
-my @seen_id_types;
-
-# Read Solids
-my @solids;
-if ($counts[$solids_rk] > 0) {
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# Solids");
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# ID Name lambda rho cp power");
- for (my $s = 0; $s < $counts[$solids_rk] ; $s++) {
- $line = <STDIN>;
- chomp $line;
- my @tmp = split("\t", $line);
- die "Wrong number of values!" unless (scalar(@tmp) == 6);
- die "Wrong ID" unless ($tmp[0] >= 0) && (! defined $seen_id_types[$tmp[0]]);
- $seen_id_types[$tmp[0]] = 'S';
- my %new_elt = (
- ID=>$tmp[0],
- TEMP=>-1, # Imposed temperature for solids not yet in Stardis-app
- NAME=>$tmp[1],
- LAMBDA=>$tmp[2],
- RHO=>$tmp[3],
- CP=>$tmp[4],
- POWER=>$tmp[5]
- );
- push @solids, \%new_elt;
- }
-}
-
-# Read Fluids
-my @fluids;
-if ($counts[$fluids_rk] > 0) {
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# Fluids");
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# ID Name rho cp");
- for (my $f = 0; $f < $counts[$fluids_rk] ; $f++) {
- $line = <STDIN>;
- chomp $line;
- my @tmp = split("\t", $line);
- die "Wrong number of values!" unless (scalar(@tmp) == 4);
- die "Wrong ID" unless ($tmp[0] >= 0) && (! defined $seen_id_types[$tmp[0]]);
- $seen_id_types[$tmp[0]] = 'F';
- my %new_elt = (
- ID=>$tmp[0],
- TEMP=>-1, # Imposed temperature for fluids not yet in Stardis-app
- NAME=>$tmp[1],
- RHO=>$tmp[2],
- CP=>$tmp[3],
- POWER=>0 # Volumic Power for fluids not yet in Stardis-app
- );
- push @fluids, \%new_elt;
- }
-}
-
-# Read T Boundaries
-my @t_boundaries;
-if ($counts[$tbounds_rk] > 0) {
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# T Boundaries");
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# ID Name temperature");
- for (my $b = 0; $b < $counts[$tbounds_rk] ; $b++) {
- $line = <STDIN>;
- chomp $line;
- my @tmp = split("\t", $line);
- die "Wrong number of values!" unless (scalar(@tmp) == 3);
- die "Wrong ID" unless ($tmp[0] >= 0) && (! defined $seen_id_types[$tmp[0]]);
- $seen_id_types[$tmp[0]] = 'T';
- my %new_elt = (
- ID=>$tmp[0],
- NAME=>$tmp[1],
- TEMP=>$tmp[2]
- );
- push @t_boundaries, \%new_elt;
- }
-}
-
-# Read H Boundaries
-my @h_boundaries;
-if ($counts[$hbounds_rk] > 0) {
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# H Boundaries");
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# ID Name emissivity specular_fraction hc hc_max T_env");
- for (my $b = 0; $b < $counts[$hbounds_rk] ; $b++) {
- $line = <STDIN>;
- chomp $line;
- my @tmp = split("\t", $line);
- die "Wrong number of values!" unless (scalar(@tmp) == 7);
- die "Wrong ID" unless ($tmp[0] >= 0) && (! defined $seen_id_types[$tmp[0]]);
- $seen_id_types[$tmp[0]] = 'H';
- my %new_elt = (
- ID=>$tmp[0],
- NAME=>$tmp[1],
- EMISSIVITY=>$tmp[2],
- SPEC_FRACTION=>$tmp[3],
- HC=>$tmp[4],
- HC_MAX=>$tmp[5],
- T_ENV=>$tmp[6]
- );
- push @h_boundaries, \%new_elt;
- }
-}
-
-# Read F Boundaries
-my @f_boundaries;
-if ($counts[$fbounds_rk] > 0) {
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# F Boundaries");
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# ID Name flux");
- for (my $b = 0; $b < $counts[$fbounds_rk] ; $b++) {
- $line = <STDIN>;
- chomp $line;
- my @tmp = split("\t", $line);
- die "Wrong number of values!" unless (scalar(@tmp) == 3);
- die "Wrong ID" unless ($tmp[0] >= 0) && (! defined $seen_id_types[$tmp[0]]);
- $seen_id_types[$tmp[0]] = 'X';
- my %new_elt = (
- ID=>$tmp[0],
- NAME=>$tmp[1],
- FLUX=>$tmp[2]
- );
- push @f_boundaries, \%new_elt;
- }
-}
-
-# Read Radiative Temperatures
-my $radiative_temp;
-my $linear_temp;
-my $rad_temp_id;
-{
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# Radiative Temperatures");
- $line = <STDIN>;
- chomp $line;
- die "Unexpected content found ($line)!" unless ($line eq "# ID Rad_Temp Lin_Temp");
- $line = <STDIN>;
- chomp $line;
- my @tmp = split("\t", $line);
- die "Wrong number of values!" unless (scalar(@tmp) == 3);
- die "Wrong ID" unless ($tmp[0] >= 0) && (! defined $seen_id_types[$tmp[0]]);
- $seen_id_types[$tmp[0]] = 'R';
- die "Wrong Temperature!" unless ($tmp[1] >= 0) && ($tmp[2] >= 0);
- $radiative_temp = $tmp[1];
- $linear_temp = $tmp[2];
- $rad_temp_id = $tmp[0];
-}
-
-# Read Samples headers
-$line = <STDIN>;
-chomp $line;
-die "Unexpected content found ($line)!" unless ($line eq "# Samples");
-$line = <STDIN>;
-chomp $line;
-die "Unexpected content found ($line)!" unless ($line eq "# end #power_terms #flux_terms power_term_1 ... power_term_n flux_term_1 ... flux_term_n");
-$line = <STDIN>;
-chomp $line;
-die "Unexpected content found ($line)!" unless ($line eq "# end = end_type end_id; end_type = T | H | X | R | F | S");
-$line = <STDIN>;
-chomp $line;
-die "Unexpected content found ($line)!" unless ($line eq "# power_term_i = power_type_i power_id_i factor_i");
-$line = <STDIN>;
-chomp $line;
-die "Unexpected content found ($line)!" unless ($line eq "# flux_term_i = flux_id_i factor_i");
-
-# read samples
-my @samples;
-for(my $s=0; $s < $counts[$ok_rk]; $s++) {
- $line = <STDIN>;
- die "Unexpected end of data ($s samples read)" unless $line;
- chomp $line;
- my @tmp = split("\t", $line);
-
- # Check read data
- die "Wrong number of values!" unless (scalar(@tmp) >= 4);
- my $pw_count = $tmp[2];
- my $fx_count = $tmp[3];
- die "Wrong power_terms count!" unless ($pw_count >=0);
- die "Wrong flux_terms count!" unless ($fx_count >= 0);
- die "Wrong number of terms!" unless (scalar(@tmp) == 4 + $pw_count * 3 + $fx_count * 2);
-
- my @pw_types;
- my @pw_ids;
- my @pw_factors;
- my @fx_ids;
- my @fx_factors;
- for(my $n = 0; $n < $pw_count; $n++) {
- my $ty = $tmp[4+3*$n] ;
- my $id = $tmp[5+3*$n];
- my $fc =$tmp[6+3*$n];
- die "Wrong ID" unless ($id >= 0) && (defined $seen_id_types[$id]) && ($seen_id_types[$id] eq $ty);
- push @pw_types, $ty;
- push @pw_ids, $id;
- push @pw_factors, $fc;
- }
- for(my $n = 0; $n < $fx_count; $n++) {
- my $id = $tmp[4+$pw_count*3+2*$n];
- my $fc =$tmp[5+$pw_count*3+2*$n];
- die "Wrong ID" unless ($id >= 0) && (defined $seen_id_types[$id]) && ($seen_id_types[$id] eq 'X');
- push @fx_ids, $id;
- push @fx_factors, $fc;
- }
- my %new_elt = (
- END_TYPE=>$tmp[0],
- END_ID=>$tmp[1],
- PW_COUNT=>$tmp[2],
- FX_COUNT=>$tmp[3],
- PW_TYPES=>\@pw_types,
- PW_IDS=>\@pw_ids,
- PW_FACTORS=>\@pw_factors,
- FX_IDS=>\@fx_ids,
- FX_FACTORS=>\@fx_factors
- );
- push @samples, \%new_elt;
-}
-
-# Check end of data
-$line = <STDIN>;
-chomp $line;
-die "Unexpected content in file" unless ($line eq "---END GREEN---");
-
-#
-# Export data into an xlsx file
-#
-
-binmode( STDOUT );
-my $workbook = Excel::Writer::XLSX->new( \*STDOUT );
-$workbook->set_properties(
- title => 'Green calculator',
- author => 'meso-star.com',
- comments => 'Created with Perl, Excel::Writer::XLSX and Stardis-app Post-Process',
- );
-
-my $title = $workbook->add_format();
-$title->set_locked(1);
-$title->set_bold();
-$title->set_top(5);
-$title->set_bottom(2);
-my $locked = $workbook->add_format();
-$locked->set_locked(1);
-$locked->set_bg_color('#F2F2F2'); #light gray
-my $unlocked = $workbook->add_format();
-$unlocked->set_locked(0);
-
-my $model_current_line = 0;
-my @id_end_cells;
-my @id_pw_cells;
-my @id_fx_cells;
-my @cpt_end_cells;
-my @cpt_pw_cells;
-my @cpt_fx_cells;
-my @name_end_cells;
-my @name_pw_cells;
-my @name_fx_cells;
-
-# One sheet for all but samples
-my $model = $workbook->add_worksheet('Model');
-$model->protect(); # Cannot edit cells unless unlocked
-$model->set_column(0, 5, 24); # Column 0 to 5: width = 24
-
-# Create table for solids
-if($counts[$solids_rk] > 0) {
- my @solid_colnames = ('Solid Name', 'Imposed Temperature', 'Lambda', 'Rho', 'Cp', 'Volumic Power');
- $model->write($model_current_line, 0, \@solid_colnames, $title);
- $model_current_line++;
- for(my $s = 0; $s < $counts[$solids_rk] ; $s++) {
- my $solid = $solids[$s];
- my $temp = $solid->{TEMP};
- my $pw = $solid->{POWER};
- my $name = $solid->{NAME};
- $model->write($model_current_line, 0, $name, $locked);
- $model->write_number($model_current_line, 1, $temp, ($temp >= 0) ? $unlocked : $locked);
- $model->write_number($model_current_line, 2, $solid->{LAMBDA}, $locked);
- $model->write_number($model_current_line, 3, $solid->{RHO}, $locked);
- $model->write_number($model_current_line, 4, $solid->{CP}, $locked);
- $model->write_number($model_current_line, 5, $pw, ($pw > 0) ? $unlocked : $locked);
- $model_current_line++;
- my $id = $solid->{ID};
- die "Inconsistency!" unless ($id >= 0) && (defined $seen_id_types[$id]) && ($seen_id_types[$id] eq 'S');
- die "Inconsistency!" unless (!defined $id_end_cells[$id]) && (!defined $cpt_end_cells[$id]) && (!defined $name_end_cells[$id]);
- die "Inconsistency!" unless (!defined $id_pw_cells[$id]) && (!defined $cpt_pw_cells[$id]) && (!defined $name_pw_cells[$id]);
- if ($temp >= 0) {
- $id_end_cells[$id] = "B$model_current_line";
- $cpt_end_cells[$id] = 0;
- $name_end_cells[$id] = "$name.Imposed_Temperature";
- }
- $id_pw_cells[$id] = "F$model_current_line";
- $cpt_pw_cells[$id] = 0;
- $name_pw_cells[$id] = "$name.Volumic_Power";
- }
- $model_current_line++;
-}
-
-# Create table for fluids
-if($counts[$fluids_rk] > 0) {
- my @fluid_colnames = ('Fluid Name', 'Imposed Temperature', 'Rho', 'Cp', 'Volumic Power');
- $model->write($model_current_line, 0, \@fluid_colnames, $title);
- $model_current_line++;
- for(my $f = 0; $f < $counts[$fluids_rk] ; $f++) {
- my $fluid = $fluids[$f];
- my $temp = $fluid->{TEMP};
- my $pw = $fluid->{POWER};
- my $name = $fluid->{NAME};
- $model->write($model_current_line, 0, $name, $locked);
- $model->write_number($model_current_line, 1, $temp, ($temp >= 0) ? $unlocked : $locked);
- $model->write_number($model_current_line, 2, $fluid->{RHO}, $locked);
- $model->write_number($model_current_line, 3, $fluid->{CP}, $locked);
- $model->write_number($model_current_line, 4, $pw, ($pw > 0) ? $unlocked : $locked);
- $model_current_line++;
- my $id = $fluid->{ID};
- die "Inconsistency!" unless ($id >= 0) && (defined $seen_id_types[$id]) && ($seen_id_types[$id] eq 'F');
- die "Inconsistency!" unless (!defined $id_end_cells[$id]) && (!defined $cpt_end_cells[$id]) && (!defined $name_end_cells[$id]);
- die "Inconsistency!" unless (!defined $id_pw_cells[$id]) && (!defined $cpt_pw_cells[$id]) && (!defined $name_pw_cells[$id]);
- if ($temp >= 0) {
- $id_end_cells[$id] = "B$model_current_line";
- $cpt_end_cells[$id] = 0;
- $name_end_cells[$id] = "$name.Imposed_Temperature";
- }
- $id_pw_cells[$id] = "E$model_current_line";
- $cpt_pw_cells[$id] = 0;
- $name_pw_cells[$id] = "$name.Volumic_Power";
- }
- $model_current_line++;
-}
-
-# Create table for T boundaries
-if($counts[$tbounds_rk] > 0) {
- my @tbound_colnames = ('T Boundary Name', 'Temperature');
- $model->write($model_current_line, 0, \@tbound_colnames, $title);
- $model_current_line++;
- for(my $b = 0; $b < $counts[$tbounds_rk] ; $b++) {
- my $tbound = $t_boundaries[$b];
- my $temp = $tbound->{TEMP};
- my $name = $tbound->{NAME};
- $model->write($model_current_line, 0, $name, $locked);
- $model->write_number($model_current_line, 1, $temp, ($temp >= 0) ? $unlocked : $locked);
- $model_current_line++;
- my $id = $tbound->{ID};
- die "Inconsistency!" unless ($id >= 0) && (defined $seen_id_types[$id]) && ($seen_id_types[$id] eq 'T');
- die "Inconsistency!" unless (!defined $id_end_cells[$id]) && (!defined $cpt_end_cells[$id]) && (!defined $name_end_cells[$id]);
- $id_end_cells[$id] = "B$model_current_line";
- $cpt_end_cells[$id] = 0;
- $name_end_cells[$id] = "$name.Temperature";
- # No volumic power at boundaries!
- }
- $model_current_line++;
-}
-
-# Create table for H boundaries
-if($counts[$hbounds_rk] > 0) {
- my @hbound_colnames = ('H Boundary Name', 'Emissivity', 'Specular Fraction', 'Hc', 'Hc Max', 'Environment Temperature');
- $model->write($model_current_line, 0, \@hbound_colnames, $title);
- $model_current_line++;
- for(my $b = 0; $b < $counts[$hbounds_rk] ; $b++) {
- my $hbound = $h_boundaries[$b];
- my $name = $hbound->{NAME};
- $model->write($model_current_line, 0, $name, $locked);
- $model->write_number($model_current_line, 1, $hbound->{EMISSIVITY}, $locked);
- $model->write_number($model_current_line, 2, $hbound->{SPEC_FRACTION}, $locked);
- $model->write_number($model_current_line, 3, $hbound->{HC}, $locked);
- $model->write_number($model_current_line, 4, $hbound->{HC_MAX}, $locked);
- $model->write_number($model_current_line, 5, $hbound->{T_ENV}, $unlocked);
- $model_current_line++;
- my $id = $hbound->{ID};
- die "Inconsistency!" unless ($id >= 0) && (defined $seen_id_types[$id]) && ($seen_id_types[$id] eq 'H');
- die "Inconsistency!" unless (!defined $id_end_cells[$id]) && (!defined $cpt_end_cells[$id]) && (!defined $name_end_cells[$id]);
- $id_end_cells[$id] = "F$model_current_line";
- $cpt_end_cells[$id] = 0;
- $name_end_cells[$id] = "$name.Environment_Temperature";
- # No volumic power at boundaries!
- }
- $model_current_line++;
-}
-
-# Create table for F boundaries
-if($counts[$fbounds_rk] > 0) {
- my @fbound_colnames = ('F Boundary Name', 'Flux');
- $model->write($model_current_line, 0, \@fbound_colnames, $title);
- $model_current_line++;
- for(my $b = 0; $b < $counts[$fbounds_rk] ; $b++) {
- my $fbound = $f_boundaries[$b];
- my $name = $fbound->{NAME};
- $model->write($model_current_line, 0, $name, $locked);
- $model->write_number($model_current_line, 1, $fbound->{FLUX}, $unlocked);
- $model_current_line++;
- my $id = $fbound->{ID};
- die "Inconsistency!" unless ($id >= 0) && (defined $seen_id_types[$id]) && ($seen_id_types[$id] eq 'X');
- die "Inconsistency!" unless (!defined $id_fx_cells[$id]) && (!defined $cpt_fx_cells[$id]) && (!defined $name_fx_cells[$id]);
- $id_fx_cells[$id] = "B$model_current_line";
- $cpt_fx_cells[$id] = 0;
- $name_fx_cells[$id] = "$name.Flux";
- # Cannot end at F boundaries
- # No volumic power at boundaries!
- }
- $model_current_line++;
-}
-
-# Create table for Radiative Temperatures
-{
- my @radiative_colnames = ('Radiative Temperature', 'Linearisation Temperature');
- $model->write($model_current_line, 0, \@radiative_colnames, $title);
- $model_current_line++;
- $model->write_number($model_current_line, 0, $radiative_temp, $unlocked);
- $model->write_number($model_current_line, 1, $linear_temp, $locked);
- $model_current_line++;
- die "Inconsistency!" unless (!defined $id_end_cells[$rad_temp_id]) && (!defined $cpt_end_cells[$rad_temp_id]) && (!defined $name_end_cells[$rad_temp_id]);
- $id_end_cells[$rad_temp_id]= "A$model_current_line";
- $cpt_end_cells[$rad_temp_id] = 0;
- $name_end_cells[$rad_temp_id] = 'Radiative_Temperature';
- # No volumic power at infinity!
- $model_current_line++;
-}
-
-# One sheet for samples
-my $green = $workbook->add_worksheet('Samples');
-$green->protect();
-
-my $samples_current_line = 0;
-if($counts[$ok_rk] > 0) {
- my @samples_colnames = ('T(sample)');
- $green->write($samples_current_line, 0, \@samples_colnames, $title);
- $samples_current_line++;
- for(my $s = 0; $s < $counts[$ok_rk] ; $s++) {
- my $sample = $samples[$s];
- my $end_type = $sample->{END_TYPE};
- my $end_id = $sample->{END_ID};
- my $pw_count = $sample->{PW_COUNT};
- my $pw_types = $sample->{PW_TYPES};
- my $pw_ids = $sample->{PW_IDS};
- my $pw_factors = $sample->{PW_FACTORS};
- my $fx_count = $sample->{FX_COUNT};
- my $fx_ids = $sample->{FX_IDS};
- my $fx_factors = $sample->{FX_FACTORS};
-
- # create a cell with formula:
- # end_temp + pw_term_1 + ... + pw_term_n + fx_term_1 + ... + fx_term_n
- unless (($end_id >= 0) && ($end_id < scalar(@seen_id_types)) && (defined $seen_id_types[$end_id]) && ($seen_id_types[$end_id] eq $end_type)) {
- print STDERR "Sample: ", Dumper($sample);
- print STDERR "END ID: $end_id\n";
- print STDERR "END TYPE: $seen_id_types[$end_id] VS $end_type\n";
- print STDERR "END TYPES: @seen_id_types\n";
- die "Invalid ID!" ;
- }
- unless ((defined $id_end_cells[$end_id]) && (defined $cpt_end_cells[$end_id])) {
- print STDERR "Sample: ", Dumper($sample);
- print STDERR "END ID: $end_id\n";
- print STDERR "END END CELLS: @id_end_cells\n";
- die "Inconsistency!"
- }
- my $formula = "=Model!$id_end_cells[$end_id]";
- $cpt_end_cells[$end_id] += 1;
- for(my $n = 0; $n < $pw_count; $n++) {
- my $ty = @$pw_types[$n];
- my $id = @$pw_ids[$n];
- my $fc = @$pw_factors[$n];
- die "Invalid ID ($id)!" unless ($id >= 0) && (defined $seen_id_types[$id]) && ($seen_id_types[$id] eq $ty);
- die "Inconsistency!" unless (defined $id_pw_cells[$id]);
- $formula = $formula."+Model!$id_pw_cells[$id]*$fc";
- $cpt_pw_cells[$id] += $fc;
- }
- for(my $n = 0; $n < $fx_count; $n++) {
- my $id = @$fx_ids[$n];
- my $fc = @$fx_factors[$n];
- die "Invalid ID ($id)!" unless ($id >= 0) && (defined $seen_id_types[$id]);
- die "Inconsistency!" unless (defined $id_fx_cells[$id]);
- $formula = $formula."+Model!$id_fx_cells[$id]*$fc";
- $cpt_fx_cells[$id] += $fc;
- }
- $green->write_formula($samples_current_line, 0, $formula, $locked);
- $samples_current_line++;
- }
-}
-
-# Create the polynom for temperature
-my $sep ="";
-my $poly = "T = ";
-for (my $n = 0; $n < scalar(@cpt_end_cells); $n++) {
- if (defined $cpt_end_cells[$n] && $cpt_end_cells[$n] != 0) {
- my $w = $cpt_end_cells[$n] / $counts[$ok_rk];
- $poly .= $sep . $w . " * $name_end_cells[$n]";
- $sep = " + ";
- }
-}
-for(my $n = 0; $n < scalar(@cpt_pw_cells); $n++) {
- if (defined $cpt_pw_cells[$n] && $cpt_pw_cells[$n] != 0) {
- my $w = $cpt_pw_cells[$n] / $counts[$ok_rk];
- $poly .= $sep . $w . " * $name_pw_cells[$n]";
- $sep = " + ";
- }
-}
-for(my $n = 0; $n < scalar(@cpt_fx_cells); $n++) {
- if (defined $cpt_fx_cells[$n] && $cpt_fx_cells[$n] != 0) {
- my $w = $cpt_fx_cells[$n] / $counts[$ok_rk];
- $poly .= $sep . $w . " * $name_fx_cells[$n]";
- $sep = " + ";
- }
-}
-
-
-# The MC Estimator and STDERR table
-my @result_colnames = ('Estimate', 'Sigma');
-$model->write($model_current_line, 0, \@result_colnames, $title);
-$model_current_line++;
-$model->write_formula($model_current_line, 0, "AVERAGE(Samples!A2:A$samples_current_line)", $locked);
-# Here we use STDEVP (Standard dev for a entire population) to remains coherent with stardis-solver
-# One could consider STDEV (Standard dev for a sampling of the population) more appropriate
-$model->write_formula($model_current_line, 1, "STDEVP(Samples!A2:A$samples_current_line)/SQRT($samples_current_line-1)", $locked);
-$model_current_line++;
-
-$model->write_string($model_current_line, 0, $poly);
-$workbook->close();
-
-print STDERR "$poly\n";
