/************
*
*   This file is part of a tool for reading 3D content in the PRC format.
*   Copyright (C) 2008 Orest Shardt <shardtor (at) gmail dot com>
*
*   This program is free software: you can redistribute it and/or modify
*   it under the terms of the GNU Lesser General Public License as published by
*   the Free Software Foundation, either version 3 of the License, or
*   (at your option) any later version.
*
*   This program is distributed in the hope that it will be useful,
*   but WITHOUT ANY WARRANTY; without even the implied warranty of
*   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
*   GNU Lesser General Public License for more details.
*
*   You should have received a copy of the GNU Lesser General Public License
*   along with this program.  If not, see <http://www.gnu.org/licenses/>.
*
*************/

#include <iostream>
#include <fstream>
#include <iomanip>
#include <string>
#include <sstream>

#include "../include/prc/PRC.h"
#include "describePRC.h"

using std::ostringstream; using std::cout; using std::endl;
using std::hex; using std::dec; using std::string; using std::setw;
using std::setfill;

// describe sections

void describeGlobals(BitByBitData &mData)
{
  mData.tellPosition();
  cout << getIndent() << "--Globals--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_ASM_FileStructureGlobals))
    return;
  indent();

  describeContentPRCBase(mData,false);
  unsigned int numberOfReferencedFileStructures = mData.readUnsignedInt();
  cout << getIndent() << "numberOfReferencedFileStructures "
      << numberOfReferencedFileStructures << endl;
  indent();
  for(unsigned int i = 0; i < numberOfReferencedFileStructures; ++i)
  {
    describeCompressedUniqueID(mData);
  }
  dedent();

  double tessellation_chord_height_ratio = mData.readDouble();
  cout << getIndent() << "tessellation_chord_height_ratio "
      << tessellation_chord_height_ratio << endl;

  double tessellation_angle_degree = mData.readDouble();
  cout << getIndent() << "tessellation_angle_degree "
      << tessellation_angle_degree << endl;

  string default_font_family_name = mData.readString();
  cout << getIndent() << "default_font_family_name \""
      << default_font_family_name << '\"' << endl;

  unsigned int number_of_fonts = mData.readUnsignedInt();
  cout << getIndent() << "number_of_fonts " << number_of_fonts << endl;

  indent();
  for(unsigned int q = 0; q < number_of_fonts; ++q)
  {
    string font_name = mData.readString();
    cout << getIndent() << "font_name \"" << font_name << '\"' << endl;
    unsigned int char_set = mData.readUnsignedInt();
    cout << getIndent() << "char_set " << char_set << endl;
    unsigned int number_of_font_keys = mData.readUnsignedInt();
    cout << getIndent() << "number_of_font_keys " << number_of_font_keys
        << endl;
    indent();
    for(unsigned int i = 0; i < number_of_font_keys; i++)
    {
      unsigned int font_size = mData.readUnsignedInt() - 1;
      cout << getIndent() << "font_size " << font_size << endl;
      unsigned char attributes = mData.readChar();
      cout << getIndent() << "attributes "
          << static_cast<unsigned int>(attributes) << endl;
    }
    dedent();
  }
  dedent();

  unsigned int number_of_colours = mData.readUnsignedInt();
  cout << getIndent() << "number_of_colours " << number_of_colours << endl;
  indent();
  for(unsigned int i = 0; i < number_of_colours; ++i)
    describeRGBColour(mData);
  dedent();

  unsigned int number_of_pictures = mData.readUnsignedInt();
  cout << getIndent() << "number_of_pictures " << number_of_pictures << endl;
  indent();
  for(unsigned int i=0;i<number_of_pictures;i++)
    describePicture(mData);
  dedent();

  unsigned int number_of_texture_definitions = mData.readUnsignedInt();
  cout << getIndent() << "number_of_texture_definitions "
      << number_of_texture_definitions << endl;
  indent();
  for(unsigned int i=0;i<number_of_texture_definitions;i++)
    describeTextureDefinition(mData);
  dedent();


  unsigned int number_of_materials = mData.readUnsignedInt();
  cout << getIndent() << "number_of_materials " << number_of_materials << endl;
  indent();
  for(unsigned int i=0;i<number_of_materials;i++)
    describeMaterial(mData);
  dedent();

  unsigned int number_of_line_patterns = mData.readUnsignedInt();
  cout << getIndent() << "number_of_line_patterns "
      << number_of_line_patterns << endl;
  indent();
  for(unsigned int i=0;i<number_of_line_patterns;i++)
    describeLinePattern(mData);
  dedent();

  unsigned int number_of_styles = mData.readUnsignedInt();
  cout << getIndent() << "number_of_styles " << number_of_styles << endl;
  indent();
  for(unsigned int i=0;i<number_of_styles;i++)
    describeCategory1LineStyle(mData);
  dedent();

  unsigned int number_of_fill_patterns = mData.readUnsignedInt();
  cout << getIndent() << "number_of_fill_patterns "
      << number_of_fill_patterns << endl;
  indent();
  for(unsigned int i=0;i<number_of_fill_patterns;i++)
    describeFillPattern(mData);
  dedent();

  unsigned int number_of_reference_coordinate_systems = mData.readUnsignedInt();
  cout << getIndent() << "number_of_reference_coordinate_systems "
      << number_of_reference_coordinate_systems << endl;
  indent();
  for(unsigned int i=0;i<number_of_reference_coordinate_systems;i++)
    //NOTE: must be PRC_TYPE_RI_CoordinateSystem
    describeRepresentationItem(mData);
  dedent();

  describeUserData(mData);
  dedent();
  mData.tellPosition();
}

void describeTree(BitByBitData &mData)
{
  mData.tellPosition();
  cout << getIndent() << "--Tree--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_ASM_FileStructureTree))
    return;
  indent();
  describeContentPRCBase(mData,false);

  unsigned int number_of_part_definitions = mData.readUnsignedInt();
  cout << getIndent() << "number_of_part_definitions "
      << number_of_part_definitions << endl;
  indent();
  for(unsigned int i = 0; i < number_of_part_definitions; ++i)
  {
    describePartDefinition(mData);
  }
  dedent();

  unsigned int number_of_product_occurrences = mData.readUnsignedInt();
  cout << getIndent() << "number_of_product_occurrences "
      << number_of_product_occurrences << endl;
  indent();
  for(unsigned int i = 0; i < number_of_product_occurrences; ++i)
  {
    describeProductOccurrence(mData);
  }
  dedent();

  describeFileStructureInternalData(mData);

  describeUserData(mData);
  dedent();
  mData.tellPosition();
}

void describeTessellation(BitByBitData &mData)
{
  mData.tellPosition();
  cout << getIndent() << "--Tessellation--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_ASM_FileStructureTessellation))
    return;
  indent();

  describeContentPRCBase(mData,false);

  unsigned int number_of_tessellations = mData.readUnsignedInt();
  cout << getIndent() << "number_of_tessellations "
      << number_of_tessellations << endl;
  indent();
  for(unsigned int i = 0; i < number_of_tessellations; ++i)
  {
    unsigned int type = mData.readUnsignedInt();
    cout << getIndent() << "tessellation type " << type << endl;
    switch(type)
    {
      case PRC_TYPE_TESS_3D:
        describe3DTess(mData);
        break;
      case PRC_TYPE_TESS_3D_Wire:
        describe3DWireTess(mData);
        break;
      case PRC_TYPE_TESS_Markup:
        describe3DMarkupTess(mData);
        break;
      case PRC_TYPE_TESS_3D_Compressed:
        describeHighlyCompressed3DTess(mData);
        break;
      default:
        cout << getIndent() << "Unrecognized tessellation data type "
            << type << endl;
        break;
    }
  }
  dedent();

  describeUserData(mData);

  dedent();
  mData.tellPosition();
}

void describeGeometry(BitByBitData &mData)
{
  mData.tellPosition();
  cout << getIndent() << "--Geometry--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_ASM_FileStructureGeometry))
    return;
  indent();

  describeContentPRCBase(mData,false);

  unsigned int number_of_topological_contexts = mData.readUnsignedInt();
  cout << getIndent() << "number_of_topological_contexts "
      << number_of_topological_contexts << endl;
  indent();
  for(unsigned int i = 0; i < number_of_topological_contexts; ++i)
  {
    describeTopoContext(mData);
    unsigned int number_of_bodies = mData.readUnsignedInt();
    cout << getIndent() << "number_of_bodies " << number_of_bodies << endl;
    for(unsigned int i = 0; i < number_of_bodies; ++i)
    {
      describeBody(mData);
    }
  }
  dedent();

  describeUserData(mData);
  dedent();
  mData.tellPosition();
}

void describeExtraGeometry(BitByBitData &mData)
{
  mData.tellPosition();
  cout << getIndent() << "--Extra Geometry--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_ASM_FileStructureExtraGeometry))
    return;
  indent();

  describeContentPRCBase(mData,false);

  unsigned int number_of_contexts = mData.readUnsignedInt();
  cout << getIndent() << "number_of_contexts " << number_of_contexts << endl;
  indent();
  for(unsigned int i = 0; i < number_of_contexts; ++i)
  {
    // geometry summary
    unsigned int number_of_bodies = mData.readUnsignedInt();
    cout << getIndent() << "number_of_bodies " << number_of_bodies << endl;
    indent();
    for(unsigned int j = 0; j < number_of_bodies; ++j)
    {
      unsigned int serial_type = mData.readUnsignedInt();
      cout << getIndent() << "serial_type " << serial_type << endl;
      indent();
      if(isCompressedSerialType(serial_type))
        cout << getIndent() << "serialTolerance " << mData.readDouble();
      dedent();
    }
    dedent();
    // context graphics
    resetCurrentGraphics();
    unsigned int number_of_treat_types = mData.readUnsignedInt();
    cout << getIndent() << "number_of_treat_types "
        << number_of_treat_types << endl;
    indent();
    for(unsigned int i = 0; i < number_of_treat_types; ++i)
    {
      cout << getIndent() << "element_type " << mData.readUnsignedInt() << endl;
      unsigned int number_of_elements = mData.readUnsignedInt();
      cout << getIndent() << "number_of_elements "
          << number_of_elements << endl;
      indent();
      for(unsigned int j = 0; j < number_of_elements; ++j)
      {
        if(mData.readBit())
        {
          describeGraphics(mData);
        }
        else
        {
          cout << getIndent() << "Element has no graphics" << endl;
        }
      }
      dedent();
    }
    dedent();
  }
  dedent();

  describeUserData(mData);
  dedent();
  mData.tellPosition();
}


void describeModelFileData(BitByBitData &mData,
                           unsigned int numberOfFileStructures)
{
  mData.tellPosition();
  cout << getIndent() << "--Model File--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_ASM_ModelFile))
    return;
  indent();

  describeContentPRCBase(mData,false);

  describeUnit(mData);
  unsigned int numberOfProductOccurrences = mData.readUnsignedInt();
  cout << getIndent() << "Number of Product Occurrences "
      << numberOfProductOccurrences << endl;
  indent();
  for(unsigned int i = 0; i < numberOfProductOccurrences; ++i)
  {
    describeCompressedUniqueID(mData);
    cout << getIndent() << "index_position + 1 = "
        << mData.readUnsignedInt() << endl;
    cout << getIndent() << "active? " << (mData.readBit()?"yes":"no")
        << endl << endl;
  }
  dedent();
  for(unsigned int i = 0; i < numberOfFileStructures; ++i)
  {
    cout << getIndent() << "File Structure Index in Model File "
        << mData.readUnsignedInt() << endl;
  }

  describeUserData(mData);
  dedent();
  mData.tellPosition();
}

// subsections
// void describe(BitByBitData &mData)
// {
//
// }

void describeLight(BitByBitData &mData)
{
  unsigned int ID = mData.readUnsignedInt();
  if(ID == PRC_TYPE_GRAPH_AmbientLight)
    cout << getIndent() << "--Ambient Light--" << endl;
  else
    return;

  indent();

  describeContentPRCBase(mData,true);

  cout << getIndent() << "ambient colour index: " << mData.readUnsignedInt()-1 << endl
      << getIndent() << "diffuse colour index: " << mData.readUnsignedInt()-1 << endl
      << getIndent() << "specular colour index: " << mData.readUnsignedInt()-1 << endl;

  describeUserData(mData);
  describeUserData(mData); // why?
  dedent();
}

void describeCamera(BitByBitData &mData)
{
  cout << getIndent() << "--Camera--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_GRAPH_Camera))
    return;
  indent();

  describeContentPRCBase(mData,true);

  cout << getIndent() << (mData.readBit()?"orthographic":"perspective") << endl;
  cout << getIndent() << "Camera Position" << endl;
  describeVector3d(mData);
  cout << getIndent() << "Look At Point" << endl;
  describeVector3d(mData);
  cout << getIndent() << "Up" << endl;
  describeVector3d(mData);
  cout << getIndent() << "X field of view angle (perspective) || X scale (orthographic) "
      << mData.readDouble() << endl;
  cout << getIndent() << "Y field of view angle (perspective) || Y scale (orthographic) "
      << mData.readDouble() << endl;
  cout << getIndent() << "aspect ratio x/y " << mData.readDouble() << endl;
  cout << getIndent() << "near z clipping plane distance from viewer " << mData.readDouble() << endl;
  cout << getIndent() << "far z clipping plane distance from viewer " << mData.readDouble() << endl;
  cout << getIndent() << "zoom factor " << mData.readDouble() << endl;
  dedent();
}

bool describeContentCurve(BitByBitData &mData)
{
  describeBaseGeometry(mData);
  cout << getIndent() << "extend_info " << mData.readUnsignedInt() << endl;
  bool is_3d = mData.readBit();
  cout << getIndent() << "is_3d " << (is_3d?"yes":"no") << endl;
  return is_3d;
}

void describeParameterization(BitByBitData &mData)
{
  cout << getIndent() << "--Parameterization--" << endl;
  indent();
  describeExtent1d(mData);
  cout << getIndent() << "parameterization_coeff_a " << mData.readDouble()
      << endl;
  cout << getIndent() << "parameterization_coeff_b " << mData.readDouble()
      << endl;
  dedent();
}

void describeCurvCircle(BitByBitData &mData)
{
  cout << getIndent() << "--Circle--" << endl;
  indent();

  if(describeContentCurve(mData))
    describeTransformation3d(mData);
  else
    describeTransformation2d(mData);

  describeParameterization(mData);

  cout << getIndent() << "radius " << mData.readDouble() << endl;

  dedent();
}

void describeCurvLine(BitByBitData &mData)
{
  cout << getIndent() << "--Line--" << endl;
  indent();

  if(describeContentCurve(mData))
    describeTransformation3d(mData);
  else
    describeTransformation2d(mData);

  describeParameterization(mData);

  dedent();
}

void describeContentWireEdge(BitByBitData &mData)
{
  cout << getIndent() << "--WireEdge--" << endl;
  indent();

  describeBaseTopology(mData);

  describeObject(mData); //3d_curve

  bool curve_trim_interval = mData.readBit();
  cout << getIndent() << "curve_trim_interval "
      << (curve_trim_interval?"yes":"no") << endl;
  if(curve_trim_interval)
  {
    describeExtent1d(mData);
  }

  dedent();
}

void describeUVParametrization(BitByBitData &mData)
{
  cout << getIndent() << "--UV Parameterization--" << endl;
  indent();

  cout << getIndent() << "swap_uv " << (mData.readBit()?"yes":"no") << endl;
  cout << getIndent() << "Domain" << endl;
  indent(); describeExtent2d(mData); dedent();
  cout << getIndent() << "parameterization_on_u_coeff_a "
      << mData.readDouble() << endl;
  cout << getIndent() << "parameterization_on_v_coeff_a "
      << mData.readDouble() << endl;
  cout << getIndent() << "parameterization_on_u_coeff_b "
      << mData.readDouble() << endl;
  cout << getIndent() << "parameterization_on_v_coeff_b "
      << mData.readDouble() << endl;

  dedent();
}

bool isCompressedSerialType(unsigned int type)
{
  return false; // TODO: actually check the type!!!
}

void describeSurfNURBS(BitByBitData &mData)
{
  cout << getIndent() << "--NURBS surface--" << endl;
  indent();

  describeContentSurface(mData);

  bool is_rational = mData.readBit();
  cout << getIndent() << "is_rational " << (is_rational?"yes":"no") << endl;

  unsigned int degree_in_u = mData.readUnsignedInt();
  cout << getIndent() << "degree_in_u " << degree_in_u << endl;
  unsigned int degree_in_v = mData.readUnsignedInt();
  cout << getIndent() << "degree_in_v " << degree_in_v << endl;

  unsigned int number_of_control_points_in_u = mData.readUnsignedInt()+1;
  cout << getIndent() << "number_of_control_points_in_u "
      << number_of_control_points_in_u << endl;
  unsigned int number_of_control_points_in_v = mData.readUnsignedInt()+1;
  cout << getIndent() << "number_of_control_points_in_v "
      << number_of_control_points_in_v << endl;

  unsigned int number_of_knots_in_u = mData.readUnsignedInt()+1;
  cout << getIndent() << "number_of_knots_in_u "
      << number_of_knots_in_u << endl;
  unsigned int number_of_knots_in_v = mData.readUnsignedInt()+1;
  cout << getIndent() << "number_of_knots_in_v "
      << number_of_knots_in_v << endl;

  indent();
  for(unsigned int i = 0; i < number_of_control_points_in_u; ++i)
  {
    for(unsigned int j = 0; j < number_of_control_points_in_v; ++j)
    {
      double x = mData.readDouble();
      double y = mData.readDouble();
      double z = mData.readDouble();
      cout << getIndent() << "control point " << i << ' ' << j << ": ("
          << x << ',' << y << ',' << z;

      if(is_rational)
        cout << ',' << mData.readDouble();

      cout << ')' << endl;
    }
  }
  dedent();

  cout << getIndent() << "knots in u" << endl;
  indent();
  for(unsigned int i = 0; i < number_of_knots_in_u; ++i)
    cout << getIndent() << mData.readDouble() << endl;
  dedent();

  cout << getIndent() << "knots in v" << endl;
  indent();
  for(unsigned int i = 0; i < number_of_knots_in_v; ++i)
    cout << getIndent() << mData.readDouble() << endl;
  dedent();

  cout << getIndent() << "knot_type " << mData.readUnsignedInt() << endl;
  cout << getIndent() << "surface_form " << mData.readUnsignedInt() << endl;

  dedent();
}

void describeCurvNURBS(BitByBitData &mData)
{
  cout << getIndent() << "--NURBS curve--" << endl;
  indent();
  describeContentCurve(mData);

  bool is_rational = mData.readBit();
  cout << getIndent() << "is_rational " << (is_rational?"yes":"no") << endl;

  unsigned int degree = mData.readUnsignedInt();
  cout << getIndent() << "degree " << degree << endl;

  unsigned int number_of_control_points = mData.readUnsignedInt()+1;
  cout << getIndent() << "number_of_control_points "
      << number_of_control_points << endl;

  unsigned int number_of_knots = mData.readUnsignedInt()+1;
  cout << getIndent() << "number_of_knots " << number_of_knots << endl;

  indent();
  for(unsigned int i = 0; i < number_of_control_points; ++i)
  {
    double x = mData.readDouble();
    double y = mData.readDouble();
    double z = mData.readDouble();
    cout << getIndent() << "control point " << i << ": ("
        << x << ',' << y << ',' << z;

    if(is_rational)
      cout << ',' << mData.readDouble();

    cout << ')' << endl;
  }
  dedent();

  cout << getIndent() << "knots" << endl;
  indent();
  for(unsigned int i = 0; i < number_of_knots; ++i)
    cout << getIndent() << mData.readDouble() << endl;
  dedent();

  cout << getIndent() << "knot_type " << mData.readUnsignedInt() << endl;
  cout << getIndent() << "surface_form " << mData.readUnsignedInt() << endl;

  dedent();
}

void describeCurvPolyLine(BitByBitData &mData)
{
  cout << getIndent() << "--PolyLine--" << endl;
  indent();
  describeContentCurve(mData);
  describeTransformation3d(mData);
  describeParameterization(mData);

  unsigned int number_of_points = mData.readUnsignedInt();
  cout << getIndent() << "number_of_points " << number_of_points << endl;

  indent();
  for(unsigned int i = 0; i < number_of_points; ++i)
  {
    describeVector3d(mData);
  }
  dedent();

  dedent();
}

void describeSurfCylinder(BitByBitData &mData)
{
  cout << getIndent() << "--Cylinder surface--" << endl;
  indent();

  describeContentSurface(mData);
  describeTransformation3d(mData);
  describeUVParametrization(mData);
  cout << getIndent() << "radius " << mData.readDouble() << endl;

  dedent();
}

void describeSurfPlane(BitByBitData &mData)
{
  cout << getIndent() << "--Plane surface--" << endl;
  indent();
  mData.setShowBits(true);
  describeContentSurface(mData);
  mData.tellPosition();
  //TODO: something is wrong, very wrong!!!
  // For now, all this does is search until the end of the data block,
  // assuming that the default parameterization [-inf,inf]x[-inf,inf] was used
  BitPosition bp = mData.getPosition();
  double nInf1 = mData.readDouble();
  double nInf2 = mData.readDouble();
  double inf1 = mData.readDouble();
  double inf2 = mData.readDouble();
  double one1 = mData.readDouble();
  double one2 = mData.readDouble();
  double zero1 = mData.readDouble();
  double zero2 = mData.readDouble();
  while(!(nInf1 == -12345.0 && nInf2 == -12345.0 && inf1 == 12345.0
          && inf2 == 12345.0 && one1 == 1.0 && one2 == 1.0 && zero1 == 0.0
          && zero2 == 0.0))
  {
    mData.setPosition(bp);
    cout << mData.readBit();
    bp = mData.getPosition();

    nInf1 = mData.readDouble();
    nInf2 = mData.readDouble();
    inf1 = mData.readDouble();
    inf2 = mData.readDouble();
    one1 = mData.readDouble();
    one2 = mData.readDouble();
    zero1 = mData.readDouble();
    zero2 = mData.readDouble();
  }
  cout << endl;
  if(bp.bitIndex == 0)
  {
    bp.bitIndex = 7;
    --bp.byteIndex;
  }
  else
  {
    --bp.bitIndex;
  }
  mData.setPosition(bp);

/*
  // this is what the 8137 docs say it should be
  describeTransformation3d(mData);
  cout << getIndent() << "UV domain" << endl;
  indent();
  cout << getIndent() << "Min: " << endl;
  describeVector2d(mData);
  cout << getIndent() << "Max: " << endl;
  describeVector2d(mData);
  dedent();

  cout << getIndent() << "u coef. a = " << mData.readDouble() << endl;
  cout << getIndent() << "v coef. a = " << mData.readDouble() << endl;
  cout << getIndent() << "u coef. b = " << mData.readDouble() << endl;
  cout << getIndent() << "v coef. b = " << mData.readDouble() << endl;
*/
  mData.setShowBits(false);
  dedent();
}

void describeTopoFace(BitByBitData &mData)
{
  cout << getIndent() << "--Face--" << endl;
  indent();

  describeBaseTopology(mData);

  cout << getIndent() << "base_surface" << endl;
  indent();
  describeObject(mData);
  dedent();

  bool surface_trim_domain = mData.readBit();
  cout << getIndent() << "surface_trim_domain "
      << (surface_trim_domain?"yes":"no") << endl;
  if(surface_trim_domain)
  {
    indent();
    describeExtent2d(mData);
    dedent();
  }

  bool have_tolerance = mData.readBit();
  cout << getIndent() << "have_tolerance "
      << (have_tolerance?"yes":"no") << endl;
  if(have_tolerance)
    cout << getIndent() << "tolerance " << mData.readDouble() << endl;

  unsigned int number_of_loops = mData.readUnsignedInt();
  cout << getIndent() << "number_of_loops " << number_of_loops << endl;
  cout << getIndent() << "outer_loop_index " << mData.readInt() << endl;
  indent();
  for(unsigned int i = 0; i < number_of_loops; ++i)
  {
    describeObject(mData);
  }
  dedent();

  dedent();
}

void describeTopoLoop(BitByBitData &mData)
{
  cout << getIndent() << "--Loop--" << endl;
  indent();

  describeBaseTopology(mData);

  cout << getIndent() << "orientation_with_surface "
      << static_cast<unsigned int>(mData.readChar()) << endl;
  unsigned int number_of_coedge = mData.readUnsignedInt();
  cout << getIndent() << "number_of_coedge " << number_of_coedge << endl;
  indent();
  for(unsigned int i = 0; i < number_of_coedge; ++i)
  {
    describeObject(mData);
    cout << getIndent() << "neigh_serial_index "
        << mData.readUnsignedInt() << endl;
  }
  dedent();

  dedent();
}

void describeTopoCoEdge(BitByBitData &mData)
{
  cout << getIndent() << "--CoEdge--" << endl;
  indent();

  describeBaseTopology(mData);

  describeObject(mData); // edge
  describeObject(mData); // uv_curve
  cout << getIndent() << "orientation_with_loop "
      << static_cast<unsigned int>(mData.readChar()) << endl;
  cout << getIndent() << "orientation_uv_with_loop "
      << static_cast<unsigned int>(mData.readChar()) << endl;
  dedent();
}

void describeTopoEdge(BitByBitData &mData)
{
  cout << getIndent() << "--Edge--" << endl;
  indent();

  describeContentWireEdge(mData);

  describeObject(mData); // vertex_start
  describeObject(mData); // vertex_end

  bool have_tolerance = mData.readBit();
  cout << getIndent() << "have_tolerance "
      << (have_tolerance?"yes":"no") << endl;
  if(have_tolerance)
    cout << getIndent() << "tolerance " << mData.readDouble() << endl;
  dedent();
}

void describeTopoUniqueVertex(BitByBitData &mData)
{
  cout << getIndent() << "--Unique Vertex--" << endl;
  indent();

  describeBaseTopology(mData);
  describeVector3d(mData);

  bool have_tolerance = mData.readBit();
  cout << getIndent() << "have_tolerance "
      << (have_tolerance?"yes":"no") << endl;
  if(have_tolerance)
    cout << getIndent() << "tolerance " << mData.readDouble() << endl;

  dedent();
}

void describeTopoConnex(BitByBitData &mData)
{
  cout << getIndent() << "--Connex--" << endl;
  indent();
  describeBaseTopology(mData);
  unsigned int number_of_shells = mData.readUnsignedInt();
  cout << getIndent() << "number_of_shells "
      << number_of_shells << endl;
  indent();
  for(unsigned int i = 0; i < number_of_shells; ++i)
  {
    // NOTE: this does not check if the objects are actually shells!
    describeObject(mData);
  }
  dedent();

  dedent();
}

void describeTopoShell(BitByBitData &mData)
{
  cout << getIndent() << "--Shell--" << endl;
  indent();

  describeBaseTopology(mData);

  cout << getIndent() << "shell_is_closed "
      << (mData.readBit()?"yes":"no") << endl;

  unsigned int number_of_faces = mData.readUnsignedInt();
  cout << getIndent() << "number_of_faces " << number_of_faces << endl;
  for(unsigned int i = 0; i < number_of_faces; ++i)
  {
    // NOTE: this does not check if the objects are actually faces!
    describeObject(mData);
    unsigned char orientation = mData.readChar();
    cout << getIndent() << "orientation_surface_with_shell "
        << static_cast<unsigned int>(orientation) << endl;
  }

  dedent();
}

void describeObject(BitByBitData &mData)
{
  cout << getIndent() << "--Object--" << endl;
  bool already_stored = mData.readBit();
  cout << getIndent() << "already_stored "
      << (already_stored?"yes":"no") << endl;
  if(already_stored) // reverse of documentation?
  {
    cout << getIndent() << "index of stored item "
        << mData.readUnsignedInt() << endl;
  }
  else
  {
    unsigned int type = mData.readUnsignedInt();
    switch(type)
    {
      case PRC_TYPE_ROOT:
        cout << getIndent() << "NULL Object" << endl;
        break;
      // topological items
      case PRC_TYPE_TOPO_Connex:
        describeTopoConnex(mData);
        break;
      case PRC_TYPE_TOPO_Shell:
        describeTopoShell(mData);
        break;
      case PRC_TYPE_TOPO_Face:
        describeTopoFace(mData);
        break;
      case PRC_TYPE_TOPO_Loop:
        describeTopoLoop(mData);
        break;
      case PRC_TYPE_TOPO_CoEdge:
        describeTopoCoEdge(mData);
        break;
      case PRC_TYPE_TOPO_Edge:
        describeTopoEdge(mData);
        break;
      case PRC_TYPE_TOPO_UniqueVertex:
        describeTopoUniqueVertex(mData);
        break;
      case PRC_TYPE_TOPO_WireEdge:
        describeContentWireEdge(mData);
        break;

      // curves
      case PRC_TYPE_CRV_Circle:
        describeCurvCircle(mData);
        break;
      case PRC_TYPE_CRV_NURBS:
        describeCurvNURBS(mData);
        break;
      case PRC_TYPE_CRV_PolyLine:
        describeCurvPolyLine(mData);
        break;
      case PRC_TYPE_CRV_Line:
        describeCurvLine(mData);
        break;
      // surfaces
      case PRC_TYPE_SURF_NURBS:
        describeSurfNURBS(mData);
        break;
      case PRC_TYPE_SURF_Cylinder:
        describeSurfCylinder(mData);
        break;
      case PRC_TYPE_SURF_Plane:
        describeSurfPlane(mData);
        break;

      // topological items
      case PRC_TYPE_TOPO_Item:
      case PRC_TYPE_TOPO_MultipleVertex:


      // curves
      case PRC_TYPE_CRV_Base:
      case PRC_TYPE_CRV_Blend02Boundary:
      case PRC_TYPE_CRV_Composite:
      case PRC_TYPE_CRV_OnSurf:
      case PRC_TYPE_CRV_Ellipse:
      case PRC_TYPE_CRV_Equation:
      case PRC_TYPE_CRV_Helix:
      case PRC_TYPE_CRV_Hyperbola:
      case PRC_TYPE_CRV_Intersection:
      case PRC_TYPE_CRV_Offset:
      case PRC_TYPE_CRV_Parabola:
      case PRC_TYPE_CRV_Transform:

      // surfaces
      case PRC_TYPE_SURF_Base:
      case PRC_TYPE_SURF_Blend01:
      case PRC_TYPE_SURF_Blend02:
      case PRC_TYPE_SURF_Blend03:
      case PRC_TYPE_SURF_Cone:
      case PRC_TYPE_SURF_Cylindrical:
      case PRC_TYPE_SURF_Offset:
      case PRC_TYPE_SURF_Pipe:
      case PRC_TYPE_SURF_Ruled:
      case PRC_TYPE_SURF_Sphere:
      case PRC_TYPE_SURF_Revolution:
      case PRC_TYPE_SURF_Extrusion:
      case PRC_TYPE_SURF_FromCurves:
      case PRC_TYPE_SURF_Torus:
      case PRC_TYPE_SURF_Transform:
      case PRC_TYPE_SURF_Blend04:
        cout << getIndent() << "TODO: Unhandled object of type "
            << type << endl;
        break;
      default:
        cout << getIndent() << "Invalid object of type " << type << endl;
        break;
    }
  }
}

void describeBaseTopology(BitByBitData &mData)
{
  bool base_information = mData.readBit();
  cout << getIndent() << "base_information " <<
      (base_information?"yes":"no") << endl;
  if(base_information)
  {
    describeAttributes(mData);
    describeName(mData);
    cout << getIndent() << "identifier " << mData.readUnsignedInt() << endl;
  }
}

void describeBaseGeometry(BitByBitData &mData)
{
  bool base_information = mData.readBit();
  cout << getIndent() << "base_information " <<
      (base_information?"yes":"no") << endl;
  if(base_information)
  {
    describeAttributes(mData);
    describeName(mData);
    cout << getIndent() << "identifier " << mData.readUnsignedInt() << endl;
  }
}

unsigned int describeContentBody(BitByBitData &mData)
{
  describeBaseTopology(mData);
  unsigned int behaviour = static_cast<unsigned int>(mData.readChar());
  cout << getIndent() << "behaviour " << behaviour << endl;
  return behaviour;
}

void describeContentSurface(BitByBitData &mData)
{
  describeBaseGeometry(mData);
  cout << getIndent() << "extend_info " << mData.readUnsignedInt() << endl;
}

void describeBody(BitByBitData &mData)
{
  cout << getIndent() << "--Body--" << endl;
  unsigned int type = mData.readUnsignedInt();
  switch(type)
  {
    case PRC_TYPE_TOPO_BrepData:
    {
      cout << getIndent() << "--PRC_TYPE_TOPO_BrepData--" << endl;
      unsigned int behaviour = describeContentBody(mData);

      unsigned int number_of_connex = mData.readUnsignedInt();
      cout << getIndent() << "number_of_connex " << number_of_connex << endl;
      indent();
      for(unsigned int i = 0; i < number_of_connex; ++i)
      {
        describeObject(mData);
      }
      dedent();
      if(behaviour != 0)
      {
        cout << getIndent() << "bbox " << endl;
        indent();
        describeExtent3d(mData);
        dedent();
      }
      break;
    }
    case PRC_TYPE_TOPO_SingleWireBody:
    {
      cout << getIndent() << "--PRC_TYPE_TOPO_SingleWireBody--" << endl;
      // unsigned int behaviour = describeContentBody(mData);
      // TODO: is behaviour needed to get data about how to describe?
      describeContentBody(mData);

      describeObject(mData); // curve

      break;
    }
    case PRC_TYPE_TOPO_BrepDataCompress:
    case PRC_TYPE_TOPO_SingleWireBodyCompress:
      cout << getIndent() << "TODO: Unhandled body type " << type << endl;
      break;
    default:
      cout << getIndent() << "Invalid body type " << type << endl;
      break;
  }
}

void describeTopoContext(BitByBitData &mData)
{
  cout << getIndent() << "--Topological Context--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_TOPO_Context))
    return;
  indent();

  describeContentPRCBase(mData,false);

  cout << getIndent() << "behaviour "
      << static_cast<unsigned int>(mData.readChar()) << endl;
  cout << getIndent() << "granularity " << mData.readDouble() << endl;
  cout << getIndent() << "tolerance " << mData.readDouble() << endl;

  bool have_smallest_face_thickness = mData.readBit();
  cout << getIndent() << "have_smallest_face_thickness "
      << (have_smallest_face_thickness?"yes":"no") << endl;
  if(have_smallest_face_thickness)
    cout << getIndent() << "smallest_thickness " << mData.readDouble() << endl;

  bool have_scale = mData.readBit();
  cout << getIndent() << "have_scale " << (have_scale?"yes":"no") << endl;
  if(have_scale)
    cout << getIndent() << "scale " << mData.readDouble() << endl;

  dedent();
}

void describeLineAttr(BitByBitData& mData)
{
  cout << getIndent() << "index_of_line_style "
      << mData.readUnsignedInt()-1 << endl;
}

void describeArrayRGBA(BitByBitData& mData, int number_of_colours,
                       int number_by_vector)
{
  // bool new_colour = true; // not currently used
  for(int i = 0; i < number_by_vector; ++i)
  {
    cout << getIndent() << static_cast<unsigned int>(mData.readChar()) << ' ';
    cout << static_cast<unsigned int>(mData.readChar()) << ' ';
    cout << static_cast<unsigned int>(mData.readChar()) << endl;
    //TODO: finish this
  }
}

void describeContentBaseTessData(BitByBitData &mData)
{
  cout << getIndent() << "is_calculated "
      << (mData.readBit()?"yes":"no") << endl;
  unsigned int number_of_coordinates = mData.readUnsignedInt();
  cout << getIndent() << "number_of_coordinates "
      << number_of_coordinates << endl;
  indent();
  for(unsigned int i = 0; i < number_of_coordinates; ++i)
  {
    cout << getIndent() << mData.readDouble() << endl;
  }
  dedent();
}

void describeTessFace(BitByBitData &mData)
{
  cout << getIndent() << "--Tessellation Face--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_TESS_Face))
    return;
  indent();

  unsigned int size_of_line_attributes = mData.readUnsignedInt();
  cout << getIndent() << "size_of_line_attributes "
      << size_of_line_attributes << endl;
  indent();
  for(unsigned int i = 0; i < size_of_line_attributes; ++i)
  {
    describeLineAttr(mData);
  }
  dedent();

  unsigned int start_wire = mData.readUnsignedInt();
  cout << getIndent() << "start_wire " << start_wire << endl;
  unsigned int size_of_sizes_wire = mData.readUnsignedInt();
  cout << getIndent() << "size_of_sizes_wire " << size_of_sizes_wire << endl;
  indent();
  for(unsigned int i = 0; i < size_of_sizes_wire; ++i)
  {
    cout << getIndent() << mData.readUnsignedInt() << endl;
  }
  dedent();


  unsigned int used_entities_flag = mData.readUnsignedInt();
  cout << getIndent() << "used_entities_flag " << used_entities_flag << endl;

  unsigned int start_triangulated = mData.readUnsignedInt();
  cout << getIndent() << "start_triangulated " << start_triangulated << endl;
  unsigned int size_of_sizes_triangulated = mData.readUnsignedInt();
  cout << getIndent() << "size_of_sizes_triangulated "
      << size_of_sizes_triangulated << endl;
  indent();
  for(unsigned int i = 0; i < size_of_sizes_triangulated; ++i)
  {
    cout << getIndent() << mData.readUnsignedInt() << endl;
  }
  dedent();

  cout << getIndent() << "number_of_texture_coordinate_indexes "
      << mData.readUnsignedInt() << endl;

  bool has_vertex_colors = mData.readBit();
  cout << getIndent() << "has_vertex_colors "
      << (has_vertex_colors?"yes":"no") << endl;
  indent();
  if(has_vertex_colors)
  {
    bool is_rgba = mData.readBit();
    cout << getIndent() << "is_rgba " << (is_rgba?"yes":"no") << endl;

    bool b_optimised = mData.readBit();
    cout << getIndent() << "b_optimised " << (b_optimised?"yes":"no") << endl;
    if(!b_optimised)
    {
      indent();
      //TODO: compute size of Array and pass it instead of 0
      describeArrayRGBA(mData,0,(is_rgba ? 4 : 3));
      dedent();
    }
    else
    {
      // not described
      // what does this mean? that this should not happen?
      // or is omitted from the docs?
    }
  }
  dedent();

  if(size_of_line_attributes)
  {
    cout << getIndent() << "behaviour " << mData.readUnsignedInt() << endl;
  }

  dedent();
}

void describe3DTess(BitByBitData &mData)
{
  cout << getIndent() << "--3D Tessellation--" << endl;
  indent();

  describeContentBaseTessData(mData);

  cout << getIndent() << "has_faces " << (mData.readBit()?"yes":"no") << endl;
  cout << getIndent() << "has_loops " << (mData.readBit()?"yes":"no") << endl;

  bool must_recalculate_normals = mData.readBit();
  cout << getIndent() << "must_recalculate_normals "
      << (must_recalculate_normals?"yes":"no") << endl;
  indent();
  if(must_recalculate_normals)
  {
    cout << getIndent()
        << "Docs were wrong: must_recalculate_normals is true." << endl;
    cout << getIndent() << "normals_recalculation_flags "
        << static_cast<unsigned int>(mData.readChar()) << endl;
    cout << getIndent() << "crease_angle " << mData.readDouble() << endl;
  }
  dedent();

  unsigned int number_of_normal_coordinates = mData.readUnsignedInt();
  cout << getIndent() << "number_of_normal_coordinates "
      << number_of_normal_coordinates << endl;
  indent();
  for(unsigned int i = 0; i < number_of_normal_coordinates; ++i)
  {
    cout << getIndent() << mData.readDouble() << endl;
  }
  dedent();

  unsigned int number_of_wire_indices = mData.readUnsignedInt();
  cout << getIndent() << "number_of_wire_indices "
      << number_of_wire_indices << endl;
  indent();
  for(unsigned int i = 0; i < number_of_wire_indices; ++i)
  {
    cout << getIndent() << mData.readUnsignedInt() << endl;
  }
  dedent();

  unsigned int number_of_triangulated_indices = mData.readUnsignedInt();
  cout << getIndent() << "number_of_triangulated_indices "
      << number_of_triangulated_indices << endl;
  indent();
  for(unsigned int i = 0; i < number_of_triangulated_indices; ++i)
  {
    cout << getIndent() << mData.readUnsignedInt() << endl;
  }
  dedent();

  unsigned int number_of_face_tessellation = mData.readUnsignedInt();
  cout << getIndent() << "number_of_face_tessellation "
      << number_of_face_tessellation << endl;
  indent();
  for(unsigned int i = 0; i < number_of_face_tessellation; ++i)
  {
    describeTessFace(mData);
  }
  dedent();

  unsigned int number_of_texture_coordinates = mData.readUnsignedInt();
  cout << getIndent() << "number_of_texture_coordinates "
      << number_of_texture_coordinates << endl;
  indent();
  for(unsigned int i = 0; i < number_of_texture_coordinates; ++i)
  {
    cout << getIndent() << mData.readDouble() << endl;
  }
  dedent();

  dedent();
}

void describe3DWireTess(BitByBitData &mData)
{
  //TODO
}

void describe3DMarkupTess(BitByBitData &mData)
{
  //TODO
}

void describeHighlyCompressed3DTess(BitByBitData &mData)
{
  //TODO
}

void describeSceneDisplayParameters(BitByBitData &mData)
{
  cout << getIndent() << "--Scene Display Parameters--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_GRAPH_SceneDisplayParameters))
    return;
  indent();
  describeContentPRCBase(mData,true);

  cout << getIndent() << "is active? " << (mData.readBit()?"yes":"no") << endl;

  unsigned int number_of_lights = mData.readUnsignedInt();
  cout << getIndent() << "number of lights " << number_of_lights << endl;
  indent();
  for(unsigned int i = 0; i < number_of_lights; ++i)
  {
    describeLight(mData);
  }
  dedent();

  bool camera = mData.readBit();
  cout << getIndent() << "camera? " << (camera?"yes":"no") << endl;
  if(camera)
    describeCamera(mData);

  bool rotation_centre = mData.readBit();
  cout << getIndent() << "rotation centre? " << (rotation_centre?"yes":"no") << endl;
  if(rotation_centre)
    describeVector3d(mData);

  unsigned int number_of_clipping_planes = mData.readUnsignedInt();
  cout << getIndent() << "number of clipping planes " << number_of_clipping_planes << endl;
  indent();
  for(unsigned int i = 0; i < number_of_clipping_planes; ++i)
  {
    cout << "Can't describe planes!!!" << endl;
    //describePlane(mData);
  }
  dedent();

  cout << getIndent() << "Background line style index: " << mData.readUnsignedInt()-1 << endl;
  cout << getIndent() << "Default line style index: " << mData.readUnsignedInt()-1 << endl;

  unsigned int number_of_default_styles_per_type = mData.readUnsignedInt();
  cout << getIndent() << "number_of_default_styles_per_type " << number_of_default_styles_per_type << endl;
  indent();
  for(unsigned int i = 0; i < number_of_default_styles_per_type; ++i)
  {
    cout << getIndent() << "type " << mData.readUnsignedInt() << endl;
    cout << getIndent() << "line style index: " << mData.readUnsignedInt()-1 << endl;
  }
  dedent();

  dedent();
}

void describeCartesionTransformation3d(BitByBitData& mData)
{
  cout << getIndent() << "--3d Cartesian Transformation--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_MISC_CartesianTransformation))
    return;
  indent();
  unsigned char behaviour = mData.readChar();
  cout << getIndent() << "behaviour "
      << static_cast<unsigned int>(behaviour) << endl;
  if((behaviour & PRC_TRANSFORMATION_Translate) != 0)
  {
    cout << getIndent() << "Translation" << endl;
    describeVector3d(mData);
  }

  if((behaviour & PRC_TRANSFORMATION_NonOrtho) != 0)
  {
    cout << getIndent() << "Non orthogonal transformation" << endl;
    cout << getIndent() << "X" << endl; describeVector3d(mData);
    cout << getIndent() << "Y" << endl; describeVector3d(mData);
    cout << getIndent() << "Z" << endl; describeVector3d(mData);
  }
  else if((behaviour & PRC_TRANSFORMATION_Rotate) != 0)
  {
    cout << getIndent() << "Rotation" << endl;
    cout << getIndent() << "X" << endl; describeVector3d(mData);
    cout << getIndent() << "Y" << endl; describeVector3d(mData);
  }

  // this is different from the docs!!! but it works...
  if ((behaviour & PRC_TRANSFORMATION_NonUniformScale) != 0)
  {
    cout << getIndent() << "Non-uniform scale by " << endl;
    describeVector3d(mData);
  }

  // this is different from the docs!!! but it works...
  if((behaviour & PRC_TRANSFORMATION_Scale) != 0)
  {
    cout << getIndent() << "Uniform Scale by " << mData.readDouble() << endl;
  }

  if((behaviour & PRC_TRANSFORMATION_Homogeneous) != 0)
  {
    cout << getIndent() << "transformation has homogenous values" << endl;
    cout << getIndent() << "x = " << mData.readDouble() << endl;
    cout << getIndent() << "y = " << mData.readDouble() << endl;
    cout << getIndent() << "z = " << mData.readDouble() << endl;
    cout << getIndent() << "w = " << mData.readDouble() << endl;
  }
  dedent();
}

void describeTransformation3d(BitByBitData& mData)
{
  cout << getIndent() << "--3d Transformation--" << endl;
  indent();
  bool has_transformation = mData.readBit();
  cout << getIndent() << "has_transformation "
      << (has_transformation?"yes":"no") << endl;
  if(has_transformation)
  {
    unsigned char behaviour = mData.readChar();
    cout << getIndent() << "behaviour "
        << static_cast<unsigned int>(behaviour) << endl;
    if((behaviour & PRC_TRANSFORMATION_Translate) != 0)
    {
      cout << getIndent() << "Translation" << endl;
      describeVector3d(mData);
    }
    if((behaviour & PRC_TRANSFORMATION_Rotate) != 0)
    {
      cout << getIndent() << "Rotation" << endl;
      cout << getIndent() << "X" << endl; describeVector3d(mData);
      cout << getIndent() << "Y" << endl; describeVector3d(mData);
    }

    if((behaviour & PRC_TRANSFORMATION_Scale) != 0)
    {
      cout << getIndent() << "Uniform Scale by " << mData.readDouble() << endl;
    }
  }
  dedent();
}

void describeTransformation2d(BitByBitData& mData)
{
  cout << getIndent() << "--2d Transformation--" << endl;
  indent();
  bool has_transformation = mData.readBit();
  cout << "has_transformation " << (has_transformation?"yes":"no") << endl;
  if(has_transformation)
  {
    unsigned char behaviour = mData.readChar();
    cout << getIndent() << "behaviour "
        << static_cast<unsigned int>(behaviour) << endl;
    if((behaviour & PRC_TRANSFORMATION_Translate) != 0)
    {
      cout << getIndent() << "Translation" << endl;
      describeVector2d(mData);
    }
    if((behaviour & PRC_TRANSFORMATION_Rotate) != 0)
    {
      cout << getIndent() << "Rotation" << endl;
      cout << getIndent() << "X" << endl; describeVector2d(mData);
      cout << getIndent() << "Y" << endl; describeVector2d(mData);
    }

    if((behaviour & PRC_TRANSFORMATION_Scale) != 0)
    {
      cout << getIndent() << "Uniform Scale by " << mData.readDouble() << endl;
    }
  }
  dedent();
}

void describeFileStructureInternalData(BitByBitData &mData)
{
  cout << getIndent() << "--File Structure Internal Data--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_ASM_FileStructure))
    return;
  indent();
  describeContentPRCBase(mData,false);
  cout << getIndent() << "next_available_index "
      << mData.readUnsignedInt() << endl;
  cout << getIndent() << "index_product_occurence "
      << mData.readUnsignedInt() << endl;
  dedent();
}

void describeProductOccurrence(BitByBitData &mData)
{
  cout << getIndent() << "--Product Occurrence--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_ASM_ProductOccurence))
    return;
  indent();

  describeContentPRCBaseWithGraphics(mData,true);

  cout << getIndent() << "index_part "
      << static_cast<int>(mData.readUnsignedInt()-1) << endl;
  unsigned int index_prototype = mData.readUnsignedInt()-1;
  cout << getIndent() << "index_prototype "
      << static_cast<int>(index_prototype) << endl;
  if(index_prototype+1 != 0)
  {
    bool prototype_in_same_file_structure = mData.readBit();
    cout << getIndent() << "prototype_in_same_file_structure "
        << (prototype_in_same_file_structure?"yes":"no") << endl;
    if(!prototype_in_same_file_structure)
      describeCompressedUniqueID(mData);
  }

  unsigned int index_external_data = mData.readUnsignedInt()-1;
  cout << getIndent() << "index_external_data "
      << static_cast<int>(index_external_data) << endl;
  if(index_external_data+1 != 0)
  {
    bool external_data_in_same_file_structure = mData.readBit();
    cout << getIndent() << "external_data_in_same_file_structure "
        << (external_data_in_same_file_structure?"yes":"no") << endl;
    if(!external_data_in_same_file_structure)
      describeCompressedUniqueID(mData);
  }

  unsigned int number_of_son_product_occurences = mData.readUnsignedInt();
  cout << getIndent() << "number_of_son_product_occurences "
      << number_of_son_product_occurences << endl;
  indent();
  for(unsigned int i = 0; i < number_of_son_product_occurences; ++i)
    cout << getIndent() << mData.readUnsignedInt() << endl;
  dedent();

  cout << getIndent() << "product_behaviour "
      << static_cast<unsigned int>(mData.readChar()) << endl;

  describeUnit(mData);
  cout << getIndent() << "Product information flags "
      << static_cast<unsigned int>(mData.readChar()) << endl;
  cout << getIndent() << "product_load_status "
      << mData.readUnsignedInt() << endl;

  bool has_location = mData.readBit();
  cout << getIndent() << "has_location " << has_location << endl;
  if(has_location)
  {
    describeCartesionTransformation3d(mData);
  }

  unsigned int number_of_references = mData.readUnsignedInt();
  cout << getIndent() << "number_of_references "
      << number_of_references << endl;
  indent();
  for(unsigned int i = 0; i < number_of_references; ++i)
  {
    //TODO: describeReferenceToPRCBase(mData);
  }
  dedent();

  describeMarkups(mData);

  unsigned int number_of_views = mData.readUnsignedInt();
  cout << getIndent() << "number_of_views " << number_of_views << endl;
  indent();
  for(unsigned int i = 0; i < number_of_views; ++i)
  {
    describeAnnotationView(mData);
  }
  dedent();

  bool has_entity_filter = mData.readBit();
  cout << getIndent() << "has_entity_filter "
      << (has_entity_filter?"yes":"no") << endl;
  if(has_entity_filter)
  {
    //TODO: describeEntityFilter(mData);
  }

  unsigned int number_of_display_filters = mData.readUnsignedInt();
  cout << getIndent() << "number_of_display_filters "
      << number_of_display_filters << endl;
  indent();
  for(unsigned int i = 0; i < number_of_display_filters; ++i)
  {
    //TODO: describeFilter(mData);
  }
  dedent();

  unsigned int number_of_scene_display_parameters = mData.readUnsignedInt();
  cout << getIndent() << "number_of_scene_display_parameters "
      << number_of_scene_display_parameters << endl;
  indent();
  for(unsigned int i = 0; i < number_of_scene_display_parameters; ++i)
  {
    describeSceneDisplayParameters(mData);
  }
  dedent();

  describeUserData(mData);
  dedent();
}

void describeGraphics(BitByBitData &mData)
{
  bool sameGraphicsAsCurrent = mData.readBit();
  cout << getIndent() << "Same graphics as current graphics? "
      << (sameGraphicsAsCurrent?"yes":"no") << endl;
  if(!sameGraphicsAsCurrent)
  {
    layer_index = mData.readUnsignedInt()-1;
    cout << getIndent() << "layer_index " << layer_index << endl;
    index_of_line_style = mData.readUnsignedInt()-1;
    cout << getIndent() << "index_of_line_style "
        << index_of_line_style << endl;
    unsigned char c1 = mData.readChar();
    unsigned char c2 = mData.readChar();
    behaviour_bit_field = c1 | (static_cast<unsigned short>(c2) << 8);
    cout << getIndent() << "behaviour_bit_field "
        << behaviour_bit_field << endl;
  }
}

void describeContentPRCBaseWithGraphics(BitByBitData &mData, bool efr)
{
  describeContentPRCBase(mData,efr);
  describeGraphics(mData);
}

void describePartDefinition(BitByBitData &mData)
{
  cout << getIndent() << "--Part Definition--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_ASM_PartDefinition))
    return;
  indent();

  describeContentPRCBaseWithGraphics(mData,true);
  describeExtent3d(mData);

  unsigned int number_of_representation_items = mData.readUnsignedInt();
  cout << getIndent() << "number_of_representation_items "
      << number_of_representation_items << endl;
  indent();
  for(unsigned int i = 0; i < number_of_representation_items; ++i)
  {
    describeRepresentationItem(mData);
  }
  dedent();

  describeMarkups(mData);

  unsigned int number_of_views = mData.readUnsignedInt();
  cout << getIndent() << "number_of_views " << number_of_views << endl;
  indent();
  for(unsigned int i = 0; i < number_of_views; ++i)
  {
    describeAnnotationView(mData);
  }
  dedent();

  describeUserData(mData);
  dedent();
}

void describeMarkups(BitByBitData& mData)
{
  cout << getIndent() << "--Markups--" << endl;
  indent();

  unsigned int number_of_linked_items = mData.readUnsignedInt();
  cout << getIndent() << "number_of_linked_items "
      << number_of_linked_items << endl;
  for(unsigned int i = 0; i < number_of_linked_items; ++i)
  {
    cout << "describe linked item!" << endl;
  }

  unsigned int number_of_leaders = mData.readUnsignedInt();
  cout << getIndent() << "number_of_leaders " << number_of_leaders << endl;
  for(unsigned int i = 0; i < number_of_leaders; ++i)
  {
    cout << "describe leader!" << endl;
  }

  unsigned int number_of_markups = mData.readUnsignedInt();
  cout << getIndent() << "number_of_markups " << number_of_markups << endl;
  for(unsigned int i=0; i < number_of_markups; ++i)
  {
    cout << "describe markup!" << endl;
  }

  unsigned int number_of_annotation_entities = mData.readUnsignedInt();
  cout << getIndent() << "number_of_annotation_entities "
      << number_of_annotation_entities << endl;
  for(unsigned int i=0; i < number_of_annotation_entities; ++i)
  {
    cout << "describe annotation entity!" << endl;
  }

  dedent();
}

void describeAnnotationView(BitByBitData &mData)
{
  cout << getIndent() << "--Annotation View--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_MKP_View))
    return;
  indent();
  describeContentPRCBaseWithGraphics(mData,true);
  unsigned int number_of_annotations = mData.readUnsignedInt();
  for(unsigned int i = 0; i < number_of_annotations; ++i)
  {
    //TODO: describeReferenceUniqueIdentifier(mData);
  }
  //TODO: describePlane(mData);
  bool scene_display_parameters = mData.readBit();
  if(scene_display_parameters)
  {
    describeSceneDisplayParameters(mData);
  }
  describeUserData(mData);
  dedent();
}

void describeExtent3d(BitByBitData &mData)
{ // I suspect the order of min/max should be flipped
  cout << getIndent() << "Minimum" << endl;
  indent(); describeVector3d(mData); dedent();
  cout << getIndent() << "Maximum" << endl;
  indent(); describeVector3d(mData); dedent();
}

void describeExtent1d(BitByBitData &mData)
{
  cout << getIndent() << "Minimum " << mData.readDouble() << endl;
  cout << getIndent() << "Maximum " << mData.readDouble() << endl;
}

void describeExtent2d(BitByBitData &mData)
{
  cout << getIndent() << "Minimum" << endl;
  indent(); describeVector2d(mData); dedent();
  cout << getIndent() << "Maximum" << endl;
  indent(); describeVector2d(mData); dedent();
}

void describeVector3d(BitByBitData &mData)
{
  double x = mData.readDouble();
  double y = mData.readDouble();
  double z = mData.readDouble();
  cout << getIndent() << '(' << x << ',' << y << ',' << z << ')' << endl;
}

void describeVector2d(BitByBitData &mData)
{
  double x = mData.readDouble();
  double y = mData.readDouble();
  cout << getIndent() << '(' << x << ',' << y << ')' << endl;
}

void describePicture(BitByBitData &mData)
{
  cout << getIndent() << "--Picture--" << endl;
  unsigned int sectionCode = mData.readUnsignedInt();
  if(sectionCode != PRC_TYPE_GRAPH_Picture)
  {
    cout << getIndent() << "Invalid section code." << endl;
  }

  describeContentPRCBase(mData,false);

  int format = mData.readInt();
  switch(format)
  {
    case KEPRCPicture_PNG:
      cout << getIndent() << "PNG format" << endl;
      break;
    case KEPRCPicture_JPG:
      cout << getIndent() << "JPG format" << endl;
      break;
    case KEPRCPicture_BITMAP_RGB_BYTE:
      cout << getIndent() << "gzipped pixel data (see PRC base compression). Each element is a RGB triple. (3 components)" << endl;
      break;
    case KEPRCPicture_BITMAP_RGBA_BYTE:
      cout << getIndent() << "gzipped pixel data (see PRC base compression). Each element is a complete RGBA element. (4 components)" << endl;
      break;
    case KEPRCPicture_BITMAP_GREY_BYTE:
      cout << getIndent() << "gzipped pixel data (see PRC base compression). Each element is a single luminance value. (1 components)" << endl;
      break;
    case KEPRCPicture_BITMAP_GREYA_BYTE:
      cout << getIndent() << "gzipped pixel data (see PRC base compression). Each element is a luminance/alpha pair. (2 components)" << endl;
      break;
    default:
      cout << getIndent() << "Invalid picture format." << endl;
      break;
  }
  cout << getIndent() << "uncompressed_file_index "
      << mData.readUnsignedInt()-1 << endl;
  cout << getIndent() << "pixel width " << mData.readUnsignedInt() << endl;
  cout << getIndent() << "pixel height " << mData.readUnsignedInt() << endl;
}

void describeTextureDefinition(BitByBitData &mData)
{
  cout << getIndent() << "--Texture Definition--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_GRAPH_TextureDefinition))
    return;

  cout << getIndent() << "TODO: Can't describe textures yet." << endl;
}

void describeMaterial(BitByBitData &mData)
{
  cout << getIndent() << "--Material--" << endl;
  unsigned int code = mData.readUnsignedInt();
  if(code == PRC_TYPE_GRAPH_Material)
  {
    describeContentPRCBase(mData,true);
    cout << getIndent() << "index of ambient color "
        << mData.readUnsignedInt() - 1 << endl;
    cout << getIndent() << "index of diffuse color "
        << mData.readUnsignedInt() - 1 << endl;
    cout << getIndent() << "index of emissive color "
        << mData.readUnsignedInt() - 1 << endl;
    cout << getIndent() << "index of specular color "
        << mData.readUnsignedInt() - 1 << endl;
    cout << getIndent() << "shininess " << mData.readDouble() << endl;
    cout << getIndent() << "ambient_alpha " << mData.readDouble() << endl;
    cout << getIndent() << "diffuse_alpha " << mData.readDouble() << endl;
    cout << getIndent() << "emissive_alpha " << mData.readDouble() << endl;
    cout << getIndent() << "specular_alpha " << mData.readDouble() << endl;
  }
  else if(code == PRC_TYPE_GRAPH_TextureApplication)
  {
    describeContentPRCBase(mData,true);
    cout << getIndent() << "material_generic_index "
        << mData.readUnsignedInt() - 1 << endl;
    cout << getIndent() << "texture_definition_index "
        << mData.readUnsignedInt() - 1 << endl;
    cout << getIndent() << "next_texture_index "
        << mData.readUnsignedInt() - 1 << endl;
    cout << getIndent() << "UV_coordinates_index "
        << mData.readUnsignedInt() - 1 << endl;
  }
  else
  {
    cout << getIndent() << "Invalid section code in material definition."
        << endl;
  }
}

void describeLinePattern(BitByBitData &mData)
{
  cout << getIndent() << "--Line Pattern--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_GRAPH_LinePattern))
    return;
  indent();

  describeContentPRCBase(mData,true);
  unsigned int size_lengths = mData.readUnsignedInt();
  cout << getIndent() << "size_lengths " << size_lengths << endl;
  indent();
  for(unsigned int i=0;i<size_lengths;i++)
  {
    cout << getIndent() << "length " << mData.readDouble() << endl;
  }
  dedent();
  cout << getIndent() << "phase " << mData.readDouble() << endl;
  cout << getIndent() << "is real length "
      << (mData.readBit()?"yes":"no") << endl;

  dedent();
}

void describeCategory1LineStyle(BitByBitData &mData)
{
  cout << getIndent() << "--Category 1 Line Style--" << endl;
  if(!checkSectionCode(mData,PRC_TYPE_GRAPH_Style))
    return;
  indent();

  // This was missing from the docs!!!!! argh!
  describeContentPRCBase(mData,true);


  cout << getIndent() << "line_width " << mData.readDouble() << " mm" << endl;

  cout << getIndent() << "is_vpicture " << (mData.readBit()?"yes":"no") << endl;

  cout << getIndent() << "line_pattern_index/vpicture_index "
      << static_cast<int>(mData.readUnsignedInt()-1) << endl;
  cout << getIndent() << "is_material " << (mData.readBit()?"yes":"no") << endl;
  cout << getIndent() << "color_index / material_index "
      << static_cast<int>(mData.readUnsignedInt()-1) << endl;

  bool is_transparency_defined = mData.readBit();
  cout << getIndent() << "is_transparency_defined "
      << (is_transparency_defined?"yes":"no") << endl;
  if(is_transparency_defined)
  {
    indent();
    cout << getIndent() << "transparency "
        << static_cast<unsigned int>(mData.readChar()) << endl;
    dedent();
  }

  bool is_additional_1_defined = mData.readBit();
  cout << getIndent() << "is_additional_1_defined "
      << (is_additional_1_defined?"yes":"no") << endl;
  if(is_additional_1_defined)
  {
    indent();
    cout << getIndent() << "additional_1 "
        << static_cast<unsigned int>(mData.readChar()) << endl;
    dedent();
  }

  bool is_additional_2_defined = mData.readBit();
  cout << getIndent() << "is_additional_2_defined "
      << (is_additional_2_defined?"yes":"no") << endl;
  if(is_additional_2_defined)
  {
    indent();
    cout << getIndent() << "additional_2 "
        << static_cast<unsigned int>(mData.readChar()) << endl;
    dedent();
  }

  bool is_additional_3_defined = mData.readBit();
  cout << getIndent() << "is_additional_3_defined "
      << (is_additional_3_defined?"yes":"no") << endl;
  if(is_additional_3_defined)
  {
    indent();
    cout << getIndent() << "additional_3 "
        << static_cast<unsigned int>(mData.readChar()) << endl;
    dedent();
  }
  dedent();
}

void describeFillPattern(BitByBitData &mData)
{
  cout << getIndent() << "--Fill Pattern--" << endl;
  unsigned int type = mData.readUnsignedInt();
  cout << getIndent() << "type " << type << endl;
  switch(type)
  {
    //TODO: actually describe fill patterns
    default:
      cout << getIndent() << "Invalid fill pattern type " << type << endl;
  }
}

void describeRepresentationItemContent(BitByBitData &mData)
{
  describeContentPRCBaseWithGraphics(mData,true);
  unsigned int index_local_coordinate_system = mData.readUnsignedInt()-1;
  unsigned int index_tessellation = mData.readUnsignedInt()-1;
  //cast to int will not be right for big indices
  cout << getIndent() << "index_local_coordinate_system "
      << static_cast<int>(index_local_coordinate_system) << endl;
  cout << getIndent() << "index_tessellation "
      << static_cast<int>(index_tessellation) << endl;
}

void describeRepresentationItem(BitByBitData &mData)
{
  cout << getIndent() << "--Representation Item--" << endl;
  unsigned int type = mData.readUnsignedInt();
  switch(type)
  {
    case PRC_TYPE_RI_Curve:
    {
      cout << getIndent() << "--PRC_TYPE_RI_Curve--" << endl;
      describeRepresentationItemContent(mData);
      bool has_wire_body = mData.readBit();
      if(has_wire_body)
      {
        cout << getIndent() << "context_id " << mData.readUnsignedInt() << endl;
        cout << getIndent() << "body_id " << mData.readUnsignedInt() << endl;
      }
      describeUserData(mData);
      break;
    }
    case PRC_TYPE_RI_PolyBrepModel:
    {
      cout << getIndent() << "--PRC_TYPE_RI_PolyBrepModel--" << endl;
      describeRepresentationItemContent(mData);
      cout << getIndent() << "is_closed "
          << (mData.readBit()?"yes":"no") << endl;
      describeUserData(mData);
      break;
    }
    case PRC_TYPE_RI_BrepModel:
    {
      cout << getIndent() << "--PRC_TYPE_RI_BrepModel--" << endl;
      describeRepresentationItemContent(mData);
      bool has_brep_data = mData.readBit();
      cout << getIndent() << "has_brep_data "
          << (has_brep_data?"yes":"no") << endl;
      if(has_brep_data)
      {
        cout << getIndent() << "context_id " << mData.readUnsignedInt() << endl;
        cout << getIndent() << "object_id " << mData.readUnsignedInt() << endl;
      }
      cout << getIndent() << "is_closed "
          << (mData.readBit()?"yes":"no") << endl;
      describeUserData(mData);
      break;
    }
    case PRC_TYPE_RI_Direction:

    case PRC_TYPE_RI_Plane:
    case PRC_TYPE_RI_CoordinateSystem:
    case PRC_TYPE_RI_PointSet:
    case PRC_TYPE_RI_Set:
    case PRC_TYPE_RI_PolyWire:
      cout << getIndent() << "TODO: Unhandled representation item "
          << type << endl;
      break;
    default:
      cout << getIndent() << "Invalid representation item type "
          << type << endl;
      break;
  }
}

void describeRGBColour(BitByBitData &mData)
{
  cout << getIndent() << "R: " << mData.readDouble();
  cout << " G: " << mData.readDouble();
  cout << " B: " << mData.readDouble() << endl;
}

void describeSchema(BitByBitData &mData)
{
  cout << getIndent() << "--Schema--" << endl;
  indent();
  unsigned int numSchemas = mData.readUnsignedInt();
  cout << getIndent() << "Number of Schemas " << numSchemas << endl;
  if(numSchemas != 0)
  {
    cout << "Error: Don't know how to handle multiple schemas." << endl;
  }
  dedent();
}

string currentName;
int layer_index;
int index_of_line_style;
unsigned short behaviour_bit_field;

void resetCurrentGraphics()
{
  layer_index = -1;
  index_of_line_style = -1;
  behaviour_bit_field = 1;
}

void unFlushSerialization()
{
  currentName = "";

  resetCurrentGraphics();
}

void describeName(BitByBitData &mData)
{
  bool sameNameAsCurrent = mData.readBit();
  cout << getIndent() << "Same name as current name? "
      << (sameNameAsCurrent?"yes":"no") << endl;
  if(!sameNameAsCurrent)
    currentName = mData.readString();
  cout << getIndent() << "Name \"" << currentName << '\"' << endl;
}

void describeUnit(BitByBitData &mData)
{
  cout << getIndent() << "Unit is from CAD file? "
      << (mData.readBit()?"yes":"no") << endl;
  cout << getIndent() << "Unit is " << mData.readDouble() << " mm" << endl;
}

void describeAttributes(BitByBitData &mData)
{
  cout << getIndent() << "--Attributes--" << endl;
  indent();

  unsigned int numAttribs = mData.readUnsignedInt();
  cout << getIndent() << "Number of Attributes " << numAttribs << endl;
  indent();
  for(unsigned int i = 0; i < numAttribs; ++i)
  {
    cout << getIndent() << "PRC_TYPE_MISC_Attribute "
        << mData.readUnsignedInt() << endl;
    bool titleIsInt = mData.readBit();
    cout << getIndent() << "Title is integer? "
        << (titleIsInt?"yes":"no") << endl;
    indent();
    if(titleIsInt)
    {
      cout << getIndent() << "Title " << mData.readUnsignedInt() << endl;
    }
    else
    {
      cout << getIndent() << "Title \"" << mData.readString() << '\"' << endl;
    }
    unsigned int sizeOfAttributeKeys = mData.readUnsignedInt();
    cout << getIndent() << "Size of Attribute Keys "
        << sizeOfAttributeKeys << endl;
    for(unsigned int a = 0; a < sizeOfAttributeKeys; ++a)
    {
      bool titleIsInt = mData.readBit();
      cout << getIndent() << "Title is integer? "
          << (titleIsInt?"yes":"no") << endl;
      indent();
      if(titleIsInt)
      {
        cout << getIndent() << "Title " << mData.readUnsignedInt() << endl;
      }
      else
      {
        cout << getIndent() << "Title \"" << mData.readString() << '\"' << endl;
      }
      dedent();
      unsigned int attributeType = mData.readUnsignedInt();
      cout << getIndent() << "Attribute Type " << attributeType << endl;
      switch(attributeType)
      {
        case KEPRCModellerAttributeTypeInt:
          cout << getIndent() << "Attribute Value (int) "
              << mData.readInt() << endl;
          break;
        case KEPRCModellerAttributeTypeReal:
          cout << getIndent() << "Attribute Value (double) "
              << mData.readDouble() << endl;
          break;
        case KEPRCModellerAttributeTypeTime:
          cout << getIndent() << "Attribute Value (time_t) "
              << mData.readUnsignedInt() << endl;
          break;
        case KEPRCModellerAttributeTypeString:
          cout << getIndent() << "Attribute Value (string) \""
              << mData.readString() << '\"' << endl;
          break;
        default:
          break;
      }
    }
    dedent();

    cout << endl;
  }
  dedent();


  dedent();
}

void describeContentPRCBase(BitByBitData &mData, bool typeEligibleForReference)
{
  cout << getIndent() << "--ContentPRCBase--" << endl;
  indent();
  describeAttributes(mData);
  describeName(mData);
  if(typeEligibleForReference)
  {
    cout << getIndent() << "CAD_identifier "
        << mData.readUnsignedInt() << endl;
    cout << getIndent() << "CAD_persistent_identifier "
        << mData.readUnsignedInt() << endl;
    cout << getIndent() << "PRC_unique_identifier "
        << mData.readUnsignedInt() << endl;
  }
  dedent();
}

void describeCompressedUniqueID(BitByBitData &mData)
{
  cout << getIndent() << "UUID: " << hex << setfill('0');
  for(int i = 0; i < 4; ++i)
    cout << setw(8) << mData.readUnsignedInt() << ' ';
  cout << dec << setfill(' ') << endl;
}

void describeUserData(BitByBitData &mData)
{
  unsigned int bits = mData.readUnsignedInt();
  cout << getIndent() << bits << " bits of user data" << endl;
  indent();
  for(unsigned int i = 0; i < bits; ++i)
  {
    if(i%64 == 0)
      cout << getIndent();
    cout << mData.readBit();
    if(i%64 == 63)
      cout << endl;
  }
  if(bits%64 != 0)
    cout << endl;
  dedent();
}

bool checkSectionCode(BitByBitData &mData, unsigned int code)
{
  unsigned int num = mData.readUnsignedInt();
  if(code != num)
  {
    cout << getIndent() << "Invalid section code " << num <<
        ". Expected " << code << " at "; mData.tellPosition();
    return false;
  }
  else
  {
    cout << getIndent() << "Section code " << code << endl;
    return true;
  }
}

unsigned int currentIndent = 0;

string getIndent()
{
  ostringstream out;
  for(unsigned int i = 0; i < currentIndent; ++i)
    out << "  ";
  return out.str();
}

void indent()
{
  ++currentIndent;
}

void dedent()
{
  --currentIndent;
}