0

这是一个cgal相关的问题，但我认为这也是一个通用的C++问题，所以我在这里问它。

我正在尝试使用Alpha_shape_2该类，并将其分配给AlphaShapeCg名为GetAlphaShalCg. 问题是其中的某些函数Alpha_shape_2没有返回正确的结果。

这是我的代码，它真的很简单，但我不太明白为什么Alpha_shape_2在子例程中分配给包装器，然后访问父例程中的成员和Alpha_shape_2直接访问之间存在差异。

如果您安装了CGAL ，这是您可以编译和使用的完整代码。

#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/algorithm.h>
#include <CGAL/Delaunay_triangulation_2.h>
#include <CGAL/Alpha_shape_2.h>

#include <iostream>
#include <fstream>
#include <vector>
#include <list>


typedef CGAL::Exact_predicates_inexact_constructions_kernel K;

typedef K::FT FT;

typedef K::Point_2  Point;
typedef K::Segment_2  Segment;


typedef CGAL::Alpha_shape_vertex_base_2<K> Vb;
typedef CGAL::Alpha_shape_face_base_2<K>  Fb;
typedef CGAL::Triangulation_data_structure_2<Vb,Fb> Tds;
typedef CGAL::Delaunay_triangulation_2<K,Tds> Triangulation_2;

typedef CGAL::Alpha_shape_2<Triangulation_2>  Alpha_shape_2;


template <class OutputIterator>
bool
file_input(OutputIterator out)
{
  std::ifstream is("./data/fin", std::ios::in);

  if(is.fail()){
    std::cerr << "unable to open file for input" << std::endl;
    return false;
  }

  int n;
  is >> n;
  std::cout << "Reading " << n << " points from file" << std::endl;
  CGAL::copy_n(std::istream_iterator<Point>(is), n, out);

  return true;
}

//------------------ main -------------------------------------------


struct AlphaShapeCg
{

    Alpha_shape_2 *AlphaShape;
};

void GetAlphaShalCg(AlphaShapeCg *ashape,  std::list<Point> points)
{

      Alpha_shape_2 A(points.begin(), points.end(),
          FT(100000),
          Alpha_shape_2::GENERAL);
    ashape->AlphaShape=&A;
}



int main()
{
  std::list<Point> points;
  if(! file_input(std::back_inserter(points))){
    return -1;
  }

   AlphaShapeCg ashape;


   GetAlphaShalCg(&ashape, points);

   Alpha_shape_2 *APtrs=(ashape.AlphaShape);
   int alphaEigenValue = APtrs->number_of_alphas(); // gives incorrect result; alphaEigenValue=0

  //Alpha_shape_2 A(points.begin(), points.end(),
  //  FT(100000),
  //  Alpha_shape_2::GENERAL);
  //   int alphaEigenValue = APtrs->number_of_alphas(); // gives correct result; alphaEigenValue!=0

}

更新：我尝试使用

Alpha_shape_2 =new A(points.begin(), points.end(), FT(100000), Alpha_shape_2::GENERAL);

但是由于这个错误，这段代码根本无法编译：

错误 C2513：“CGAL::Alpha_shape_2”：在“=”之前没有声明变量

4

1 回答 1

1

您正在为退出函数时被破坏的局部变量分配一个指针。

如果你想在函数中创建对象并返回它的地址 - 你应该使用动态分配（它，当你完成它时new不要忘记）。delete

于 2011-10-05T07:33:33.307 回答