102 lines
2.4 KiB
C++
102 lines
2.4 KiB
C++
#include "bsp_1_b.h"
|
|
|
|
#include <cassert>
|
|
#include <fstream>
|
|
#include <iostream>
|
|
#include <stdexcept>
|
|
#include <string>
|
|
#include <vector>
|
|
#include <cmath>
|
|
using namespace std;
|
|
|
|
|
|
void meandevvec(vector<short int> x, float &a, float &g, float &h, float &d)
|
|
{
|
|
a = 0;
|
|
for(size_t k=1; k<=x.size(); k=k+1)
|
|
{
|
|
a = a+x.at(k-1);
|
|
}
|
|
a = a/x.size();
|
|
|
|
g = 1;
|
|
for(size_t k=1; k<=x.size(); k=k+1)
|
|
{
|
|
g = g*pow(x.at(k-1),1.0/x.size());
|
|
}
|
|
|
|
h = 0;
|
|
for(size_t k=1; k<=x.size(); k=k+1)
|
|
{
|
|
h = h + 1.0/x.at(k-1);
|
|
}
|
|
h = x.size()/h;
|
|
|
|
d = 0;
|
|
for(size_t k=1; k<=x.size(); k=k+1)
|
|
{
|
|
d = d + pow(x.at(k-1)-a,2);
|
|
}
|
|
d = d/x.size();
|
|
d = sqrt(d);
|
|
|
|
return;
|
|
}
|
|
|
|
|
|
// for reading a file and writing into a txt-file
|
|
// [Str10, p.364]
|
|
void fill_vector(istream& istr, vector<short int>& v)
|
|
{
|
|
short int d=0;
|
|
while ( istr >> d) v.push_back(d); // Einlesen
|
|
if (!istr.eof())
|
|
{ // Fehlerbehandlung
|
|
cout << " Error handling \n";
|
|
if ( istr.bad() ) throw runtime_error("Schwerer Fehler in istr");
|
|
if ( istr.fail() ) // Versuch des Aufraeumens
|
|
{
|
|
cout << " Failed in reading all data.\n";
|
|
istr.clear();
|
|
}
|
|
}
|
|
v.shrink_to_fit(); // C++11
|
|
return;
|
|
}
|
|
|
|
|
|
void read_vector_from_file(const string& file_name, vector<short int>& v)
|
|
{
|
|
ifstream fin(file_name); // Oeffne das File im ASCII-Modus
|
|
if( fin.is_open() ) // File gefunden:
|
|
{
|
|
v.clear(); // Vektor leeren
|
|
fill_vector(fin, v);
|
|
}
|
|
else // File nicht gefunden:
|
|
{
|
|
cout << "\nFile " << file_name << " has not been found.\n\n" ;
|
|
assert( fin.is_open() && "File not found." ); // exeption handling for the poor programmer
|
|
}
|
|
|
|
return;
|
|
}
|
|
|
|
void write_vector_to_file(const string& file_name, const vector<float>& v)
|
|
{
|
|
ofstream fout(file_name); // Oeffne das File im ASCII-Modus
|
|
if( fout.is_open() )
|
|
{
|
|
for (unsigned int k=0; k<v.size(); ++k)
|
|
{
|
|
fout << v.at(k) << endl;
|
|
}
|
|
}
|
|
else
|
|
{
|
|
cout << "\nFile " << file_name << " has not been opened.\n\n" ;
|
|
assert( fout.is_open() && "File not opened." ); // exeption handling for the poor programmer
|
|
}
|
|
|
|
return;
|
|
}
|