#include <iostream>
#include <gp_Ax2.hxx>
#include <TopoDS_Shape.hxx>
#include <BRepPrimAPI_MakeTorus.hxx>
#include <StlAPI_Writer.hxx>
#include <BRepMesh_IncrementalMesh.hxx>

void Generate_Shape(TopoDS_Shape * store)
{
  gp_Ax2 axis = gp_Ax2(gp_Pnt(0,0,0), gp_Dir(0,0,1));
  Standard_Real r1 = 1.0;
  Standard_Real r2 = 0.1;
  Standard_Real a1 = -1.5708;
  Standard_Real a2 = 1.5708;
  Standard_Real angle = 2.094;
  TopoDS_Shape shape = BRepPrimAPI_MakeTorus(axis, r1, r2, a1, a2, angle).Shape();

  *store = shape;
}

void Write_STL(TopoDS_Shape * store)
{
  StlAPI_Writer stl_writer;
  BRepMesh_IncrementalMesh Mesh( *store, 0.01 );
  Mesh.Perform();

  if (stl_writer.Write(*store, "demo1.stl")) {
    std::cout << std::endl << "True ";
  } else {
    std::cout << std::endl << "False ";
  }
}

int main()
{
  TopoDS_Shape *store = new TopoDS_Shape();
  Generate_Shape(store);
  Write_STL(store);

  return 0;
}