path: root/Build/source/utils/asymptote/prc/test.cc

/************
*
*   This file is part of a tool for producing 3D content in the PRC format.
*   Copyright (C) 2008  Orest Shardt <shardtor (at) gmail dot com> and
*                       Michail Vidiassov <master@iaas.msu.ru>
*
*   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 <iomanip>
#include <fstream>
#include <cmath>
#include <string>
#include "oPRCFile.h"

using namespace std;

extern const double pi;

int main()
{
  #if 0
  // List of pictures used; keep track of memory allocated to free it in the end
  // shared pointers or garbage collector may be an alternative
  uint8_t *picture1 = NULL;
  uint8_t *picture2 = NULL;
  uint8_t *picture3 = NULL;
  uint8_t *picture4 = NULL;

  oPRCFile file("test.prc");

  const size_t N_COLOURS = 32;
  RGBAColour colours[N_COLOURS];
  for(size_t i = 0; i < N_COLOURS; ++i)
  {
    colours[i%N_COLOURS].R = 0.0;
    colours[i%N_COLOURS].G = (i%N_COLOURS)/static_cast<double>(N_COLOURS);
    colours[i%N_COLOURS].B = 0.95;
    colours[i%N_COLOURS].A = 0.75;
  }
  
  PRCmaterial materials[N_COLOURS];
  for(size_t i = 0; i < N_COLOURS; ++i)
  {
    materials[i%N_COLOURS].diffuse.R = 0.0;
    materials[i%N_COLOURS].diffuse.G = (i%N_COLOURS)/static_cast<double>(N_COLOURS);
    materials[i%N_COLOURS].diffuse.B = 0.95;
    materials[i%N_COLOURS].diffuse.A = 0.75;
    materials[i%N_COLOURS].specular.R = 0.01*0.0;
    materials[i%N_COLOURS].specular.G = 0.01*(i%N_COLOURS)/static_cast<double>(N_COLOURS);
    materials[i%N_COLOURS].specular.B = 0.01*0.95;
    materials[i%N_COLOURS].specular.A = 0.01*0.75;
    materials[i%N_COLOURS].emissive.R = 0.20*0.0;
    materials[i%N_COLOURS].emissive.G = 0.20*(i%N_COLOURS)/static_cast<double>(N_COLOURS);
    materials[i%N_COLOURS].emissive.B = 0.20*0.95;
    materials[i%N_COLOURS].emissive.A = 0.20*0.75;
    materials[i%N_COLOURS].ambient.R  = 0.05*0.0;
    materials[i%N_COLOURS].ambient.G  = 0.05*(i%N_COLOURS)/static_cast<double>(N_COLOURS);
    materials[i%N_COLOURS].ambient.B  = 0.05*0.95;
    materials[i%N_COLOURS].ambient.A  = 0.05*0.75;
    materials[i%N_COLOURS].alpha      = 0.75;
    materials[i%N_COLOURS].shininess  = 0.1;
  }
  
  if(1) {
  double knotsU[] = {1,1,1,1,2,2,2,2};
  double knotsV[] = {1,1,1,1,2,2,2,2};
  const size_t NUMBER_OF_PATCHES = 32;
  double controlPoints[NUMBER_OF_PATCHES][16][3] = 
  {
    { // Patch 0
      {1.4,0,2.4},{1.4,-0.784,2.4},{0.784,-1.4,2.4},{0,-1.4,2.4},
      {1.3375,0,2.53125},{1.3375,-0.749,2.53125},{0.749,-1.3375,2.53125},{0,-1.3375,2.53125},
      {1.4375,0,2.53125},{1.4375,-0.805,2.53125},{0.805,-1.4375,2.53125},{0,-1.4375,2.53125},
      {1.5,0,2.4},{1.5,-0.84,2.4},{0.84,-1.5,2.4},{0,-1.5,2.4},
    },
    { // Patch 1
      {0,-1.4,2.4},{-0.784,-1.4,2.4},{-1.4,-0.784,2.4},{-1.4,0,2.4},
      {0,-1.3375,2.53125},{-0.749,-1.3375,2.53125},{-1.3375,-0.749,2.53125},{-1.3375,0,2.53125},
      {0,-1.4375,2.53125},{-0.805,-1.4375,2.53125},{-1.4375,-0.805,2.53125},{-1.4375,0,2.53125},
      {0,-1.5,2.4},{-0.84,-1.5,2.4},{-1.5,-0.84,2.4},{-1.5,0,2.4},
    },
    { // Patch 2
      {-1.4,0,2.4},{-1.4,0.784,2.4},{-0.784,1.4,2.4},{0,1.4,2.4},
      {-1.3375,0,2.53125},{-1.3375,0.749,2.53125},{-0.749,1.3375,2.53125},{0,1.3375,2.53125},
      {-1.4375,0,2.53125},{-1.4375,0.805,2.53125},{-0.805,1.4375,2.53125},{0,1.4375,2.53125},
      {-1.5,0,2.4},{-1.5,0.84,2.4},{-0.84,1.5,2.4},{0,1.5,2.4},
    },
    { // Patch 3
      {0,1.4,2.4},{0.784,1.4,2.4},{1.4,0.784,2.4},{1.4,0,2.4},
      {0,1.3375,2.53125},{0.749,1.3375,2.53125},{1.3375,0.749,2.53125},{1.3375,0,2.53125},
      {0,1.4375,2.53125},{0.805,1.4375,2.53125},{1.4375,0.805,2.53125},{1.4375,0,2.53125},
      {0,1.5,2.4},{0.84,1.5,2.4},{1.5,0.84,2.4},{1.5,0,2.4},
    },
    { // Patch 4
      {1.5,0,2.4},{1.5,-0.84,2.4},{0.84,-1.5,2.4},{0,-1.5,2.4},
      {1.75,0,1.875},{1.75,-0.98,1.875},{0.98,-1.75,1.875},{0,-1.75,1.875},
      {2,0,1.35},{2,-1.12,1.35},{1.12,-2,1.35},{0,-2,1.35},
      {2,0,0.9},{2,-1.12,0.9},{1.12,-2,0.9},{0,-2,0.9},
    },
    { // Patch 5
      {0,-1.5,2.4},{-0.84,-1.5,2.4},{-1.5,-0.84,2.4},{-1.5,0,2.4},
      {0,-1.75,1.875},{-0.98,-1.75,1.875},{-1.75,-0.98,1.875},{-1.75,0,1.875},
      {0,-2,1.35},{-1.12,-2,1.35},{-2,-1.12,1.35},{-2,0,1.35},
      {0,-2,0.9},{-1.12,-2,0.9},{-2,-1.12,0.9},{-2,0,0.9},
    },
    { // Patch 6
      {-1.5,0,2.4},{-1.5,0.84,2.4},{-0.84,1.5,2.4},{0,1.5,2.4},
      {-1.75,0,1.875},{-1.75,0.98,1.875},{-0.98,1.75,1.875},{0,1.75,1.875},
      {-2,0,1.35},{-2,1.12,1.35},{-1.12,2,1.35},{0,2,1.35},
      {-2,0,0.9},{-2,1.12,0.9},{-1.12,2,0.9},{0,2,0.9},
    },
    { // Patch 7
      {0,1.5,2.4},{0.84,1.5,2.4},{1.5,0.84,2.4},{1.5,0,2.4},
      {0,1.75,1.875},{0.98,1.75,1.875},{1.75,0.98,1.875},{1.75,0,1.875},
      {0,2,1.35},{1.12,2,1.35},{2,1.12,1.35},{2,0,1.35},
      {0,2,0.9},{1.12,2,0.9},{2,1.12,0.9},{2,0,0.9},
    },
    { // Patch 8
      {2,0,0.9},{2,-1.12,0.9},{1.12,-2,0.9},{0,-2,0.9},
      {2,0,0.45},{2,-1.12,0.45},{1.12,-2,0.45},{0,-2,0.45},
      {1.5,0,0.225},{1.5,-0.84,0.225},{0.84,-1.5,0.225},{0,-1.5,0.225},
      {1.5,0,0.15},{1.5,-0.84,0.15},{0.84,-1.5,0.15},{0,-1.5,0.15},
    },
    { // Patch 9
      {0,-2,0.9},{-1.12,-2,0.9},{-2,-1.12,0.9},{-2,0,0.9},
      {0,-2,0.45},{-1.12,-2,0.45},{-2,-1.12,0.45},{-2,0,0.45},
      {0,-1.5,0.225},{-0.84,-1.5,0.225},{-1.5,-0.84,0.225},{-1.5,0,0.225},
      {0,-1.5,0.15},{-0.84,-1.5,0.15},{-1.5,-0.84,0.15},{-1.5,0,0.15},
    },
    { // Patch 10
      {-2,0,0.9},{-2,1.12,0.9},{-1.12,2,0.9},{0,2,0.9},
      {-2,0,0.45},{-2,1.12,0.45},{-1.12,2,0.45},{0,2,0.45},
      {-1.5,0,0.225},{-1.5,0.84,0.225},{-0.84,1.5,0.225},{0,1.5,0.225},
      {-1.5,0,0.15},{-1.5,0.84,0.15},{-0.84,1.5,0.15},{0,1.5,0.15},
    },
    { // Patch 11
      {0,2,0.9},{1.12,2,0.9},{2,1.12,0.9},{2,0,0.9},
      {0,2,0.45},{1.12,2,0.45},{2,1.12,0.45},{2,0,0.45},
      {0,1.5,0.225},{0.84,1.5,0.225},{1.5,0.84,0.225},{1.5,0,0.225},
      {0,1.5,0.15},{0.84,1.5,0.15},{1.5,0.84,0.15},{1.5,0,0.15},
    },
    { // Patch 12
      {-1.6,0,2.025},{-1.6,-0.3,2.025},{-1.5,-0.3,2.25},{-1.5,0,2.25},
      {-2.3,0,2.025},{-2.3,-0.3,2.025},{-2.5,-0.3,2.25},{-2.5,0,2.25},
      {-2.7,0,2.025},{-2.7,-0.3,2.025},{-3,-0.3,2.25},{-3,0,2.25},
      {-2.7,0,1.8},{-2.7,-0.3,1.8},{-3,-0.3,1.8},{-3,0,1.8},
    },
    { // Patch 13
      {-1.5,0,2.25},{-1.5,0.3,2.25},{-1.6,0.3,2.025},{-1.6,0,2.025},
      {-2.5,0,2.25},{-2.5,0.3,2.25},{-2.3,0.3,2.025},{-2.3,0,2.025},
      {-3,0,2.25},{-3,0.3,2.25},{-2.7,0.3,2.025},{-2.7,0,2.025},
      {-3,0,1.8},{-3,0.3,1.8},{-2.7,0.3,1.8},{-2.7,0,1.8},
    },
    { // Patch 14
      {-2.7,0,1.8},{-2.7,-0.3,1.8},{-3,-0.3,1.8},{-3,0,1.8},
      {-2.7,0,1.575},{-2.7,-0.3,1.575},{-3,-0.3,1.35},{-3,0,1.35},
      {-2.5,0,1.125},{-2.5,-0.3,1.125},{-2.65,-0.3,0.9375},{-2.65,0,0.9375},
      {-2,0,0.9},{-2,-0.3,0.9},{-1.9,-0.3,0.6},{-1.9,0,0.6},
    },
    { // Patch 15
      {-3,0,1.8},{-3,0.3,1.8},{-2.7,0.3,1.8},{-2.7,0,1.8},
      {-3,0,1.35},{-3,0.3,1.35},{-2.7,0.3,1.575},{-2.7,0,1.575},
      {-2.65,0,0.9375},{-2.65,0.3,0.9375},{-2.5,0.3,1.125},{-2.5,0,1.125},
      {-1.9,0,0.6},{-1.9,0.3,0.6},{-2,0.3,0.9},{-2,0,0.9},
    },
    { // Patch 16
      {1.7,0,1.425},{1.7,-0.66,1.425},{1.7,-0.66,0.6},{1.7,0,0.6},
      {2.6,0,1.425},{2.6,-0.66,1.425},{3.1,-0.66,0.825},{3.1,0,0.825},
      {2.3,0,2.1},{2.3,-0.25,2.1},{2.4,-0.25,2.025},{2.4,0,2.025},
      {2.7,0,2.4},{2.7,-0.25,2.4},{3.3,-0.25,2.4},{3.3,0,2.4},
    },
    { // Patch 17
      {1.7,0,0.6},{1.7,0.66,0.6},{1.7,0.66,1.425},{1.7,0,1.425},
      {3.1,0,0.825},{3.1,0.66,0.825},{2.6,0.66,1.425},{2.6,0,1.425},
      {2.4,0,2.025},{2.4,0.25,2.025},{2.3,0.25,2.1},{2.3,0,2.1},
      {3.3,0,2.4},{3.3,0.25,2.4},{2.7,0.25,2.4},{2.7,0,2.4},
    },
    { // Patch 18
      {2.7,0,2.4},{2.7,-0.25,2.4},{3.3,-0.25,2.4},{3.3,0,2.4},
      {2.8,0,2.475},{2.8,-0.25,2.475},{3.525,-0.25,2.49375},{3.525,0,2.49375},
      {2.9,0,2.475},{2.9,-0.15,2.475},{3.45,-0.15,2.5125},{3.45,0,2.5125},
      {2.8,0,2.4},{2.8,-0.15,2.4},{3.2,-0.15,2.4},{3.2,0,2.4},
    },
    { // Patch 19
      {3.3,0,2.4},{3.3,0.25,2.4},{2.7,0.25,2.4},{2.7,0,2.4},
      {3.525,0,2.49375},{3.525,0.25,2.49375},{2.8,0.25,2.475},{2.8,0,2.475},
      {3.45,0,2.5125},{3.45,0.15,2.5125},{2.9,0.15,2.475},{2.9,0,2.475},
      {3.2,0,2.4},{3.2,0.15,2.4},{2.8,0.15,2.4},{2.8,0,2.4},
    },
    { // Patch 20
      {0,0,3.15},{0,0,3.15},{0,0,3.15},{0,0,3.15},
      {0.8,0,3.15},{0.8,-0.45,3.15},{0.45,-0.8,3.15},{0,-0.8,3.15},
      {0,0,2.85},{0,0,2.85},{0,0,2.85},{0,0,2.85},
      {0.2,0,2.7},{0.2,-0.112,2.7},{0.112,-0.2,2.7},{0,-0.2,2.7},
    },
    { // Patch 21
      {0,0,3.15},{0,0,3.15},{0,0,3.15},{0,0,3.15},
      {0,-0.8,3.15},{-0.45,-0.8,3.15},{-0.8,-0.45,3.15},{-0.8,0,3.15},
      {0,0,2.85},{0,0,2.85},{0,0,2.85},{0,0,2.85},
      {0,-0.2,2.7},{-0.112,-0.2,2.7},{-0.2,-0.112,2.7},{-0.2,0,2.7},
    },
    { // Patch 22
      {0,0,3.15},{0,0,3.15},{0,0,3.15},{0,0,3.15},
      {-0.8,0,3.15},{-0.8,0.45,3.15},{-0.45,0.8,3.15},{0,0.8,3.15},
      {0,0,2.85},{0,0,2.85},{0,0,2.85},{0,0,2.85},
      {-0.2,0,2.7},{-0.2,0.112,2.7},{-0.112,0.2,2.7},{0,0.2,2.7},
    },
    { // Patch 23
      {0,0,3.15},{0,0,3.15},{0,0,3.15},{0,0,3.15},
      {0,0.8,3.15},{0.45,0.8,3.15},{0.8,0.45,3.15},{0.8,0,3.15},
      {0,0,2.85},{0,0,2.85},{0,0,2.85},{0,0,2.85},
      {0,0.2,2.7},{0.112,0.2,2.7},{0.2,0.112,2.7},{0.2,0,2.7},
    },
    { // Patch 24
      {0.2,0,2.7},{0.2,-0.112,2.7},{0.112,-0.2,2.7},{0,-0.2,2.7},
      {0.4,0,2.55},{0.4,-0.224,2.55},{0.224,-0.4,2.55},{0,-0.4,2.55},
      {1.3,0,2.55},{1.3,-0.728,2.55},{0.728,-1.3,2.55},{0,-1.3,2.55},
      {1.3,0,2.4},{1.3,-0.728,2.4},{0.728,-1.3,2.4},{0,-1.3,2.4},
    },
    { // Patch 25
      {0,-0.2,2.7},{-0.112,-0.2,2.7},{-0.2,-0.112,2.7},{-0.2,0,2.7},
      {0,-0.4,2.55},{-0.224,-0.4,2.55},{-0.4,-0.224,2.55},{-0.4,0,2.55},
      {0,-1.3,2.55},{-0.728,-1.3,2.55},{-1.3,-0.728,2.55},{-1.3,0,2.55},
      {0,-1.3,2.4},{-0.728,-1.3,2.4},{-1.3,-0.728,2.4},{-1.3,0,2.4},
    },
    { // Patch 26
      {-0.2,0,2.7},{-0.2,0.112,2.7},{-0.112,0.2,2.7},{0,0.2,2.7},
      {-0.4,0,2.55},{-0.4,0.224,2.55},{-0.224,0.4,2.55},{0,0.4,2.55},
      {-1.3,0,2.55},{-1.3,0.728,2.55},{-0.728,1.3,2.55},{0,1.3,2.55},
      {-1.3,0,2.4},{-1.3,0.728,2.4},{-0.728,1.3,2.4},{0,1.3,2.4},
    },
    { // Patch 27
      {0,0.2,2.7},{0.112,0.2,2.7},{0.2,0.112,2.7},{0.2,0,2.7},
      {0,0.4,2.55},{0.224,0.4,2.55},{0.4,0.224,2.55},{0.4,0,2.55},
      {0,1.3,2.55},{0.728,1.3,2.55},{1.3,0.728,2.55},{1.3,0,2.55},
      {0,1.3,2.4},{0.728,1.3,2.4},{1.3,0.728,2.4},{1.3,0,2.4},
    },
    { // Patch 28
      {0,0,0},{0,0,0},{0,0,0},{0,0,0},
      {1.425,0,0},{1.425,0.798,0},{0.798,1.425,0},{0,1.425,0},
      {1.5,0,0.075},{1.5,0.84,0.075},{0.84,1.5,0.075},{0,1.5,0.075},
      {1.5,0,0.15},{1.5,0.84,0.15},{0.84,1.5,0.15},{0,1.5,0.15},
    },
    { // Patch 29
      {0,0,0},{0,0,0},{0,0,0},{0,0,0},
      {0,1.425,0},{-0.798,1.425,0},{-1.425,0.798,0},{-1.425,0,0},
      {0,1.5,0.075},{-0.84,1.5,0.075},{-1.5,0.84,0.075},{-1.5,0,0.075},
      {0,1.5,0.15},{-0.84,1.5,0.15},{-1.5,0.84,0.15},{-1.5,0,0.15},
    },
    { // Patch 30
      {0,0,0},{0,0,0},{0,0,0},{0,0,0},
      {-1.425,0,0},{-1.425,-0.798,0},{-0.798,-1.425,0},{0,-1.425,0},
      {-1.5,0,0.075},{-1.5,-0.84,0.075},{-0.84,-1.5,0.075},{0,-1.5,0.075},
      {-1.5,0,0.15},{-1.5,-0.84,0.15},{-0.84,-1.5,0.15},{0,-1.5,0.15},
    },
    { // Patch 31
      {0,0,0},{0,0,0},{0,0,0},{0,0,0},
      {0,-1.425,0},{0.798,-1.425,0},{1.425,-0.798,0},{1.425,0,0},
      {0,-1.5,0.075},{0.84,-1.5,0.075},{1.5,-0.84,0.075},{1.5,0,0.075},
      {0,-1.5,0.15},{0.84,-1.5,0.15},{1.5,-0.84,0.15},{1.5,0,0.15},
    },
  };
    file.begingroup("Teapot");
    for(size_t i = 0; i < NUMBER_OF_PATCHES; ++i)
    {
       if(1) file.addPatch(controlPoints[i],materials[i%N_COLOURS]);
       if(0) file.addSurface(3,3,4,4,controlPoints[i],knotsU,knotsV,materials[i%N_COLOURS],NULL); // use (too) general API for the same result as above
    }
    file.endgroup();

    double t[4][4];
    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=6;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=0;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;

    PRCoptions teapotopt;
    teapotopt.no_break = true;
    teapotopt.do_break = false;
    teapotopt.compression = 0.0001;
    // force joining together of patches, damaging transparency
    file.begingroup("Teapot rendered in the way of opaque surfaces and transferred",&teapotopt,t);
    for(size_t i = 0; i < NUMBER_OF_PATCHES; ++i)
    {
       file.addPatch(controlPoints[i],materials[i%N_COLOURS]);
    }
    file.endgroup();
  }
  
  const size_t NUMBER_OF_POINTS = 31;
  double points[NUMBER_OF_POINTS][3];
  for(size_t i = 0; i < NUMBER_OF_POINTS; ++i)
  {
    points[i][0] = 3.5*cos(3.0*i/NUMBER_OF_POINTS*2.0*pi);
    points[i][1] = 3.5*sin(3.0*i/NUMBER_OF_POINTS*2.0*pi);
    points[i][2] = 5.0*i/NUMBER_OF_POINTS-1.0;
  }
  const size_t NUMBER_OF_WIRES = 2;
  double shifted_points[NUMBER_OF_WIRES][NUMBER_OF_POINTS][3];
  for(size_t wire = 0; wire < NUMBER_OF_WIRES; ++wire)
    for(size_t point = 0; point < NUMBER_OF_POINTS; ++point)
      {
        shifted_points[wire][point][0] = points[point][0];
        shifted_points[wire][point][1] = points[point][1];
        shifted_points[wire][point][2] = points[point][2]+0.1*wire+0.1;
      }
  double knots[3+NUMBER_OF_POINTS+1];
  knots[0] = 1;
  for(size_t i = 1; i < 3+NUMBER_OF_POINTS; ++i)
  {
    knots[i] = (i+2)/3; // integer division is intentional
  }
  knots[3+NUMBER_OF_POINTS] = (3+NUMBER_OF_POINTS+1)/3;

  PRCoptions grpopt;
  grpopt.no_break = true;
  grpopt.do_break = false;
  grpopt.tess = true;
  if(1){
    double point1[3] = {11,0,0};
    double point2[3] = {12,0,0};
    double points[2][3] = {{9,0,0},{10,0,0}};
    file.begingroup("points",&grpopt);
    file.addPoint(point1, RGBAColour(1.0,0.0,0.0));
    file.addPoint(point2, RGBAColour(1.0,0.0,0.0));
    file.addPoints(2, points, RGBAColour(1.0,0.0,0.0,0.5),10);
    file.endgroup();
  }

  if(1){
    PRCoptions grpopt;
    grpopt.no_break = true;
    grpopt.do_break = false;
    grpopt.tess = true;
    grpopt.closed = true;

    double t[4][4];

    const size_t nP = 5;
    double P[nP][3] = {{0,0,0},{1,0,0},{1,1,0},{0,1,0},{0,2,0}};
    const size_t nI = 3;
    uint32_t PI[nI][3] = {{0,1,3},{1,2,3},{3,2,4}};
    const size_t nM = 2;
    PRCmaterial M[nM];
    M[0] = PRCmaterial(
      RGBAColour(0.0,0.0,0.18),
      RGBAColour(0.0,0.0,0.878431),
      RGBAColour(0.0,0.0,0.32),
      RGBAColour(0.0,0.0,0.072),
      1.0,0.1);
    M[1] = PRCmaterial(
      RGBAColour(0.18,0.0,0.0),
      RGBAColour(0.878431,0.0,0.0),
      RGBAColour(0.32,0.0,0.0),
      RGBAColour(0.072,0.0,0.0),
      0.5,0.1);
    uint32_t MI[nI] = {0,1,0};
    const size_t nN = 2;
    double N[nN][3] = {{0,0,1},{0,0,-1}};
    uint32_t NI[nI][3] = {{0,0,0},{0,0,0},{1,1,1}};

    const uint32_t nC = 3;
    RGBAColour C[nC];
    uint32_t CI[nI][3] = {{0,0,0},{1,1,1},{1,1,2}};

    PRCmaterial materialGreen(
      RGBAColour(0.0,0.18,0.0),
      RGBAColour(0.0,0.878431,0.0),
      RGBAColour(0.0,0.32,0.0),
      RGBAColour(0.0,0.072,0.0),
      1.0,0.1);

    if(1){
    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=-1;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    file.begingroup("triangles_onecolor_with_normals",&grpopt,t);
    file.addTriangles(nP, P, nI, PI, materialGreen, nN, N, NI, 0, NULL, NULL, 0, NULL, NULL, 0, NULL, NULL);
    file.endgroup();
    }

    if(1){
    file.begingroup("triangles_onecolor",&grpopt);
    file.addTriangles(nP, P, nI, PI, materialGreen, 0, NULL, NULL, 0, NULL, NULL, 0, NULL, NULL, 0, NULL, NULL);
    file.endgroup();
    }

    if(1){
    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=1;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    file.begingroup("triangles_manymaterials",&grpopt,t);
    file.addTriangles(nP, P, nI, PI, materialGreen, 0, NULL, NULL, 0, NULL, NULL, 0, NULL, NULL, nM, M, MI);
    file.endgroup();
    }

    if(1){
    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=2;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    PRCmaterial materialBase(
      RGBAColour(0.1,0.1,0.1,1),
      RGBAColour(1,1,1,1),
      RGBAColour(0.1,0.1,0.1,1),
      RGBAColour(0.1,0.1,0.1,1),
      1.0,0.1);
    C[0] = RGBAColour(1,0,0,0.1);
    C[1] = RGBAColour(0,1,0,0.5);
    C[2] = RGBAColour(0,0,1,0.9);
    file.begingroup("triangles_rgba_vertexcolors_on_opaque_a_component_of_vertexcolor_ignored",&grpopt,t);
    file.addTriangles(nP, P, nI, PI, materialBase, 0, NULL, NULL, 0, NULL, NULL, nC, C, CI, 0, NULL, NULL);
    file.endgroup();

    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=3;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    PRCmaterial materialTransparent(
      RGBAColour(0.1,0.1,0.1,0.3),
      RGBAColour(1,1,1,0.3),
      RGBAColour(0.1,0.1,0.1,0.3),
      RGBAColour(0.1,0.1,0.1,0.3),
      0.3,0.1);
    C[0] = RGBAColour(1,0,0,0.1);
    C[1] = RGBAColour(0,1,0,0.5);
    C[2] = RGBAColour(0,0,1,0.9);
    file.begingroup("triangles_rgba_vertexcolors_on_transparent_a_component_of_vertexcolor_ignored",&grpopt,t);
    file.addTriangles(nP, P, nI, PI, materialTransparent, 0, NULL, NULL, 0, NULL, NULL, nC, C, CI, 0, NULL, NULL);
    file.endgroup();

    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=4;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    C[0] = RGBAColour(1,0,0,0.5);
    C[1] = RGBAColour(0,1,0,0.5);
    C[2] = RGBAColour(0,0,1,0.5);
    file.begingroup("triangles_rgb_vertexcolors_on_transparent_may_not_work_in_OpenGL",&grpopt,t);
    file.addTriangles(nP, P, nI, PI, materialTransparent, 0, NULL, NULL, 0, NULL, NULL, nC, C, CI, 0, NULL, NULL);
    file.endgroup();
    }

    if(1){
    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=5;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    const uint32_t picture_width=2;
    const uint32_t picture_height=2;
    const uint8_t picRGB[picture_width*picture_height*3] =
      {255,0,0, 0,255,0, 0,0,255, 0,0,0 };
//      {255,255,0, 255,0,0, 255,0,0, 255,0,0 };
//     { 255,0,0, 255,255,0, 255,255,0, 255,255,255 };
//      {255,0,0, 0,255,0, 0,0,255, 0,0,0 };

    const uint8_t picRGBA[picture_width*picture_height*4] =
      {255,0,0,255, 0,255,0,150, 0,0,255,150, 0,0,0,100 };
// (1,0) 2 3 (1,1)
// (0,0) 0 1 (1,0)
    uint8_t *pictureRGB = new uint8_t[picture_width*picture_height*3];
    for(size_t i=0; i<picture_width*picture_height*3; i++) pictureRGB[i]=picRGB[i];
    picture1 = pictureRGB;
    uint8_t *pictureRGBA = new uint8_t[picture_width*picture_height*4];
    for(size_t i=0; i<picture_width*picture_height*4; i++) pictureRGBA[i]=picRGBA[i];
    picture2 = pictureRGBA;
    const uint32_t nT = 4;
    const double T[nT][2] = { {0.1,0.1}, {0.9,0.1}, {0.9,0.9}, {0.9,0.1} };
    uint32_t TI[nI][3] = {{0,1,3},{1,2,3},{3,2,3}};

    PRCmaterial materialBase(
      RGBAColour(0.1,0.1,0.1,1),
      RGBAColour(1,1,1,1),
      RGBAColour(0.1,0.1,0.1,1),
      RGBAColour(0.1,0.1,0.1,1),
      1.0,0.1,
      pictureRGB, picture_width, picture_height, false, false);
    file.begingroup("triangles_rgb_texture",&grpopt,t);
    file.addTriangles(nP, P, nI, PI, materialBase, 0, NULL, NULL, nT, T, TI, 0, NULL, NULL, 0, NULL, NULL);
    file.endgroup();

    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=6;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    PRCmaterial materialTransparent(
      RGBAColour(0.1,0.1,0.1,0.3),
      RGBAColour(1,1,1,0.3),
      RGBAColour(0.1,0.1,0.1,0.3),
      RGBAColour(0.1,0.1,0.1,0.3),
      0.3,0.1,
      pictureRGBA, picture_width, picture_height, true, true);
    file.begingroup("triangles_rgba_texture_replacing",&grpopt,t);
    file.addTriangles(nP, P, nI, PI, materialTransparent, 0, NULL, NULL, nT, T, TI, 0, NULL, NULL, 0, NULL, NULL);
    file.endgroup();
    }

  }

  if(1){
    PRCoptions grpopt;
    grpopt.no_break = true;
    grpopt.do_break = false;
    grpopt.tess = true;
    grpopt.closed = true;

    double t[4][4];

    const size_t nP = 6;
    double P[nP][3] = {{3+0,0,0},{3+1,0,0},{3+0,1,0},{3+1,1,0},{3+0,2,0},{3+1,2,0}};
    const size_t nI = 2;
    uint32_t PI[nI][4] = {{0,1,3,2},{2,3,5,4}};
    const size_t nM = 2;
    PRCmaterial M[nM];
    M[0] = PRCmaterial(
      RGBAColour(0.0,0.0,0.18),
      RGBAColour(0.0,0.0,0.878431),
      RGBAColour(0.0,0.0,0.32),
      RGBAColour(0.0,0.0,0.072),
      1.0,0.1);
    M[1] = PRCmaterial(
      RGBAColour(0.18,0.0,0.0),
      RGBAColour(0.878431,0.0,0.0),
      RGBAColour(0.32,0.0,0.0),
      RGBAColour(0.072,0.0,0.0),
      0.5,0.1);
    uint32_t MI[nI] = {0,1};
    const size_t nN = 2;
    double N[nN][3] = {{0,0,1},{0,0,-1}};
    uint32_t NI[nI][4] = {{0,0,0,0},{1,1,1,1}};

    const uint32_t nC = 3;
    RGBAColour C[nC];
    uint32_t CI[nI][4] = {{0,0,1,1},{1,1,2,2}};
    C[0] = RGBAColour(1,0,0,0.5);
    C[1] = RGBAColour(0,1,0,0.5);
    C[2] = RGBAColour(0,0,1,0.5);

    PRCmaterial materialGreen(
      RGBAColour(0.0,0.18,0.0),
      RGBAColour(0.0,0.878431,0.0),
      RGBAColour(0.0,0.32,0.0),
      RGBAColour(0.0,0.072,0.0),
      1.0,0.1);

    PRCmaterial materialBase(
      RGBAColour(0.1,0.1,0.1,1),
      RGBAColour(1,1,1,1),
      RGBAColour(0.1,0.1,0.1,1),
      RGBAColour(0.1,0.1,0.1,1),
      1.0,0.1);

    PRCmaterial materialTransparent(
      RGBAColour(0.1,0.1,0.1,0.3),
      RGBAColour(1,1,1,0.3),
      RGBAColour(0.1,0.1,0.1,0.3),
      RGBAColour(0.1,0.1,0.1,0.3),
      0.3,0.1);

    if(1){
    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=-1;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    file.begingroup("quads_onecolor_with_normals",&grpopt,t);
    file.addQuads(nP, P, nI, PI, materialGreen, nN, N, NI, 0, NULL, NULL, 0, NULL, NULL, 0, NULL, NULL);
    file.endgroup();
    }

    if(1){
    file.begingroup("quads_onecolor",&grpopt);
    file.addQuads(nP, P, nI, PI, materialGreen, 0, NULL, NULL, 0, NULL, NULL, 0, NULL, NULL, 0, NULL, NULL);
    file.endgroup();
    }

    if(1){
    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=1;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    file.begingroup("quads_manymaterials",&grpopt,t);
    file.addQuads(nP, P, nI, PI, materialGreen, 0, NULL, NULL, 0, NULL, NULL, 0, NULL, NULL, nM, M, MI);
    file.endgroup();
    }

    if(1){
    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=2;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    file.begingroup("quads_rgb_vertexcolors_on_transparent_may_not_work_in_OpenGL",&grpopt,t);
    file.addQuads(nP, P, nI, PI, materialTransparent, 0, NULL, NULL, 0, NULL, NULL, nC, C, CI, 0, NULL, NULL);
    file.endgroup();
    }

    if(1){
    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=5;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    const uint32_t picture_width=2;
    const uint32_t picture_height=2;
    const uint8_t picRGB[picture_width*picture_height*3] =
      {255,0,0, 0,255,0, 0,0,255, 0,0,0 };
//      {255,255,0, 255,0,0, 255,0,0, 255,0,0 };
//     { 255,0,0, 255,255,0, 255,255,0, 255,255,255 };
//      {255,0,0, 0,255,0, 0,0,255, 0,0,0 };

    const uint8_t picRGBA[picture_width*picture_height*4] =
      {255,0,0,255, 0,255,0,150, 0,0,255,150, 0,0,0,100 };
// (1,0) 2 3 (1,1)
// (0,0) 0 1 (1,0)
    uint8_t *pictureRGB = new uint8_t[picture_width*picture_height*3];
    for(size_t i=0; i<picture_width*picture_height*3; i++) pictureRGB[i]=picRGB[i];
    picture3 = pictureRGB;
    uint8_t *pictureRGBA = new uint8_t[picture_width*picture_height*4];
    for(size_t i=0; i<picture_width*picture_height*4; i++) pictureRGBA[i]=picRGBA[i];
    picture4 = pictureRGBA;
    const uint32_t nT = 3;
    const double T[nT][2] = { {0.1,0.1}, {0.5,0.5}, {0.9,0.9} };
    uint32_t TI[nI][4] = {{0,0,1,1},{1,1,2,2}};

    PRCmaterial materialBase(
      RGBAColour(0.1,0.1,0.1,1),
      RGBAColour(1,1,1,1),
      RGBAColour(0.1,0.1,0.1,1),
      RGBAColour(0.1,0.1,0.1,1),
      1.0,0.1,
      pictureRGB, picture_width, picture_height, false, false);
    file.begingroup("quads_rgb_texture",&grpopt,t);
    file.addQuads(nP, P, nI, PI, materialBase, 0, NULL, NULL, nT, T, TI, 0, NULL, NULL, 0, NULL, NULL);
    file.endgroup();

    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=6;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    PRCmaterial materialTransparent(
      RGBAColour(0.1,0.1,0.1,0.3),
      RGBAColour(1,1,1,0.3),
      RGBAColour(0.1,0.1,0.1,0.3),
      RGBAColour(0.1,0.1,0.1,0.3),
      0.6,0.1,
      pictureRGBA, picture_width, picture_height, true, true);
    file.begingroup("quads_rgba_texture_replacing",&grpopt,t);
    file.addQuads(nP, P, nI, PI, materialTransparent, 0, NULL, NULL, nT, T, TI, 0, NULL, NULL, 0, NULL, NULL);
    file.endgroup();
    }

  }

  if(1){
    file.begingroup("polyline",&grpopt);
    file.addLine(NUMBER_OF_POINTS, points, RGBAColour(1.0,0.0,0.0));
    file.endgroup();

    file.begingroup("polylines",&grpopt);
    file.addLine(NUMBER_OF_POINTS, shifted_points[0], RGBAColour(0.0,1.0,0.0),2.0);
    file.addLine(NUMBER_OF_POINTS, shifted_points[1], RGBAColour(1.0,1.0,0.0),2.0);
    file.endgroup();
  }

  if(1){
    PRCoptions grpopt;
    grpopt.no_break = true;
    grpopt.do_break = false;
    grpopt.tess = true;
    grpopt.closed = true;

    double t[4][4];

    const size_t nP = 5;
    double P[nP][3] = {{0,0,0},{1,0,0},{1,1,0},{0,1,0},{0,2,0}};
    const size_t nI = 1+5 + 1+3;
    uint32_t PI[nI] = {5,0,1,2,3,0, 3,2,4,3 }; // square of 5 points, wire of 2 points
    const uint32_t nC = 5;  // number of colors
    RGBAColour C[nC]; // color for segments
    C[0] = RGBAColour(1,0,0);
    C[1] = RGBAColour(0,1,0);
    C[2] = RGBAColour(0,0,1);
    C[3] = RGBAColour(0,0,0);
    C[4] = RGBAColour(1,1,0);
    const uint32_t nCI = 5-1 + 3-1;  // number of segments
    uint32_t CI[nCI] = { 0,1,2,3, 4,4 };

    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=5;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    file.begingroup("polyline_indexed_one_color",&grpopt,t);
    file.addLines(nP, P, nI, PI, RGBAColour(1,0,0), 1, false, 0, NULL, 0, NULL);
    file.endgroup();

    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=6;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    file.begingroup("polyline_indexed_color_per_segment",&grpopt,t);
    file.addLines(nP, P, nI, PI, RGBAColour(1,0,0), 2, true, nC, C, nCI, CI);
    file.endgroup();

    t[0][0]=1; t[0][1]=0; t[0][2]=0; t[0][3]=0;
    t[1][0]=0; t[1][1]=1; t[1][2]=0; t[1][3]=0;
    t[2][0]=0; t[2][1]=0; t[2][2]=1; t[2][3]=7;
    t[3][0]=0; t[3][1]=0; t[3][2]=0; t[3][3]=1;
    const uint32_t nCI2 = 5 + 3;  // number of vertices
    uint32_t CI2[nCI2] = { 0,1,2,3,0, 4,4,4 };
    file.begingroup("polyline_indexed_color_per_vertex_is_broken",&grpopt,t);
    file.addLines(nP, P, nI, PI, RGBAColour(1,0,0), 2, false, nC, C, nCI2, CI2);
    file.endgroup();
  }

  if(1)
  {
    file.begingroup("bezier_wire");
    file.addBezierCurve(NUMBER_OF_POINTS,points,RGBAColour(1.0,1.0,1.0));
    file.endgroup();
  }
  if(1)
  {
    file.begingroup("NURBS_wire");
    file.addCurve(3, NUMBER_OF_POINTS, points, knots, RGBAColour(1.0,1.0,1.0), NULL); // genarl API for the above
    file.endgroup();
  }

// following box examples show a) different ways to represent a surface consisting of flat rectangles
// b) that the only way to have almost working transparency is a set of NURBS bodies.
// (Or may be other topology types like plane also work
// demonstration how non-transparent materials work the same for all kinds of objects  

  if(1) { // demonstration how non-transparent materials work the same for all kinds of objects  
    const size_t NUMBER_OF_PATCHES = 6;
    double vertices[NUMBER_OF_PATCHES][4][3] = 
    {
      { // Patch 0
       {-1,-1,-1},
       { 1,-1,-1},
       {-1, 1,-1},
       { 1, 1,-1}
      },
      { // Patch 1
       {-1,-1, 1},
       { 1,-1, 1},
       {-1, 1, 1},
       { 1, 1, 1}
      },
      { // Patch 2
       {-1,-1,-1},
       { 1,-1,-1},
       {-1,-1, 1},
       { 1,-1, 1}
      },
      { // Patch 3
       {-1, 1,-1},
       { 1, 1,-1},
       {-1, 1, 1},
       { 1, 1, 1}
      },
      { // Patch 4
       {-1,-1,-1},
       {-1, 1,-1},
       {-1,-1, 1},
       {-1, 1, 1}
      },
      { // Patch 5
       { 1,-1,-1},
       { 1, 1,-1},
       { 1,-1, 1},
       { 1, 1, 1}
      }
    };
    const size_t NUMBER_OF_BOXES = 6;
    double shifted_vertices[NUMBER_OF_BOXES][NUMBER_OF_PATCHES][4][3];
    for(size_t box = 0; box < NUMBER_OF_BOXES; ++box)
      for(size_t patch = 0; patch < NUMBER_OF_PATCHES; ++patch)
        for(size_t i = 0; i < 4; ++i)
        {
          shifted_vertices[box][patch][i][0] = vertices[patch][i][0]+3*box;
          shifted_vertices[box][patch][i][1] = vertices[patch][i][1];
          shifted_vertices[box][patch][i][2] = vertices[patch][i][2]-2;
        }
    PRCmaterial materialGreen(
      RGBAColour(0.0,0.18,0.0),
      RGBAColour(0.0,0.878431,0.0),
      RGBAColour(0.0,0.32,0.0),
      RGBAColour(0.0,0.072,0.0),
      1.0,0.1);

    PRCoptions grpopt;
    file.begingroup("TransparentBox");
    grpopt.no_break = false;
    grpopt.do_break = true;
    grpopt.tess = false;
    file.begingroup("TransparentBoxSetOfNURBSBodies",&grpopt);
    for(size_t patch = 0; patch < NUMBER_OF_PATCHES; ++patch)
    {
      file.addRectangle(shifted_vertices[0][patch], materials[(patch*5)%N_COLOURS]);
    }
    file.endgroup();
    grpopt.no_break = true;
    grpopt.do_break = false;
    grpopt.tess = false;
    file.begingroup("TransparentBoxNURBSFaces",&grpopt);
    for(size_t patch = 0; patch < NUMBER_OF_PATCHES; ++patch)
    {
      file.addRectangle(shifted_vertices[1][patch], materials[(patch*5)%N_COLOURS]);
    }
    file.endgroup();
    grpopt.no_break = true;
    grpopt.do_break = false;
    grpopt.tess = true;
    file.begingroup("TransparentBoxTessellated",&grpopt);
    for(size_t patch = 0; patch < NUMBER_OF_PATCHES; ++patch)
    {
      file.addRectangle(shifted_vertices[2][patch], materials[(patch*5)%N_COLOURS]);
    }
    file.endgroup();
    file.endgroup();

    file.begingroup("Box");
    grpopt.no_break = false;
    grpopt.do_break = false;
    grpopt.tess = true;
    file.begingroup("BoxTessellated",&grpopt);
    for(size_t patch = 0; patch < NUMBER_OF_PATCHES; ++patch)
    {
      file.addRectangle(shifted_vertices[3][patch], materialGreen);
    }
    file.endgroup();
    grpopt.no_break = true;
    grpopt.do_break = false;
    grpopt.tess = false;
    file.begingroup("BoxNURBSFaces",&grpopt);
    for(size_t patch = 0; patch < NUMBER_OF_PATCHES; ++patch)
    {
      file.addRectangle(shifted_vertices[4][patch], materialGreen);
    }
    file.endgroup();
    grpopt.no_break = false;
    grpopt.do_break = true;
    grpopt.tess = false;
    file.begingroup("BoxSetOfNURBSBodies",&grpopt);
    for(size_t patch = 0; patch < NUMBER_OF_PATCHES; ++patch)
    {
      file.addRectangle(shifted_vertices[5][patch], materialGreen);
    }
    file.endgroup();
    file.endgroup();
  }

  if(1) {
    PRCmaterial materialGreen(
      RGBAColour(0.0,0.18,0.0),
      RGBAColour(0.0,0.878431,0.0),
      RGBAColour(0.0,0.32,0.0),
      RGBAColour(0.0,0.072,0.0),
      1.0,0.1);
    PRCoptions grpopt;
    grpopt.closed = true;

     file.begingroup("Complex elements, one-sided");

     const double disk_origin[3] = {11,0,2};
     const double disk_x_axis[3] = {1,0,0};
     const double disk_y_axis[3] = {0,-1,0};
     const double disk_scale = 2;
     file.begingroup("disk",&grpopt);
     file.addDisk(1,materialGreen,disk_origin,disk_x_axis,disk_y_axis,disk_scale);
     file.endgroup();
     const double hs_origin[3] = {11,0,2};
     const double hs_x_axis[3] = {1,0,0};
     const double hs_y_axis[3] = {0,1,0};
     const double hs_scale = 2;
     file.begingroup("hemisphere",&grpopt);
     file.addHemisphere(1,materialGreen,hs_origin,hs_x_axis,hs_y_axis,hs_scale);
     file.endgroup();
     const double cyl_origin[3] = {11,0,1};
     const double cyl_x_axis[3] = {1,0,0};
     const double cyl_y_axis[3] = {0,1,0};
     const double cyl_scale = 1;
     file.begingroup("cylinder",&grpopt);
     file.addCylinder(1,1,materialGreen,cyl_origin,cyl_x_axis,cyl_y_axis,cyl_scale);
     file.endgroup();
     const double cone_origin[3] = {11,0,1};
     const double cone_x_axis[3] = {1,0,0};
     const double cone_y_axis[3] = {0,-1,0};
     const double cone_scale = 1;
     file.begingroup("cone",&grpopt);
     file.addCone(1,1,materialGreen,cone_origin,cone_x_axis,cone_y_axis,cone_scale);
     file.endgroup();
     const double sp_origin[3] = {11,0,-1};
     const double sp_x_axis[3] = {1,0,0};
     const double sp_y_axis[3] = {0,1,0};
     const double sp_scale = 1;
     file.begingroup("sphere",&grpopt);
     file.addSphere(0.5,materialGreen,sp_origin,sp_x_axis,sp_y_axis,sp_scale);
     file.endgroup();
     const double tor_origin[3] = {11,0,0};
     const double tor_x_axis[3] = {1,0,0};
     const double tor_y_axis[3] = {0,1,0};
     const double tor_scale = 1;
     file.begingroup("torus",&grpopt);
     file.addTorus(0.5,0.1,0,360,materialGreen,tor_origin,tor_x_axis,tor_y_axis,tor_scale);
     file.endgroup();

     const size_t NUMBER_OF_POINTS = 4;
     double cpoints[NUMBER_OF_POINTS][3];
     cpoints[0][0] = 1; cpoints[0][1] = 0;                 cpoints[0][2] = 0;
     cpoints[1][0] = 1; cpoints[1][1] = 0.552284749830793; cpoints[1][2] = 0;
     cpoints[2][0] = 0.552284749830793; cpoints[2][1] = 1; cpoints[2][2] = 0;
     cpoints[3][0] = 0;                 cpoints[3][1] = 1; cpoints[3][2] = 0;
     double opoints[NUMBER_OF_POINTS][3];
     opoints[0][0] = 1.1*1; opoints[0][1] = 1.1*0;                 opoints[0][2] = 0;
     opoints[1][0] = 1.1*1; opoints[1][1] = 1.1*0.552284749830793; opoints[1][2] = 0;
     opoints[2][0] = 1.1*0.552284749830793; opoints[2][1] = 1.1*1; opoints[2][2] = 0;
     opoints[3][0] = 1.1*0;                 opoints[3][1] = 1.1*1; opoints[3][2] = 0;
     const double tube_origin[3] = {11,0,0};
     const double tube_x_axis[3] = {1,0,0};
     const double tube_y_axis[3] = {0,1,0};
     const double tube_scale = 1;
     file.begingroup("tube",&grpopt);
     file.addTube(NUMBER_OF_POINTS,cpoints,opoints,false,materialGreen,tube_origin,tube_x_axis,tube_y_axis,tube_scale);
     file.endgroup();

     file.endgroup();
  }

  file.finish();

  delete[] picture1;
  delete[] picture2;
  delete[] picture3;
  delete[] picture4;

#endif
  return 0;
}