privateのメンバ変数の値が勝手に変わる

前提

c++でclassのprivateなメンバ変数が意図せず変わります。
関数を使用するとその後に変数の値が0になります。

実現したいこと

ソースコード中の、printRate関数を使用した後もメンバ変数の値を変えないようにする。
それと、原因が分かる人は教えてほしいです。

発生している問題

printRate関数の前後で変数の値が変わっている。

0.0010653
|100.00%|[==================================================>] : solveK
0.00

該当のソースコード

C++
1#include <iostream>
2#include <fstream>
3#include <iomanip>
4#include <cmath>
5#include <stdlib.h>
6#include <string>
7#include <time.h>
8#include <random>
9#define pi 3.1415926f
10using namespace std;
11using std::ofstream;
12
13void printRate(float rate,string text=""){
14    fflush(stdout);
15    cout << "\r";
16    int width = 50;
17    int before = width*rate;
18    int after = width-before-1;
19    cout << "|" << std::setw(6) << fixed << setprecision(2)  << rate*100.0 << "%" << "|" << "[";//6文字分のスペースを確保、fixedで小数であることを伝える、小数点以下2桁表示
20    //cout << "|" << setw(6)  << setprecision(2)  << rate*100.0 << "%" << "|" << "[";
21    for(int i=0;i<before;i++){
22        cout << "=";
23    }
24    cout << ">";
25    for(int i=0;i<after;i++){
26        cout << ".";
27    }
28    cout << "] : " << text;
29}
30
31float solveK_HERE(float* m,float** v,int N){
32    float K=0;
33    int n=N/64;
34    for(int i=0;i<N;i++){
35        float v2 = v[i][0]*v[i][0] + v[i][1]*v[i][1] + v[i][2]*v[i][2];
36        K += m[i]*v2/2.0;
37        if((i+1)%n==0){
38            printRate((float)(i+1)/(float)N,"solveK");
39        }
40    }
41    cout << endl;
42    return K;
43}
44
45class Parameter {
46    public:
47    double R,T,K;
48    double factor_L,factor_T,factor_M;
49
50    void setKbar(float value){
51        k_bar = value;
52    }
53
54    float readKbar(){
55        return k_bar;
56    }
57    
58    //初期化
59    Parameter(){
60        factor_L = 1.5e10;
61        factor_T = 3.1e7;
62        factor_M = 2e30;
63        R = 8.314f;
64        T = 30.;
65    }
66
67    //クラス関数
68    void initialize();
69    void printInfo();
70
71    private:
72    float k_bar;
73};
74
75void Parameter::initialize(){
76    setKbar(R*T /pow(factor_L/factor_T,2));
77}
78
79int main(void){
80    Parameter para;
81    para.initialize();
82    
83    cout << para.readKbar() << endl;
84    
85    float* m=new float[100];
86    float** v=new float*[100];
87    for(int i=0;i<100;i++){
88        v[i] = new float[3];
89    }
90    solveK_HERE(m,v,100);
91
92    //この次にk_barを表示すると中身が0になってる
93    cout << para.readKbar();
94    
95    
96}
97
98