OpenFOAM开发编程基础01 输入输出

// terminalcd /home/aerosand/aerosand/ofspmkdir 01_IOcd 01_IO

#include <iostream> // IO标准库#include <string> // 字符串库#include <fstream> // 文件流库#include <cassert> // 异常判断库
int main(int argc, char *argv[]){    std::string name;    int year;    std::string var[3];
    std::fstream infile; // 定义一个文件流    infile.open("input.dat",std::fstream::in); // 打开文件并准备读取    if (!infile) { // 异常判断，如果打开失败，则执行下面        std::cout << "# WARNING: NO input!" << std::endl;        assert(infile); // 终止程序并抛出错误信息（频繁调用影响性能）    }    infile >> var[0] >> name;    infile >> var[1] >> year;    infile.close(); // 必须关闭文件
    std::cout << var[0] << "\t\t@\t" << var[1] << std::endl;    std::cout << name << "\t@\t" << year << std::endl;

    std::fstream outfile;    outfile.open("output.dat",std::fstream::out); // 打开文件并准备写入    outfile << var[0] <<"\t" << name << std::endl;    outfile << var[1] << "\t" << year << std::endl;    outfile.close(); // 必须关闭文件
	// 现代 C++ 可以不写 return 0;}

arg0    Aerosandarg1    2023

// terminalmake run

g++ -std=c++17 -Wall -Wextra -g -Iinclude -c -MMD src/main.cpp  -o src/main.osrc/main.cpp: In function ‘int main(int, char**)’:src/main.cpp:6:14: warning: unused parameter ‘argc’ [-Wunused-parameter]    6 | int main(int argc, char *argv[])      |          ~~~~^~~~src/main.cpp:6:26: warning: unused parameter ‘argv’ [-Wunused-parameter]    6 | int main(int argc, char *argv[])      |                    ~~~~~~^~~~~~g++ -std=c++17 -Wall -Wextra -g -Iinclude -o output/main src/main.o  -LlibExecuting all complete!./output/mainarg0		@	arg1Aerosand	@	2023Executing run: all complete!

arg0		@	arg1Aerosand	@	2023

arg0	Aerosandarg1	2023

// terminalcd /home/aerosand/aerosand/ofsp/01_IO/foamNewApp 01_02_IOcd 01_02_IOcp -r $FOAM_TUTORIALS/incompressible/icoFoam/cavity/cavity debug_casecode .

|- 01_02_IO/	|- debug_case/		|- 0/		|- constant/		|- system/	|- Make/ 		|- files		|- options	|- 01_02_IO.C

// terminalcode _appmake.sh _appclean.sh _caserun.sh _caseclean.sh README.md

#!/bin/bash
wmake

#!/bin/bash
wclean

#!/bin/bash
blockMesh -case debug_case | tee debug_case/log.meshecho "Meshing done."
01_02_IO -case debug_case | tee debug_case/log.run

#!/bin/bash
wcleanrm -rf debug_case/log.*foamCleanTutorials debug_caseecho "Cleaning done."

#include "fvCFD.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
int main(int argc, char *argv[]){    #include "setRootCase.H"    #include "createTime.H"
    #include "createMesh.H"

    /*     * Reading from the dictionary    */    const word dictName("customProperties"); // 创建word类型变量保存字典名称    IOobject dictIO // 创建IOobject类型变量，从参数列表初始化    (        dictName, // 字典名称        runTime.constant(), // 字典位置        mesh, // 和mesh相关        IOobject::MUST_READ, // 必须从字典读取        IOobject::NO_WRITE // 不向字典写入    );
    if (!dictIO.typeHeaderOk<dictionary>(true))    {        FatalErrorIn(args.executable()) << "Cannot open specified dictionary"            << dictName << exit(FatalError);    }    // 如果字典文件在文件头中指定的不是 dictionary 类型，则报错
    dictionary myDictionary;    myDictionary = IOdictionary(dictIO);    // 从IOobject变量中创建字典对象
    // 一般使用下面这种紧凑写法    Info<< "Reading myProperties\n" << endl;     IOdictionary myProperties // 字典变量名和字典文件名取相同    (        IOobject        (            "myProperties",            runTime.constant(),            mesh,            IOobject::MUST_READ,            IOobject::NO_WRITE        )    );
    word solver; // 创建word类型变量    myProperties.lookup("application") >> solver;    // 从myProperties文件中查找到关键词，并取值赋给solver
    word format(myProperties.lookup("writeFormat"));    // 或者写成更紧凑的形式
    scalar timeStep(myProperties.lookupOrDefault("deltaT", scalar(0.01)));    // 也可以写成    // scalar timeStep(myProperties.lookupOrDefault<scalar>("deltaT", 0.01));    // 如果字典中没有提供这一关键词，则使用此句提供的默认值
    bool ifPurgeWrite(myProperties.lookupOrDefault<Switch>("purgeWrite",0));    // bool类型也可以按关键词查找并读取
    List<scalar> pointList(myProperties.lookup("point"));    // 列表也可以读取
    HashTable<vector,word> sourceField(myProperties.lookup("source"));    // 哈希表也可以读取（无所谓什么是哈希表，没有关系）
	// 输出读取的内容    Info<< nl        << "application: " << solver << nl << nl        << "writeFormat: " << format << nl << nl        << "deltaT: " << timeStep << nl << nl        << "purgeWrite: " << ifPurgeWrite << nl << nl        << "point: " << pointList << nl << nl        << "source: " << sourceField << nl << nl        << endl;

    /*     * Writing to files    */    fileName outputDir = runTime.path()/"processing"; // 创建outputDir变量并赋值路径    mkDir(outputDir); // 创建上句路径的文件夹
    autoPtr<OFstream> outputFilePtr; // 输出文件流的指针    outputFilePtr.reset(new OFstream(outputDir/"myOutPut.dat")); // 给指针定向
	// 通过指针给输出文件写入信息    outputFilePtr() << "processing/myOutPut.dat" << endl;    outputFilePtr() << "0 1 2 3 ..." << endl;    sourceField.insert("U3", vector(1, 0.0, 0.0));    outputFilePtr() << sourceField << endl;
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    Info<< nl;    runTime.printExecutionTime(Info);
    Info<< "End\n" << endl;
    return 0;}

// ************************************************************************* //

FoamFile{    version     2.0;    format      ascii;    class       dictionary;    location    "constant";    object      customProperties;}

FoamFile{    version     2.0;    format      ascii;    class       dictionary;    location    "constant";    object      myProperties;}
application     icoFoam;
writeFormat     ascii;
purgeWrite      1;
point(    0    1    2);
source(    U1 (0 0 0)    U2 (1 0 0));

// terminalsh _appmake.shsh _caserun.sh

// terminal
Create time
Create mesh for time = 0
Reading myProperties

application: icoFoam
writeFormat: ascii
deltaT: 0.01
purgeWrite: 1
point: 3(0 1 2)
source: 2(U1 (0 0 0)U2 (1 0 0))


ExecutionTime = 0.01 s  ClockTime = 0 s
End

processing/myOutPut.dat0 1 2 3 ...
3(U1 (0 0 0)U3 (1 0 0)U2 (1 0 0))