vector_valued_function.h

/* Copyright (C) 2004-2014 MBSim Development Team
 *
 * This library is free software; you can redistribute it and/or 
 * modify it under the terms of the GNU Lesser General Public 
 * License as published by the Free Software Foundation; either 
 * version 2.1 of the License, or (at your option) any later version. 
 *  
 * This library 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 
 * Lesser General Public License for more details. 
 *  
 * You should have received a copy of the GNU Lesser General Public 
 * License along with this library; if not, write to the Free Software 
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301 USA
 *
 * Contact: martin.o.foerg@gmail.com
 */
 
#ifndef _VECTOR_VALUED_FUNCTIONS_H_
#define _VECTOR_VALUED_FUNCTIONS_H_
 
#include "mbsim/functions/function.h"
#include "mbsim/utils/utils.h"
 
namespace MBSim {
 
  template<typename Sig> class VectorValuedFunction;
 
  template<typename Ret, typename Arg>
  class VectorValuedFunction<Ret(Arg)> : public Function<Ret(Arg)> {
    using B = fmatvec::Function<Ret(Arg)>; 
 
    public:
      VectorValuedFunction() = default;
      VectorValuedFunction(std::vector<Function<double(Arg)> *> component_) : component(std::move(component_)) {
        for(auto it=component.begin(); it!=component.end(); ++it)
          (*it)->setParent(this);
      }
      ~VectorValuedFunction() override {
        for (unsigned int i=0; i<component.size(); i++)
          delete component[i];
      }
      void addComponent(Function<double(Arg)> *function) {
        component.push_back(function);
        function->setParent(this);
      }
      int getArgSize() const override { return 1; }
      std::pair<int, int> getRetSize() const override { return std::make_pair(component.size(),1); }
      Ret operator()(const Arg &x) override {
        for(unsigned int i=0; i<component.size(); i++)
          y[i]=(*component[i])(x);
        return Ret(y);
      }
      typename B::DRetDArg parDer(const Arg &x) override {
        for(unsigned int i=0; i<component.size(); i++)
          dy[i]=component[i]->parDer(x);
        return FromStdvec<typename fmatvec::Der<double, Arg>::type>::cast(dy);
      }
      typename B::DRetDArg parDerDirDer(const Arg &xDir, const Arg &x) override {
        for(unsigned int i=0; i<component.size(); i++)
          ddy[i]=component[i]->parDerDirDer(xDir,x);
        return FromStdvec<typename fmatvec::Der<double, Arg>::type>::cast(ddy);
      }
      void initializeUsingXML(xercesc::DOMElement *element) override {
        xercesc::DOMElement *e=MBXMLUtils::E(element)->getFirstElementChildNamed(MBSIM%"components")->getFirstElementChild();
        while (e) {
          addComponent(ObjectFactory::createAndInit<Function<double(Arg)>>(e));
          e=e->getNextElementSibling();
        }
      }
      void init(Element::InitStage stage, const InitConfigSet &config) override {
        Function<Ret(Arg)>::init(stage, config);
        for(auto it=component.begin(); it!=component.end(); ++it)
          (*it)->init(stage, config);
        if(stage==Element::preInit) {
          y.resize(component.size());
          dy.resize(component.size());
          ddy.resize(component.size());
        }
      }
    private:
      std::vector<Function<double(Arg)> *> component;
      std::vector<double> y;
      std::vector<typename fmatvec::Der<double, Arg>::type> dy ,ddy;
  };
 
}
 
#endif