蘇軍偉博文openfoam中矢量類(lèi)vector使用解析

1、OpenFOAM中向量類(lèi)vector和C+中的vector雖然同名，但功能可謂大相徑庭。C+中vector是一個(gè)模板類(lèi)，類(lèi)本身是一個(gè)數(shù)組，用來(lái)連續(xù)存放一些變量。而OpenFOAM中的vector是向量類(lèi)，而非模板(Vector為該類(lèi)的模板，vector = Vector)，vector對(duì)象有三個(gè)分量，也就是三維空間中的xyz。openfoam對(duì)幾何的點(diǎn)（point）定義就是采用向量形式，其實(shí)point為vector的一個(gè)別名，采用如下定義typedef vector point; （1） vector里面有兩個(gè)靜態(tài)變量one(1,1,1),zero(0, 0 ,0），通?？梢允褂眠@兩個(gè)值進(jìn)行初

2、始化，調(diào)用這樣個(gè)值可以采用vector:one; vector:zero. （2）vector計(jì)算相關(guān)（a b為兩個(gè)矢量） 和: a+b 差： a-b外積 a*b /得到是一個(gè)張量，很少用 內(nèi)積 a&b這里的&本身是按位操作符，OpenFOAM對(duì)其進(jìn)行重載，變成內(nèi)集 叉積 ab 這里OpenFOAM重載了“” 模 mag(a) 模的平方 magSqr(a) (3) test/ OpenFOAM給了一個(gè)該類(lèi)的test，在applicationstestvector；該test并沒(méi)有將向量的所有功進(jìn)行測(cè)試#include vector.H /向量頭文件，別將其寫(xiě)成vector.h了，不然就變成了c

3、+的了#include IOstreams.H /輸入using namespace Foam;int main() Info vector:zero endl /打印zero （0 0 0） vector:one endl /打印 one （1 1 1） vector:dim endl /打印維數(shù),返回3 vector:rank endl; /打印階數(shù)，返回1， 標(biāo)量：0階，向量 1階，張量2階 vector d(0.5, 0.5, 0.5); /聲明一個(gè)向量 d /= mag(d); /求向量方向。 /下面主要測(cè)試小向量的誤差。 vector dSmall = (1e-100)*d; dSmall /= mag(dSmall); Info (dSmall - d) endl; /向量和標(biāo)量相乘，各分量都相乘。 d *= 4.0;/輸出d Info d endl; /測(cè)試求和 Info d + d endl; /模的平方 Info magSqr(d) endl; / vector d2(0.5, 0.51, -0.5); /生命另外一個(gè)變量，注意個(gè)數(shù)的差別，看看輸出結(jié)果是啥吧。 Info cmptMax(d2) cmptSum(d2) cmptMag(d2) endl; Info min(d, d2) endl; return 0;上面列出的其他功能自己向里面寫(xiě)吧，如果有