ycwang125-CSDN博客

原创在C++中混编调用CHEMKIN

填周六的坑，CHEMKIN是一个燃烧经常使用的Fortran代码包。笔者当前的代码为C++的MPI并行代码，为了在代码中使用CHEMKIN的模块，需要进行C++和Fortran的混编。基础的混编测试可以查看前一篇https://blog.csdn.net/qq_40583925/article/details/123024724调试花了几天的时间，这里先讲需要注意的点，然后直接贴代码~编译选项的使用前面已经遇到过C++代码混编和C代码混编不同的问题，解决方案便是使用extern "C"功能，指示编

2022-02-21 17:52:16 1256

原创 Linux下C++和Fortran的混编

事情的起因是我需要计算多组分燃烧问题中的粘性项，实现起来非常麻烦。并且CHEMKIN中有现成的代码，但问题在于我的程序为.cpp的MPI并行代码。将Fortran代码翻译成C++代码再用也可，但是仅仅计算粘性的单独模块TRANLIB.F大概在2000行左右，改写起来很麻烦。因此！干脆将C++和Fortran混编起来吧！为此专门研究了一下怎么弄C语言和Fortran混编C语言和Fortran混编的博客还是挺多的，虽然网上好像在互相抄同一段代码。似乎最早的源头是这篇https://www.cnblogs.

2022-02-19 22:21:58 1926

原创 Shell学习笔记

通过几个例子，快速的学习或者复习Shell脚本的使用方法

2021-12-19 23:36:48 681

原创 Tecplot进阶——如何用Tecplot制作一张满足论文投稿要求的图片

如题，一两年前刚刚做Tecplot入门，满足自己基础的作图需求。之后投了自己的第一篇文章，Tecplot作为作图主力被大量使用。除了基础功能外，还有很多额外的操作被用到，这里进行一个简单的总结。基础操作https://blog.csdn.net/qq_40583925/article/details/101852339?spm=1001.2014.3001.5502嗯，之前写过这个，这里主要做二维的图。我们以如下一个图为例讲解。调整XY的取值范围双击坐标轴，然后可以得到一系列的调整，选到Rang

2021-04-14 15:08:18 29156 2

原创杂谈——并行计算的四种道路

马上要找工作了呀…虽然鄙人方向是计算流体力学，但是这个方向找工作，呃…微妙…所以…我只好碰瓷人家并行方向。嘛，不过虽然不是专门做并行的，不过并行大规模的使用是没有瞎说的。这里在正式的面试之前，进行一个简单的总结。目前笔者掌握的并行计算的方法有四种，分别是openMP并行，多线程并行MPI并行，以及CUDA并行。从前往后，从串行代码改写到并行代码的麻烦程度递增，不过加速比也是越来越好。这里要分别进行介绍(但愿面试第一轮面试前我能写完这个)openMP并行OK！我们先从最简单的开始，openMP并行其实是

2021-03-31 22:06:28 862 3

原创 win下的Qt Creator如何在C++代码使用MPI 和 openMP

如题，我们希望在Windows下的Qt Creator中运行如下的一段并行代码：#include <math.h>#include "mpi.h"#include <iostream>using namespace std;int main(int argc, char *argv[]){ /*myid 和 numprocs分别记录某一个并行执行进程的标识和所有参加计算的进程的个数*/ int myid, numprocs; MPI_Init

2020-11-24 20:24:23 2011 5

原创 OpenFOAM中的Chemistry类

化学反应相关的类，以热学相关的类为基础，它位于路径src/chemistryModel/。它的类间关系如下：接下来我们逐个介绍这些类的内容basicChemistryModel最基础的类，头文件内容如下：class basicChemistryModel: public IOdictionary{protected: //- Reference to the mesh database const fvMesh& mesh_;

2020-08-06 21:35:58 1701 1

原创 openFOAM中文件流的获取——lookup函数

openFOAM中的几乎所有类都提供了文件流操作的重定义，即>> <<即可实现文件流操作。例如对于标量来说，读取标量时：scalar readScalar(Istream& is){ scalar rs; is >> rs; return rs;}我们只需要提供文件流is，就可以通过调用函数返回一个标量值。那么这样的文件流如何获取呢？这里讨论这部分内容这里介绍的是lookup函数，他位于路径src/OpenFOAM/db/d

2020-08-05 17:29:02 1895

原创 openFOAM的基础数据结构汇总

openFOAM将数组链表等数据结构也进行了封装，这里进行一个汇总。可能陆续也会更新标量scalar其实就是浮点数，不过浮点数有多重类型精度，比如float double longdouble，这里将类型统一为scalar这个类型使用，具体见如下链接：https://blog.csdn.net/qq_40583925/article/details/107735351向量vector...

2020-08-05 15:02:24 2183 1

原创 openFOAM中的label

其实就是指i,j,k这类浮标使用的类型。我们通常就使用int就可以，但是这里也进行了封装。其头文件的源码如下：#include "int.H"#define INT_ADD_SIZE(x,s,y) x ## s ## y#define INT_ADD_DEF_SIZE(x,s,y) INT_ADD_SIZE(x,s,y)#define INT_SIZE(x,y) INT_ADD_DEF_SIZE(x,WM_LABEL_SIZE,y)#if WM_LABEL_SIZE != 32 &&

2020-08-05 15:01:56 1462

原创 openFOAM中的Field

就是域的创建，速度域等等量直接在这里存储他分别位于若干个位置，所以下面会单独罗列，其类间关系如下这里对其中的Fields DimensionedFields GeometricFields进行解读Field其路径为src/OpenFOAM/fields/Fields/Field他包含若干个头文件：FieldFunctions.HFieldFunctionsM.HFieldM.HField.CFieldFunctions.CFieldFunctionsM.C我们从头文件开始解读，代

2020-08-05 11:56:08 2323 5

原创 openFOAM中的tensor

是指3*3的张量，其路径为src/OpenFOAM/primitives/Tensor对应的文件关系如下：基本结构和vector非常相似MatrixSpace他用VectorSpace创建，头文件如下：template<class Form, class Cmpt, direction Mrows, direction Ncols>class MatrixSpace: public VectorSpace<Form, Cmpt, Mrows*Ncols>{publ

2020-08-01 22:30:54 680

原创 openFOAM中的vector

首先需要声明这里的vector和我们平时C++中说的数组，或者STL中的vector不是一个。它是openFOAM从九十年代改动至今留下来的一个习惯，这里的vector特质(x,y,z)这样的三维向量，用来表示坐标速度等等。他位于路径OpenFOAM/primitives/Vector/complexVector我们首先来看一下文件的依赖关系VectorSpace头文件内容如下：template<class Form, class Cmpt, direction Ncmpts>cl

2020-08-01 20:41:25 1485 6

原创 openFOAM中的Scalar

位于路径src/OpenFOAM/primitives/Scalar标量，其实就是浮点数。但是浮点数具有不同的精度，最终文件将根据宏的设置决定程序使用哪种精度文件依赖关系如下：注意这里大小写的scalar是两个文件，大写的Scalar是各个精度的浮点数的依赖文件，而这些具体的精度的实现又是小写的scalar的实现floatScalar这里给出其中一个文件的具体实现方式，其余的类似头文件内容如下：typedef float floatScalar;定义最大最小的浮点数floatScala

2020-08-01 19:50:01 1867 2

原创 openFOAM学习笔记(五)——chemFoam的运行过程

在前面的帖子中已经大概给出了chemFoam主程序的结构，这里给出一个比较全面的总结首先程序结构如下：添加头文件//*****************************//int main(int argc, char *argv[]){ 初始化 while(runTime.run()) { 计算化学反应输出中间结果 } return 0;}我们接下来介绍各个部分的具体内容while循环中的燃烧计算其源码如下：while (run

2020-07-31 15:38:45 3843

原创 openFOAM中时间步长计算参数是怎样读进程序的——Time类

接着前面的内容，在createTime.H中创建了Time类的对象runTimeFoam::Info<< "Create time\n" << Foam::endl;Foam::Time runTime(Foam::Time::controlDictName, args);这里，我们尝试解读这第二行代码中更加具体的内容。Time类的内容源码有500+行，这里不粘贴了，只给出大概的内容，首先文件使用了一系列的头文件，后面会陆续解读：#include "TimePaths.

2020-07-15 16:55:28 1998 1

原创 openFOAM的基础类型汇总

openFOAM代码非常非常多，但是让人比较郁闷的是很多是基础类的实现，openFOAM中很多基础的类，包括char,string,int等等，都添加功能封装在了Foam这个名字域中。好处就是同样使用int，可以比在C++使用更多的功能，坏处就是一个很简单的功能，依赖的文件数量惊人。把这么多的封装全部读完总结完可能头发都白了，但是只把常用的拿出来总结一下是可以的，这篇帖子就是做这个工作，当然这里是针对chemFoam这个求解器。这个帖子会在编写其他帖子的过程中持续更新基础类的源码组织特点通常为一个x

2020-07-15 15:16:18 3769

原创 openFOAM中时间步长计算参数是怎样读进程序的——初步理解框架

内容如题，在阅读完大量的基础代码之后，我们现在可以尝试解读一下chemFOAM中的程序流程。因为内容很多，所以会分好几个子标题讲完。我们使用openFOAM的过程中，通常是从tutorials中找到对应算例的配置文件，复制到run的文件夹中，然后使用编译好程序运行，而控制则通过配置文件中的参数设置。那么，首先我们需要理解这些控制文件中的参数是怎样读入到程序中。配置文件的内容本文以chemFOAM中的氢气燃烧算例为例，对应的文件夹为tutorials/combustion/chemFoam/h2，内容有

2020-07-10 21:21:10 3053 1

原创 openFOAM学习笔记（四）—— openFOAM中的List

又是一个很底层的部分，但是也非常重要

2020-06-27 23:13:55 3269 1

原创 openFOAM中的forAll

这是用户二次开发代码中非常常用的一个语法，我们这里看一下它的底层实现。当然不同的class会各自实现一个forall，这里给出的是UList的实现，具体路径位于src/OpenFOAM/container/Lists/UList/文件夹中定义了UList类，使用一个连续存储，通过制定类型T的指针给定首地址template<class T>class UList { // Private Data //- Number of elements in UList

2020-06-27 22:20:40 2196

原创 openFOAM学习笔记（三）—— char和string相关的类

openFOAM中的char和string并不是直接使用C++提供的string，而是进行了一系列的封装，并添加了更多的功能。这里进行一个总结。char类型的更多功能openFOAM中并没有对char进行封装，而是直接在上面添加功能。主要依赖的是如下路径的几个文件：src/OpenFOAM/primitives/chars/charchar.H charIO.Csrc/OpenFOAM/primitives/chars/wcharwchar.H wcharIO.C添加了readCha

2020-06-25 22:38:15 907

原创 openFOAM代码阅读——基础代码中的string

代码仍然是src/OpenFOAM/primitives/路径中的，为文件夹strings，其中有内容fileName keyType lists string stringOps word wordRe这里的多个文件夹都使用到了这里的word文件夹中的文件，而word用到了string中的内容，所以我们先读这个吧。这里包含了如下几个文件：string.C string.H stringI.H stringIO.C stringIOList.C stringIOList.H从string.H开始

2020-06-11 23:52:43 929

原创 openFOAM代码阅读——基础代码中的wchar

路径src/OpenFOAM/primitives/chars中还有另外一个文件夹wchar，我们这里来看看这里面有什么。头文件wchar.H的内容如下：#include <cwchar>#include <string>// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //namespace Foam{class Istream;class Ostream;

2020-06-11 17:18:56 519

原创 openFOAM代码阅读——基础代码中的char

openFOAM的源码文件夹src中的内容如下：dyfluid@dyfluid:~/OpenFOAM/OpenFOAM-7/src$ lsAllwmake genericPatchFields rigidBodyDynamicsatmosphericModels lagrangian rigidBodyMeshMotioncombustionModels mesh rigidBodyState

2020-06-11 16:41:58 1671

原创 openFOAM代码阅读——综述和版本说明

这里阅读的是openFOAM7的源码。openFOAM的教程要多少有多少，各种书籍线上教程，甚至还有培训班。但是openFOAM的代码的真正解读的资料缺并不多，这里笔者只好自行补齐部分基础的代码阅读。首先说明一下我们阅读的主要内容。openFOAM中的代码主要分了三个部分：其一，是用来进行计算的参数设置网格设置文件，在文件夹tutorials中其二，是用来支持各个类型的算例运行的具体求解器的源码，在文件夹application中其三，是application的支持性文件，也就是真正的源码所在的文

2020-06-11 16:13:33 3228

原创 chemFoam的源码提取

ok！经过一系列的铺垫，我们来尝试一下提取出openFOAM中的chemFoam程序的源码，然后用Makefile编译运行一下，这样可以帮助我们进行二次开发。呃。。。这个怎么说呢。。。openFOAM一个单独算例的程序依赖深度超过我的预期了。首先程序编译的过程是有记录的，保存如下路径下的文件：dyfluid@dyfluid:~/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/applications/solvers/combustion/chemFo

2020-06-09 15:49:18 852

空空如也

WinEdt7.0使用dvi2pdf生成pdf文件时eps图片无法显示完全如何解决