| | | | |

2017年2月4日

GFORTRAN在64位時候設置線程數目使用OMP_SET_NUM_THREADS無法生效，查看glibc代碼，可以看出OMP_SET_NUM_THREADS是用C語言實現的，fortran參數傳過來默認為地址，而OMP_SET_NUM_THREADS把傳過來地址默認處理為32位地址，所以調用失敗，64位下要使用專門處理64位地址的對應函數OMP_SET_NUM_THREADS_8_，上述兩個函數輸入參數不需要確認輸入的integer是32位還是64位，這個與Intel Fortran不一樣，Intel Fortran只是需要確認OMP_SET_NUM_THREADS的參數是32位還是64位

posted @ 2017-02-04 10:59 隔壁老劉閱讀(632) | 評論 (0) | 編輯收藏

2016年4月12日

C/C++位運算的陷阱

昨天用C/C++寫Base64編碼，結果老是不對，
比如

unsigned char A =;

unsigned char B = ;

unsigned char C = ;
unsigned char a, b, c, d;

a = A >> 2;

b=(A << 6) >> 2|(B>>4);

c = ((B << 4)) >> 2) | (C >> 6);

d = (C << 2) >> 2;
結果發現bcd值都不對
但是寫成

b = A<< 6;

b = b >> 2;
b=b|B>>4;
就可以，查看b計算兩種方式匯編
第一種

00BACB66 movzx eax,byte ptr [ebp-1F9h]

00BACB6D shl eax,6

00BACB70 sar eax,2

00BACB73 mov byte ptr [ebp-1C9h],al

第二種

00BACB79 movzx eax,byte ptr [ebp-1F9h]

00BACB80 shl eax,6

00BACB83 mov byte ptr [ebp-1EDh],al

00BACB89 movzx eax,byte ptr [ebp-1EDh]

00BACB90 sar eax,2

00BACB93 mov byte ptr [ebp-1EDh],al
問題就在指令sar，這是算數右移位指令，會保留最高位，而第二種可以是因為movzx 將高位已經置零
找到了原因，可以寫為

a = A >> 2;

b=((unsigned char)(A << 6)) >> 2|(B>>4);

c = (((unsigned char)(B << 4)) >> 2) | (C >> 6);

d = ((unsigned char)(C << 2)) >> 2;
只要涉及到右移，保證操作數是無符號就可以了。

posted @ 2016-04-12 10:41 隔壁老劉閱讀(406) | 評論 (2) | 編輯收藏

2016年2月27日

Intel OpenMP運行時錯誤 OMP: Error #133 解決辦法

OMP: Error #133: Inconsistent THREADPRIVATE common block declarations are non-conforming and are unsupported. Either all threadprivate common blocks must be declared identically, or the largest instance of each threadprivate common block must be referenced first during the run.

發生以上錯誤確實很麻煩，而且可參考文檔基本沒有。Fortran代碼中的Common塊數據大小可以不一樣大（相同名稱的Common），但是在OpenMP中確要求，要么大小一樣，要么占用內存最大的先調用，變量如果是隱式聲明，一定要注意，變量名稱不同，可能會編譯為不同長度的內存

posted @ 2016-02-27 12:13 隔壁老劉閱讀(955) | 評論 (0) | 編輯收藏

2015年11月12日

OpenCASCADE 6.9.1 Bug

OpenCASCADE 6.9.1（估計從6.9.0開始），使用新的屏幕選擇算法，我使用6.8開發的程序，升級到6.9.1之后，使用IVTK部分模型無法選中幾何元素，經過調試發現6.9.1代碼有BUG，導致數據丟失，無法判斷選中的元素是否在選擇范圍，BUG出現的文件為SelectMgr_RectangularFrustum.cxx中的函數ScaleAndTransform

NCollection_Handle<SelectMgr_BaseFrustum> SelectMgr_RectangularFrustum::ScaleAndTransform (const Standard_Integer theScaleFactor,
const gp_Trsf& theTrsf)
{
Standard_ASSERT_RAISE (theScaleFactor > 0,
"Error! Pixel tolerance for selection should be greater than zero");

SelectMgr_RectangularFrustum* aRes = new SelectMgr_RectangularFrustum();
const Standard_Boolean isToScale = theScaleFactor != 1;
const Standard_Boolean isToTrsf = theTrsf.Form() != gp_Identity;

if (!isToScale && !isToTrsf)
return aRes ;//這里，當部分模型滿足該條件，返回一個空的數據，使用該數據無法判斷模型是否在區域

aRes->myIsOrthographic = myIsOrthographic;
SelectMgr_RectangularFrustum* aRef = this;

if (isToScale)
{
    aRes->myNearPickedPnt = myNearPickedPnt;
    aRes->myFarPickedPnt = myFarPickedPnt;
    aRes->myViewRayDir    = myViewRayDir;

    const gp_Pnt2d aMinPnt (myMousePos.X() - theScaleFactor * 0.5,
                            myMousePos.Y() - theScaleFactor * 0.5);
    const gp_Pnt2d aMaxPnt (myMousePos.X() + theScaleFactor * 0.5,
                            myMousePos.Y() + theScaleFactor * 0.5);

// recompute base frustum characteristics from scratch
computeFrustum (aMinPnt, aMaxPnt, myBuilder, aRes->myVertices, aRes->myEdgeDirs);

aRef = aRes;
}

if (isToTrsf)
{
    aRes->myNearPickedPnt = aRef->myNearPickedPnt.Transformed (theTrsf);
    aRes->myFarPickedPnt = aRef->myFarPickedPnt.Transformed (theTrsf);
    aRes->myViewRayDir    = aRes->myFarPickedPnt.XYZ() - aRes->myNearPickedPnt.XYZ();

      // LeftTopNear
    aRes->myVertices[0] = aRef->myVertices[0].Transformed (theTrsf);
    // LeftTopFar
    aRes->myVertices[1] = aRef->myVertices[1].Transformed (theTrsf);
    // LeftBottomNear
    aRes->myVertices[2] = aRef->myVertices[2].Transformed (theTrsf);
    // LeftBottomFar
    aRes->myVertices[3] = aRef->myVertices[3].Transformed (theTrsf);
    // RightTopNear
    aRes->myVertices[4] = aRef->myVertices[4].Transformed (theTrsf);
    // RightTopFar
    aRes->myVertices[5] = aRef->myVertices[5].Transformed (theTrsf);
    // RightBottomNear
    aRes->myVertices[6] = aRef->myVertices[6].Transformed (theTrsf);
    // RightBottomFar
    aRes->myVertices[7] = aRef->myVertices[7].Transformed (theTrsf);

    // Horizontal
    aRes->myEdgeDirs[0] = aRes->myVertices[4].XYZ() - aRes->myVertices[0].XYZ();
    // Vertical
    aRes->myEdgeDirs[1] = aRes->myVertices[2].XYZ() - aRes->myVertices[0].XYZ();
    // LeftLower
    aRes->myEdgeDirs[2] = aRes->myVertices[2].XYZ() - aRes->myVertices[3].XYZ();
    // RightLower
    aRes->myEdgeDirs[3] = aRes->myVertices[6].XYZ() - aRes->myVertices[7].XYZ();
    // LeftUpper
    aRes->myEdgeDirs[4] = aRes->myVertices[0].XYZ() - aRes->myVertices[1].XYZ();
    // RightUpper
    aRes->myEdgeDirs[5] = aRes->myVertices[4].XYZ() - aRes->myVertices[5].XYZ();
}

// compute frustum normals
computeNormals (aRes->myEdgeDirs, aRes->myPlanes);

cacheVertexProjections (aRes);

return NCollection_Handle<SelectMgr_BaseFrustum> (aRes);
}

修改辦法：

if (!isToScale && !isToTrsf)
{
aRes->myPixelTolerance=myPixelTolerance;
aRes->myIsOrthographic=myIsOrthographic;

aRes->myBuilder= myBuilder;

   aRes->myIsOrthographic = myIsOrthographic;
   aRes->myNearPickedPnt = myNearPickedPnt;
   aRes->myFarPickedPnt = myFarPickedPnt;
   aRes->myViewRayDir = myViewRayDir;
   aRes->myMousePos = myMousePos;

   for (int i = 0; i <6;++i)
   {
    aRes->myPlanes[i] = myPlanes[i];
   }
   for (int i = 0; i < 6; ++i)
   {
    aRes->myVertices[i] = myVertices[i];
   }
   for (int i = 0; i < 6; ++i)
   {
    aRes->myMaxVertsProjections[i] = myMaxVertsProjections[i];
   }
   for (int i = 0; i < 6; ++i)
   {
    aRes->myMinVertsProjections[i] = myMinVertsProjections[i];
   }
   for (int i = 0; i <3; ++i)
   {
    aRes->myMaxOrthoVertsProjections[i] = myMaxOrthoVertsProjections[i];
   }
   for (int i = 0; i <3; ++i)
   {
    aRes->myMinOrthoVertsProjections[i] = myMinOrthoVertsProjections[i];
   }
   for (int i = 0; i < 6; ++i)
   {
    aRes->myEdgeDirs[i] = myEdgeDirs[i];
   }
   return NCollection_Handle<SelectMgr_BaseFrustum>(aRes);
}

posted @ 2015-11-12 15:20 隔壁老劉閱讀(583) | 評論 (0) | 編輯收藏

2015年11月7日

C++智能指針類用不好害死人

C++智能指針類使用起來方便，但是不維護好引用計數，不能隨時控制對象的釋放，反而造成內存耗盡，這樣說來，反而不能隨時控制對象的生命周期，所以使用要謹慎啊。

posted @ 2015-11-07 21:43 隔壁老劉閱讀(613) | 評論 (0) | 編輯收藏

僅列出標題

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

常用鏈接

留言簿

隨筆檔案

文章檔案

搜索

最新評論

閱讀排行榜

評論排行榜

2017年2月4日

2016年4月12日

2016年2月27日

2015年11月12日

2015年11月7日