1

我正在使用 libMesh FEM 库并尝试开发一个继承自 libMesh 的类（EqCore）。这个类将提供一些额外的特性，这些特性被我想要实际使用的类（MainEq）再次继承。

这两个函数 set_constant 和 get_constant 导致以下错误。它们使用不同的继承方案如图所示工作（请参阅模板类的继承与 C++ 中的模板成员函数）。这个问题的不同之处在于，现在模板参数（Type）实际上是一个被继承的类。这是一种危险的做法吗？

我将不胜感激使此代码正常工作或找到替代方法的任何帮助。

错误信息：

在成员函数 'void EqCore::set_constant(std::string, ParamType)': test_libmesh.cpp:26:57: error: expected primary-expression before '>' token

在成员函数“ParamType EqCore::get_constant(std::string)”中：/home/slaughter/Documents/programs/source/test_libmesh.cpp:31:76: 错误：“>”标记之前的预期主表达式

程序：

//! \example test_libmesh.cpp

#include <string>
using std::string;

// libMesh includes
#include <libmesh.h>
#include <libmesh_common.h> 
#include <equation_systems.h>
#include <transient_system.h>
#include <explicit_system.h>
#include <parameters.h>
#include <mesh.h>
using namespace libMesh;

// Fundamental behavior that will be used among many classes
template <typename Type> class EqCore : Type{
    public:

        // Class constructor
        EqCore(EquationSystems& sys, string name) : Type(sys, name, 1){}

        // A function for storing a constant value (causes error)
        template<typename ParamType> void set_constant(std::string name, ParamType var){  
            Type::get_equation_systems().parameters.set<ParamType>(name) = var;
        }

        // A function for retrieving a constant value (causes error)
        template<typename ParamType> ParamType get_constant(std::string name){  
            ParamType output = Type::get_equation_systems().parameters.get<ParamType>(name);
            return output;
        } 
};

// A test class derived
class MainEq : public EqCore<ExplicitSystem>{
    public: 

        // Constructor
        MainEq(EquationSystems& sys) : EqCore(sys, "main"){ }   

};  


// Begin main function
int main (int argc, char** argv){

    // Initialize libMesh and create an empty mesh
    LibMeshInit init (argc, argv);
    Mesh mesh;

    // Test w/o any of the above classes
    EquationSystems eq_sys(mesh);
    eq_sys.parameters.set<double>("test1") = 1;
    printf("Test 1: %f\n", eq_sys.parameters.get<double>("test1"));

    // Test my class set/get functions
    MainEq eq(eq_sys);
    eq.set_constant<double>("test2", 2);
    printf("Test 2: %f\n", eq.get_constant<double>("test2"));   
}
4

2 回答 2

6

因为您在模板中，所以编译器无法set在解析期间自动确定这是一个模板，并且假设它set是一个非模板，因此解析失败。

解决方案是显式通知编译器set是成员模板，因此。

Type::get_equation_systems().parameters.template set<ParamType>(name) = var

于 2012-08-13T20:07:09.793 回答
2

在C++ Template Metaprogramming: Concepts, Tools, and Techniques from Boost and Beyond中，David Abrahams 和 Aleksey Gurtovoy ( Amazon ) 解释如下：

double const pi = 3.14159265359;

template <class T>
int f(T& x)
{
    return x.convert<3>(pi);
}

T::convert可能是成员函数模板，在这种情况下，突出显示的代码传递pi给convert<3>. 它也可能是一个数据成员，在这种情况下f返回 (x.convert < 3 ) > pi。这不是一个非常有用的计算，但编译器不知道。

关键字告诉编译器template依赖名称是成员模板：

template <class T>
int f(T& x)
{
    return x.template convert<3>(pi);
}

如果我们省略template，编译器假定x.convert没有命名模板，并且它后面的 < 被解析为小于运算符。

于 2012-08-13T20:20:22.020 回答