求改进：计算曲线选择集交点

Love-Lisp · 发表于 2013-8-20 14:23:46

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？立即注册

×

初学arx,很多方面知识欠缺！自己尝试写了个曲线选择集计算交点的函数，代码如下：

void CurveInters(ads_name ss,AcGePoint3dArray & pPoints)
{
long len,total;
acedSSLength(ss,&len);
if (len < 2)
{
return ;
}
AcDbVoidPtrArray entArray;
ads_name ent;
AcDbObjectId entId;
for (int i = 0;i < len; i++)
{
acedSSName(ss, i, ent);
acdbGetObjectId(entId,ent);
AcDbEntity *pEnt;
acdbOpenObject(pEnt,entId, AcDb::kForRead);
entArray.append(pEnt);
}
AcDbEntity *pEnt1,*pEnt2;
for (int i = 0;i < (len - 1);i++)
{
pEnt1 = AcDbEntity::cast((AcRxObject*)entArray);
for (int j = i + 1;j < len;j++)
{
pEnt2 = AcDbEntity::cast((AcRxObject*)entArray[j]);
AcGePoint3dArray pPtArray;
pEnt1->intersectWith(pEnt2,AcDb::kOnBothOperands,pPtArray);
pPoints.append(pPtArray);
}
}
for (int i = 0;i < len;i++)
{
pEnt1 = AcDbEntity::cast((AcRxObject*)entArray);
pEnt1->close();
}
return ;
}

void CurveInters(ads_name ss,AcGePoint3dArray & pPoints) 
{ 
  long len,total; 
  acedSSLength(ss,&len); 
  if (len < 2) 
  { 
    return ; 
  } 
   AcDbVoidPtrArray entArray; 
  ads_name ent; 
  AcDbObjectId entId; 
  for (int i = 0;i < len; i++) 
  { 
    acedSSName(ss, i, ent); 
    acdbGetObjectId(entId,ent); 
    AcDbEntity *pEnt; 
    acdbOpenObject(pEnt,entId, AcDb::kForRead); 
    entArray.append(pEnt); 
  } 
  AcDbEntity *pEnt1,*pEnt2; 
  for (int i = 0;i < (len - 1);i++) 
  { 
    pEnt1 = AcDbEntity::cast((AcRxObject*)entArray); 
    for (int j = i + 1;j < len;j++) 
    { 
 
      pEnt2 = AcDbEntity::cast((AcRxObject*)entArray[j]); 
      AcGePoint3dArray pPtArray; 
      pEnt1->intersectWith(pEnt2,AcDb::kOnBothOperands,pPtArray); 
      pPoints.append(pPtArray); 
 
    } 
 
  } 
  for (int i = 0;i < len;i++) 
  { 
    pEnt1 = AcDbEntity::cast((AcRxObject*)entArray); 
    pEnt1->close(); 
  } 
 
  return ; 
}

问题：该函数计算4000个实体交点的用时是晓东API函数 XDRX_GETSSINTERS 用时的3倍！
我相信，肯定是我的代码有问题，望高手们指点指点！不胜感激！

Free-Lancer · 发表于 2013-8-20 14:38:28

不会ARX，只是听说过有AcDb和AcGe区别，用AcDb和ActivX的 IntersectWith 类似吧，效率高不上去

Love-Lisp · 发表于 2013-8-20 14:43:00

Free-Lancer 发表于 2013-8-20 14:38
不会ARX，只是听说过有AcDb和AcGe区别，用AcDb和ActivX的 IntersectWith 类似吧，效率高不上去

还没开始学习几何类的用法！

Free-Lancer · 发表于 2013-8-20 14:49:57

Love-Lisp 发表于 2013-8-20 14:43
还没开始学习几何类的用法！

呵呵，借用LoveARX一句话“不会AcGe库就是不会ARX”，努力啊

XDSoft · 发表于 2013-8-20 15:17:20

楼主，确信用的是 xdrx_getSSInters 还是 xdrx_curve_getinters ?

如果是xdrx_getSSInters ，你用的什么时候版本的API？现在新的API已经没这函数了。

Love-Lisp · 发表于 2013-8-20 15:37:05

XDSoft 发表于 2013-8-20 15:17
楼主，确信用的是 xdrx_getSSInters 还是 xdrx_curve_getinters ?

如果是xdrx_getSSInters ，你用的什么 ...

之前测试的xdrx_getSSInters 函数是8月5号未更新的版本！
更新8月20号的API,用xdrx_curve_getinters计算交点，好像效率低了很多啊！不过4000个实体求交点还是比我定义的arx快近1倍！
老大给解释解释原因啊！

XDSoft · 发表于 2013-8-20 15:40:10

10多年前写的代码了，和你的比较，

没把选择集先放到AcDbVoidPtrArray里面，直接选择集提出实体，你可能这个地方慢了,这个地方先放到数组里面多余了。

Adesk::Boolean RetIntersSec(ads_name ss ,AcGePoint3dArray& nPts,int IntExtMode)
{
long len;
ads_name nEname1,nEname2;
AcDbObjectId id1,id2;
AcDbCurve *pCurve1,*pCurve2;;
AcDbEntity *pEnt1,*pEnt2;
Acad::ErrorStatus es;
if ((ads_sslength(ss,&len)!=RTNORM)&&(len<2))
return Adesk::kFalse;
if (!(nPts.isEmpty())){
nPts.setPhysicalLength(0);
nPts.setLogicalLength(0);
}
for (int i=0;i<len;i++){
ads_ssname(ss,i,nEname1);
acdbGetObjectId(id1,nEname1);
es=acdbOpenAcDbEntity(pEnt1,id1,AcDb::kForRead);
if (es!=Acad::eOk)
continue;
if ((pCurve1=AcDbCurve::cast(pEnt1))==NULL){
pEnt1->close();
continue;
}
for (int j=i+1;j<len;j++){
ads_ssname(ss,j,nEname2);
acdbGetObjectId(id2,nEname2);
es=acdbOpenAcDbEntity(pEnt2,id2,AcDb::kForRead);
if (es!=Acad::eOk)
continue;
if ((pCurve2=AcDbCurve::cast(pEnt2))==NULL){
pEnt2->close();
continue;
}
if (IntExtMode==1){
pEnt1->intersectWith(pEnt2,AcDb::kExtendBoth,nPts);
}
else if (IntExtMode==0){
pEnt1->intersectWith(pEnt2,AcDb::kOnBothOperands,nPts);
}
pEnt2->close();
}
pEnt1->close();
}
return ((nPts.isEmpty())?Adesk::kFalse:Adesk::kTrue);
}

Adesk::Boolean RetIntersSec(ads_name ss ,AcGePoint3dArray& nPts,int IntExtMode) 
{ 
        long len; 
        ads_name nEname1,nEname2; 
        AcDbObjectId id1,id2; 
        AcDbCurve *pCurve1,*pCurve2;; 
        AcDbEntity *pEnt1,*pEnt2; 
        Acad::ErrorStatus es; 
        if ((ads_sslength(ss,&len)!=RTNORM)&&(len<2)) 
                return Adesk::kFalse; 
        if (!(nPts.isEmpty())){ 
                nPts.setPhysicalLength(0); 
                nPts.setLogicalLength(0); 
        } 
        for (int i=0;i<len;i++){ 
                ads_ssname(ss,i,nEname1); 
                acdbGetObjectId(id1,nEname1); 
                es=acdbOpenAcDbEntity(pEnt1,id1,AcDb::kForRead); 
                if (es!=Acad::eOk) 
                        continue; 
                if ((pCurve1=AcDbCurve::cast(pEnt1))==NULL){ 
                        pEnt1->close(); 
                        continue; 
                } 
                for (int j=i+1;j<len;j++){ 
                        ads_ssname(ss,j,nEname2); 
                        acdbGetObjectId(id2,nEname2); 
                        es=acdbOpenAcDbEntity(pEnt2,id2,AcDb::kForRead); 
                        if (es!=Acad::eOk) 
                                continue; 
                        if ((pCurve2=AcDbCurve::cast(pEnt2))==NULL){ 
                                pEnt2->close(); 
                                continue; 
                        } 
                        if (IntExtMode==1){ 
                                pEnt1->intersectWith(pEnt2,AcDb::kExtendBoth,nPts); 
                        } 
                        else if (IntExtMode==0){ 
                                pEnt1->intersectWith(pEnt2,AcDb::kOnBothOperands,nPts); 
                        } 
                        pEnt2->close(); 
                } 
                pEnt1->close(); 
        } 
        return ((nPts.isEmpty())?Adesk::kFalse:Adesk::kTrue); 
}

XDSoft · 发表于 2013-8-20 15:42:28

Love-Lisp 发表于 2013-8-20 15:37
之前测试的xdrx_getSSInters 函数是8月5号未更新的版本！
更新8月20号的API,用xdrx_curve_getinters计算 ...

上面的代码是 xdrx_getSSInters实现的核心部分，你放你的代码里面看看，速度提高不？

Love-Lisp · 发表于 2013-8-20 16:39:40

本帖最后由 Love-Lisp 于 2013-8-20 16:40 编辑

XDSoft 发表于 2013-8-20 15:42
上面的代码是 xdrx_getSSInters实现的核心部分，你放你的代码里面看看，速度提高不？

通过看你的代码，我找到慢的原因了！问题在这里：          AcGePoint3dArray pPtArray;
      pEnt1->intersectWith(pEnt2,AcDb::kOnBothOperands,pPtArray);
      pPoints.append(pPtArray);
      改成这样即可：
      pEnt1->intersectWith(pEnt2,AcDb::kOnBothOperands,pPoints);

   还是对指针理解不透彻！
   现在修改后的4000个物体的计算速度已经比xdrx_curve_getinters函数要快30%！

st788796 · 发表于 2013-8-20 16:49:48

Love-Lisp 发表于 9 分钟前

 登录/注册后可看大图

通过看你的代码，我找到慢的原因了！问题在这里： AcGePoint3dArray pPtArray;
pEnt1->intersectWith(pEnt2

就是说xdrx_curve_getinters还有优化空间？来自: Android客户端

XDSoft · 发表于 2013-8-20 18:10:59

Love-Lisp 发表于 2013-8-20 16:39
通过看你的代码，我找到慢的原因了！问题在这里： AcGePoint3dArray pPtArray;
pE ...

你的函数传递进来的是引用，所以直接往里面赋值就传回去了，你慢的原因是在循环里面，创建数组对象，没循环次都要创建，析构。

另外，修改后，你比较了和 xdrx_getSSInters 的快慢了吗？原先你说慢3倍，现在比curve_getinters快30%，好像还是比xdrx_getSSInters慢啊。

Love-Lisp · 发表于 2013-8-20 19:11:53

本帖最后由 Love-Lisp 于 2013-8-20 19:18 编辑

XDSoft 发表于 2013-8-20 18:10
你的函数传递进来的是引用，所以直接往里面赋值就传回去了，你慢的原因是在循环里面，创建数组对象，没循 ...

经测试，4000个实体，我修改后的函数运行效率要比xdrx_getSSInters高！计算出的交点个数同为64419个！
xdrx_getSSInters :
用时 94.646 秒

CURVEINTERS :
用时 87.158 秒
我想7楼给出的代码里，每次计算交点，实体都要先打开，再关闭，这样打开关闭的次数是 n+ (n-1) + ... 1次，我的代码是先将所有的对象一次先打开，全部交点计算完成后再一次性关闭，是不是这儿也会有不少时间差呢？

xdrx_curve_getinters函数效率比xdrx_getSSInters低，是不是考虑了容差导致的？

Love-Lisp · 发表于 2013-8-20 19:32:35

同样测试4000条直线，交点个数为12800，
xdrx_getSSInters :
用时 1.67 秒
CURVEINTERS :
用时 1.248 秒
要是实体包含其他各种曲线，正如12楼的用时，和纯直线比起来，计算速度差距也太大了！

XDSoft · 发表于 2013-8-20 19:41:26

Love-Lisp 发表于 2013-8-20 19:11
经测试，4000个实体，我修改后的函数运行效率要比xdrx_getSSInters高！计算出的交点个数同为64419个！
x ...

看代码你的就应该快，我的程序里面加了曲线的判断，也加了两种扩展方式的判断，在大量循环的时候肯定要浪费不少时间，一次性都打开提高不了太多的性能，以前测试过，打开对象写表格不关闭，和打开关闭时间差不多多少。

不推荐你一次性都打开对象，然后等着操作，万一有故障什么的程序不正常退出了呢，那些对象都打开，CAD就要出错了。

Love-Lisp · 发表于 2013-8-20 20:24:22

XDSoft 发表于 2013-8-20 19:41
看代码你的就应该快，我的程序里面加了曲线的判断，也加了两种扩展方式的判断，在大量循环的时候肯定要浪 ...

明白了！一次性打开万万使不得！程序中途意外中断就麻烦了！修改后的代码用时基本和xdrx_getssinters一样了！

void SSInters(ads_name ss,AcGePoint3dArray & pPoints)
{
Acad::ErrorStatus es;
long len;
acedSSLength(ss,&len);
if (len < 2)
{
return ;
}
ads_name e1,e2;
AcDbEntity *pEnt1,*pEnt2;
AcDbObjectId Id1,Id2;
for (int i = 0;i < (len - 1);i++)
{
acedSSName(ss, i, e1);
acdbGetObjectId(Id1,e1);
es = acdbOpenObject(pEnt1,Id1, AcDb::kForRead);
if (es != Acad::eOk)
{
pEnt1->close();
continue;
}/*如果只计算曲线交点，请使用如下语句
if (AcDbCurve::cast(pEnt1) == NULL)
{
pEnt1->close();
continue;
}*/
for (int j = i + 1;j < len;j++)
{
acedSSName(ss, j, e2);
acdbGetObjectId(Id2,e2);
es = acdbOpenObject(pEnt2,Id2, AcDb::kForRead);
if (es != Acad::eOk)
{
pEnt2->close();
continue;
}/*如果只计算曲线交点，请使用如下语句
if (AcDbCurve::cast(pEnt2) == NULL)
{
pEnt2->close();
continue;
}*/
pEnt1->intersectWith(pEnt2,AcDb::kOnBothOperands,pPoints);
pEnt2->close();
}
pEnt1->close();
}
return ;
}

void SSInters(ads_name ss,AcGePoint3dArray & pPoints) 
{ 
  Acad::ErrorStatus es; 
  long len; 
  acedSSLength(ss,&len); 
  if (len < 2) 
  { 
    return ; 
  } 
 
  ads_name e1,e2; 
  AcDbEntity *pEnt1,*pEnt2; 
  AcDbObjectId Id1,Id2; 
  for (int i = 0;i < (len - 1);i++) 
  { 
    acedSSName(ss, i, e1); 
    acdbGetObjectId(Id1,e1); 
        es = acdbOpenObject(pEnt1,Id1, AcDb::kForRead); 
    if (es != Acad::eOk) 
    { 
      pEnt1->close(); 
      continue; 
    }/*如果只计算曲线交点，请使用如下语句 
    if (AcDbCurve::cast(pEnt1) == NULL) 
    { 
      pEnt1->close(); 
      continue; 
    }*/ 
    for (int j = i + 1;j < len;j++) 
    { 
      acedSSName(ss, j, e2); 
      acdbGetObjectId(Id2,e2); 
      es = acdbOpenObject(pEnt2,Id2, AcDb::kForRead); 
      if (es != Acad::eOk) 
      { 
        pEnt2->close(); 
        continue; 
      }/*如果只计算曲线交点，请使用如下语句 
      if (AcDbCurve::cast(pEnt2) == NULL) 
      { 
        pEnt2->close(); 
        continue; 
      }*/ 
 
      pEnt1->intersectWith(pEnt2,AcDb::kOnBothOperands,pPoints); 
      pEnt2->close(); 
 
    } 
    pEnt1->close(); 
 
  } 
 
  return ; 
}

账号		自动登录	找回密码
密码			立即注册

[求助] 求改进：计算曲线选择集交点

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

点评

点评

点评

点评

点评

点评

点评

点评

点评