Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. 自定义类中调用PrimitivePatch类中的faceCentre的bug

自定义类中调用PrimitivePatch类中的faceCentre的bug

已定时 已固定 已锁定 已移动 OpenFOAM
2 帖子 2 发布者 1.1k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • L 离线
    L 离线
    lmessi1011
    写于 最后由 编辑
    #1

    各位大神,我最近在写一个类的时候遇到一个问题,想请假一下。首先我在这个spanwiseAverage类的patchAddressing方法中需要用到某一patch的点坐标和面心坐标,我分别通过patch的localPoints()方法和facecentres()方法获取了两个pointField的引用。随后我需要通过点坐标来计算一个变量geoMin,问题就出在了此处。如果我在这个句代码前面通过Info输出一下patchCentre的大小或者是值是可以得到正确的结果的,后面的程序也可以正常的运行。但是如果不在计算geoMin变量之前对patchCentre进行操作,在计算完geoMin后再输出patchCentre的大小就会变为0。我在计算geoMin的时候并没有对patchCentre这个变量进行任何改变。想请问一下,有没有大神遇到过类似的问题或者这个bug产生的原因。

    const Foam::labelList& Foam::spanwiseAverage::patchAddressing(const polyPatch& patch) const
    {
        const edgeList& patchEdge = patch.edges();
        const pointField& patchPoint = patch.localPoints();
        const pointField& patchCentre = patch.faceCentres();
    
        Info << patchCentre.size() << endl;  // 若在此处输出patchCentre的size是可以获取正确的值的
    
        scalar geoMin = min(patchPoint & averagePatchsDir_);
    
        Info << patchCentre.size() << endl;  // 若注释掉上一条Info语句,在这里输出patchCentre的size得到的值变为了0
    
        ......
    }
    
    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #2
    Info << patchCentre.size() << endl;  // 若注释掉上一条Info语句,在这里输出patchCentre的size得到的值变为了0
    
    

    换成这个呢?

    Info << patch.faceCentres();
    

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]