fe_vectors.h

Go to the documentation of this file.

// ----------------------------------------------------------------------------
//
// FCST: Fuel Cell Simulation Toolbox
//
// Copyright (C) 2006-2009 by Guido Kanschat
// Copyright (C) 2006-2014 by Energy Systems Design Laboratory, University of Alberta
//
// This software is distributed under the MIT License
// For more information, see the README file in /doc/LICENSE
//
// - Class: fe_vectors.h
// - Description: This class implements data type
//                used in function calls of applications
// - Developers: Guido Kanschat,     Texas A&M University
//               Valentin N. Zingan, University of Alberta
//               Marc Secanell,      University of Alberta
// - Id: $Id: fe_vectors.h 2605 2014-08-15 03:36:44Z secanell $
//
// ----------------------------------------------------------------------------

#ifndef _FUEL_CELL_APPLICATION_CORE_FE_VECTORS_H_
#define _FUEL_CELL_APPLICATION_CORE_FE_VECTORS_H_

#include <deal.II/lac/vector_memory.h>
#include <deal.II/lac/block_vector.h>

#include <algorithm>
#include <fstream>

using namespace dealii;

namespace FuelCell
{
namespace ApplicationCore
{

typedef BlockVector<double> FEVector;

class FEVectors : public Subscriptor
{
public:

  FEVectors();



  void add_scalar(double&            s,
                  const std::string& name);

  void add_vector(FEVector&          v,
                  const std::string& name);

  void add_vector(const FEVector&    v,
                  const std::string& name);

  void merge(FEVectors& other);

  void merge(const FEVectors& other);




  unsigned int n_scalars() const;

  unsigned int n_vectors() const;

  double& scalar(unsigned int i);

  const double& scalar(unsigned int i) const;

  FEVector& vector(unsigned int i);

  const FEVector& vector(unsigned int i) const;

  const std::string& scalar_name(unsigned int i) const;

  const std::string& vector_name(unsigned int i) const;

  unsigned int find_scalar(const std::string& name) const;

  unsigned int find_vector(const std::string& name) const;

  unsigned int count_vector(const std::string& name) const;

  template<typename OUT>
  void list(OUT&         out,
            unsigned int verbosity = 1) const;


  DeclException2(ExcNameMismatch,
                 int,
                 std::string,
                 << "Name at position " << arg1 << " is not equal to " << arg2);

private:

  bool is_constant;

  std::vector<double*> scalars;

  std::vector<std::string> scalar_names;

  std::vector<FEVector*> vectors;

  std::vector<std::string> vector_names;
};

} // ApplicationCore

} // FuelCell

#endif