composite_function.h

浏览该文件的文档.

// The libMesh Finite Element Library.
// Copyright (C) 2002-2023 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner

// 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., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

#ifndef LIBMESH_COMPOSITE_FUNCTION_H
#define LIBMESH_COMPOSITE_FUNCTION_H

// libMesh includes
#include "libmesh/dense_vector.h"
#include "libmesh/function_base.h"
#include "libmesh/int_range.h"
#include "libmesh/libmesh.h"
#include "libmesh/point.h"

// C++ includes
#include <algorithm>
#include <utility>
#include <vector>

namespace libMesh
{

template <typename Output=Number>
class CompositeFunction : public FunctionBase<Output>
{
public:
  explicit
  CompositeFunction() = default;

  CompositeFunction(CompositeFunction &&) = default;
  CompositeFunction & operator= (CompositeFunction &&) = default;

  CompositeFunction(const CompositeFunction &) = delete;
  CompositeFunction & operator= (const CompositeFunction &) = delete;

  virtual ~CompositeFunction() = default;

  void attach_subfunction(const FunctionBase<Output> & f,
                           std::vector<unsigned int> index_map)
  {
    const unsigned int subfunction_index =
      cast_int<unsigned int>(subfunctions.size());
    libmesh_assert_equal_to(subfunctions.size(), index_maps.size());

    subfunctions.push_back(f.clone());

    unsigned int max_index =
      *std::max_element(index_map.begin(), index_map.end());

    if (max_index >= reverse_index_map.size())
      reverse_index_map.resize
        (max_index+1, std::make_pair(libMesh::invalid_uint,
                                     libMesh::invalid_uint));

    for (auto j : index_range(index_map))
      {
        libmesh_assert_less(index_map[j], reverse_index_map.size());
        libmesh_assert_equal_to(reverse_index_map[index_map[j]].first,
                                libMesh::invalid_uint);
        libmesh_assert_equal_to(reverse_index_map[index_map[j]].second,
                                libMesh::invalid_uint);
        reverse_index_map[index_map[j]] =
          std::make_pair(subfunction_index, j);
      }

    // 检查是否有时间依赖性
    // 我们只检查刚刚添加的函数，而不是重新研究所有子函数
    // 如果这是第一个子函数，则确定时间依赖性。
    if (subfunctions.size() == 1)
      this->_is_time_dependent = subfunctions[0]->is_time_dependent();

    // 否则，我们已经有了多于1个函数。
    // 如果 _is_time_dependent 为 true，则之前的某个子函数是时间依赖的，因此此 CompositeFunction 是时间依赖的。
    // 如果 _is_time_dependent 为 false，则刚刚添加的子函数确定时间依赖性。
    else if (!this->_is_time_dependent)
      this->_is_time_dependent = (subfunctions.back())->is_time_dependent();

    index_maps.push_back(std::move(index_map));
  }

  virtual Output operator()(const Point & p,
                             const Real time = 0) override
  {
    return this->component(0,p,time);
  }

  virtual void operator()(const Point & p,
                           const Real time,
                           DenseVector<Output> & output) override
  {
    libmesh_assert_greater_equal (output.size(),
                                  reverse_index_map.size());

    // 必要时，清零未被任何子函数覆盖的输出分量
    output.zero();

    DenseVector<Output> temp;
    for (auto i : index_range(subfunctions))
      {
        temp.resize(cast_int<unsigned int>(index_maps[i].size()));
        (*subfunctions[i])(p, time, temp);
        for (auto j : index_range(temp))
          output(index_maps[i][j]) = temp(j);
      }
  }

  virtual Output component(unsigned int i,
                            const Point & p,
                            Real time) override
  {
    if (i >= reverse_index_map.size() ||
        reverse_index_map[i].first == libMesh::invalid_uint)
      return 0;

    libmesh_assert_less(reverse_index_map[i].first,
                        subfunctions.size());
    libmesh_assert_not_equal_to(reverse_index_map[i].second,
                                libMesh::invalid_uint);
    return subfunctions[reverse_index_map[i].first]->
      component(reverse_index_map[i].second,p,time);
  }

  virtual std::unique_ptr<FunctionBase<Output>> clone() const override
  {
    CompositeFunction * returnval = new CompositeFunction();
    for (auto i : index_range(subfunctions))
      returnval->attach_subfunction(*subfunctions[i], index_maps[i]);
    return std::unique_ptr<FunctionBase<Output>> (returnval);
  }

  unsigned int n_subfunctions() const
  {
    return subfunctions.size();
  }

  unsigned int n_components() const
  {
    return reverse_index_map.size();
  }

private:
  // 填充值的函数列表
  std::vector<std::unique_ptr<FunctionBase<Output>>> subfunctions;

  std::vector<std::vector<unsigned int>> index_maps;


  std::vector<std::pair<unsigned int, unsigned int>> reverse_index_map;
};


} // namespace libMesh

#endif // LIBMESH_COMPOSITE_FUNCTION_H