Step Write Sample

 
 
#ifdef WIN32
#include <windows.h>
#endif  
 
#include "testwriters.h"
//DATAKIT headers needed
#include "datakit.h"
#include "stepw/stepw.hpp"
 
 
//  For tessellation
#include "tess/tess.h"
 
/******************************************/
 
/****
Creation of data : topological data (a Body and its content), FDT, etc
****/
 
#include "SampleWriter/testcreatecylfdt.hpp"
#include "SampleWriter/testcreatefdt.hpp"
 
/*********************************************************/
/****
Writing a Product : a Part or an Assembly
 
Also with example of writing ATTRIBUTES attached at level of the Product , Part or Assembly
In particular, VALIDATION PROPERTIES
****/
 
extern int init_file_to_read( char *filename );
extern void close_file_to_read();
 
/****
Creation of a Sample View
****/
Dtk_ModelDisplayPtr CreateView()
{
     Dtk_ModelDisplayPtr view;
 
     Dtk_string name( "SampleView:1" );
     Dtk_transfo axes;
     Dtk_pnt projpoint( 0., 0., 50. );
     Dtk_pnt location( 0., 0., 0. );
     Dtk_dir AY = axes.getYdir();
     Dtk_dir target = axes.getZdir();
     double scale = 1.;
     double angle = 45.;
     double viewdist = 50.;
     //  creation
     Dtk_Camera::ProjectionTypeEnum ProjType = Dtk_Camera::ProjectionTypeParallel;
     Dtk_CameraPtr Camera = Dtk_Camera::Create( projpoint, location, AY, angle, viewdist, scale, ProjType );
     Dtk_PlaneSurfacePtr section_plane = NULL;
     /*  Section Plane ?
     fscanf(F,"%s",ligne);
     if (ligne[0] == 'S') {  // Section (1 plan)
     Dtk_transfo trsf_plan = Dtk_StepFdt_TestIn_TRSF(F,ligne);
     section_plane = Dtk_PlaneSurface::Create(trsf_plan.getOrigin(),trsf_plan.getZdir(),trsf_plan.getXdir());
 
     fscanf(F,"%s",ligne);
     }
     */
     Dtk_EntityPtr section_ent = Dtk_EntityPtr::DtkDynamicCast( section_plane );
     view = Dtk_ModelDisplay::Create( Camera, section_ent, DTK_TRUE );
     //  Dtk_FdtCapturePtr capture = Dtk_FdtCapture::Create(Camera, section_plane);
     Dtk_InfoPtr& viewinfo = view->info();
     if( viewinfo.IsNULL() ) viewinfo = Dtk_Info::create();
     viewinfo->SetName( name );
 
     //  Links to FDTs : when writing the STEP file
 
     return view;
}
 
/****
Creation of an Axis System (on the cube)
****/
Dtk_AxisSystemPtr CreateAxisSystem()
{
     Dtk_AxisSystemPtr axis = Dtk_AxisSystem::create();
     Dtk_transfo trsf;
     trsf.setOrigin( Dtk_pnt( 100., 100., 100. ) );
     axis->SetMatrix( trsf );
     axis->SetName( "Axis_System" );
     return axis;
}
 
//start_of_mainstepwrite
//This sample show you how to write a body with an unbounded surface and wireframe entity.
// IFDT = 1 to also write FDT, View, etc - 0 none
// IFDT = 2 to use this Part by a FDT in Assembly (needs to declare ANCHOR entities)
 
int WritePart( const Dtk_string& inReferenceName, const int ifdt )
{
     // We construct a body 
     Dtk_BodyPtr body = sampleWriter::CreateCube_2();
 
     //  TEST : adding wireframe
     Dtk_BodyPtr sbody = sampleWriter::CreateCurves_2();
 
     //You init the Part writer 
     int PartID;
     if( stepw_InitProduct( inReferenceName, PartID ) )
          return PartID;
 
     //  Following is required for writing of FDT with geometric links - either inside the Part or external designating anchors
     //  To be called before starting the part
     stepw_SetAnchorProduct( PartID, PartID );
 
     //  IFDT = 2 : declaring some entities as ANCHOR for a FDT in Assembly
     //  This is needed where the Assembly is written in several files
     if( ifdt == 2 )
     {
          int geomid1 = 73;
          Dtk_string anchor1( "73" );
          stepw_AddAnchorItem( geomid1, PartID, anchor1 );
          int geomid2 = 33;
          Dtk_string anchor2( "33" );
          stepw_AddAnchorItem( geomid2, PartID, anchor2 );
     }
 
     // test : properties defined as Metadata (at part level)
     stepw_SetModeProp( 2 );
     Dtk_MetaDataPtr metadata_1 = Dtk_MetaData::CreateMetaData( Dtk_MetaData::TypeConfigurationProperty, "Revision", "test_revision_4" );
 
     Dtk_MetaDataPtr metadata_2 = Dtk_MetaData::CreateMetaData( Dtk_MetaData::TypeProperty, "Cost", "12500 $" );
 
     stepw_Init3DPart( PartID );
     stepw_Add3DPartProperty( metadata_1 );
     stepw_Add3DPartProperty( metadata_2 );
 
     //  test : properties at part level - on the product_definition_shape
     // test : validation properties - AP242 form (no intermediate name)
     Dtk_InfoPtr part_info = Dtk_Info::create();
     Dtk_Val valprop_area( "1234 mm2" );
     part_info->AddAttribute( "geometric validation property :  : surface area measure", valprop_area );
     Dtk_Val valprop_volume( "5678 mm3" );
     part_info->AddAttribute( "geometric validation property :  : volume measure", valprop_volume );
     Dtk_Val valprop_CG( Dtk_pnt( 12., 34., 56. ) );
     part_info->AddAttribute( "geometric validation property :  : centre point", valprop_CG );
     stepw_SetPartProperties( part_info, 1 );
     // test : properties at part level - on the product_definition
     Dtk_InfoPtr prod_info = Dtk_Info::create();
     Dtk_Val mat_name( "Steel YC38" );
     prod_info->AddAttribute( "material property : material name : Steel", mat_name );
     Dtk_Val mat_density( "7.89 g/cm3" );
     prod_info->AddAttribute( "material property : density : density measure", mat_density );
     stepw_SetPartProperties( prod_info, 2 );
 
     // And we write body in the step file
     stepw_Write3DPartBody( body );
 
     stepw_Write3DPartBody( sbody );  // test wire frame
 
                                              //  IFDT = 1 : writing a FDT within the Part (here, a Datum)
     if( ifdt == 1 )
     {
 
          //  Test Write FDT
          Dtk_FdtPtr Fdt1 = sampleWriter::CreateFdtDatum();
          int nodeid_Fdt1 = 1001;
          stepw_InitNodeContext( nodeid_Fdt1 );
          stepw_Add3DPartFDT( Fdt1 );
          // Then write the geometrical links : within the same Part
          int geomid1 = 73;
          stepw_ER ElementReference_1;
          stepw_CreateReference( ElementReference_1, geomid1 );
          stepw_AddReference( ElementReference_1 );
          stepw_EndNodeContext();
 
          //  Test Write View
          Dtk_ModelDisplayPtr View1 = CreateView();
          int nodeid_View1 = 2001;
          stepw_InitNodeContext( nodeid_View1 );
          stepw_Add3DModelDisplay( View1, 0 );
          stepw_ER ElementReference_Fdt1;
          stepw_CreateReferenceToNode( ElementReference_Fdt1, nodeid_Fdt1, PartID );
          stepw_AddReference( ElementReference_Fdt1 );
          stepw_EndNodeContext();
 
          //  Test Construction Geometries & Axis System
          Dtk_AxisSystemPtr AxisSystem = CreateAxisSystem();
          stepw_Write3DAxisSystem( AxisSystem );
          Dtk_BodyPtr constr_plane = sampleWriter::CreateConstructionPlane();
          stepw_Write3DConstructionGeometry( constr_plane );
 
     }
 
     // Close the current part
     stepw_End3DPart();
 
     stepw_EndProduct( PartID );
 
     return PartID;
}
 
//  Another Part, to test other data - with colors and tessellation
// itessel = 0 or 1 : body only
// itessel = 2 : tessellation only
// itessel = 3 : body + tessellation
 
int WritePartCyl( const Dtk_string& inReferenceName, const int itessel )
{
     // We construct a body : a cylinder
     Dtk_BodyPtr body = sampleWriter::CreateCylinder();
 
 
     //You init the Part writer 
     int PartID;
     if( stepw_InitProduct( inReferenceName, PartID ) )
          return PartID;
 
     //  Following is required for writing of FDT with geometric links
     //  To be called before starting the part
     stepw_SetAnchorProduct( PartID, PartID );
 
     stepw_Init3DPart( PartID );
 
     //  test : properties at part level - on the product_definition_shape
     // test : validation properties - AP242 form (no intermediate name)
     Dtk_InfoPtr part_info = Dtk_Info::create();
     Dtk_Val valprop_area( "1234 mm2" );
     part_info->AddAttribute( "geometric validation property :  : surface area measure", valprop_area );
     Dtk_Val valprop_volume( "5678 mm3" );
     part_info->AddAttribute( "geometric validation property :  : volume measure", valprop_volume );
     Dtk_Val valprop_CG( Dtk_pnt( 12., 34., 56. ) );
     part_info->AddAttribute( "geometric validation property :  : centre point", valprop_CG );
     stepw_SetPartProperties( part_info, 1 );
     // test : properties at part level - on the product_definition
     Dtk_InfoPtr prod_info = Dtk_Info::create();
 
     stepw_SetPartProperties( prod_info, 2 );
     //  Standard use : write a Body : the cylinder
     stepw_Write3DPartBody( body );
     /*
     // Test of Tessellation - NOT RECOMMENDED ON VERSIONS STEP BEFORE  AP242
     //  To get the functions of tessellation : tess/tess.h
     Dtk_MeshPtr mesh;
     Dtk_tab<Dtk_MeshPtr> listMesh;
     Dtk_tab<Dtk_Int32> IsSolid;
     Dtk_bool TessWireframe=false;
     tess_InitTesselation("../../SampleFiles/dtk/",0.2);
     tess_BodyToMeshes (body,listMesh,IsSolid,TessWireframe);
     int nbmesh = listMesh.size();
     if (nbmesh >= 1)
     mesh = listMesh[0];
     // Writing Tessellation
     //  Option 1 : Tessellation "alone" i.e. not linked to the body
     // we write body in the step file, then the tessellation separately
     if (itessel <= 1) stepw_Write3DPartBody (body);
     if (itessel == 2) stepw_Write3DPartMesh (mesh);
     if (itessel == 3) stepw_Write3DPartBodyWithMesh(body,mesh);
     */
     // Close the current part
     stepw_End3DPart();
 
     stepw_EndProduct( PartID );
 
     return PartID;
}
 
/*
Variants of writing an assembly
All write an assembly the same way, but the variants focus on a specific capability (not all together for more clarity)
WriteAsm is the basic variant, completed with validation properties and user attributes on the assembly or items of it
*/
 
void WriteAsm()
{
     int status = 0;
     Dtk_transfo TransformationMatrix0, TransformationMatrix1, TransformationMatrix2;
     Dtk_pnt orig1;  orig1[ 0 ] = orig1[ 1 ] = 0.;  orig1[ 2 ] = 200.;
     Dtk_pnt orig2;  orig2[ 0 ] = orig2[ 2 ] = 0.;  orig2[ 1 ] = 300.;
     TransformationMatrix1.setOrigin( orig1 );
     TransformationMatrix2.setOrigin( orig2 );
 
     //The ProcessPart is an user implemented function to create a part
     int PartID = WritePart( L"essaiPart", 0 );
 
     int ProdID, ProdID2;
 
     //we init a product context in STEP : sub-assembly ...
     status = stepw_InitProduct( L"SubProduct", ProdID2 );
 
     //...to insert a Part instance of "essaiPart"
     status = stepw_AddInstance( ProdID2, PartID, TransformationMatrix0, "PartInstance" );
 
     // test : validation properties - AP242 form (no intermediate name)
     Dtk_InfoPtr part_info = Dtk_Info::create();
     stepw_SetPartProperties( part_info, 1 );
     Dtk_Val valprop_notional_CG( Dtk_pnt( 10., 20., 30. ) );
     part_info->AddAttribute( "assembly validation property :  : centre point", valprop_notional_CG );
     stepw_SetPartProperties( part_info, 1 );
     Dtk_InfoPtr prod_info = Dtk_Info::create();
     Dtk_Val valprop_nbchildren( 3 );
     prod_info->AddAttribute( "assembly validation property :  : number of children", valprop_nbchildren );
     stepw_SetPartProperties( prod_info, 2 );
 
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID2 );
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID2 );
 
 
     //we create another assembly node (main assembly)
     status = stepw_InitProduct( L"RootProduct", ProdID );
     //and we insert another time the a part instance - multinstancing management -
     // test : instance with property
     Dtk_InfoPtr instance_info = Dtk_Info::create();
     Dtk_Val instance_kind( "first" );
     instance_info->AddAttribute( "user defined attribute : : instance kind", instance_kind );
     Dtk_Val instance_valprop( 1 );
     instance_info->AddAttribute( "attribute validation property :  : text user attributes", instance_valprop );
     //   status = stepw_AddInstance(ProdID, PartID, TransformationMatrix0,"PartInstance");
     status = stepw_AddInstanceWithInfo( ProdID, PartID, TransformationMatrix0, "PartInstance", instance_info );
 
     //and two instances to the first sub-assembly
     status = stepw_AddInstance( ProdID, ProdID2, TransformationMatrix1, "PartInstance1" );
     status = stepw_AddInstance( ProdID, ProdID2, TransformationMatrix2, "PartInstance2" );
 
     //finally we close the second assembly (root) ...
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID );
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID );
 
}
 
/****
Variant to write the different products in several files
Similar to WriteAsm (just duplicated), with differences of external files : calls to stepw_AddExternalReference
level = 0 : monolithic file (so, like WriteAsm)
level = 1 : basic external references, only the leaf parts are written separately
level = 2 : nested external references, each component is written separately
Warning : these are the only allowed combinations, others may be not or bad supported by readers
 
asm_only = 0 : complete mode, all components (assemblies and leaf part) are written - monolithic writing is always complete
asm_only = 1 : only assemblies are written, part are declared as external products but not written
****/
void WriteAsmExt( const int level, const int asm_only )
{
     int status = 0;
     Dtk_transfo TransformationMatrix0, TransformationMatrix1, TransformationMatrix2;
     Dtk_pnt orig1;  orig1[ 0 ] = orig1[ 1 ] = 0.;  orig1[ 2 ] = 200.;
     Dtk_pnt orig2;  orig2[ 0 ] = orig2[ 2 ] = 0.;  orig2[ 1 ] = 300.;
     TransformationMatrix1.setOrigin( orig1 );
     TransformationMatrix2.setOrigin( orig2 );
 
     //The ProcessPart is an user implemented function to create a part
     //  Option to write the Part as external file
     int PartID;
     if( level ) stepw_AddExternalReference( L"essaiPart", L"RootPart.step", PartID );
 
     //with external references, it is possible to declare parts as external products but not write them
     //  Each product is identified by its NAME - there is a unique equivalence between Product Name and returned vale of ProdID
     if( !asm_only || level == 0 )
          PartID = WritePart( L"essaiPart", 0 );
 
     int ProdID, ProdID2;
 
     //we init a product context in STEP : sub-assembly ...
     //  Option to write the sub-assembly as external file
     if( level > 1 ) stepw_AddExternalReference( L"SubProduct", L"RootSubAssy.step", ProdID2 );
     status = stepw_InitProduct( L"SubProduct", ProdID2 );
 
     //...to insert a Part instance of "essaiPart"
     status = stepw_AddInstance( ProdID2, PartID, TransformationMatrix0, "PartInstance" );
 
     // test : validation properties - AP242 form (no intermediate name)
     Dtk_InfoPtr part_info = Dtk_Info::create();
     stepw_SetPartProperties( part_info, 1 );
     Dtk_Val valprop_notional_CG( Dtk_pnt( 10., 20., 30. ) );
     part_info->AddAttribute( "assembly validation property :  : centre point", valprop_notional_CG );
     stepw_SetPartProperties( part_info, 1 );
     Dtk_InfoPtr prod_info = Dtk_Info::create();
     Dtk_Val valprop_nbchildren( 3 );
     prod_info->AddAttribute( "assembly validation property :  : number of children", valprop_nbchildren );
     stepw_SetPartProperties( prod_info, 2 );
 
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID2 );
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID2 );
 
 
     //we create another assembly node (main assembly)
     status = stepw_InitProduct( L"RootProduct", ProdID );
     //and we insert another time the a part instance - multinstancing management -
     // test : instance with property
     Dtk_InfoPtr instance_info = Dtk_Info::create();
     Dtk_Val instance_kind( "first" );
     instance_info->AddAttribute( "user defined attribute : : instance kind", instance_kind );
     Dtk_Val instance_valprop( 1 );
     instance_info->AddAttribute( "attribute validation property : : text user attributes", instance_valprop );
     //   status = stepw_AddInstance(ProdID, PartID, TransformationMatrix0,"PartInstance");
     status = stepw_AddInstanceWithInfo( ProdID, PartID, TransformationMatrix0, "PartInstance", instance_info );
 
     //and two instances to the first sub-assembly
     status = stepw_AddInstance( ProdID, ProdID2, TransformationMatrix1, "PartInstance1" );
     status = stepw_AddInstance( ProdID, ProdID2, TransformationMatrix2, "PartInstance2" );
 
     //finally we close the second assembly (root) ...
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID );
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID );
 
}
 
 
/****
Variant to write an assembly with a FDT at assembly level - monolithic or with external references
Similar to WriteAsmExt (just duplicated),including external files, but with a FDT in more, and validation properties in less
level = 0 : monolithic file (so, like WriteAsm)
level = 1 : basic external references, only the leaf parts are written separately
level = 2 : nested external references, each component is written separately
Warning : these are the only allowed combinations, others may be not or bad supported by readers
****/
void WriteAsmFdt( const int level )
{
     int status = 0;
     Dtk_transfo TransformationMatrix0, TransformationMatrix1, TransformationMatrix2;
     Dtk_pnt orig1;  orig1[ 0 ] = orig1[ 1 ] = 0.;  orig1[ 2 ] = 200.;
     Dtk_pnt orig2;  orig2[ 0 ] = orig2[ 2 ] = 0.;  orig2[ 1 ] = 300.;
     TransformationMatrix1.setOrigin( orig1 );
     TransformationMatrix2.setOrigin( orig2 );
 
     //The ProcessPart is an user implemented function to create a part
     //  Option to write the Part as external file, with anchors entities for the FDT
     int PartID;
     if( level ) stepw_AddExternalReference( L"essaiPart", L"RootPart.step", PartID );
     PartID = WritePart( L"essaiPart", 2 );
 
     int ProdID, ProdID2;
 
     //we init a product context in STEP : sub-assembly ...
     //  Option to write the sub-assembly as external file
 
     // To write FDT, Instances must have a user identifier
     Dtk_InfoPtr instance_info = Dtk_Info::create();
 
     if( level > 1 ) stepw_AddExternalReference( L"SubProduct", L"RootSubAssy.step", ProdID2 );
     status = stepw_InitProduct( L"SubProduct", ProdID2 );
 
     //...to insert a Part instance of "essaiPart"
     int InstID1, InstID2, InstID3, InstID4;
     int UserInstID1 = 1, UserInstID2 = 2, UserInstID3 = 3, UserInstID4 = 4;
     status = stepw_AddInstanceWithInfo( ProdID2, PartID, TransformationMatrix0, "PartInstance", instance_info, InstID1, UserInstID1, "RD1" );
 
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID2 );
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID2 );
 
     //we create another assembly node (main assembly)
     status = stepw_InitProduct( L"RootProduct", ProdID );
     //and we insert another time the a part instance - multinstancing management -
 
     //   status = stepw_AddInstance(ProdID, PartID, TransformationMatrix0,"PartInstance");
     status = stepw_AddInstanceWithInfo( ProdID, PartID, TransformationMatrix0, "PartInstance", instance_info, InstID2, UserInstID2, "RD2" );
 
     //and two instances to the first sub-assembly
     status = stepw_AddInstanceWithInfo( ProdID, ProdID2, TransformationMatrix1, "PartInstance1", instance_info, InstID3, UserInstID3, "RD3" );
     status = stepw_AddInstanceWithInfo( ProdID, ProdID2, TransformationMatrix2, "PartInstance2", instance_info, InstID4, UserInstID4, "RD4" );
 
     //finally we close the second assembly (root) ...
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID );
 
     // Here is the specifics of this variant : writing a FDT at assembly level
     Dtk_FdtPtr fdtasm = sampleWriter::CreateFdtDimAssy();
     int nodeid_fdtasm = 1002;
     stepw_InitNodeContext( nodeid_fdtasm );
 
     stepw_Add3DPartFDT( fdtasm );
     // Then write the geometrical links : within the same Part
     int geomid1 = 73;
     stepw_ER ElementReference_1;
     stepw_ERP ElementReferencePath_1;
     stepw_CreateInstancePath( ElementReferencePath_1 );
     stepw_AddInstanceToPath( ElementReferencePath_1, UserInstID2 );
     stepw_CreateReference( ElementReference_1, geomid1, PartID );
     stepw_SetReferencePath( ElementReference_1, ElementReferencePath_1 );
     stepw_AddReference( ElementReference_1 );
 
     int geomid2 = 33;
     stepw_ER ElementReference_2;
     stepw_ERP ElementReferencePath_2;
     stepw_CreateInstancePath( ElementReferencePath_2 );
     stepw_AddInstanceToPath( ElementReferencePath_2, UserInstID4 );
     stepw_AddInstanceToPath( ElementReferencePath_2, UserInstID1 );
     stepw_CreateReference( ElementReference_2, geomid2, PartID );
     stepw_SetReferencePath( ElementReference_2, ElementReferencePath_2 );
     stepw_AddReference( ElementReference_2 );
 
     stepw_EndNodeContext();
 
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID );
 
}
 
 
 
/****
Main routine of writing a STEP file
- Opening and Closing of the File
- (writing of content between the 2)
- Global Settings
****/
 
/****    mainstep0 : writing all in one file    ****/
int
StepWriteSample( const Dtk_string &inResultDirectory )
{
     Dtk_string outputFileName, outputDirectory;
 
     cout << endl << "----------------------------------------------" << endl;
     cout << "Step Write start" << endl;
 
     // Choosing output directory and file name
     outputDirectory = inResultDirectory + L"dtk/Step/";
     outputDirectory.FixPathSeparator();
     outputDirectory.mkdir();
     outputFileName = outputDirectory + L"RootProduct.step";
 
     // First we initialize writing with name of files and protection function
     // First you have to initialize the Step writer
     // StepVersion 4 is for AP242
     int status = stepw_InitFile( outputFileName, " sample ", 4 );
     if( status )
     {
          cout << "erreur : " << dtkTypeError( status ).c_str() << endl;
          return status;
     }
 
     stepw_SetModeProp( 1 );
 
     stepw_SetModeFdt( 1 );
 
     //    Test :
     //           WritePartCyl for just one part (cylinder) with tessellation
     //           WritePart for just one part (cube) with properties, with or without FDT (2nd argument = 1 to have FDT, 0 none)
     //           WriteAsm  for an assembly (all in one file)
     //           WriteAsmExt for an assembly in several files - complete or assembly only
     //           WriteAsmFdt for an assembly with FDT at assembly level - monolithic or in several files
     //   WritePartCyl("CylPart",2);
     //   WritePart("onlyonePart",1);  // 1 : with FDT
     //   WriteAsm();
     WriteAsmExt( 2, 0 );  // level=2 : nested external references - asm_only=0 : complete writing including leaf parts
                                // WriteAsmFdt(2);
 
                                //...and the Main writer
     stepw_EndFile();
     cout << "=> " << outputFileName.c_str() << endl;
     cout << "Step Write end" << endl;
 
     return 0;
}
 
//end_of_mainstepwrite
 

//start_of_mainstepwrite
//This sample show you how to write a body with an unbounded surface and wireframe entity.
// IFDT = 1 to also write FDT, View, etc - 0 none
// IFDT = 2 to use this Part by a FDT in Assembly (needs to declare ANCHOR entities)
 
int WritePart( const Dtk_string& inReferenceName, const int ifdt )
{
     // We construct a body 
     Dtk_BodyPtr body = sampleWriter::CreateCube_2();
 
     //  TEST : adding wireframe
     Dtk_BodyPtr sbody = sampleWriter::CreateCurves_2();
 
     //You init the Part writer 
     int PartID;
     if( stepw_InitProduct( inReferenceName, PartID ) )
          return PartID;
 
     //  Following is required for writing of FDT with geometric links - either inside the Part or external designating anchors
     //  To be called before starting the part
     stepw_SetAnchorProduct( PartID, PartID );
 
     //  IFDT = 2 : declaring some entities as ANCHOR for a FDT in Assembly
     //  This is needed where the Assembly is written in several files
     if( ifdt == 2 )
     {
          int geomid1 = 73;
          Dtk_string anchor1( "73" );
          stepw_AddAnchorItem( geomid1, PartID, anchor1 );
          int geomid2 = 33;
          Dtk_string anchor2( "33" );
          stepw_AddAnchorItem( geomid2, PartID, anchor2 );
     }
 
     // test : properties defined as Metadata (at part level)
     stepw_SetModeProp( 2 );
     Dtk_MetaDataPtr metadata_1 = Dtk_MetaData::CreateMetaData( Dtk_MetaData::TypeConfigurationProperty, "Revision", "test_revision_4" );
 
     Dtk_MetaDataPtr metadata_2 = Dtk_MetaData::CreateMetaData( Dtk_MetaData::TypeProperty, "Cost", "12500 $" );
 
     stepw_Init3DPart( PartID );
     stepw_Add3DPartProperty( metadata_1 );
     stepw_Add3DPartProperty( metadata_2 );
 
     //  test : properties at part level - on the product_definition_shape
     // test : validation properties - AP242 form (no intermediate name)
     Dtk_InfoPtr part_info = Dtk_Info::create();
     Dtk_Val valprop_area( "1234 mm2" );
     part_info->AddAttribute( "geometric validation property :  : surface area measure", valprop_area );
     Dtk_Val valprop_volume( "5678 mm3" );
     part_info->AddAttribute( "geometric validation property :  : volume measure", valprop_volume );
     Dtk_Val valprop_CG( Dtk_pnt( 12., 34., 56. ) );
     part_info->AddAttribute( "geometric validation property :  : centre point", valprop_CG );
     stepw_SetPartProperties( part_info, 1 );
     // test : properties at part level - on the product_definition
     Dtk_InfoPtr prod_info = Dtk_Info::create();
     Dtk_Val mat_name( "Steel YC38" );
     prod_info->AddAttribute( "material property : material name : Steel", mat_name );
     Dtk_Val mat_density( "7.89 g/cm3" );
     prod_info->AddAttribute( "material property : density : density measure", mat_density );
     stepw_SetPartProperties( prod_info, 2 );
 
     // And we write body in the step file
     stepw_Write3DPartBody( body );
 
     stepw_Write3DPartBody( sbody );  // test wire frame
 
                                              //  IFDT = 1 : writing a FDT within the Part (here, a Datum)
     if( ifdt == 1 )
     {
 
          //  Test Write FDT
          Dtk_FdtPtr Fdt1 = sampleWriter::CreateFdtDatum();
          int nodeid_Fdt1 = 1001;
          stepw_InitNodeContext( nodeid_Fdt1 );
          stepw_Add3DPartFDT( Fdt1 );
          // Then write the geometrical links : within the same Part
          int geomid1 = 73;
          stepw_ER ElementReference_1;
          stepw_CreateReference( ElementReference_1, geomid1 );
          stepw_AddReference( ElementReference_1 );
          stepw_EndNodeContext();
 
          //  Test Write View
          Dtk_ModelDisplayPtr View1 = CreateView();
          int nodeid_View1 = 2001;
          stepw_InitNodeContext( nodeid_View1 );
          stepw_Add3DModelDisplay( View1, 0 );
          stepw_ER ElementReference_Fdt1;
          stepw_CreateReferenceToNode( ElementReference_Fdt1, nodeid_Fdt1, PartID );
          stepw_AddReference( ElementReference_Fdt1 );
          stepw_EndNodeContext();
 
          //  Test Construction Geometries & Axis System
          Dtk_AxisSystemPtr AxisSystem = CreateAxisSystem();
          stepw_Write3DAxisSystem( AxisSystem );
          Dtk_BodyPtr constr_plane = sampleWriter::CreateConstructionPlane();
          stepw_Write3DConstructionGeometry( constr_plane );
 
     }
 
     // Close the current part
     stepw_End3DPart();
 
     stepw_EndProduct( PartID );
 
     return PartID;
}
 
//  Another Part, to test other data - with colors and tessellation
// itessel = 0 or 1 : body only
// itessel = 2 : tessellation only
// itessel = 3 : body + tessellation
 
int WritePartCyl( const Dtk_string& inReferenceName, const int itessel )
{
     // We construct a body : a cylinder
     Dtk_BodyPtr body = sampleWriter::CreateCylinder();
 
 
     //You init the Part writer 
     int PartID;
     if( stepw_InitProduct( inReferenceName, PartID ) )
          return PartID;
 
     //  Following is required for writing of FDT with geometric links
     //  To be called before starting the part
     stepw_SetAnchorProduct( PartID, PartID );
 
     stepw_Init3DPart( PartID );
 
     //  test : properties at part level - on the product_definition_shape
     // test : validation properties - AP242 form (no intermediate name)
     Dtk_InfoPtr part_info = Dtk_Info::create();
     Dtk_Val valprop_area( "1234 mm2" );
     part_info->AddAttribute( "geometric validation property :  : surface area measure", valprop_area );
     Dtk_Val valprop_volume( "5678 mm3" );
     part_info->AddAttribute( "geometric validation property :  : volume measure", valprop_volume );
     Dtk_Val valprop_CG( Dtk_pnt( 12., 34., 56. ) );
     part_info->AddAttribute( "geometric validation property :  : centre point", valprop_CG );
     stepw_SetPartProperties( part_info, 1 );
     // test : properties at part level - on the product_definition
     Dtk_InfoPtr prod_info = Dtk_Info::create();
 
     stepw_SetPartProperties( prod_info, 2 );
     //  Standard use : write a Body : the cylinder
     stepw_Write3DPartBody( body );
     /*
     // Test of Tessellation - NOT RECOMMENDED ON VERSIONS STEP BEFORE  AP242
     //  To get the functions of tessellation : tess/tess.h
     Dtk_MeshPtr mesh;
     Dtk_tab<Dtk_MeshPtr> listMesh;
     Dtk_tab<Dtk_Int32> IsSolid;
     Dtk_bool TessWireframe=false;
     tess_InitTesselation("../../SampleFiles/dtk/",0.2);
     tess_BodyToMeshes (body,listMesh,IsSolid,TessWireframe);
     int nbmesh = listMesh.size();
     if (nbmesh >= 1)
     mesh = listMesh[0];
     // Writing Tessellation
     //  Option 1 : Tessellation "alone" i.e. not linked to the body
     // we write body in the step file, then the tessellation separately
     if (itessel <= 1) stepw_Write3DPartBody (body);
     if (itessel == 2) stepw_Write3DPartMesh (mesh);
     if (itessel == 3) stepw_Write3DPartBodyWithMesh(body,mesh);
     */
     // Close the current part
     stepw_End3DPart();
 
     stepw_EndProduct( PartID );
 
     return PartID;
}
 
/*
Variants of writing an assembly
All write an assembly the same way, but the variants focus on a specific capability (not all together for more clarity)
WriteAsm is the basic variant, completed with validation properties and user attributes on the assembly or items of it
*/
 
void WriteAsm()
{
     int status = 0;
     Dtk_transfo TransformationMatrix0, TransformationMatrix1, TransformationMatrix2;
     Dtk_pnt orig1;  orig1[ 0 ] = orig1[ 1 ] = 0.;  orig1[ 2 ] = 200.;
     Dtk_pnt orig2;  orig2[ 0 ] = orig2[ 2 ] = 0.;  orig2[ 1 ] = 300.;
     TransformationMatrix1.setOrigin( orig1 );
     TransformationMatrix2.setOrigin( orig2 );
 
     //The ProcessPart is an user implemented function to create a part
     int PartID = WritePart( L"essaiPart", 0 );
 
     int ProdID, ProdID2;
 
     //we init a product context in STEP : sub-assembly ...
     status = stepw_InitProduct( L"SubProduct", ProdID2 );
 
     //...to insert a Part instance of "essaiPart"
     status = stepw_AddInstance( ProdID2, PartID, TransformationMatrix0, "PartInstance" );
 
     // test : validation properties - AP242 form (no intermediate name)
     Dtk_InfoPtr part_info = Dtk_Info::create();
     stepw_SetPartProperties( part_info, 1 );
     Dtk_Val valprop_notional_CG( Dtk_pnt( 10., 20., 30. ) );
     part_info->AddAttribute( "assembly validation property :  : centre point", valprop_notional_CG );
     stepw_SetPartProperties( part_info, 1 );
     Dtk_InfoPtr prod_info = Dtk_Info::create();
     Dtk_Val valprop_nbchildren( 3 );
     prod_info->AddAttribute( "assembly validation property :  : number of children", valprop_nbchildren );
     stepw_SetPartProperties( prod_info, 2 );
 
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID2 );
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID2 );
 
 
     //we create another assembly node (main assembly)
     status = stepw_InitProduct( L"RootProduct", ProdID );
     //and we insert another time the a part instance - multinstancing management -
     // test : instance with property
     Dtk_InfoPtr instance_info = Dtk_Info::create();
     Dtk_Val instance_kind( "first" );
     instance_info->AddAttribute( "user defined attribute : : instance kind", instance_kind );
     Dtk_Val instance_valprop( 1 );
     instance_info->AddAttribute( "attribute validation property :  : text user attributes", instance_valprop );
     //   status = stepw_AddInstance(ProdID, PartID, TransformationMatrix0,"PartInstance");
     status = stepw_AddInstanceWithInfo( ProdID, PartID, TransformationMatrix0, "PartInstance", instance_info );
 
     //and two instances to the first sub-assembly
     status = stepw_AddInstance( ProdID, ProdID2, TransformationMatrix1, "PartInstance1" );
     status = stepw_AddInstance( ProdID, ProdID2, TransformationMatrix2, "PartInstance2" );
 
     //finally we close the second assembly (root) ...
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID );
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID );
 
}
 
/****
Variant to write the different products in several files
Similar to WriteAsm (just duplicated), with differences of external files : calls to stepw_AddExternalReference
level = 0 : monolithic file (so, like WriteAsm)
level = 1 : basic external references, only the leaf parts are written separately
level = 2 : nested external references, each component is written separately
Warning : these are the only allowed combinations, others may be not or bad supported by readers
 
asm_only = 0 : complete mode, all components (assemblies and leaf part) are written - monolithic writing is always complete
asm_only = 1 : only assemblies are written, part are declared as external products but not written
****/
void WriteAsmExt( const int level, const int asm_only )
{
     int status = 0;
     Dtk_transfo TransformationMatrix0, TransformationMatrix1, TransformationMatrix2;
     Dtk_pnt orig1;  orig1[ 0 ] = orig1[ 1 ] = 0.;  orig1[ 2 ] = 200.;
     Dtk_pnt orig2;  orig2[ 0 ] = orig2[ 2 ] = 0.;  orig2[ 1 ] = 300.;
     TransformationMatrix1.setOrigin( orig1 );
     TransformationMatrix2.setOrigin( orig2 );
 
     //The ProcessPart is an user implemented function to create a part
     //  Option to write the Part as external file
     int PartID;
     if( level ) stepw_AddExternalReference( L"essaiPart", L"RootPart.step", PartID );
 
     //with external references, it is possible to declare parts as external products but not write them
     //  Each product is identified by its NAME - there is a unique equivalence between Product Name and returned vale of ProdID
     if( !asm_only || level == 0 )
          PartID = WritePart( L"essaiPart", 0 );
 
     int ProdID, ProdID2;
 
     //we init a product context in STEP : sub-assembly ...
     //  Option to write the sub-assembly as external file
     if( level > 1 ) stepw_AddExternalReference( L"SubProduct", L"RootSubAssy.step", ProdID2 );
     status = stepw_InitProduct( L"SubProduct", ProdID2 );
 
     //...to insert a Part instance of "essaiPart"
     status = stepw_AddInstance( ProdID2, PartID, TransformationMatrix0, "PartInstance" );
 
     // test : validation properties - AP242 form (no intermediate name)
     Dtk_InfoPtr part_info = Dtk_Info::create();
     stepw_SetPartProperties( part_info, 1 );
     Dtk_Val valprop_notional_CG( Dtk_pnt( 10., 20., 30. ) );
     part_info->AddAttribute( "assembly validation property :  : centre point", valprop_notional_CG );
     stepw_SetPartProperties( part_info, 1 );
     Dtk_InfoPtr prod_info = Dtk_Info::create();
     Dtk_Val valprop_nbchildren( 3 );
     prod_info->AddAttribute( "assembly validation property :  : number of children", valprop_nbchildren );
     stepw_SetPartProperties( prod_info, 2 );
 
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID2 );
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID2 );
 
 
     //we create another assembly node (main assembly)
     status = stepw_InitProduct( L"RootProduct", ProdID );
     //and we insert another time the a part instance - multinstancing management -
     // test : instance with property
     Dtk_InfoPtr instance_info = Dtk_Info::create();
     Dtk_Val instance_kind( "first" );
     instance_info->AddAttribute( "user defined attribute : : instance kind", instance_kind );
     Dtk_Val instance_valprop( 1 );
     instance_info->AddAttribute( "attribute validation property : : text user attributes", instance_valprop );
     //   status = stepw_AddInstance(ProdID, PartID, TransformationMatrix0,"PartInstance");
     status = stepw_AddInstanceWithInfo( ProdID, PartID, TransformationMatrix0, "PartInstance", instance_info );
 
     //and two instances to the first sub-assembly
     status = stepw_AddInstance( ProdID, ProdID2, TransformationMatrix1, "PartInstance1" );
     status = stepw_AddInstance( ProdID, ProdID2, TransformationMatrix2, "PartInstance2" );
 
     //finally we close the second assembly (root) ...
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID );
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID );
 
}
 
 
/****
Variant to write an assembly with a FDT at assembly level - monolithic or with external references
Similar to WriteAsmExt (just duplicated),including external files, but with a FDT in more, and validation properties in less
level = 0 : monolithic file (so, like WriteAsm)
level = 1 : basic external references, only the leaf parts are written separately
level = 2 : nested external references, each component is written separately
Warning : these are the only allowed combinations, others may be not or bad supported by readers
****/
void WriteAsmFdt( const int level )
{
     int status = 0;
     Dtk_transfo TransformationMatrix0, TransformationMatrix1, TransformationMatrix2;
     Dtk_pnt orig1;  orig1[ 0 ] = orig1[ 1 ] = 0.;  orig1[ 2 ] = 200.;
     Dtk_pnt orig2;  orig2[ 0 ] = orig2[ 2 ] = 0.;  orig2[ 1 ] = 300.;
     TransformationMatrix1.setOrigin( orig1 );
     TransformationMatrix2.setOrigin( orig2 );
 
     //The ProcessPart is an user implemented function to create a part
     //  Option to write the Part as external file, with anchors entities for the FDT
     int PartID;
     if( level ) stepw_AddExternalReference( L"essaiPart", L"RootPart.step", PartID );
     PartID = WritePart( L"essaiPart", 2 );
 
     int ProdID, ProdID2;
 
     //we init a product context in STEP : sub-assembly ...
     //  Option to write the sub-assembly as external file
 
     // To write FDT, Instances must have a user identifier
     Dtk_InfoPtr instance_info = Dtk_Info::create();
 
     if( level > 1 ) stepw_AddExternalReference( L"SubProduct", L"RootSubAssy.step", ProdID2 );
     status = stepw_InitProduct( L"SubProduct", ProdID2 );
 
     //...to insert a Part instance of "essaiPart"
     int InstID1, InstID2, InstID3, InstID4;
     int UserInstID1 = 1, UserInstID2 = 2, UserInstID3 = 3, UserInstID4 = 4;
     status = stepw_AddInstanceWithInfo( ProdID2, PartID, TransformationMatrix0, "PartInstance", instance_info, InstID1, UserInstID1, "RD1" );
 
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID2 );
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID2 );
 
     //we create another assembly node (main assembly)
     status = stepw_InitProduct( L"RootProduct", ProdID );
     //and we insert another time the a part instance - multinstancing management -
 
     //   status = stepw_AddInstance(ProdID, PartID, TransformationMatrix0,"PartInstance");
     status = stepw_AddInstanceWithInfo( ProdID, PartID, TransformationMatrix0, "PartInstance", instance_info, InstID2, UserInstID2, "RD2" );
 
     //and two instances to the first sub-assembly
     status = stepw_AddInstanceWithInfo( ProdID, ProdID2, TransformationMatrix1, "PartInstance1", instance_info, InstID3, UserInstID3, "RD3" );
     status = stepw_AddInstanceWithInfo( ProdID, ProdID2, TransformationMatrix2, "PartInstance2", instance_info, InstID4, UserInstID4, "RD4" );
 
     //finally we close the second assembly (root) ...
     // the part for the sub-assembly is REQUIRED - assembly node : created empty
     status = stepw_Init3DPart( ProdID );
 
     // Here is the specifics of this variant : writing a FDT at assembly level
     Dtk_FdtPtr fdtasm = sampleWriter::CreateFdtDimAssy();
     int nodeid_fdtasm = 1002;
     stepw_InitNodeContext( nodeid_fdtasm );
 
     stepw_Add3DPartFDT( fdtasm );
     // Then write the geometrical links : within the same Part
     int geomid1 = 73;
     stepw_ER ElementReference_1;
     stepw_ERP ElementReferencePath_1;
     stepw_CreateInstancePath( ElementReferencePath_1 );
     stepw_AddInstanceToPath( ElementReferencePath_1, UserInstID2 );
     stepw_CreateReference( ElementReference_1, geomid1, PartID );
     stepw_SetReferencePath( ElementReference_1, ElementReferencePath_1 );
     stepw_AddReference( ElementReference_1 );
 
     int geomid2 = 33;
     stepw_ER ElementReference_2;
     stepw_ERP ElementReferencePath_2;
     stepw_CreateInstancePath( ElementReferencePath_2 );
     stepw_AddInstanceToPath( ElementReferencePath_2, UserInstID4 );
     stepw_AddInstanceToPath( ElementReferencePath_2, UserInstID1 );
     stepw_CreateReference( ElementReference_2, geomid2, PartID );
     stepw_SetReferencePath( ElementReference_2, ElementReferencePath_2 );
     stepw_AddReference( ElementReference_2 );
 
     stepw_EndNodeContext();
 
     status = stepw_End3DPart();
 
     //we close the CATProduct context and retrieve the resulting DocID
     status = stepw_EndProduct( ProdID );
 
}
 
 
 
/****
Main routine of writing a STEP file
- Opening and Closing of the File
- (writing of content between the 2)
- Global Settings
****/
 
/****    mainstep0 : writing all in one file    ****/
int
StepWriteSample( const Dtk_string &inResultDirectory )
{
     Dtk_string outputFileName, outputDirectory;
 
     cout << endl << "----------------------------------------------" << endl;
     cout << "Step Write start" << endl;
 
     // Choosing output directory and file name
     outputDirectory = inResultDirectory + L"dtk/Step/";
     outputDirectory.FixPathSeparator();
     outputDirectory.mkdir();
     outputFileName = outputDirectory + L"RootProduct.step";
 
     // First we initialize writing with name of files and protection function
     // First you have to initialize the Step writer
     // StepVersion 4 is for AP242
     int status = stepw_InitFile( outputFileName, " sample ", 4 );
     if( status )
     {
          cout << "erreur : " << dtkTypeError( status ).c_str() << endl;
          return status;
     }
 
     stepw_SetModeProp( 1 );
 
     stepw_SetModeFdt( 1 );
 
     //    Test :
     //           WritePartCyl for just one part (cylinder) with tessellation
     //           WritePart for just one part (cube) with properties, with or without FDT (2nd argument = 1 to have FDT, 0 none)
     //           WriteAsm  for an assembly (all in one file)
     //           WriteAsmExt for an assembly in several files - complete or assembly only
     //           WriteAsmFdt for an assembly with FDT at assembly level - monolithic or in several files
     //   WritePartCyl("CylPart",2);
     //   WritePart("onlyonePart",1);  // 1 : with FDT
     //   WriteAsm();
     WriteAsmExt( 2, 0 );  // level=2 : nested external references - asm_only=0 : complete writing including leaf parts
                                // WriteAsmFdt(2);
 
                                //...and the Main writer
     stepw_EndFile();
     cout << "=> " << outputFileName.c_str() << endl;
     cout << "Step Write end" << endl;
 
     return 0;
}
 
//end_of_mainstepwrite