experimental_data.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//
// C++ Implementation: experimental_data
//
// Description: Used to read in an array of experimental data with column headers.
//              For use with parameter estimation or generating polarization curves.
//
// Author: Marc Secanell <secanell@ualberta.ca>, Peter Dobson <pdobson@ualberta.ca> (C) 2009,2010
//
// Copyright: See COPYING file that comes with this distribution
//
// Created by: Marc Secanell, Peter Dobson (December 2009)
// Last updated: January, 2010
//
//---------------------------------------------------------------------------

#ifndef _FUELCELL__EXPERIMENTAL_DATA_H
#define _FUELCELL__EXPERIMENTAL_DATA_H

//appframe
#include <application_core/fe_vectors.h>
#include <application_core/application_data.h>
#include <application_core/event.h>
#include <application_core/application_wrapper.h>

//STL libraries:
#include<string>
#include<fstream>
#include<iostream>
#include<vector>

using namespace dealii;

namespace SIM
{
  class ExperimentalData
  {
  public:

    ExperimentalData(std::string filename);
    ~ExperimentalData()
        {};

    void parse_data_names() const;

    void print_data() const;
    void extract_vector(std::string label, std::vector<double>& data_vector);

    void get_experimental_values(std::vector<std::string>& name_OC, std::vector< std::vector<double> >& value_OC) const;

//     /**
//      * function that gets the string vector and value array
//      */
//     void store_OCV_values(std::vector<std::string>& name_OC, std::vector< std::vector<double> >& value_OC, std::vector<double>& current);

//     /**
//      * function that gets the string vector and value array
//      */
//     bool find_OCV(std::vector<double>& conditions, double& OCV);

   // Inline Functions
    inline int get_num_rows() const
        {
        return experimental_values.size();
        }
    inline int get_num_columns() const
        {
        return experimental_values[0].size();
        }


  private:
    void read_data_file();

    std::vector< std::vector<double> > experimental_values;
    std::vector<std::string> data_names;
    const std::string data_file;
    std::vector<double> OCV_list;

    std::vector <std::vector<double> > OCV_conditions;

  };

}

#endif