stencilet.cpp

来自「c++经典教材 Blitz++ v0.8」· C++ 代码 · 共 22 行

#include <blitz/array.h>
#include <blitz/array/stencil-et.h>

using namespace blitz;

int main()
{
    Array<float,2> A(Range(-7,7),Range(-7,7)), B(Range(-7,7),Range(-7,7));
    A = 0.0;
    A(0,0) = 1;
    B = 0.0;

    Range I(-5,5), J(-5,5);

    Array<float,2> Ai = A(I,J), Bi = B(I,J);

    Bi = Laplacian2D4n(Ai) * 0.1;

    cout << Bi << endl;
}