在C中,有更好的方法来计算不确定性值吗?

时间:2013-05-22 16:57:33

标签: c macros c-preprocessor conditional-compilation

我一直在做的是使用#define UNC(不确定性)来切换计算x(值)和dx(不确定性)的功能。它工作得很好,但它不容易读或写。有一个更好的方法吗?在C ++中,我可能使用std::pair,或者某种元组。我也可以使用同时包含xdx的结构,有时会保留dx未定义,但我计划在非常大的范围内运行此代码并且不要UNC后,通过处理所有这些不必要的项目,不希望它减慢编译器的速度。

你们有没有遇到类似的问题?

以下是代码,供参考:

#include <stdio.h>

#define UNC

typedef struct{
    double x;
#ifdef UNC
    double dx;
#endif
    } Foo;

int force(Foo* m, Foo* a, Foo* f){
    f->x = m->x * a->x;
#ifdef UNC
    f->dx = (f->x)*(m->dx/m->x + a->dx/a->x);
#endif
    return 0;
}

int main(){
    Foo m; m.x = 3.0;
    Foo a; a.x = 2.0;
#ifdef UNC
    m.dx = 0.3;
    a.dx = 0.2;
#endif
    Foo f;
    force(&m,&a,&f);
    printf(" f is %f \n", f.x);
#ifdef UNC
    printf("df is %f \n", f.dx);
#endif
    return 0;
}

1 个答案:

答案 0 :(得分:2)

您可以按如下方式定义一些宏,这样可以解决问题......

#ifdef UNC
    #define SetDx(f, v) (f)->dx = (v)
    #define PRINT_UNC(x) printf x
#else
    #define SetDx(f, v)
    #define PRINT_UNC(x)
#endif

然后将您的代码编写为

int force(Foo* m, Foo* a, Foo* f){
    f->x = m->x * a->x;
    SetDx(f, (f->x)*(m->dx/m->x + a->dx/a->x));
    return 0;
}

int main(){
    Foo m; m.x = 3.0;
    Foo a; a.x = 2.0;
    SetDx(&m, 0.3);
    SetDx(&a, 0.2);

    Foo f;
    force(&m,&a,&f);
    printf(" f is %f \n", f.x);
    PRINT_UNC(("df is %f \n", f.dx));
    return 0;
}

这样,您可以最小化代码中#ifdef块的数量,这使得难以阅读。当您转动主要定义时,SetDxPRINT_UNC宏将不执行任何操作。 (注意需要在PRINT_UNC宏中加倍括号。)

只是一个关于你的问题部分的快速问题:

  

...通过处理所有这些不必要的项来减慢编译器的速度   一旦UNC关闭......

您是否担心预处理速度?我听起来好像你可能会担心编译器花费时间编译 ifdef作为一种切换......如果是这种情况,不要担心,预处理器会杀死未定义的块,所以编译器永远不会看到它们。如果您担心预处理速度,那么您可以做的事情并不多,我不会想到,但我不认为预处理开销会很大吗?