user_guide:tutorials:latest:optimization

Differences

This shows you the differences between two versions of the page.


user_guide:tutorials:latest:optimization [2023/11/06 10:57] (current) – created - external edit 127.0.0.1
Line 1: Line 1:
 +====== polymake and Optimization ======
 +
 +By Sebastian Pokutta and Marc E. Pfetsch.
 +
 +===== Introduction =====
 +
 +Polymake offers many interesting features that can help (discrete) optimizers to analyze optimization problems. For example
 +
 +  * linear optimization can be performed exactly and visualized in small dimensions
 +  * the convex hull of feasible points of an integer program can be computed and analyzed
 +  * Hilbert bases can be computed
 +
 +There are several other tutorials that cover similar topics:
 +
 +  * [[ilp_and_hilbertbases|ILP and Hilbert Bases]]
 +  * [[branchnbound|Gomory Cuts]]
 +  * [[gomory|Branch and Bound]]
 +  * [[lattice_polytopes_tutorial|Tutorial for Lattice Polytopes]]
 +
 +This tutorial is targeted towards the optimization community, since, surprisingly, polymake does not seem to be well known in this community. In particular, the community still tends to use the quite old program ''%%porta%%'' to compute convex hulls or inequality descriptions. While ''%%porta%%'' still does a decent job here, ''%%polymake%%'' offers a much broader feature set. Polymake supports several convex hull algorithms which might be better suited depending on the data. Moreover it offers many visualization tools that can help to better //understand// a given polytope. We think that polymake has many advantages for discrete optimizers and hope that this tutorial will help to spread the usage of polymake.
 +
 +You can find files of the example LPs in the folder demo/files/optimization in your polymake directory.
 +
 +===== Input: lp2poly =====
 +
 +The first important step is to get the desired input into ''%%polymake%%''. While it is possible to define polytopes in ''%%polymake%%'' right away, a more convenient methods, especially for people working in discrete optimization, is the possibility to read an optimization problem from a file. This has the advantage that input that is usually generated outside ''%%polymake%%'' can simply be loaded. A standard file format for linear (integer) programs is the so-called CPLEX LP-format, which is well known in the optimization community and can be read into ''%%polymake%%'' in a very straightforward way as shown in the following example:
 +
 +<code>
 +Minimize
 + obj:  x1 + x2 + x3
 +Subject to
 + C1: x1 + x2 + x3 <= 2
 +Bounds
 + 0 <= x1 <= 1
 + 0 <= x2 <= 1
 + 0 <= x3 <= 1
 +End
 +</code>
 +Thus, the file describes a 0/1-cube in three dimensions. It should be easy to adapt this format to other cases (If for example ''%%x1%%'' does not have any bounds you can write ''%%x1 free%%'' instead).
 +
 +Now assume that this example is contained in file ''%%c3t.lp%%''. We create a polytope from the file via:
 +
 +<code perl>
 +> $f=lp2poly('files/optimization/c3t.lp');
 +</code>
 +The polytope ''%%$f%%'' is coded via floating point numbers:
 +
 +<code perl>
 +> print $f->type->full_name;
 +Polytope<Rational>
 +</code>
 +We convert it to a rational polytope via:
 +
 +<code perl>
 +> $p = new Polytope<Rational>($f);
 +</code>
 +Now, ''%%$p%%'' contains data for the linear optimization problem in the file ''%%ct3.lp%%''. The corresponding polyhedron is stored directly in ''%%$p%%''. The information about the linear objective is contained in ''%%$p->LP%%''.
 +
 +===== Linear Optimization =====
 +
 +Polymake can be used to perform several actions related to linear optimization (linear programming - LP). For instance, one can exactly solve a linear program (via lrs or cdd). Before we explain the corresponding usage, we first need to have a linear optimization problem at hand.
 +
 +Assuming that we are given the above example in variable ''%%$p%%'', the linear objective can be accessed via
 +
 +<code perl>
 +> print $p->LP->LINEAR_OBJECTIVE;
 +0 1 1 1
 +</code>
 +Thus - as described in the file - the objective function coefficients are 1 for all three variables (and there is an offset of 0).
 +
 +Now, we can solve the corresponding linear program via
 +
 +<code perl>
 +> print $p->LP->MAXIMAL_VALUE;
 +2
 +</code>
 +Thus, the maximal value that we can obtain via the above linear objective function is 2. We can also get an optimal vertex via
 +
 +<code perl>
 +> print $p->LP->MAXIMAL_VERTEX;
 +1 0 1 1
 +</code>
 +This vertex corresponds to setting ''%%x1=0, x2=1, x3=1%%''. The optimal face can also be computed:
 +
 +<code perl>
 +> print $p->LP->MAXIMAL_FACE;
 +{4 5 6}
 +</code>
 +This means that the optimal face is the convex hull of three vertices (with indices 4, 5, 6).
 +
 +Of course, by replacing ''%%MAXIMAL%%'' by ''%%MINIMAL%%'', we obtain the corresponding values for minimization.
 +
 +The directed graph obtained by directing the graph of the polytope in the direction of increasing objective function can be obtained via
 +
 +<code perl>
 +> $p->VISUAL->DIRECTED_GRAPH;
 +</code>
 +<HTML>
 +<!--
 +polymake for knusper
 +Thu Mar  3 00:38:23 2022
 +c3t
 +-->
 +
 +
 +<html>
 +   <head>
 +      <meta charset=utf-8>
 +      <title>c3t</title>
 +      <style>
 +/*
 +// COMMON_CODE_BLOCK_BEGIN
 +*/
 +         html {overflow: scroll;}
 +         strong{font-size: 18px;}
 +         canvas { z-index: 8; }
 +         input[type='radio'] {margin-left:0;}
 +         input[type='checkbox'] {margin-right:7px; margin-left: 0px; padding-left:0px;}
 +         .group{padding-bottom: 15px;}
 +         .settings * {z-index: 11; }
 +         .settings{z-index: 10; font-family: Arial, Helvetica, sans-serif; margin-left: 30px; visibility: hidden; width: 14em; height: 96%; border: solid 1px silver; padding: 2px; overflow-y: scroll; box-sizing: border-box; background-color: white; position: absolute;}
 +         .indented{margin-left: 20px; margin-top: 10px; padding-bottom: 0px;} 
 +         .shownObjectsList{overflow: auto; max-width: 150px; max-height: 150px;}
 +         .showSettingsButton{visibility: visible; z-index: 12; position: absolute }
 +         .hideSettingsButton{visibility: hidden; z-index: 12; position: absolute; opacity: 0.5}
 +         button{margin-left: 0; margin-top: 10px}
 +         img{cursor: pointer;}
 +         .suboption{padding-top: 15px;}
 +         #model84526085210 { width: 100%; height: 100%; }
 +         .threejs_container { width: 100%; height: 75vh;}
 +         .settings{max-height: 74vh} 
 +         input[type=range] {
 +           -webkit-appearance: none;
 +           padding:0; 
 +           width:90%; 
 +           margin-left: auto;
 +           margin-right: auto;
 +           margin-top: 15px;
 +           margin-bottom: 15px;
 +           display: block;
 +         }
 +         input[type=range]:focus {
 +           outline: none;
 +         }
 +         input[type=range]::-webkit-slider-runnable-track {
 +           height: 4px;
 +           cursor: pointer;
 +           animate: 0.2s;
 +           box-shadow: 0px 0px 0px #000000;
 +           background: #E3E3E3;
 +           border-radius: 0px;
 +           border: 0px solid #000000;
 +         }
 +         input[type=range]::-webkit-slider-thumb {
 +           box-shadow: 1px 1px 2px #B8B8B8;
 +           border: 1px solid #ABABAB;
 +           height: 13px;
 +           width: 25px;
 +           border-radius: 20px;
 +           background: #E0E0E0;
 +           cursor: pointer;
 +           -webkit-appearance: none;
 +           margin-top: -5px;
 +         }
 +         input[type=range]:focus::-webkit-slider-runnable-track {
 +           background: #E3E3E3;
 +         }
 +         input[type=range]::-moz-range-track {
 +           height: 4px;
 +           cursor: pointer;
 +           animate: 0.2s;
 +           box-shadow: 0px 0px 0px #000000;
 +           background: #E3E3E3;
 +           border-radius: 0px;
 +           border: 0px solid #000000;
 +         }
 +         input[type=range]::-moz-range-thumb {
 +           box-shadow: 1px 1px 2px #B8B8B8;
 +           border: 1px solid #ABABAB;
 +           height: 13px;
 +           width: 25px;
 +           border-radius: 20px;
 +           background: #E0E0E0;
 +           cursor: pointer;
 +         }
 +         input[type=range]::-ms-track {
 +           height: 4px;
 +           cursor: pointer;
 +           animate: 0.2s;
 +           background: transparent;
 +           border-color: transparent;
 +           color: transparent;
 +         }
 +         input[type=range]::-ms-fill-lower {
 +           background: #E3E3E3;
 +           border: 0px solid #000000;
 +           border-radius: 0px;
 +           box-shadow: 0px 0px 0px #000000;
 +         }
 +         input[type=range]::-ms-fill-upper {
 +           background: #E3E3E3;
 +           border: 0px solid #000000;
 +           border-radius: 0px;
 +           box-shadow: 0px 0px 0px #000000;
 +         }
 +         input[type=range]::-ms-thumb {
 +           box-shadow: 1px 1px 2px #B8B8B8;
 +           border: 1px solid #ABABAB;
 +           height: 13px;
 +           width: 25px;
 +           border-radius: 20px;
 +           background: #E0E0E0;
 +           cursor: pointer;
 +         }
 +         input[type=range]:focus::-ms-fill-lower {
 +           background: #E3E3E3;
 +         }
 +         input[type=range]:focus::-ms-fill-upper {
 +           background: #E3E3E3;
 +         }
 +/*
 +// COMMON_CODE_BLOCK_END
 +*/
 + </style>
 +   </head>
 +<body>
 +   <div class='threejs_container'>
 + <div id='settings_0' class='settings'>
 + <div class=group id='explode_0'>
 + <strong>Explode</strong>
 + <input id='explodeRange_0' type='range' min='0.00001' max=6 step=0.01 value=0.00001>
 + <div class=indented><input id='explodeCheckbox_0' type='checkbox'>Automatic explosion</div>
 + <div class=suboption>Exploding speed</div>
 + <input id='explodingSpeedRange_0' type='range' min=0 max=0.5 step=0.001 value=0.05>
 + </div>
 +
 + <div class=group id='transparency_0' class='transparency'>
 + <strong>Transparency</strong>
 + <input id='transparencyRange_0' type='range' min=0 max=1 step=0.01 value=0>
 +            <div class=indented><input id='depthWriteCheckbox_0' type='checkbox'>depthWrite</div>
 + </div>
 +
 + <div class=group id='rotation_0'>
 + <strong>Rotation</strong>
 + <div class=indented>
 + <div><input type='checkbox' id='changeRotationX_0'> x-axis</div>
 + <div><input type='checkbox' id='changeRotationY_0'> y-axis</div>
 + <div><input type='checkbox' id='changeRotationZ_0'> z-axis</div>
 + <button id='resetButton_0'>Reset</button>
 + </div>
 +
 + <div class=suboption>Rotation speed</div>
 + <input id='rotationSpeedRange_0' type='range' min=0 max=5 step=0.01 value=2>
 + </div>
 +
 +
 + <div class=group id='display_0'>
 + <strong>Display</strong>
 + <div class=indented>
 + <div id='shownObjectTypesList_0' class='shownObjectsList'></div>
 + </div>
 + <div class=suboption>Objects</div>
 + <div class=indented>
 +    <div id='shownObjectsList_0' class='shownObjectsList'></div>
 + </div>
 + </div>
 +         
 +         <div class=group id='camera_0'>
 +            <strong>Camera</strong>
 +            <div class=indented>
 +               <form>
 +                  <select id="cameraType_0">
 +                     <option value='perspective' selected> Perspective<br></option>
 +                     <option value='orthographic' > Orthographic<br></option>
 +                  </select>
 +               </form>
 +            </div>
 +         </div>
 +
 + <div class=group id='svg_0'>
 + <strong>SVG</strong>
 + <div class=indented>
 + <form>
 + <input type="radio" name='screenshotMode' value='download' id='download_0' checked> Download<br>
 + <input type="radio" name='screenshotMode' value='tab' id='tab_0' > New tab<br>
 + </form>
 + <button id='takeScreenshot_0'>Screenshot</button>
 + </div>
 + </div>
 +
 + </div> <!-- end of settings -->
 + <img id='hideSettingsButton_0' class='hideSettingsButton' src='/kernelspecs/r118/polymake/close.svg' width=20px">
 + <img id='showSettingsButton_0' class='showSettingsButton' src='/kernelspecs/r118/polymake/menu.svg' width=20px">
 +<div id="model84526085210"></div>
 +</div>
 +   <script>
 +    requirejs.config({
 +      paths: {
 +        three: '/kernelspecs/r118/polymake/three',
 +        TrackballControls: '/kernelspecs/r118/polymake/TrackballControls',
 +        OrbitControls: '/kernelspecs/r118/polymake/OrbitControls',
 +        Projector: '/kernelspecs/r118/polymake/Projector',
 +        SVGRenderer: '/kernelspecs/r118/polymake/SVGRenderer',
 +        WEBGL: '/kernelspecs/r118/polymake/WebGL',
 +      },
 +      shim: {
 +        'three': { exports: 'THREE'},
 +        'SVGRenderer': { deps: [ 'three' ], exports: 'THREE.SVGRenderer' },
 +        'WEBGL': { deps: [ 'three' ], exports: 'THREE.WEBGL' },
 +        'Projector': { deps: [ 'three' ], exports: 'THREE.Projector' },
 +        'TrackballControls': { deps: [ 'three' ], exports: 'THREE.TrackballControls' },
 +        'OrbitControls': { deps: [ 'three' ], exports: 'THREE.OrbitControls' },
 +      }
 +    });
 +    
 +    require(['three'],function(THREE){
 +        window.THREE = THREE;
 +      require(['TrackballControls', 'OrbitControls', 'Projector', 'SVGRenderer', 'WEBGL'],
 +               function(TrackballControls, OrbitControls, Projector, SVGRenderer, WEBGL) {
 +    THREE.TrackballControls = TrackballControls;
 +    THREE.OrbitControls = OrbitControls;
 +    THREE.Projector = Projector;
 +    THREE.SVGRenderer = SVGRenderer;
 +    THREE.WEBGL = WEBGL;
 +
 +// COMMON_CODE_BLOCK_BEGIN
 +
 +const intervalLength = 25; // for automatic animations
 +const explodableModel = true; 
 +const modelContains = { points: false, pointlabels: false, lines: false, edgelabels: false, faces: false, arrowheads: false };
 +const foldables = [];
 +
 +var three = document.getElementById("model84526085210");
 +var scene = new THREE.Scene();
 +var renderer = new THREE.WebGLRenderer( { antialias: true } );
 +var svgRenderer = new THREE.SVGRenderer( { antialias: true } );
 +renderer.setPixelRatio( window.devicePixelRatio );
 +renderer.setClearColor(0xFFFFFF, 1);
 +svgRenderer.setClearColor(0xFFFFFF, 1);
 +three.appendChild(renderer.domElement);
 +
 +var frustumSize = 4;
 +var cameras = [new THREE.PerspectiveCamera(75, 1, 0.1, 1000), new THREE.OrthographicCamera()];
 +cameras.forEach(function(cam) {
 +    cam.position.set(0, 0, 5);
 +    cam.lookAt(0, 0, 0);  
 +    cam.up.set(0, 1, 0);         
 +});
 +var controls = [new THREE.TrackballControls(cameras[0], three), new THREE.OrbitControls(cameras[1], three)];
 +var camera, control;
 +
 +controls[0].zoomSpeed = 0.2;
 +controls[0].rotateSpeed = 4;
 +
 +
 +// class to allow move points together with labels and spheres
 +var PMPoint = function (x,y,z) {
 +   this.vector = new THREE.Vector3(x,y,z);
 +   this.sprite = null;
 +   this.sphere = null;
 +}
 +PMPoint.prototype.addLabel = function(labelsprite) {
 +   this.sprite = labelsprite;
 +   this.sprite.position.copy(this.vector);
 +}
 +PMPoint.prototype.addSphere = function(spheremesh) {
 +   this.sphere = spheremesh;
 +   this.sphere.position.copy(this.vector);
 +}
 +PMPoint.prototype.set = function(x,y,z) {
 +   this.vector.set(x,y,z);
 +   if (this.sprite) {
 +      this.sprite.position.copy(this.vector);
 +   }
 +   if (this.sphere) {
 +      this.sphere.position.copy(this.vector);
 +   }
 +}
 +PMPoint.prototype.radius = function() {
 +   if (this.sphere) {
 +      return this.sphere.geometry.parameters.radius;
 +   } else {
 +      return 0;
 +   }
 +};
 +// select the target node
 +var target = document.querySelector('#model84526085210');
 +
 +// create an observer instance
 +var observer = new MutationObserver(function(mutations) {
 +   mutations.forEach(function(mutation) {
 +      if (mutation.removedNodes && mutation.removedNodes.length > 0) {
 +         cancelAnimationFrame(renderId);
 +         observer.disconnect();
 +         console.log("cancelled frame "+renderId);
 +      }
 +   });
 +});
 +
 +// configuration of the observer:
 +var config = { childList: true, characterData: true }
 +
 +// pass in the target node, as well as the observer options
 +while (target) {
 +   if (target.className=="output") {
 +      observer.observe(target, config);
 +      break;
 +   }
 +   target = target.parentNode;
 +}
 +
 +// COMMON_CODE_BLOCK_END
 +
 +var obj0 = new THREE.Object3D();
 +obj0.name = "c3t";
 +obj0.userData.explodable = 1;
 +obj0.userData.points = [];
 +obj0.userData.points.push(new PMPoint(0, 0, 1));
 +obj0.userData.points.push(new PMPoint(1, 0, 0));
 +obj0.userData.points.push(new PMPoint(0, 0, 0));
 +obj0.userData.points.push(new PMPoint(0, 1, 0));
 +obj0.userData.points.push(new PMPoint(1, 1, 0));
 +obj0.userData.points.push(new PMPoint(1, 0, 1));
 +obj0.userData.points.push(new PMPoint(0, 1, 1));
 +
 +obj0.userData.pointradii = 0.02;
 +   <!-- Vertex style -->
 +obj0.userData.pointmaterial = new THREE.MeshBasicMaterial( { color: 0xFF0000, side: THREE.DoubleSide, transparent: false } );
 +obj0.userData.pointlabels = ["0", "1", "2", "3", "4", "5", "6"];
 +obj0.userData.edgeindices = [0, 2, 1, 2, 2, 3, 1, 4, 3, 4, 0, 5, 1, 5, 4, 5, 0, 6, 3, 6, 4, 6, 5, 6];
 +obj0.userData.facets = [[5, 4, 6], [6, 3, 2, 0], [1, 4, 5], [5, 0, 2, 1], [6, 4, 3], [3, 4, 1, 2], [5, 6, 0]];
 +   <!-- Facet style -->
 +obj0.userData.facetmaterial = new THREE.MeshBasicMaterial( { color: 0x77EC9E, depthFunc: THREE.LessDepth, depthWrite: false, opacity: 1, polygonOffset: true, polygonOffsetFactor: 1, polygonOffsetUnits: 0.5, side: THREE.DoubleSide, transparent: true } );
 +init_object(obj0);
 +scene.add(obj0);
 +
 +var obj1 = new THREE.Object3D();
 +obj1.name = "GRAPH directed with LP unnamed#0";
 +obj1.userData.explodable = 1;
 +obj1.userData.points = [];
 +obj1.userData.points.push(new PMPoint(0, 0, 1));
 +obj1.userData.points.push(new PMPoint(1, 0, 0));
 +obj1.userData.points.push(new PMPoint(0, 0, 0));
 +obj1.userData.points.push(new PMPoint(0, 1, 0));
 +obj1.userData.points.push(new PMPoint(1, 1, 0));
 +obj1.userData.points.push(new PMPoint(1, 0, 1));
 +obj1.userData.points.push(new PMPoint(0, 1, 1));
 +
 +obj1.userData.edgeindices = [0, 5, 0, 6, 1, 4, 1, 5, 2, 0, 2, 1, 2, 3, 3, 4, 3, 6];
 +   <!-- Edge style -->
 +obj1.userData.edgematerial = new THREE.LineBasicMaterial( { color: 0x000000, linewidth: 1.5, transparent: false } );
 +obj1.userData.arrowstyle = 1;
 +init_object(obj1);
 +scene.add(obj1);
 +
 +// COMMON_CODE_BLOCK_BEGIN
 +function textSpriteMaterial(message, parameters) {
 +    if ( parameters === undefined ) parameters = {};
 +    var fontface = "Helvetica";
 +    var fontsize = parameters.hasOwnProperty("fontsize") ? parameters["fontsize"] : 15;
 +    fontsize = fontsize*10;
 +    var lines = message.split('\\n');
 +    var size = 512;
 +    for(var i = 0; i<lines.length; i++){
 +        var tmp = lines[i].length;
 +        while(tmp*fontsize > size){
 +           fontsize--;
 +        }
 +    }
 +    
 +    var canvas = document.createElement('canvas');
 +    canvas.width = size;
 +    canvas.height = size;
 +    var context = canvas.getContext('2d');
 +    context.fillStyle = "rgba(255, 255, 255, 0)";
 +    context.fill();
 +    context.font = fontsize + "px " + fontface;
 +    
 +    // text color
 +    context.fillStyle = "rgba(0, 0, 0, 1.0)";
 +     for(var i = 0; i<lines.length; i++){
 +        context.fillText(lines[i], size/2, size/2+i*fontsize);
 +     }
 +    
 +    // canvas contents will be used for a texture
 +    var texture = new THREE.Texture(canvas);
 +    texture.needsUpdate = true;
 +    
 +    var spriteMaterial = new THREE.SpriteMaterial({map: texture, depthTest: true, depthWrite: false, polygonOffset: true, polygonOffsetFactor: -1, polygonOffsetUnits: 1 });
 +    return spriteMaterial;
 +}
 +
 +
 +// ---------------------- INITIALIZING OBJECTS--------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +function init_object(obj) {
 +    if (obj.userData.hasOwnProperty("pointmaterial")) {
 +        init_points(obj);
 +        modelContains.points = true;
 +    }
 +    if (obj.userData.hasOwnProperty("pointlabels")) {
 +        init_pointlabels(obj);
 +        modelContains.pointlabels = true;
 +    }
 +    if (obj.userData.hasOwnProperty("edgematerial")) {
 +        init_lines(obj);
 +        modelContains.lines = true;
 +    }
 +    if (obj.userData.hasOwnProperty("edgelabels")) {
 +        init_edgelabels(obj);
 +        modelContains.edgelabels = true;
 +    }
 +    if (obj.userData.hasOwnProperty("arrowstyle")) {
 +        init_arrowheads(obj);
 +        modelContains.arrowheads = true;
 +    }
 +    if (obj.userData.hasOwnProperty("facetmaterial")) {
 +        init_faces(obj);
 +        modelContains.faces = true;
 +    }
 +}
 +
 +function init_points(obj) {
 +    var pointgroup = new THREE.Group();
 +    pointgroup.name = "points";
 +    var points = obj.userData.points;
 +    var radii = obj.userData.pointradii;
 +    var materials = obj.userData.pointmaterial;
 +    var geometry,material;
 +    if (!Array.isArray(radii)) {
 +        geometry = new THREE.SphereBufferGeometry(radii);  
 +    }
 +    if (!Array.isArray(materials)) {
 +        material = materials;
 +    }
 +    for (var i=0; i<points.length; i++) {
 +        var point = points[i];
 +        if (Array.isArray(radii)) {
 +            if (radii[i] == 0) {
 +                continue;
 +            }
 +            geometry = new THREE.SphereBufferGeometry(radii[i]);  
 +        } 
 +        if (Array.isArray(materials)) {
 +            material = materials[i];     
 +        } 
 +        var sphere = new THREE.Mesh(geometry, material);
 +        point.addSphere(sphere);
 +        pointgroup.add(sphere);
 +    }
 +    obj.add(pointgroup);
 +}
 +
 +function init_pointlabels(obj) {
 +    var points = obj.userData.points;
 +    var labels = obj.userData.pointlabels;
 +    var pointlabels = new THREE.Group();
 +    pointlabels.name = "pointlabels";
 +    if (Array.isArray(labels)) {
 +        for (var i=0; i<points.length; i++) {
 +            var point = points[i];
 +            var spriteMaterial = textSpriteMaterial( labels[i] );
 +         var sprite = new THREE.Sprite(spriteMaterial);
 +            point.addLabel(sprite);
 +            pointlabels.add(sprite);
 +        }
 +    } else {
 +        var spriteMaterial = textSpriteMaterial( labels );
 +        for (var i=0; i<points.length; i++) {
 +            var point = points[i];
 +         var sprite = new THREE.Sprite(spriteMaterial);
 +            point.addLabel(sprite);
 +            pointlabels.add(sprite);
 +        }
 +    }
 +    obj.add(pointlabels);
 +}
 +
 +function init_lines(obj) {
 +    var edgeindices = obj.userData.edgeindices;
 +    var points = obj.userData.points;
 +    var materials = obj.userData.edgematerial;
 +    var geometry = new THREE.BufferGeometry();
 +    var bufarr = new Float32Array( obj.userData.edgeindices.length * 3 );
 +    var bufattr = new THREE.Float32BufferAttribute( bufarr, 3 );
 +    var geometry = new THREE.BufferGeometry();
 +    geometry.setAttribute('position', bufattr);
 +    if (Array.isArray(materials)) {     
 +        for (var i=0; i<materials.length; i++) {
 +            geometry.addGroup(2*i,2,i);
 +        }
 +    }
 +    var lines = new THREE.LineSegments(geometry, materials);
 +    lines.name = "lines";
 +    obj.add(lines);
 +    updateEdgesPosition(obj);
 +}
 +
 +function init_edgelabels(obj) {
 +    var points = obj.userData.points;
 +    var edgeindices = obj.userData.edgeindices;
 +    var labels = obj.userData.edgelabels;
 +    var edgelabels = new THREE.Group();
 +    edgelabels.name = "edgelabels";
 +    if (Array.isArray(labels)) {
 +        for (var i=0; i<edgeindices.length/2; i++) {
 +            var spriteMaterial = textSpriteMaterial( labels[i] );
 +            var sprite = new THREE.Sprite(spriteMaterial);
 +            sprite.position.copy(new THREE.Vector3().addVectors(points[edgeindices[2*i]].vector,points[edgeindices[2*i+1]].vector).multiplyScalar(0.5));
 +            edgelabels.add(sprite);
 +        }
 +    } else {
 +        var spriteMaterial = textSpriteMaterial( labels );
 +        for (var i=0; i<edgeindices.length/2; i++) {
 +            var sprite = new THREE.Sprite(spriteMaterial);
 +            sprite.position.copy(new THREE.Vector3().addVectors(points[edgeindices[2*i]].vector,points[edgeindices[2*i+1]].vector).multiplyScalar(0.5));
 +            edgelabels.add(sprite);
 +        }
 +    }
 +    obj.add(edgelabels);
 +}
 +
 +function init_arrowheads(obj) {
 +    var arrowheads = new THREE.Group();
 +    arrowheads.name = "arrowheads";
 +    var arrowstyle = obj.userData.arrowstyle;
 +    var edgeindices = obj.userData.edgeindices;
 +    var edgematerials = obj.userData.edgematerial;
 +    var points = obj.userData.points;
 +    var material;
 +    if (!Array.isArray(edgematerials)) {
 +        material = new THREE.MeshBasicMaterial( {color: edgematerials.color} );
 +    }
 +
 +    for (var i=0; i<edgeindices.length; i=i+2) {
 +        var start = points[edgeindices[i]];
 +        var end = points[edgeindices[i+1]];
 +        var dist = start.vector.distanceTo( end.vector ) - start.radius() - end.radius();
 +        if (dist <= 0) {
 +            continue;
 +        }
 +        var dir = new THREE.Vector3().subVectors(end.vector,start.vector);
 +        dir.normalize();
 +        var axis = new THREE.Vector3().set(dir.z,0,-dir.x);
 +        axis.normalize();
 +        var radians = Math.acos( dir.y );
 +        var radius = dist/25;
 +        var height = dist/5;
 +        var geometry = new THREE.ConeBufferGeometry(radius,height);
 +        var position = new THREE.Vector3().addVectors(start.vector,dir.clone().multiplyScalar(start.radius()+dist-height/2));
 +        if (Array.isArray(edgematerials)) {
 +            material = new THREE.MeshBasicMaterial( {color: edgematerials[i].color} );
 +        }
 +        var cone = new THREE.Mesh( geometry, material );
 +        cone.quaternion.setFromAxisAngle(axis,radians);;
 +        cone.position.copy(position);;
 +        arrowheads.add(cone);
 +    }
 +    obj.add(arrowheads);
 +}
 +
 +function init_faces(obj) {
 +    var points = obj.userData.points;
 +    var facets = obj.userData.facets;
 +    obj.userData.triangleindices = [];
 +    for (var i=0; i<facets.length; i++) {
 +        facet = facets[i];
 +        for (var t=0; t<facet.length-2; t++) {
 +            obj.userData.triangleindices.push(facet[0],facet[t+1],facet[t+2]);  
 +        }
 +    }
 +    var bufarr = new Float32Array( obj.userData.triangleindices.length * 3 );
 +    var bufattr = new THREE.Float32BufferAttribute(bufarr,3);
 +    
 +    var materials = obj.userData.facetmaterial;
 +    var geometry = new THREE.BufferGeometry();
 +    var frontmaterials = [];
 +    var backmaterials = [];
 +    geometry.setAttribute('position',bufattr);
 +    if (Array.isArray(materials)) {
 +        var tricount = 0;
 +        var facet;
 +        for (var i=0; i<facets.length; i++) {
 +            facet = facets[i];
 +            geometry.addGroup(tricount,(facet.length-2)*3,i);
 +            tricount += (facet.length-2)*3;
 +        }
 +        for (var j=0; j<materials.length; j++) {
 +            var fmat = materials[j].clone()
 +            fmat.side = THREE.FrontSide;
 +            frontmaterials.push(fmat);
 +            var bmat = materials[j].clone()
 +            bmat.side = THREE.BackSide;
 +            backmaterials.push(bmat);
 +            obj.userData.facetmaterial = frontmaterials.concat(backmaterials);
 +        }
 +    } else if (materials instanceof THREE.Material) {
 +        frontmaterials = materials.clone()
 +        frontmaterials.side = THREE.FrontSide;
 +        backmaterials = materials.clone()
 +        backmaterials.side = THREE.BackSide;
 +        obj.userData.facetmaterial = [frontmaterials, backmaterials];
 +    }
 +    // duplicating the object with front and back should avoid transparency issues
 +    var backmesh = new THREE.Mesh(geometry, backmaterials);
 +    // meshname is used to show/hide objects
 +    backmesh.name = "backfaces";
 +    obj.add(backmesh);
 +    var frontmesh = new THREE.Mesh(geometry, frontmaterials);
 +    frontmesh.name = "frontfaces";
 +    obj.add(frontmesh);
 +    updateFacesPosition(obj);
 +}
 +// //INITIALIZING
 +
 +
 +function updateFacesPosition(obj) {
 +    var points = obj.userData.points;
 +    var indices = obj.userData.triangleindices;
 +    var faces = obj.getObjectByName("frontfaces");
 +    var ba = faces.geometry.getAttribute("position");
 +    for (var i=0; i<indices.length; i++) {
 +        ba.setXYZ(i, points[indices[i]].vector.x, points[indices[i]].vector.y ,points[indices[i]].vector.z); 
 +    }
 +    faces.geometry.attributes.position.needsUpdate = true;
 +    
 +}
 +
 +function updateEdgesPosition(obj) {
 +    var points = obj.userData.points;
 +    var indices = obj.userData.edgeindices;
 +    var lines = obj.getObjectByName("lines");
 +    var ba = lines.geometry.getAttribute("position"); 
 +    for (var i=0; i<indices.length; i++) {
 +        ba.setXYZ(i, points[indices[i]].vector.x, points[indices[i]].vector.y ,points[indices[i]].vector.z); 
 +    }
 +    lines.geometry.attributes.position.needsUpdate = true;
 +}
 +
 +function onWindowResize() {
 +    renderer.setSize( three.clientWidth, three.clientHeight );
 +    svgRenderer.setSize( three.clientWidth, three.clientHeight );
 +    updateCamera();
 +}
 +
 +function updateCamera() {
 +    var width = three.clientWidth;
 +    var height = three.clientHeight;
 +    var aspect = width / height;
 +    if (camera.type == "OrthographicCamera") {
 +        camera.left = frustumSize * aspect / - 2;
 +        camera.right = frustumSize * aspect / 2;
 +        camera.top = frustumSize / 2;
 +        camera.bottom = - frustumSize / 2;
 +    } else if (camera.type == "PerspectiveCamera") {
 +        camera.aspect = aspect;
 +    }
 +    camera.updateProjectionMatrix();
 +}
 +
 +function changeCamera(event) {
 +    var selindex = event.currentTarget.selectedIndex;
 +    camera = cameras[selindex];
 +    control = controls[selindex];
 +    control.enabled = true; 
 +    for (var i=0; i<controls.length; i++) {
 +        if (i!=selindex) {
 +            controls[i].enabled = false;
 +        }
 +    }
 +    updateCamera();
 +}
 +
 +var camtypenode = document.getElementById('cameraType_0');
 +camtypenode.onchange = changeCamera;
 +camtypenode.dispatchEvent(new Event('change'));
 +
 +onWindowResize();
 +window.addEventListener('resize', onWindowResize);
 +
 +
 +var xRotationEnabled = false;
 +var yRotationEnabled = false;
 +var zRotationEnabled = false;
 +var rotationSpeedFactor = 1;
 +var settingsShown = false;
 +var labelsShown = true;
 +var intervals = [];
 +var timeouts = [];
 +var explodingSpeed = 0.05;
 +var explodeScale = 0.000001;
 +var XMLS = new XMLSerializer();
 +var svgElement;
 +var renderId;
 +
 +var render = function () {
 +
 + renderId = requestAnimationFrame(render);
 +
 +// comment in for automatic explosion
 +// explode(updateFactor());
 +
 +    var phi = 0.02 * rotationSpeedFactor;
 +
 +    if (xRotationEnabled) {
 +        scene.rotation.x += phi;
 +    }
 +    if (yRotationEnabled) {
 +        scene.rotation.y += phi;
 +    }
 +    if (zRotationEnabled) {
 +        scene.rotation.z += phi;
 +    }
 +
 +    control.update();
 +    renderer.render(scene, camera);
 +};
 +
 +if ( THREE.WEBGL.isWebGLAvailable() ) {
 + render();
 +} else {
 + var warning = WEBGL.getWebGLErrorMessage();
 + three.appendChild( warning );
 +}
 +    
 +function changeTransparency() {
 +    var opacity = 1-Number(event.currentTarget.value);
 +    for (var i=0; i<scene.children.length; i++) {
 +        child = scene.children[i];
 +        if ( child.userData.hasOwnProperty("facetmaterial") ) {
 +            if (Array.isArray(child.userData.facetmaterial)) {
 +                for (var j=0; j<child.userData.facetmaterial.length; j++) {
 +                    child.userData.facetmaterial[j].opacity = opacity;
 +                }
 +            } else {
 +                child.userData.facetmaterial.opacity = opacity;
 +            }    
 +        }
 +    }
 +}
 +
 +function toggleDepthWrite(event) {
 +    depthwrite = event.currentTarget.checked;
 +    for (var i=0; i<scene.children.length; i++) {
 +        child = scene.children[i];
 +        if ( child.userData.hasOwnProperty("facetmaterial") ) {
 +            if (Array.isArray(child.userData.facetmaterial)) {
 +                for (var j=0; j<child.userData.facetmaterial.length; j++) {
 +                    child.userData.facetmaterial[j].depthWrite = depthwrite;
 +                }
 +            } else {
 +                child.userData.facetmaterial.depthWrite = depthWrite;
 +            }    
 +        }
 +    }
 +}
 +
 +function changeRotationX(event){
 +    xRotationEnabled = event.currentTarget.checked;
 +}
 +
 +function changeRotationY(event){
 +    yRotationEnabled = event.currentTarget.checked;
 +}
 +
 +function changeRotationZ(event){
 +    zRotationEnabled = event.currentTarget.checked;
 +}
 +
 +
 +function changeRotationSpeedFactor(event){
 +    rotationSpeedFactor = Number(event.currentTarget.value);
 +}
 +
 +function resetScene(){
 +    scene.rotation.set(0,0,0);
 +    camera.position.set(0,0,5);
 +    camera.up.set(0,1,0);
 +}
 +
 +function showSettings(event){
 +    document.getElementById('settings_0').style.visibility = 'visible';
 +    document.getElementById('showSettingsButton_0').style.visibility = 'hidden';
 +    document.getElementById('hideSettingsButton_0').style.visibility = 'visible';
 +    settingsShown = true;
 +}
 +
 +function hideSettings(event){
 +    document.getElementById('settings_0').style.visibility = 'hidden';
 +    document.getElementById('showSettingsButton_0').style.visibility = 'visible';
 +    document.getElementById('hideSettingsButton_0').style.visibility = 'hidden';
 +    settingsShown = false;
 +}
 +
 +
 +
 +var pos = 150* Math.PI;
 +
 +function updateFactor() {
 +    pos++;
 +    return Math.sin(.01*pos)+1;
 +}
 +
 +// ------------------------ FOLDING ------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +// rotate point p around axis defined by points p1 and p2 by given angle
 +function rotate(p, p1, p2, angle ){   
 +    angle = -angle;
 +    var x = p.x, y = p.y, z = p.z, 
 +    a = p1.x, b = p1.y, c = p1.z, 
 +    u = p2.x-p1.x, v = p2.y-p1.y, w = p2.z-p1.z;
 +    var result = [];
 +    var L = u*u + v*v + w*w;
 +    var sqrt = Math.sqrt;
 +    var cos = Math.cos;
 +    var sin = Math.sin;
 +
 +    result[0] = ((a*(v*v+w*w)-u*(b*v+c*w-u*x-v*y-w*z))*(1-cos(angle))+L*x*cos(angle)+sqrt(L)*(-c*v+b*w-w*y+v*z)*sin(angle))/L;
 +    result[1] = ((b*(u*u+w*w)-v*(a*u+c*w-u*x-v*y-w*z))*(1-cos(angle))+L*y*cos(angle)+sqrt(L)*(c*u-a*w+w*x-u*z)*sin(angle))/L;
 +    result[2] = ((c*(u*u+v*v)-w*(a*u+b*v-u*x-v*y-w*z))*(1-cos(angle))+L*z*cos(angle)+sqrt(L)*(-b*u+a*v-v*x+u*y)*sin(angle))/L;
 +
 +    return result;
 +}
 +
 +var fold = function(event){
 +    var obj = foldables[Number(event.currentTarget.name)];
 +    var foldvalue = Number(event.currentTarget.value);
 +    var scale = foldvalue - obj.userData.oldscale;
 +
 +    for (var j=0; j<obj.userData.axes.length; j++) {
 +        rotateVertices(obj, j, scale);
 +    }
 +    update(obj);
 +    obj.userData.oldscale += scale;
 +    lookAtBarycenter(obj);
 +}
 +
 +function lookAtBarycenter(obj){
 +    control.target = barycenter(obj);
 +}
 +
 +function barycenter(obj) {
 +    var center = new THREE.Vector3(0,0,0);
 +    var points = obj.userData.points;
 +    for (var i=0; i<points.length; i++){
 +        center.add(points[i].vector);
 +    }
 +    center.divideScalar(points.length);
 +    return center;
 +}
 +
 +function rotateVertices(obj, edge, scale) {
 +    var axes = obj.userData.axes;
 +    var subtrees = obj.userData.subtrees;
 +    var points = obj.userData.points;
 +    var angles = obj.userData.angles;
 +    if (edge < axes.length){
 +        for (var j=0; j<subtrees[edge].length; j++){
 +            var rotP = rotate(points[subtrees[edge][j]].vector, points[axes[edge][0]].vector,points[axes[edge][1]].vector, scale * (Math.PI - angles[edge]));
 +            points[subtrees[edge][j]].set(rotP[0],rotP[1],rotP[2]);
 +        }
 +    }
 +}
 +
 +function update(obj) {
 +   updateFacesPosition(obj);
 +   updateEdgesPosition(obj);
 +}
 +
 +if (foldables.length) {
 +    var settings = document.getElementById('settings_0');
 +    var foldDiv = document.createElement('div');
 +    foldDiv.id = 'fold_0';
 +    var title = document.createElement('strong');
 +    title.innerHTML = 'Fold';
 +    foldDiv.appendChild(title);
 +    foldDiv.className = 'group';
 +    for (var i=0; i<foldables.length; i++) {
 +        var range = document.createElement('input');
 +        range.type = 'range';
 +        range.min = 0;
 +        range.max = 1;
 +        range.value = 0;
 +        range.step = 0.001;
 +        range.name = String(i);
 +        range.oninput = fold;
 +        foldDiv.appendChild(range);
 +    }
 +    lookAtBarycenter(foldables[0]);
 +    settings.insertBefore(foldDiv,settings.childNodes[0]);
 +}
 +
 +    
 +// ---------------------- EXPLOSION ------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +if (explodableModel) {
 +    for (var i=0; i<scene.children.length; i++) {
 +        obj = scene.children[i];
 +        if ( obj.userData.explodable ) {
 +            computeCentroid(obj);
 +        }
 +    }
 +    document.getElementById('explodeRange_0').oninput = triggerExplode;
 +    document.getElementById('explodeCheckbox_0').onchange = triggerAutomaticExplode;
 +    document.getElementById('explodingSpeedRange_0').oninput = setExplodingSpeed;
 +    explode(0.000001);
 +}
 +
 +function computeCentroid(obj) {
 +    centroid = new THREE.Vector3();
 +    obj.userData.points.forEach(function(pmpoint) {
 +        centroid.add(pmpoint.vector);
 +    });
 +    centroid.divideScalar(obj.userData.points.length);
 +    obj.userData.centroid = centroid;
 +}
 +
 +function explode(factor) {
 +    for (var i=0; i<scene.children.length; i++) {
 +        var obj = scene.children[i];
 +        if (obj.userData.hasOwnProperty("centroid")) { 
 +            var c = obj.userData.centroid;
 +            obj.position.set(c.x*factor, c.y*factor, c.z*factor);
 +        }
 +    }
 +}
 +
 +function triggerExplode(event){
 +    explodeScale = Number(event.currentTarget.value);
 +    explode(explodeScale);
 +}
 +
 +function setExplodingSpeed(event){
 +    explodingSpeed = Number(event.currentTarget.value);
 +}
 +
 +function triggerAutomaticExplode(event){
 +    if (event.currentTarget.checked){
 +        startExploding();
 +    } else {
 +        clearIntervals();
 +    }
 +}
 +
 +function startExploding(){
 +    intervals.push(setInterval(explodingInterval, 25));
 +}
 +
 +
 +function explodingInterval(){
 +    explodeScale += explodingSpeed;
 +    if (explodeScale <= 6){ 
 +        explode(explodeScale);
 +    }
 +    else{
 +        explode(6);
 +        explodeScale = 6;
 +        clearIntervals();
 +        timeouts.push(setTimeout(startUnexploding, 3000));
 +    }
 +    document.getElementById('explodeRange_0').value = explodeScale;
 +}
 +
 +
 +function startUnexploding(){
 +    intervals.push(setInterval(unexplodingInterval, 25));
 +}
 +
 +function unexplodingInterval(){
 +    explodeScale -= explodingSpeed;
 +    if (explodeScale >= 0){
 +        explode(explodeScale);
 +    }
 +    else {
 +        explode(0);
 +        explodeScale = 0;
 +        clearIntervals();
 +        timeouts.push(setTimeout(startExploding, 3000));
 +    }
 +    document.getElementById('explodeRange_0').value = explodeScale;
 +}
 +
 +function clearIntervals(){
 +    intervals.forEach(function(interval){
 +        clearInterval(interval);
 +    });
 +    intervals = [];
 +    timeouts.forEach(function(timeout){
 +        clearTimeout(timeout);
 +    });
 +    timeouts = [];
 +}
 +
 +// ---------------------- DISPLAY --------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +const objectTypeInnerHTMLs = { points: "Points", pointlabels: "Point labels", lines: "Edges", edgelabels: "Edge labels", faces: "Faces", arrowheads: "Arrow heads" };
 +const objectTypeVisible = {};
 +Object.assign(objectTypeVisible,modelContains);
 +const sortedObjectTypeKeys = Object.keys(objectTypeInnerHTMLs).sort();
 +const shownObjectTypesList = document.getElementById('shownObjectTypesList_0');
 +
 +function setVisibility(bool,objname) {
 +    for (var i=0; i<scene.children.length; i++){
 +        var obj = scene.children[i].getObjectByName(objname);
 +        if (obj) {
 +            obj.visible = bool;
 +        }
 +    }
 +}
 +
 +function toggleObjectTypeVisibility(event){
 +    var name = event.currentTarget.name;
 +    var checked = event.currentTarget.checked;
 +    objectTypeVisible[name] = checked;
 +    if (name == "faces") {
 +        setVisibility(checked,"frontfaces");
 +        setVisibility(checked,"backfaces");
 +    } else {
 +        setVisibility(checked,name);
 +    }
 +}
 +
 +for (var i=0; i<sortedObjectTypeKeys.length; i++){
 +    var key = sortedObjectTypeKeys[i];
 +    if (modelContains[key]) {
 +        var objTypeNode = document.createElement('span');
 +        objTypeNode.innerHTML = objectTypeInnerHTMLs[key] + '<br>';
 +        var checkbox = document.createElement('input');
 +        checkbox.type = 'checkbox';
 +        checkbox.checked = true;
 +        checkbox.name = key;
 +        checkbox.onchange = toggleObjectTypeVisibility;
 +        shownObjectTypesList.appendChild(checkbox);
 +        shownObjectTypesList.appendChild(objTypeNode);
 +    }
 +}
 +
 +// ------------------------------------------------------
 +
 +function toggleObjectVisibility(event){
 +    var nr = Number(event.currentTarget.name);
 +    scene.children[nr].visible = event.currentTarget.checked;
 +}
 +
 +// append checkboxes for displaying or hiding objects
 +var shownObjectsList = document.getElementById('shownObjectsList_0');
 +for (var i=0; i<scene.children.length; i++){
 +    obj = scene.children[i];
 +    var objNode = document.createElement('span');
 +    objNode.innerHTML = obj.name + '<br>';
 +    var checkbox = document.createElement('input');
 +    checkbox.type = 'checkbox';
 +    checkbox.checked = true;
 +    checkbox.name = String(i);
 +    checkbox.onchange = toggleObjectVisibility;
 +    shownObjectsList.appendChild(checkbox);
 +    shownObjectsList.appendChild(objNode);
 +}
 +
 +// ---------------------- SVG ------------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +function takeSvgScreenshot() {
 +    if (objectTypeVisible["pointlabels"]) {
 +        setVisibility(false,"pointlabels");
 +    }
 +    if (objectTypeVisible["edgelabels"]) {
 +        setVisibility(false,"edgelabels");
 +    }
 +    svgRenderer.render(scene,camera);
 +    svgElement = XMLS.serializeToString(svgRenderer.domElement);
 +    
 +    if (objectTypeVisible["pointlabels"]) {
 +        setVisibility(true,"pointlabels");
 +    }
 +    if (objectTypeVisible["edgelabels"]) {
 +        setVisibility(true,"edgelabels");
 +    }
 +
 +    if (document.getElementById('tab_0').checked){
 +        //show in new tab
 +        var myWindow = window.open("","");
 +        myWindow.document.body.innerHTML = svgElement;
 +    } else{
 +        // download svg file 
 +        download("screenshot.svg", svgElement);
 +    }
 +}
 +
 +function download(filename, text) {
 +    var element = document.createElement('a');
 +    element.setAttribute('href', 'data:text/plain;charset=utf-8,' + encodeURIComponent(text));
 +    element.setAttribute('download', filename);
 +
 +    element.style.display = 'none';
 +    document.body.appendChild(element);
 +
 +    element.click();
 +
 +    document.body.removeChild(element);
 +}
 +
 +
 +document.getElementById('transparencyRange_0').oninput = changeTransparency;
 +document.getElementById('depthWriteCheckbox_0').onchange = toggleDepthWrite;
 +document.getElementById('changeRotationX_0').onchange = changeRotationX;
 +document.getElementById('changeRotationY_0').onchange = changeRotationY;
 +document.getElementById('changeRotationZ_0').onchange = changeRotationZ;
 +document.getElementById('resetButton_0').onclick = resetScene;
 +document.getElementById('rotationSpeedRange_0').oninput = changeRotationSpeedFactor;
 +document.getElementById('takeScreenshot_0').onclick = takeSvgScreenshot;
 +document.getElementById('showSettingsButton_0').onclick = showSettings;
 +document.getElementById('hideSettingsButton_0').onclick = hideSettings;
 +
 +
 +// ------------------ SHORTCUTS --------------------------------------------
 +// -------------------------------------------------------------------------
 +
 +/**
 + * http://www.openjs.com/scripts/events/keyboard_shortcuts/
 + * Version : 2.01.B
 + * By Binny V A
 + * License : BSD
 + */
 +shortcut = {
 + 'all_shortcuts':{},//All the shortcuts are stored in this array
 + 'add': function(shortcut_combination,callback,opt) {
 + //Provide a set of default options
 + var default_options = {
 + 'type':'keydown',
 + 'propagate':false,
 + 'disable_in_input':false,
 + 'target':document,
 + 'keycode':false
 + }
 + if(!opt) opt = default_options;
 + else {
 + for(var dfo in default_options) {
 + if(typeof opt[dfo] == 'undefined') opt[dfo] = default_options[dfo];
 + }
 + }
 +
 + var ele = opt.target;
 + if(typeof opt.target == 'string') ele = document.getElementById(opt.target);
 + var ths = this;
 + shortcut_combination = shortcut_combination.toLowerCase();
 +
 + //The function to be called at keypress
 + var func = function(e) {
 + e = e || window.event;
 +
 + if(opt['disable_in_input']) { //Don't enable shortcut keys in Input, Textarea fields
 + var element;
 + if(e.target) element=e.target;
 + else if(e.srcElement) element=e.srcElement;
 + if(element.nodeType==3) element=element.parentNode;
 +
 + if(element.tagName == 'INPUT' || element.tagName == 'TEXTAREA') return;
 + }
 +
 + //Find Which key is pressed
 + if (e.keyCode) code = e.keyCode;
 + else if (e.which) code = e.which;
 + var character = String.fromCharCode(code).toLowerCase();
 +
 + if(code == 188) character=","; //If the user presses , when the type is onkeydown
 + if(code == 190) character="."; //If the user presses , when the type is onkeydown
 +
 + var keys = shortcut_combination.split("+");
 + //Key Pressed - counts the number of valid keypresses - if it is same as the number of keys, the shortcut function is invoked
 + var kp = 0;
 +
 + //Work around for stupid Shift key bug created by using lowercase - as a result the shift+num combination was broken
 + var shift_nums = {
 + "`":"~",
 + "1":"!",
 + "2":"@",
 + "3":"#",
 + "4":"$",
 + "5":"%",
 + "6":"^",
 + "7":"&",
 + "8":"*",
 + "9":"(",
 + "0":")",
 + "-":"_",
 + "=":"+",
 + ";":":",
 + "'":"\"",
 + ",":"<",
 + ".":">",
 + "/":"?",
 + "\\":"|"
 + }
 + //Special Keys - and their codes
 + var special_keys = {
 + 'esc':27,
 + 'escape':27,
 + 'tab':9,
 + 'space':32,
 + 'return':13,
 + 'enter':13,
 + 'backspace':8,
 +
 + 'scrolllock':145,
 + 'scroll_lock':145,
 + 'scroll':145,
 + 'capslock':20,
 + 'caps_lock':20,
 + 'caps':20,
 + 'numlock':144,
 + 'num_lock':144,
 + 'num':144,
 +
 + 'pause':19,
 + 'break':19,
 +
 + 'insert':45,
 + 'home':36,
 + 'delete':46,
 + 'end':35,
 +
 + 'pageup':33,
 + 'page_up':33,
 + 'pu':33,
 +
 + 'pagedown':34,
 + 'page_down':34,
 + 'pd':34,
 +
 + 'left':37,
 + 'up':38,
 + 'right':39,
 + 'down':40,
 +
 + 'f1':112,
 + 'f2':113,
 + 'f3':114,
 + 'f4':115,
 + 'f5':116,
 + 'f6':117,
 + 'f7':118,
 + 'f8':119,
 + 'f9':120,
 + 'f10':121,
 + 'f11':122,
 + 'f12':123
 + }
 +
 + var modifiers = { 
 + shift: { wanted:false, pressed:false},
 + ctrl : { wanted:false, pressed:false},
 + alt  : { wanted:false, pressed:false},
 + meta : { wanted:false, pressed:false} //Meta is Mac specific
 + };
 +                        
 + if(e.ctrlKey) modifiers.ctrl.pressed = true;
 + if(e.shiftKey) modifiers.shift.pressed = true;
 + if(e.altKey) modifiers.alt.pressed = true;
 + if(e.metaKey)   modifiers.meta.pressed = true;
 +                        
 + for(var i=0; k=keys[i],i<keys.length; i++) {
 + //Modifiers
 + if(k == 'ctrl' || k == 'control') {
 + kp++;
 + modifiers.ctrl.wanted = true;
 +
 + } else if(k == 'shift') {
 + kp++;
 + modifiers.shift.wanted = true;
 +
 + } else if(k == 'alt') {
 + kp++;
 + modifiers.alt.wanted = true;
 + } else if(k == 'meta') {
 + kp++;
 + modifiers.meta.wanted = true;
 + } else if(k.length > 1) { //If it is a special key
 + if(special_keys[k] == code) kp++;
 +
 + } else if(opt['keycode']) {
 + if(opt['keycode'] == code) kp++;
 +
 + } else { //The special keys did not match
 + if(character == k) kp++;
 + else {
 + if(shift_nums[character] && e.shiftKey) { //Stupid Shift key bug created by using lowercase
 + character = shift_nums[character]; 
 + if(character == k) kp++;
 + }
 + }
 + }
 + }
 +
 + if(kp == keys.length && 
 + modifiers.ctrl.pressed == modifiers.ctrl.wanted &&
 + modifiers.shift.pressed == modifiers.shift.wanted &&
 + modifiers.alt.pressed == modifiers.alt.wanted &&
 + modifiers.meta.pressed == modifiers.meta.wanted) {
 + callback(e);
 +
 + if(!opt['propagate']) { //Stop the event
 + //e.cancelBubble is supported by IE - this will kill the bubbling process.
 + e.cancelBubble = true;
 + e.returnValue = false;
 +
 + //e.stopPropagation works in Firefox.
 + if (e.stopPropagation) {
 + e.stopPropagation();
 + e.preventDefault();
 + }
 + return false;
 + }
 + }
 + }
 + this.all_shortcuts[shortcut_combination] = {
 + 'callback':func, 
 + 'target':ele, 
 + 'event': opt['type']
 + };
 + //Attach the function with the event
 + if(ele.addEventListener) ele.addEventListener(opt['type'], func, false);
 + else if(ele.attachEvent) ele.attachEvent('on'+opt['type'], func);
 + else ele['on'+opt['type']] = func;
 + },
 +
 + //Remove the shortcut - just specify the shortcut and I will remove the binding
 + 'remove':function(shortcut_combination) {
 + shortcut_combination = shortcut_combination.toLowerCase();
 + var binding = this.all_shortcuts[shortcut_combination];
 + delete(this.all_shortcuts[shortcut_combination])
 + if(!binding) return;
 + var type = binding['event'];
 + var ele = binding['target'];
 + var callback = binding['callback'];
 +
 + if(ele.detachEvent) ele.detachEvent('on'+type, callback);
 + else if(ele.removeEventListener) ele.removeEventListener(type, callback, false);
 + else ele['on'+type] = false;
 + }
 +}
 +
 +shortcut.add("Alt+Left",function() {
 + var event = new Event('click');
 + if (settingsShown){
 + document.getElementById('hideSettingsButton_0').dispatchEvent(event);
 + } else {
 + document.getElementById('showSettingsButton_0').dispatchEvent(event);
 + }
 +});
 +
 +
 +// COMMON_CODE_BLOCK_END
 +
 +});});
 +      </script>
 +   </body>
 +</html>
 +</HTML>
 +The minimal and maximal faces can be visualized via
 +
 +<code perl>
 +> $p->VISUAL->MIN_MAX_FACE;
 +</code>
 +<HTML>
 +<!--
 +polymake for knusper
 +Thu Mar  3 00:38:23 2022
 +c3t
 +-->
 +
 +
 +<html>
 +   <head>
 +      <meta charset=utf-8>
 +      <title>c3t</title>
 +      <style>
 +/*
 +// COMMON_CODE_BLOCK_BEGIN
 +*/
 +         html {overflow: scroll;}
 +         strong{font-size: 18px;}
 +         canvas { z-index: 8; }
 +         input[type='radio'] {margin-left:0;}
 +         input[type='checkbox'] {margin-right:7px; margin-left: 0px; padding-left:0px;}
 +         .group{padding-bottom: 15px;}
 +         .settings * {z-index: 11; }
 +         .settings{z-index: 10; font-family: Arial, Helvetica, sans-serif; margin-left: 30px; visibility: hidden; width: 14em; height: 96%; border: solid 1px silver; padding: 2px; overflow-y: scroll; box-sizing: border-box; background-color: white; position: absolute;}
 +         .indented{margin-left: 20px; margin-top: 10px; padding-bottom: 0px;} 
 +         .shownObjectsList{overflow: auto; max-width: 150px; max-height: 150px;}
 +         .showSettingsButton{visibility: visible; z-index: 12; position: absolute }
 +         .hideSettingsButton{visibility: hidden; z-index: 12; position: absolute; opacity: 0.5}
 +         button{margin-left: 0; margin-top: 10px}
 +         img{cursor: pointer;}
 +         .suboption{padding-top: 15px;}
 +         #model67891861409 { width: 100%; height: 100%; }
 +         .threejs_container { width: 100%; height: 75vh;}
 +         .settings{max-height: 74vh} 
 +         input[type=range] {
 +           -webkit-appearance: none;
 +           padding:0; 
 +           width:90%; 
 +           margin-left: auto;
 +           margin-right: auto;
 +           margin-top: 15px;
 +           margin-bottom: 15px;
 +           display: block;
 +         }
 +         input[type=range]:focus {
 +           outline: none;
 +         }
 +         input[type=range]::-webkit-slider-runnable-track {
 +           height: 4px;
 +           cursor: pointer;
 +           animate: 0.2s;
 +           box-shadow: 0px 0px 0px #000000;
 +           background: #E3E3E3;
 +           border-radius: 0px;
 +           border: 0px solid #000000;
 +         }
 +         input[type=range]::-webkit-slider-thumb {
 +           box-shadow: 1px 1px 2px #B8B8B8;
 +           border: 1px solid #ABABAB;
 +           height: 13px;
 +           width: 25px;
 +           border-radius: 20px;
 +           background: #E0E0E0;
 +           cursor: pointer;
 +           -webkit-appearance: none;
 +           margin-top: -5px;
 +         }
 +         input[type=range]:focus::-webkit-slider-runnable-track {
 +           background: #E3E3E3;
 +         }
 +         input[type=range]::-moz-range-track {
 +           height: 4px;
 +           cursor: pointer;
 +           animate: 0.2s;
 +           box-shadow: 0px 0px 0px #000000;
 +           background: #E3E3E3;
 +           border-radius: 0px;
 +           border: 0px solid #000000;
 +         }
 +         input[type=range]::-moz-range-thumb {
 +           box-shadow: 1px 1px 2px #B8B8B8;
 +           border: 1px solid #ABABAB;
 +           height: 13px;
 +           width: 25px;
 +           border-radius: 20px;
 +           background: #E0E0E0;
 +           cursor: pointer;
 +         }
 +         input[type=range]::-ms-track {
 +           height: 4px;
 +           cursor: pointer;
 +           animate: 0.2s;
 +           background: transparent;
 +           border-color: transparent;
 +           color: transparent;
 +         }
 +         input[type=range]::-ms-fill-lower {
 +           background: #E3E3E3;
 +           border: 0px solid #000000;
 +           border-radius: 0px;
 +           box-shadow: 0px 0px 0px #000000;
 +         }
 +         input[type=range]::-ms-fill-upper {
 +           background: #E3E3E3;
 +           border: 0px solid #000000;
 +           border-radius: 0px;
 +           box-shadow: 0px 0px 0px #000000;
 +         }
 +         input[type=range]::-ms-thumb {
 +           box-shadow: 1px 1px 2px #B8B8B8;
 +           border: 1px solid #ABABAB;
 +           height: 13px;
 +           width: 25px;
 +           border-radius: 20px;
 +           background: #E0E0E0;
 +           cursor: pointer;
 +         }
 +         input[type=range]:focus::-ms-fill-lower {
 +           background: #E3E3E3;
 +         }
 +         input[type=range]:focus::-ms-fill-upper {
 +           background: #E3E3E3;
 +         }
 +/*
 +// COMMON_CODE_BLOCK_END
 +*/
 + </style>
 +   </head>
 +<body>
 +   <div class='threejs_container'>
 + <div id='settings_1' class='settings'>
 + <div class=group id='transparency_1' class='transparency'>
 + <strong>Transparency</strong>
 + <input id='transparencyRange_1' type='range' min=0 max=1 step=0.01 value=0>
 +            <div class=indented><input id='depthWriteCheckbox_1' type='checkbox'>depthWrite</div>
 + </div>
 +
 + <div class=group id='rotation_1'>
 + <strong>Rotation</strong>
 + <div class=indented>
 + <div><input type='checkbox' id='changeRotationX_1'> x-axis</div>
 + <div><input type='checkbox' id='changeRotationY_1'> y-axis</div>
 + <div><input type='checkbox' id='changeRotationZ_1'> z-axis</div>
 + <button id='resetButton_1'>Reset</button>
 + </div>
 +
 + <div class=suboption>Rotation speed</div>
 + <input id='rotationSpeedRange_1' type='range' min=0 max=5 step=0.01 value=2>
 + </div>
 +
 +
 + <div class=group id='display_1'>
 + <strong>Display</strong>
 + <div class=indented>
 + <div id='shownObjectTypesList_1' class='shownObjectsList'></div>
 + </div>
 + <div class=suboption>Objects</div>
 + <div class=indented>
 +    <div id='shownObjectsList_1' class='shownObjectsList'></div>
 + </div>
 + </div>
 +         
 +         <div class=group id='camera_1'>
 +            <strong>Camera</strong>
 +            <div class=indented>
 +               <form>
 +                  <select id="cameraType_1">
 +                     <option value='perspective' selected> Perspective<br></option>
 +                     <option value='orthographic' > Orthographic<br></option>
 +                  </select>
 +               </form>
 +            </div>
 +         </div>
 +
 + <div class=group id='svg_1'>
 + <strong>SVG</strong>
 + <div class=indented>
 + <form>
 + <input type="radio" name='screenshotMode' value='download' id='download_1' checked> Download<br>
 + <input type="radio" name='screenshotMode' value='tab' id='tab_1' > New tab<br>
 + </form>
 + <button id='takeScreenshot_1'>Screenshot</button>
 + </div>
 + </div>
 +
 + </div> <!-- end of settings -->
 + <img id='hideSettingsButton_1' class='hideSettingsButton' src='/kernelspecs/r118/polymake/close.svg' width=20px">
 + <img id='showSettingsButton_1' class='showSettingsButton' src='/kernelspecs/r118/polymake/menu.svg' width=20px">
 +<div id="model67891861409"></div>
 +</div>
 +   <script>
 +    requirejs.config({
 +      paths: {
 +        three: '/kernelspecs/r118/polymake/three',
 +        TrackballControls: '/kernelspecs/r118/polymake/TrackballControls',
 +        OrbitControls: '/kernelspecs/r118/polymake/OrbitControls',
 +        Projector: '/kernelspecs/r118/polymake/Projector',
 +        SVGRenderer: '/kernelspecs/r118/polymake/SVGRenderer',
 +        WEBGL: '/kernelspecs/r118/polymake/WebGL',
 +      },
 +      shim: {
 +        'three': { exports: 'THREE'},
 +        'SVGRenderer': { deps: [ 'three' ], exports: 'THREE.SVGRenderer' },
 +        'WEBGL': { deps: [ 'three' ], exports: 'THREE.WEBGL' },
 +        'Projector': { deps: [ 'three' ], exports: 'THREE.Projector' },
 +        'TrackballControls': { deps: [ 'three' ], exports: 'THREE.TrackballControls' },
 +        'OrbitControls': { deps: [ 'three' ], exports: 'THREE.OrbitControls' },
 +      }
 +    });
 +    
 +    require(['three'],function(THREE){
 +        window.THREE = THREE;
 +      require(['TrackballControls', 'OrbitControls', 'Projector', 'SVGRenderer', 'WEBGL'],
 +               function(TrackballControls, OrbitControls, Projector, SVGRenderer, WEBGL) {
 +    THREE.TrackballControls = TrackballControls;
 +    THREE.OrbitControls = OrbitControls;
 +    THREE.Projector = Projector;
 +    THREE.SVGRenderer = SVGRenderer;
 +    THREE.WEBGL = WEBGL;
 +
 +// COMMON_CODE_BLOCK_BEGIN
 +
 +const intervalLength = 25; // for automatic animations
 +const explodableModel = false; 
 +const modelContains = { points: false, pointlabels: false, lines: false, edgelabels: false, faces: false, arrowheads: false };
 +const foldables = [];
 +
 +var three = document.getElementById("model67891861409");
 +var scene = new THREE.Scene();
 +var renderer = new THREE.WebGLRenderer( { antialias: true } );
 +var svgRenderer = new THREE.SVGRenderer( { antialias: true } );
 +renderer.setPixelRatio( window.devicePixelRatio );
 +renderer.setClearColor(0xFFFFFF, 1);
 +svgRenderer.setClearColor(0xFFFFFF, 1);
 +three.appendChild(renderer.domElement);
 +
 +var frustumSize = 4;
 +var cameras = [new THREE.PerspectiveCamera(75, 1, 0.1, 1000), new THREE.OrthographicCamera()];
 +cameras.forEach(function(cam) {
 +    cam.position.set(0, 0, 5);
 +    cam.lookAt(0, 0, 0);  
 +    cam.up.set(0, 1, 0);         
 +});
 +var controls = [new THREE.TrackballControls(cameras[0], three), new THREE.OrbitControls(cameras[1], three)];
 +var camera, control;
 +
 +controls[0].zoomSpeed = 0.2;
 +controls[0].rotateSpeed = 4;
 +
 +
 +// class to allow move points together with labels and spheres
 +var PMPoint = function (x,y,z) {
 +   this.vector = new THREE.Vector3(x,y,z);
 +   this.sprite = null;
 +   this.sphere = null;
 +}
 +PMPoint.prototype.addLabel = function(labelsprite) {
 +   this.sprite = labelsprite;
 +   this.sprite.position.copy(this.vector);
 +}
 +PMPoint.prototype.addSphere = function(spheremesh) {
 +   this.sphere = spheremesh;
 +   this.sphere.position.copy(this.vector);
 +}
 +PMPoint.prototype.set = function(x,y,z) {
 +   this.vector.set(x,y,z);
 +   if (this.sprite) {
 +      this.sprite.position.copy(this.vector);
 +   }
 +   if (this.sphere) {
 +      this.sphere.position.copy(this.vector);
 +   }
 +}
 +PMPoint.prototype.radius = function() {
 +   if (this.sphere) {
 +      return this.sphere.geometry.parameters.radius;
 +   } else {
 +      return 0;
 +   }
 +};
 +// select the target node
 +var target = document.querySelector('#model67891861409');
 +
 +// create an observer instance
 +var observer = new MutationObserver(function(mutations) {
 +   mutations.forEach(function(mutation) {
 +      if (mutation.removedNodes && mutation.removedNodes.length > 0) {
 +         cancelAnimationFrame(renderId);
 +         observer.disconnect();
 +         console.log("cancelled frame "+renderId);
 +      }
 +   });
 +});
 +
 +// configuration of the observer:
 +var config = { childList: true, characterData: true }
 +
 +// pass in the target node, as well as the observer options
 +while (target) {
 +   if (target.className=="output") {
 +      observer.observe(target, config);
 +      break;
 +   }
 +   target = target.parentNode;
 +}
 +
 +// COMMON_CODE_BLOCK_END
 +
 +var obj0 = new THREE.Object3D();
 +obj0.name = "c3t";
 +obj0.userData.explodable = 1;
 +obj0.userData.points = [];
 +obj0.userData.points.push(new PMPoint(0, 0, 1));
 +obj0.userData.points.push(new PMPoint(1, 0, 0));
 +obj0.userData.points.push(new PMPoint(0, 0, 0));
 +obj0.userData.points.push(new PMPoint(0, 1, 0));
 +obj0.userData.points.push(new PMPoint(1, 1, 0));
 +obj0.userData.points.push(new PMPoint(1, 0, 1));
 +obj0.userData.points.push(new PMPoint(0, 1, 1));
 +
 +obj0.userData.pointradii = 0.02;
 +   <!-- Vertex style -->
 +obj0.userData.pointmaterial = [new THREE.MeshBasicMaterial( { color: 0xFF0000, side: THREE.DoubleSide, transparent: false } ),
 +new THREE.MeshBasicMaterial( { color: 0xFF0000, side: THREE.DoubleSide, transparent: false } ),
 +new THREE.MeshBasicMaterial( { color: 0xFFFF00, side: THREE.DoubleSide, transparent: false } ),
 +new THREE.MeshBasicMaterial( { color: 0xFF0000, side: THREE.DoubleSide, transparent: false } ),
 +new THREE.MeshBasicMaterial( { color: 0xFF0000, side: THREE.DoubleSide, transparent: false } ),
 +new THREE.MeshBasicMaterial( { color: 0xFF0000, side: THREE.DoubleSide, transparent: false } ),
 +new THREE.MeshBasicMaterial( { color: 0xFF0000, side: THREE.DoubleSide, transparent: false } )];
 +obj0.userData.pointlabels = ["0", "1", "2", "3", "4", "5", "6"];
 +obj0.userData.edgeindices = [0, 2, 1, 2, 2, 3, 1, 4, 3, 4, 0, 5, 1, 5, 4, 5, 0, 6, 3, 6, 4, 6, 5, 6];
 +   <!-- Edge style -->
 +obj0.userData.edgematerial = new THREE.LineBasicMaterial( { color: 0x000000, linewidth: 1.5, transparent: false } );
 +obj0.userData.facets = [[5, 4, 6], [6, 3, 2, 0], [1, 4, 5], [5, 0, 2, 1], [6, 4, 3], [3, 4, 1, 2], [5, 6, 0]];
 +   <!-- Facet style -->
 +obj0.userData.facetmaterial = [new THREE.MeshBasicMaterial( { color: 0xFF0000, depthFunc: THREE.LessDepth, depthWrite: false, opacity: 1, polygonOffset: true, polygonOffsetFactor: 1, polygonOffsetUnits: 0.5, side: THREE.DoubleSide, transparent: true } ),
 +new THREE.MeshBasicMaterial( { color: 0x77EC9E, depthFunc: THREE.LessDepth, depthWrite: false, opacity: 1, polygonOffset: true, polygonOffsetFactor: 1, polygonOffsetUnits: 0.5, side: THREE.DoubleSide, transparent: true } ),
 +new THREE.MeshBasicMaterial( { color: 0x77EC9E, depthFunc: THREE.LessDepth, depthWrite: false, opacity: 1, polygonOffset: true, polygonOffsetFactor: 1, polygonOffsetUnits: 0.5, side: THREE.DoubleSide, transparent: true } ),
 +new THREE.MeshBasicMaterial( { color: 0x77EC9E, depthFunc: THREE.LessDepth, depthWrite: false, opacity: 1, polygonOffset: true, polygonOffsetFactor: 1, polygonOffsetUnits: 0.5, side: THREE.DoubleSide, transparent: true } ),
 +new THREE.MeshBasicMaterial( { color: 0x77EC9E, depthFunc: THREE.LessDepth, depthWrite: false, opacity: 1, polygonOffset: true, polygonOffsetFactor: 1, polygonOffsetUnits: 0.5, side: THREE.DoubleSide, transparent: true } ),
 +new THREE.MeshBasicMaterial( { color: 0x77EC9E, depthFunc: THREE.LessDepth, depthWrite: false, opacity: 1, polygonOffset: true, polygonOffsetFactor: 1, polygonOffsetUnits: 0.5, side: THREE.DoubleSide, transparent: true } ),
 +new THREE.MeshBasicMaterial( { color: 0x77EC9E, depthFunc: THREE.LessDepth, depthWrite: false, opacity: 1, polygonOffset: true, polygonOffsetFactor: 1, polygonOffsetUnits: 0.5, side: THREE.DoubleSide, transparent: true } )];
 +init_object(obj0);
 +scene.add(obj0);
 +
 +// COMMON_CODE_BLOCK_BEGIN
 +function textSpriteMaterial(message, parameters) {
 +    if ( parameters === undefined ) parameters = {};
 +    var fontface = "Helvetica";
 +    var fontsize = parameters.hasOwnProperty("fontsize") ? parameters["fontsize"] : 15;
 +    fontsize = fontsize*10;
 +    var lines = message.split('\\n');
 +    var size = 512;
 +    for(var i = 0; i<lines.length; i++){
 +        var tmp = lines[i].length;
 +        while(tmp*fontsize > size){
 +           fontsize--;
 +        }
 +    }
 +    
 +    var canvas = document.createElement('canvas');
 +    canvas.width = size;
 +    canvas.height = size;
 +    var context = canvas.getContext('2d');
 +    context.fillStyle = "rgba(255, 255, 255, 0)";
 +    context.fill();
 +    context.font = fontsize + "px " + fontface;
 +    
 +    // text color
 +    context.fillStyle = "rgba(0, 0, 0, 1.0)";
 +     for(var i = 0; i<lines.length; i++){
 +        context.fillText(lines[i], size/2, size/2+i*fontsize);
 +     }
 +    
 +    // canvas contents will be used for a texture
 +    var texture = new THREE.Texture(canvas);
 +    texture.needsUpdate = true;
 +    
 +    var spriteMaterial = new THREE.SpriteMaterial({map: texture, depthTest: true, depthWrite: false, polygonOffset: true, polygonOffsetFactor: -1, polygonOffsetUnits: 1 });
 +    return spriteMaterial;
 +}
 +
 +
 +// ---------------------- INITIALIZING OBJECTS--------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +function init_object(obj) {
 +    if (obj.userData.hasOwnProperty("pointmaterial")) {
 +        init_points(obj);
 +        modelContains.points = true;
 +    }
 +    if (obj.userData.hasOwnProperty("pointlabels")) {
 +        init_pointlabels(obj);
 +        modelContains.pointlabels = true;
 +    }
 +    if (obj.userData.hasOwnProperty("edgematerial")) {
 +        init_lines(obj);
 +        modelContains.lines = true;
 +    }
 +    if (obj.userData.hasOwnProperty("edgelabels")) {
 +        init_edgelabels(obj);
 +        modelContains.edgelabels = true;
 +    }
 +    if (obj.userData.hasOwnProperty("arrowstyle")) {
 +        init_arrowheads(obj);
 +        modelContains.arrowheads = true;
 +    }
 +    if (obj.userData.hasOwnProperty("facetmaterial")) {
 +        init_faces(obj);
 +        modelContains.faces = true;
 +    }
 +}
 +
 +function init_points(obj) {
 +    var pointgroup = new THREE.Group();
 +    pointgroup.name = "points";
 +    var points = obj.userData.points;
 +    var radii = obj.userData.pointradii;
 +    var materials = obj.userData.pointmaterial;
 +    var geometry,material;
 +    if (!Array.isArray(radii)) {
 +        geometry = new THREE.SphereBufferGeometry(radii);  
 +    }
 +    if (!Array.isArray(materials)) {
 +        material = materials;
 +    }
 +    for (var i=0; i<points.length; i++) {
 +        var point = points[i];
 +        if (Array.isArray(radii)) {
 +            if (radii[i] == 0) {
 +                continue;
 +            }
 +            geometry = new THREE.SphereBufferGeometry(radii[i]);  
 +        } 
 +        if (Array.isArray(materials)) {
 +            material = materials[i];     
 +        } 
 +        var sphere = new THREE.Mesh(geometry, material);
 +        point.addSphere(sphere);
 +        pointgroup.add(sphere);
 +    }
 +    obj.add(pointgroup);
 +}
 +
 +function init_pointlabels(obj) {
 +    var points = obj.userData.points;
 +    var labels = obj.userData.pointlabels;
 +    var pointlabels = new THREE.Group();
 +    pointlabels.name = "pointlabels";
 +    if (Array.isArray(labels)) {
 +        for (var i=0; i<points.length; i++) {
 +            var point = points[i];
 +            var spriteMaterial = textSpriteMaterial( labels[i] );
 +         var sprite = new THREE.Sprite(spriteMaterial);
 +            point.addLabel(sprite);
 +            pointlabels.add(sprite);
 +        }
 +    } else {
 +        var spriteMaterial = textSpriteMaterial( labels );
 +        for (var i=0; i<points.length; i++) {
 +            var point = points[i];
 +         var sprite = new THREE.Sprite(spriteMaterial);
 +            point.addLabel(sprite);
 +            pointlabels.add(sprite);
 +        }
 +    }
 +    obj.add(pointlabels);
 +}
 +
 +function init_lines(obj) {
 +    var edgeindices = obj.userData.edgeindices;
 +    var points = obj.userData.points;
 +    var materials = obj.userData.edgematerial;
 +    var geometry = new THREE.BufferGeometry();
 +    var bufarr = new Float32Array( obj.userData.edgeindices.length * 3 );
 +    var bufattr = new THREE.Float32BufferAttribute( bufarr, 3 );
 +    var geometry = new THREE.BufferGeometry();
 +    geometry.setAttribute('position', bufattr);
 +    if (Array.isArray(materials)) {     
 +        for (var i=0; i<materials.length; i++) {
 +            geometry.addGroup(2*i,2,i);
 +        }
 +    }
 +    var lines = new THREE.LineSegments(geometry, materials);
 +    lines.name = "lines";
 +    obj.add(lines);
 +    updateEdgesPosition(obj);
 +}
 +
 +function init_edgelabels(obj) {
 +    var points = obj.userData.points;
 +    var edgeindices = obj.userData.edgeindices;
 +    var labels = obj.userData.edgelabels;
 +    var edgelabels = new THREE.Group();
 +    edgelabels.name = "edgelabels";
 +    if (Array.isArray(labels)) {
 +        for (var i=0; i<edgeindices.length/2; i++) {
 +            var spriteMaterial = textSpriteMaterial( labels[i] );
 +            var sprite = new THREE.Sprite(spriteMaterial);
 +            sprite.position.copy(new THREE.Vector3().addVectors(points[edgeindices[2*i]].vector,points[edgeindices[2*i+1]].vector).multiplyScalar(0.5));
 +            edgelabels.add(sprite);
 +        }
 +    } else {
 +        var spriteMaterial = textSpriteMaterial( labels );
 +        for (var i=0; i<edgeindices.length/2; i++) {
 +            var sprite = new THREE.Sprite(spriteMaterial);
 +            sprite.position.copy(new THREE.Vector3().addVectors(points[edgeindices[2*i]].vector,points[edgeindices[2*i+1]].vector).multiplyScalar(0.5));
 +            edgelabels.add(sprite);
 +        }
 +    }
 +    obj.add(edgelabels);
 +}
 +
 +function init_arrowheads(obj) {
 +    var arrowheads = new THREE.Group();
 +    arrowheads.name = "arrowheads";
 +    var arrowstyle = obj.userData.arrowstyle;
 +    var edgeindices = obj.userData.edgeindices;
 +    var edgematerials = obj.userData.edgematerial;
 +    var points = obj.userData.points;
 +    var material;
 +    if (!Array.isArray(edgematerials)) {
 +        material = new THREE.MeshBasicMaterial( {color: edgematerials.color} );
 +    }
 +
 +    for (var i=0; i<edgeindices.length; i=i+2) {
 +        var start = points[edgeindices[i]];
 +        var end = points[edgeindices[i+1]];
 +        var dist = start.vector.distanceTo( end.vector ) - start.radius() - end.radius();
 +        if (dist <= 0) {
 +            continue;
 +        }
 +        var dir = new THREE.Vector3().subVectors(end.vector,start.vector);
 +        dir.normalize();
 +        var axis = new THREE.Vector3().set(dir.z,0,-dir.x);
 +        axis.normalize();
 +        var radians = Math.acos( dir.y );
 +        var radius = dist/25;
 +        var height = dist/5;
 +        var geometry = new THREE.ConeBufferGeometry(radius,height);
 +        var position = new THREE.Vector3().addVectors(start.vector,dir.clone().multiplyScalar(start.radius()+dist-height/2));
 +        if (Array.isArray(edgematerials)) {
 +            material = new THREE.MeshBasicMaterial( {color: edgematerials[i].color} );
 +        }
 +        var cone = new THREE.Mesh( geometry, material );
 +        cone.quaternion.setFromAxisAngle(axis,radians);;
 +        cone.position.copy(position);;
 +        arrowheads.add(cone);
 +    }
 +    obj.add(arrowheads);
 +}
 +
 +function init_faces(obj) {
 +    var points = obj.userData.points;
 +    var facets = obj.userData.facets;
 +    obj.userData.triangleindices = [];
 +    for (var i=0; i<facets.length; i++) {
 +        facet = facets[i];
 +        for (var t=0; t<facet.length-2; t++) {
 +            obj.userData.triangleindices.push(facet[0],facet[t+1],facet[t+2]);  
 +        }
 +    }
 +    var bufarr = new Float32Array( obj.userData.triangleindices.length * 3 );
 +    var bufattr = new THREE.Float32BufferAttribute(bufarr,3);
 +    
 +    var materials = obj.userData.facetmaterial;
 +    var geometry = new THREE.BufferGeometry();
 +    var frontmaterials = [];
 +    var backmaterials = [];
 +    geometry.setAttribute('position',bufattr);
 +    if (Array.isArray(materials)) {
 +        var tricount = 0;
 +        var facet;
 +        for (var i=0; i<facets.length; i++) {
 +            facet = facets[i];
 +            geometry.addGroup(tricount,(facet.length-2)*3,i);
 +            tricount += (facet.length-2)*3;
 +        }
 +        for (var j=0; j<materials.length; j++) {
 +            var fmat = materials[j].clone()
 +            fmat.side = THREE.FrontSide;
 +            frontmaterials.push(fmat);
 +            var bmat = materials[j].clone()
 +            bmat.side = THREE.BackSide;
 +            backmaterials.push(bmat);
 +            obj.userData.facetmaterial = frontmaterials.concat(backmaterials);
 +        }
 +    } else if (materials instanceof THREE.Material) {
 +        frontmaterials = materials.clone()
 +        frontmaterials.side = THREE.FrontSide;
 +        backmaterials = materials.clone()
 +        backmaterials.side = THREE.BackSide;
 +        obj.userData.facetmaterial = [frontmaterials, backmaterials];
 +    }
 +    // duplicating the object with front and back should avoid transparency issues
 +    var backmesh = new THREE.Mesh(geometry, backmaterials);
 +    // meshname is used to show/hide objects
 +    backmesh.name = "backfaces";
 +    obj.add(backmesh);
 +    var frontmesh = new THREE.Mesh(geometry, frontmaterials);
 +    frontmesh.name = "frontfaces";
 +    obj.add(frontmesh);
 +    updateFacesPosition(obj);
 +}
 +// //INITIALIZING
 +
 +
 +function updateFacesPosition(obj) {
 +    var points = obj.userData.points;
 +    var indices = obj.userData.triangleindices;
 +    var faces = obj.getObjectByName("frontfaces");
 +    var ba = faces.geometry.getAttribute("position");
 +    for (var i=0; i<indices.length; i++) {
 +        ba.setXYZ(i, points[indices[i]].vector.x, points[indices[i]].vector.y ,points[indices[i]].vector.z); 
 +    }
 +    faces.geometry.attributes.position.needsUpdate = true;
 +    
 +}
 +
 +function updateEdgesPosition(obj) {
 +    var points = obj.userData.points;
 +    var indices = obj.userData.edgeindices;
 +    var lines = obj.getObjectByName("lines");
 +    var ba = lines.geometry.getAttribute("position"); 
 +    for (var i=0; i<indices.length; i++) {
 +        ba.setXYZ(i, points[indices[i]].vector.x, points[indices[i]].vector.y ,points[indices[i]].vector.z); 
 +    }
 +    lines.geometry.attributes.position.needsUpdate = true;
 +}
 +
 +function onWindowResize() {
 +    renderer.setSize( three.clientWidth, three.clientHeight );
 +    svgRenderer.setSize( three.clientWidth, three.clientHeight );
 +    updateCamera();
 +}
 +
 +function updateCamera() {
 +    var width = three.clientWidth;
 +    var height = three.clientHeight;
 +    var aspect = width / height;
 +    if (camera.type == "OrthographicCamera") {
 +        camera.left = frustumSize * aspect / - 2;
 +        camera.right = frustumSize * aspect / 2;
 +        camera.top = frustumSize / 2;
 +        camera.bottom = - frustumSize / 2;
 +    } else if (camera.type == "PerspectiveCamera") {
 +        camera.aspect = aspect;
 +    }
 +    camera.updateProjectionMatrix();
 +}
 +
 +function changeCamera(event) {
 +    var selindex = event.currentTarget.selectedIndex;
 +    camera = cameras[selindex];
 +    control = controls[selindex];
 +    control.enabled = true; 
 +    for (var i=0; i<controls.length; i++) {
 +        if (i!=selindex) {
 +            controls[i].enabled = false;
 +        }
 +    }
 +    updateCamera();
 +}
 +
 +var camtypenode = document.getElementById('cameraType_1');
 +camtypenode.onchange = changeCamera;
 +camtypenode.dispatchEvent(new Event('change'));
 +
 +onWindowResize();
 +window.addEventListener('resize', onWindowResize);
 +
 +
 +var xRotationEnabled = false;
 +var yRotationEnabled = false;
 +var zRotationEnabled = false;
 +var rotationSpeedFactor = 1;
 +var settingsShown = false;
 +var labelsShown = true;
 +var intervals = [];
 +var timeouts = [];
 +var explodingSpeed = 0.05;
 +var explodeScale = 0.000001;
 +var XMLS = new XMLSerializer();
 +var svgElement;
 +var renderId;
 +
 +var render = function () {
 +
 + renderId = requestAnimationFrame(render);
 +
 +// comment in for automatic explosion
 +// explode(updateFactor());
 +
 +    var phi = 0.02 * rotationSpeedFactor;
 +
 +    if (xRotationEnabled) {
 +        scene.rotation.x += phi;
 +    }
 +    if (yRotationEnabled) {
 +        scene.rotation.y += phi;
 +    }
 +    if (zRotationEnabled) {
 +        scene.rotation.z += phi;
 +    }
 +
 +    control.update();
 +    renderer.render(scene, camera);
 +};
 +
 +if ( THREE.WEBGL.isWebGLAvailable() ) {
 + render();
 +} else {
 + var warning = WEBGL.getWebGLErrorMessage();
 + three.appendChild( warning );
 +}
 +    
 +function changeTransparency() {
 +    var opacity = 1-Number(event.currentTarget.value);
 +    for (var i=0; i<scene.children.length; i++) {
 +        child = scene.children[i];
 +        if ( child.userData.hasOwnProperty("facetmaterial") ) {
 +            if (Array.isArray(child.userData.facetmaterial)) {
 +                for (var j=0; j<child.userData.facetmaterial.length; j++) {
 +                    child.userData.facetmaterial[j].opacity = opacity;
 +                }
 +            } else {
 +                child.userData.facetmaterial.opacity = opacity;
 +            }    
 +        }
 +    }
 +}
 +
 +function toggleDepthWrite(event) {
 +    depthwrite = event.currentTarget.checked;
 +    for (var i=0; i<scene.children.length; i++) {
 +        child = scene.children[i];
 +        if ( child.userData.hasOwnProperty("facetmaterial") ) {
 +            if (Array.isArray(child.userData.facetmaterial)) {
 +                for (var j=0; j<child.userData.facetmaterial.length; j++) {
 +                    child.userData.facetmaterial[j].depthWrite = depthwrite;
 +                }
 +            } else {
 +                child.userData.facetmaterial.depthWrite = depthWrite;
 +            }    
 +        }
 +    }
 +}
 +
 +function changeRotationX(event){
 +    xRotationEnabled = event.currentTarget.checked;
 +}
 +
 +function changeRotationY(event){
 +    yRotationEnabled = event.currentTarget.checked;
 +}
 +
 +function changeRotationZ(event){
 +    zRotationEnabled = event.currentTarget.checked;
 +}
 +
 +
 +function changeRotationSpeedFactor(event){
 +    rotationSpeedFactor = Number(event.currentTarget.value);
 +}
 +
 +function resetScene(){
 +    scene.rotation.set(0,0,0);
 +    camera.position.set(0,0,5);
 +    camera.up.set(0,1,0);
 +}
 +
 +function showSettings(event){
 +    document.getElementById('settings_1').style.visibility = 'visible';
 +    document.getElementById('showSettingsButton_1').style.visibility = 'hidden';
 +    document.getElementById('hideSettingsButton_1').style.visibility = 'visible';
 +    settingsShown = true;
 +}
 +
 +function hideSettings(event){
 +    document.getElementById('settings_1').style.visibility = 'hidden';
 +    document.getElementById('showSettingsButton_1').style.visibility = 'visible';
 +    document.getElementById('hideSettingsButton_1').style.visibility = 'hidden';
 +    settingsShown = false;
 +}
 +
 +
 +
 +var pos = 150* Math.PI;
 +
 +function updateFactor() {
 +    pos++;
 +    return Math.sin(.01*pos)+1;
 +}
 +
 +// ------------------------ FOLDING ------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +// rotate point p around axis defined by points p1 and p2 by given angle
 +function rotate(p, p1, p2, angle ){   
 +    angle = -angle;
 +    var x = p.x, y = p.y, z = p.z, 
 +    a = p1.x, b = p1.y, c = p1.z, 
 +    u = p2.x-p1.x, v = p2.y-p1.y, w = p2.z-p1.z;
 +    var result = [];
 +    var L = u*u + v*v + w*w;
 +    var sqrt = Math.sqrt;
 +    var cos = Math.cos;
 +    var sin = Math.sin;
 +
 +    result[0] = ((a*(v*v+w*w)-u*(b*v+c*w-u*x-v*y-w*z))*(1-cos(angle))+L*x*cos(angle)+sqrt(L)*(-c*v+b*w-w*y+v*z)*sin(angle))/L;
 +    result[1] = ((b*(u*u+w*w)-v*(a*u+c*w-u*x-v*y-w*z))*(1-cos(angle))+L*y*cos(angle)+sqrt(L)*(c*u-a*w+w*x-u*z)*sin(angle))/L;
 +    result[2] = ((c*(u*u+v*v)-w*(a*u+b*v-u*x-v*y-w*z))*(1-cos(angle))+L*z*cos(angle)+sqrt(L)*(-b*u+a*v-v*x+u*y)*sin(angle))/L;
 +
 +    return result;
 +}
 +
 +var fold = function(event){
 +    var obj = foldables[Number(event.currentTarget.name)];
 +    var foldvalue = Number(event.currentTarget.value);
 +    var scale = foldvalue - obj.userData.oldscale;
 +
 +    for (var j=0; j<obj.userData.axes.length; j++) {
 +        rotateVertices(obj, j, scale);
 +    }
 +    update(obj);
 +    obj.userData.oldscale += scale;
 +    lookAtBarycenter(obj);
 +}
 +
 +function lookAtBarycenter(obj){
 +    control.target = barycenter(obj);
 +}
 +
 +function barycenter(obj) {
 +    var center = new THREE.Vector3(0,0,0);
 +    var points = obj.userData.points;
 +    for (var i=0; i<points.length; i++){
 +        center.add(points[i].vector);
 +    }
 +    center.divideScalar(points.length);
 +    return center;
 +}
 +
 +function rotateVertices(obj, edge, scale) {
 +    var axes = obj.userData.axes;
 +    var subtrees = obj.userData.subtrees;
 +    var points = obj.userData.points;
 +    var angles = obj.userData.angles;
 +    if (edge < axes.length){
 +        for (var j=0; j<subtrees[edge].length; j++){
 +            var rotP = rotate(points[subtrees[edge][j]].vector, points[axes[edge][0]].vector,points[axes[edge][1]].vector, scale * (Math.PI - angles[edge]));
 +            points[subtrees[edge][j]].set(rotP[0],rotP[1],rotP[2]);
 +        }
 +    }
 +}
 +
 +function update(obj) {
 +   updateFacesPosition(obj);
 +   updateEdgesPosition(obj);
 +}
 +
 +if (foldables.length) {
 +    var settings = document.getElementById('settings_1');
 +    var foldDiv = document.createElement('div');
 +    foldDiv.id = 'fold_1';
 +    var title = document.createElement('strong');
 +    title.innerHTML = 'Fold';
 +    foldDiv.appendChild(title);
 +    foldDiv.className = 'group';
 +    for (var i=0; i<foldables.length; i++) {
 +        var range = document.createElement('input');
 +        range.type = 'range';
 +        range.min = 0;
 +        range.max = 1;
 +        range.value = 0;
 +        range.step = 0.001;
 +        range.name = String(i);
 +        range.oninput = fold;
 +        foldDiv.appendChild(range);
 +    }
 +    lookAtBarycenter(foldables[0]);
 +    settings.insertBefore(foldDiv,settings.childNodes[0]);
 +}
 +
 +    
 +// ---------------------- EXPLOSION ------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +if (explodableModel) {
 +    for (var i=0; i<scene.children.length; i++) {
 +        obj = scene.children[i];
 +        if ( obj.userData.explodable ) {
 +            computeCentroid(obj);
 +        }
 +    }
 +    document.getElementById('explodeRange_1').oninput = triggerExplode;
 +    document.getElementById('explodeCheckbox_1').onchange = triggerAutomaticExplode;
 +    document.getElementById('explodingSpeedRange_1').oninput = setExplodingSpeed;
 +    explode(0.000001);
 +}
 +
 +function computeCentroid(obj) {
 +    centroid = new THREE.Vector3();
 +    obj.userData.points.forEach(function(pmpoint) {
 +        centroid.add(pmpoint.vector);
 +    });
 +    centroid.divideScalar(obj.userData.points.length);
 +    obj.userData.centroid = centroid;
 +}
 +
 +function explode(factor) {
 +    for (var i=0; i<scene.children.length; i++) {
 +        var obj = scene.children[i];
 +        if (obj.userData.hasOwnProperty("centroid")) { 
 +            var c = obj.userData.centroid;
 +            obj.position.set(c.x*factor, c.y*factor, c.z*factor);
 +        }
 +    }
 +}
 +
 +function triggerExplode(event){
 +    explodeScale = Number(event.currentTarget.value);
 +    explode(explodeScale);
 +}
 +
 +function setExplodingSpeed(event){
 +    explodingSpeed = Number(event.currentTarget.value);
 +}
 +
 +function triggerAutomaticExplode(event){
 +    if (event.currentTarget.checked){
 +        startExploding();
 +    } else {
 +        clearIntervals();
 +    }
 +}
 +
 +function startExploding(){
 +    intervals.push(setInterval(explodingInterval, 25));
 +}
 +
 +
 +function explodingInterval(){
 +    explodeScale += explodingSpeed;
 +    if (explodeScale <= 6){ 
 +        explode(explodeScale);
 +    }
 +    else{
 +        explode(6);
 +        explodeScale = 6;
 +        clearIntervals();
 +        timeouts.push(setTimeout(startUnexploding, 3000));
 +    }
 +    document.getElementById('explodeRange_1').value = explodeScale;
 +}
 +
 +
 +function startUnexploding(){
 +    intervals.push(setInterval(unexplodingInterval, 25));
 +}
 +
 +function unexplodingInterval(){
 +    explodeScale -= explodingSpeed;
 +    if (explodeScale >= 0){
 +        explode(explodeScale);
 +    }
 +    else {
 +        explode(0);
 +        explodeScale = 0;
 +        clearIntervals();
 +        timeouts.push(setTimeout(startExploding, 3000));
 +    }
 +    document.getElementById('explodeRange_1').value = explodeScale;
 +}
 +
 +function clearIntervals(){
 +    intervals.forEach(function(interval){
 +        clearInterval(interval);
 +    });
 +    intervals = [];
 +    timeouts.forEach(function(timeout){
 +        clearTimeout(timeout);
 +    });
 +    timeouts = [];
 +}
 +
 +// ---------------------- DISPLAY --------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +const objectTypeInnerHTMLs = { points: "Points", pointlabels: "Point labels", lines: "Edges", edgelabels: "Edge labels", faces: "Faces", arrowheads: "Arrow heads" };
 +const objectTypeVisible = {};
 +Object.assign(objectTypeVisible,modelContains);
 +const sortedObjectTypeKeys = Object.keys(objectTypeInnerHTMLs).sort();
 +const shownObjectTypesList = document.getElementById('shownObjectTypesList_1');
 +
 +function setVisibility(bool,objname) {
 +    for (var i=0; i<scene.children.length; i++){
 +        var obj = scene.children[i].getObjectByName(objname);
 +        if (obj) {
 +            obj.visible = bool;
 +        }
 +    }
 +}
 +
 +function toggleObjectTypeVisibility(event){
 +    var name = event.currentTarget.name;
 +    var checked = event.currentTarget.checked;
 +    objectTypeVisible[name] = checked;
 +    if (name == "faces") {
 +        setVisibility(checked,"frontfaces");
 +        setVisibility(checked,"backfaces");
 +    } else {
 +        setVisibility(checked,name);
 +    }
 +}
 +
 +for (var i=0; i<sortedObjectTypeKeys.length; i++){
 +    var key = sortedObjectTypeKeys[i];
 +    if (modelContains[key]) {
 +        var objTypeNode = document.createElement('span');
 +        objTypeNode.innerHTML = objectTypeInnerHTMLs[key] + '<br>';
 +        var checkbox = document.createElement('input');
 +        checkbox.type = 'checkbox';
 +        checkbox.checked = true;
 +        checkbox.name = key;
 +        checkbox.onchange = toggleObjectTypeVisibility;
 +        shownObjectTypesList.appendChild(checkbox);
 +        shownObjectTypesList.appendChild(objTypeNode);
 +    }
 +}
 +
 +// ------------------------------------------------------
 +
 +function toggleObjectVisibility(event){
 +    var nr = Number(event.currentTarget.name);
 +    scene.children[nr].visible = event.currentTarget.checked;
 +}
 +
 +// append checkboxes for displaying or hiding objects
 +var shownObjectsList = document.getElementById('shownObjectsList_1');
 +for (var i=0; i<scene.children.length; i++){
 +    obj = scene.children[i];
 +    var objNode = document.createElement('span');
 +    objNode.innerHTML = obj.name + '<br>';
 +    var checkbox = document.createElement('input');
 +    checkbox.type = 'checkbox';
 +    checkbox.checked = true;
 +    checkbox.name = String(i);
 +    checkbox.onchange = toggleObjectVisibility;
 +    shownObjectsList.appendChild(checkbox);
 +    shownObjectsList.appendChild(objNode);
 +}
 +
 +// ---------------------- SVG ------------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +function takeSvgScreenshot() {
 +    if (objectTypeVisible["pointlabels"]) {
 +        setVisibility(false,"pointlabels");
 +    }
 +    if (objectTypeVisible["edgelabels"]) {
 +        setVisibility(false,"edgelabels");
 +    }
 +    svgRenderer.render(scene,camera);
 +    svgElement = XMLS.serializeToString(svgRenderer.domElement);
 +    
 +    if (objectTypeVisible["pointlabels"]) {
 +        setVisibility(true,"pointlabels");
 +    }
 +    if (objectTypeVisible["edgelabels"]) {
 +        setVisibility(true,"edgelabels");
 +    }
 +
 +    if (document.getElementById('tab_1').checked){
 +        //show in new tab
 +        var myWindow = window.open("","");
 +        myWindow.document.body.innerHTML = svgElement;
 +    } else{
 +        // download svg file 
 +        download("screenshot.svg", svgElement);
 +    }
 +}
 +
 +function download(filename, text) {
 +    var element = document.createElement('a');
 +    element.setAttribute('href', 'data:text/plain;charset=utf-8,' + encodeURIComponent(text));
 +    element.setAttribute('download', filename);
 +
 +    element.style.display = 'none';
 +    document.body.appendChild(element);
 +
 +    element.click();
 +
 +    document.body.removeChild(element);
 +}
 +
 +
 +document.getElementById('transparencyRange_1').oninput = changeTransparency;
 +document.getElementById('depthWriteCheckbox_1').onchange = toggleDepthWrite;
 +document.getElementById('changeRotationX_1').onchange = changeRotationX;
 +document.getElementById('changeRotationY_1').onchange = changeRotationY;
 +document.getElementById('changeRotationZ_1').onchange = changeRotationZ;
 +document.getElementById('resetButton_1').onclick = resetScene;
 +document.getElementById('rotationSpeedRange_1').oninput = changeRotationSpeedFactor;
 +document.getElementById('takeScreenshot_1').onclick = takeSvgScreenshot;
 +document.getElementById('showSettingsButton_1').onclick = showSettings;
 +document.getElementById('hideSettingsButton_1').onclick = hideSettings;
 +
 +
 +// ------------------ SHORTCUTS --------------------------------------------
 +// -------------------------------------------------------------------------
 +
 +/**
 + * http://www.openjs.com/scripts/events/keyboard_shortcuts/
 + * Version : 2.01.B
 + * By Binny V A
 + * License : BSD
 + */
 +shortcut = {
 + 'all_shortcuts':{},//All the shortcuts are stored in this array
 + 'add': function(shortcut_combination,callback,opt) {
 + //Provide a set of default options
 + var default_options = {
 + 'type':'keydown',
 + 'propagate':false,
 + 'disable_in_input':false,
 + 'target':document,
 + 'keycode':false
 + }
 + if(!opt) opt = default_options;
 + else {
 + for(var dfo in default_options) {
 + if(typeof opt[dfo] == 'undefined') opt[dfo] = default_options[dfo];
 + }
 + }
 +
 + var ele = opt.target;
 + if(typeof opt.target == 'string') ele = document.getElementById(opt.target);
 + var ths = this;
 + shortcut_combination = shortcut_combination.toLowerCase();
 +
 + //The function to be called at keypress
 + var func = function(e) {
 + e = e || window.event;
 +
 + if(opt['disable_in_input']) { //Don't enable shortcut keys in Input, Textarea fields
 + var element;
 + if(e.target) element=e.target;
 + else if(e.srcElement) element=e.srcElement;
 + if(element.nodeType==3) element=element.parentNode;
 +
 + if(element.tagName == 'INPUT' || element.tagName == 'TEXTAREA') return;
 + }
 +
 + //Find Which key is pressed
 + if (e.keyCode) code = e.keyCode;
 + else if (e.which) code = e.which;
 + var character = String.fromCharCode(code).toLowerCase();
 +
 + if(code == 188) character=","; //If the user presses , when the type is onkeydown
 + if(code == 190) character="."; //If the user presses , when the type is onkeydown
 +
 + var keys = shortcut_combination.split("+");
 + //Key Pressed - counts the number of valid keypresses - if it is same as the number of keys, the shortcut function is invoked
 + var kp = 0;
 +
 + //Work around for stupid Shift key bug created by using lowercase - as a result the shift+num combination was broken
 + var shift_nums = {
 + "`":"~",
 + "1":"!",
 + "2":"@",
 + "3":"#",
 + "4":"$",
 + "5":"%",
 + "6":"^",
 + "7":"&",
 + "8":"*",
 + "9":"(",
 + "0":")",
 + "-":"_",
 + "=":"+",
 + ";":":",
 + "'":"\"",
 + ",":"<",
 + ".":">",
 + "/":"?",
 + "\\":"|"
 + }
 + //Special Keys - and their codes
 + var special_keys = {
 + 'esc':27,
 + 'escape':27,
 + 'tab':9,
 + 'space':32,
 + 'return':13,
 + 'enter':13,
 + 'backspace':8,
 +
 + 'scrolllock':145,
 + 'scroll_lock':145,
 + 'scroll':145,
 + 'capslock':20,
 + 'caps_lock':20,
 + 'caps':20,
 + 'numlock':144,
 + 'num_lock':144,
 + 'num':144,
 +
 + 'pause':19,
 + 'break':19,
 +
 + 'insert':45,
 + 'home':36,
 + 'delete':46,
 + 'end':35,
 +
 + 'pageup':33,
 + 'page_up':33,
 + 'pu':33,
 +
 + 'pagedown':34,
 + 'page_down':34,
 + 'pd':34,
 +
 + 'left':37,
 + 'up':38,
 + 'right':39,
 + 'down':40,
 +
 + 'f1':112,
 + 'f2':113,
 + 'f3':114,
 + 'f4':115,
 + 'f5':116,
 + 'f6':117,
 + 'f7':118,
 + 'f8':119,
 + 'f9':120,
 + 'f10':121,
 + 'f11':122,
 + 'f12':123
 + }
 +
 + var modifiers = { 
 + shift: { wanted:false, pressed:false},
 + ctrl : { wanted:false, pressed:false},
 + alt  : { wanted:false, pressed:false},
 + meta : { wanted:false, pressed:false} //Meta is Mac specific
 + };
 +                        
 + if(e.ctrlKey) modifiers.ctrl.pressed = true;
 + if(e.shiftKey) modifiers.shift.pressed = true;
 + if(e.altKey) modifiers.alt.pressed = true;
 + if(e.metaKey)   modifiers.meta.pressed = true;
 +                        
 + for(var i=0; k=keys[i],i<keys.length; i++) {
 + //Modifiers
 + if(k == 'ctrl' || k == 'control') {
 + kp++;
 + modifiers.ctrl.wanted = true;
 +
 + } else if(k == 'shift') {
 + kp++;
 + modifiers.shift.wanted = true;
 +
 + } else if(k == 'alt') {
 + kp++;
 + modifiers.alt.wanted = true;
 + } else if(k == 'meta') {
 + kp++;
 + modifiers.meta.wanted = true;
 + } else if(k.length > 1) { //If it is a special key
 + if(special_keys[k] == code) kp++;
 +
 + } else if(opt['keycode']) {
 + if(opt['keycode'] == code) kp++;
 +
 + } else { //The special keys did not match
 + if(character == k) kp++;
 + else {
 + if(shift_nums[character] && e.shiftKey) { //Stupid Shift key bug created by using lowercase
 + character = shift_nums[character]; 
 + if(character == k) kp++;
 + }
 + }
 + }
 + }
 +
 + if(kp == keys.length && 
 + modifiers.ctrl.pressed == modifiers.ctrl.wanted &&
 + modifiers.shift.pressed == modifiers.shift.wanted &&
 + modifiers.alt.pressed == modifiers.alt.wanted &&
 + modifiers.meta.pressed == modifiers.meta.wanted) {
 + callback(e);
 +
 + if(!opt['propagate']) { //Stop the event
 + //e.cancelBubble is supported by IE - this will kill the bubbling process.
 + e.cancelBubble = true;
 + e.returnValue = false;
 +
 + //e.stopPropagation works in Firefox.
 + if (e.stopPropagation) {
 + e.stopPropagation();
 + e.preventDefault();
 + }
 + return false;
 + }
 + }
 + }
 + this.all_shortcuts[shortcut_combination] = {
 + 'callback':func, 
 + 'target':ele, 
 + 'event': opt['type']
 + };
 + //Attach the function with the event
 + if(ele.addEventListener) ele.addEventListener(opt['type'], func, false);
 + else if(ele.attachEvent) ele.attachEvent('on'+opt['type'], func);
 + else ele['on'+opt['type']] = func;
 + },
 +
 + //Remove the shortcut - just specify the shortcut and I will remove the binding
 + 'remove':function(shortcut_combination) {
 + shortcut_combination = shortcut_combination.toLowerCase();
 + var binding = this.all_shortcuts[shortcut_combination];
 + delete(this.all_shortcuts[shortcut_combination])
 + if(!binding) return;
 + var type = binding['event'];
 + var ele = binding['target'];
 + var callback = binding['callback'];
 +
 + if(ele.detachEvent) ele.detachEvent('on'+type, callback);
 + else if(ele.removeEventListener) ele.removeEventListener(type, callback, false);
 + else ele['on'+type] = false;
 + }
 +}
 +
 +shortcut.add("Alt+Left",function() {
 + var event = new Event('click');
 + if (settingsShown){
 + document.getElementById('hideSettingsButton_1').dispatchEvent(event);
 + } else {
 + document.getElementById('showSettingsButton_1').dispatchEvent(event);
 + }
 +});
 +
 +
 +// COMMON_CODE_BLOCK_END
 +
 +});});
 +      </script>
 +   </body>
 +</html>
 +</HTML>
 +===== Computing Facets =====
 +
 +An important action that is often needed to come up with new facet describing inequalities for combinatorial optimization problems is the computation of convex hulls for small examples.
 +
 +==== Pure Integer Case ====
 +
 +We begin with the case in which all variables are required to be integral, i.e., the //pure integer case//. Moreover, the approach depends on whether the polyhedron is bounded or not.
 +
 +=== Bounded Polyhedra ===
 +
 +Let us illustrate the approach via the example of the //stable set problem//: Here one is given an (undirected) Graph G = (V,E) with node set V and edges E. The goal is to find a largest subset of nodes V' such that any two nodes in V' are not connected by an edge.
 +
 +For our example consider the 5-cycle, i.e., the graph C<sub>5</sub> with five nodes {1, 2, 3, 4, 5} and edges {1,2}, {2,3}, {3,4}, {4,5}, {5,1}. A formulation of the stable set problem for this graph looks as follows:
 +
 +<code>
 +Maximize
 + obj: x#1 + x#2 + x#3 + x#4 + x#5
 +Subject to
 + edge_1: x#2 + x#1 <= 1
 + edge_2: x#3 + x#2 <= 1
 + edge_3: x#4 + x#3 <= 1
 + edge_4: x#5 + x#4 <= 1
 + edge_5: x#1 + x#5 <= 1
 +Bounds
 + 0 <= x#1 <= 1
 + 0 <= x#2 <= 1
 + 0 <= x#3 <= 1
 + 0 <= x#4 <= 1
 + 0 <= x#5 <= 1
 +General
 + x#1 x#2 x#3 x#4 x#5
 +End
 +</code>
 +Here, ''%%General%%'' encodes that the following variables should be restricted to obtain integer values. Thus, all variables ''%%x#1, x#2, x#3, x#4, x#5%%'' are restricted to binary values (0 or 1). The value 1 encodes that the corresponding node will be in an optimal stable set (and 0 otherwise). The constraints ''%%edge_*%%'' encode that a most one of the two nodes covered by an edge is selected. Clearly, the corresponding polyhedron is bounded.
 +
 +We assume that the above information is contained in the file ''%%stab.lp%%''. We now read it into polymake and convert it to rational form, as explained above:
 +
 +<code perl>
 +> $f=lp2poly('files/optimization/stab.lp');
 +> $p = new Polytope<Rational>($f);
 +</code>
 +We are now interested in all feasible solutions to the above problem, i.e., all assignments of 0 or 1 to the variables such that the above inequalities are satisfied. These feasible points can be computed via:
 +
 +<code perl>
 +> $p->LATTICE_POINTS;
 +</code>
 +To understand these points and make computational use of this information, we are interested in the convex hull of all feasible solutions. To this end, we construct a new polytope which is specified via points for which the convex hull is taken. This can be done as follows:
 +
 +<code perl>
 +> $s=new Polytope(POINTS=>$p->LATTICE_POINTS, COORDINATE_LABELS=>$p->COORDINATE_LABELS);
 +</code>
 +Here, the coordinate labels, i.e., the variable names, are copied to the new polytope.
 +
 +Now, the facets of the new polytope can be computed and listed via:
 +
 +<code perl>
 +> print_constraints($s);
 +Facets:
 +0: x#3 >= 0
 +1: -x#4 - x#5 >= -1
 +2: x#1 >= 0
 +3: -x#2 - x#3 >= -1
 +4: -x#1 - x#2 - x#3 - x#4 - x#5 >= -2
 +5: -x#1 - x#2 >= -1
 +6: x#4 >= 0
 +7: -x#1 - x#5 >= -1
 +8: -x#3 - x#4 >= -1
 +9: x#2 >= 0
 +10: x#5 >= 0
 +
 +</code>
 +The facet defining inequalities can be interpreted as follows:
 +
 +  * There are five trivial inequalities ''%%x#? >= 0%%''.
 +  * The five original 'edge' inequalities ''%%x#i + x#j <= 1%%'' define facets.
 +  * We have the so-called //odd-cycle inequality// ''%%x#1 +x#2 + x#3 + x#4 + x#5 <= 2%%'', stating that at most two nodes in an (odd) cycle of length 5 can be selected. This inequality can be generalized by taking the sum of all variables in an odd cycle and restricting the sum to be less or equal to the size of the cycle minus 1 divided by 2.
 +
 +Of course, one can also use the usual polymake output, e.g., ''%%print $s->FACETS%%''.
 +
 +This example showed one of the routine actions often performed by discrete optimizers. Of course, this action can also be performed by a script, which makes the computation a one-line command.
 +
 +Note that the size of instances that can be handled will probably be small. Usually, things become difficult from dimension 15 on, but it depends on the particular structure of your instances, i.e., on the number of facets and lattice points.
 +
 +=== Unbounded Polyhedra ===
 +
 +If the underlying polyhedron is unbounded, the approach above does not work anymore, since there are infinitely many lattice points. Arguably, this case occurs less often than the bounded case, but it is a excellent show-case for polymake's potential.
 +
 +The following mathematical insights are important to treat the unbounded case. First, we have to assume that the data, i.e., the inequality description of the polyhedron ''%%P%%'', is rational; otherwise, we cannot expect a finite description of the convex hull. Second, we write ''%%P%%'' as the sum of a finite part ''%%Q%%'' and the recession cone ''%%C%%''. If the data is rational, the recession cone of ''%%P%%'' and of the integer hull coincide. Third, it suffices to generate the integer points in ''%%Q + R%%'', where ''%%R%%'' is the parallelotope generated by the rays of ''%%C%%''. Thus, ''%%R%%'' is generated by the Minkowski sum of the interval ''%%[0,1]%%'' and the generating rays.
 +
 +To illustrate the construction, consider the following example:
 +
 +<code>
 +Minimize
 + obj:  x1 + x2
 +Subject to
 +C1: x1 + x2 >= 0.5
 +C2: x1 - 2 x2 <= 1.5
 +C3: x2 - 2 x1 <= 1.5
 +General
 + x1 x2
 +End
 +</code>
 +<html><img src="/lib/exe/fetch.php?media=tutorials/release/4.11/optimization/ip-unb.gif" class="center-block"></html><html><img></html>
 +
 +We now assume that the example is contained in the file ''%%unbounded.lp%%'' and proceed as above
 +
 +<code perl>
 +> $f = lp2poly('files/optimization/unbounded.lp');
 +> $pin = new Polytope<Rational>($f);
 +</code>
 +The visualization in the picture can be generated with ''%%$pin->VISUAL%%''. The lattice points can be shown with ''%%$pin->VISUAL->LATTICE_COLORED%%''.
 +
 +We now extract the rays of the recession cone
 +
 +<code perl>
 +> $rays = $pin->VERTICES->minor($pin->FAR_FACE, All);
 +</code>
 +This command first computes all vertices of the polyhedron (this includes unbounded vertices); note that is involves a convex hull computation. The set ''%%FAR_FACE%%'' contains the indices of all vertices that are unbounded. The result is:
 +
 +<code perl>
 +> print $rays;
 +0 1 1/2
 +0 1 2
 +</code>
 +Thus, there are two rays that are generators of the recession cone.
 +
 +We now have to construct the Minkowski hull of all intervals ''%%[0,r]%%'' for each ray ''%%r%%'' (scaled to be integral). This can be done with the following code (possibly easier):
 +
 +<code perl>
 +> $zero = unit_vector<Rational>($pin->DIM + 1, 0);
 +> $B = new Polytope<Rational>(POINTS=>$zero);
 +> foreach my $r (@$rays)
 +> {
 +>     $M = new Matrix<Rational>(primitive($r));
 +>     $M->[0]->[0] = 1;
 +>     $M = $M / $zero;
 +>     $ptemp = new Polytope<Rational>(POINTS=>$M);
 +>     $B = minkowski_sum($B, $ptemp);
 +> }
 +</code>
 +The code first generates a polytope ''%%B%%'' consisting of 0 only. It then takes each ray ''%%r%%'' in turn and creates ''%%[0,r]%%''. It then takes the Minkowski sum of this new polytope with ''%%B%%'' and stores the result in ''%%B%%''.
 +
 +The next step is to obtain the bounded part ''%%Q%%'' of ''%%P%%'', by first extracting the bounded vertices and the creating a new polytope:
 +
 +<code perl>
 +> $Qpoints = $pin->VERTICES->minor($pin->BOUNDED_VERTICES, All);
 +> $Q = new Polytope<Rational>(POINTS=>$Qpoints);
 +</code>
 +The two polytopes are now combined:
 +
 +<code perl>
 +> $p = minkowski_sum($Q, $B);
 +</code>
 +We now generate the lattice points (as in the bounded part) and add the rays from above:
 +
 +<code perl>
 +> $latticemat = new Matrix<Rational>($p->LATTICE_POINTS);
 +> $newpoints = new Matrix<Rational>($latticemat / $rays);
 +</code>
 +Here, ''%%newpoints%%'' is a matrix that contains all lattice points in ''%%Q%%'' and the rays from above.
 +
 +Finally, the polytope we are interested in is:
 +
 +<code perl>
 +> $q = new Polytope(POINTS=>$newpoints, COORDINATE_LABELS=>$pin->COORDINATE_LABELS);
 +</code>
 +The facets can be viewed as usual:
 +
 +<code perl>
 +> print_constraints($q);
 +Facets:
 +0: 2 x1 - x2 >= -1
 +1: 0 >= -1
 +2: -x1 + 2 x2 >= -1
 +3: x1 + x2 >= 1
 +
 +</code>
 +<html><img src="/lib/exe/fetch.php?media=tutorials/release/4.11/optimization/ip-unb-integerhull.gif" class="center-block"></html><html><img></html>
 +
 +Note that the upper right part (including the red vertices) arises from truncation of the polyhedron for visualization.
 +
 +==== Mixed-Integer Case ====
 +
 +Let us now briefly discuss how to proceed if there are variables that are allowed to be integral. In this case there are several different types of information that one might be interested in. Let us first consider the question of how to compute the convex hull of all feasible integral variables, i.e., we consider the projection to the integral variables and then consider the convex hull of all feasible solutions. We only consider the bounded case, i.e., the original polyhedron is bounded.
 +
 +Consider the following example:
 +
 +<code>
 +Minimize
 + obj:  x1 + x2
 +Subject to
 +C1: s1 - 10 x1 <= 0
 +C2: s2 - 10 x2 <= 0
 +C3: s1 + s2 <= 1.5
 +C4: s1 + s2 >= 0.5
 +Bounds
 + 0 <= s1
 + 0 <= s2
 + 0 <= x1 <= 1
 + 0 <= x2 <= 1
 +General
 + x1 x2
 +End
 +</code>
 +In this example there are two integral variables ''%%x1%%'' and ''%%x2%%'', while ''%%s1%%'' and ''%%s2%%'' are continuous variables. Assuming the data is contained in the file ''%%mip.lp%%'', we proceed as follows:
 +
 +<code perl>
 +> $m=lp2poly('files/optimization/mip.lp');
 +> $p = new Polytope<Rational>($m);
 +</code>
 +We project the polyhedron in ''%%$p%%'' to the third and fourth variables as follows:
 +
 +<code perl>
 +> $q=projection($p, [3,4]);
 +</code>
 +We now construct the convex hull of all feasible points as above:
 +
 +<code perl>
 +> $s=new Polytope(POINTS=>$q->LATTICE_POINTS);
 +> print_constraints($s);
 +Facets:
 +0: -x1 >= -1
 +1: -x2 >= -1
 +2: x1 + x2 >= 1
 +
 +</code>
 +Thus, as expected, the convex hull equals the triangle with vertices ''%%{(0,1),(1,0),(1,1)}%%''.
 +
 +===== Integral Polytopes and Total Unimodularity =====
 +
 +As explained in the previous example, the integral points in a polytope are of particular interest in discrete optimization. These points are called //lattice points// in polymake and the corresponding convex hull //lattice polytope//. The handling of such polytopes is explained in more detail in the [[lattice_polytopes_tutorial|Tutorial for Lattice Polytopes]].
 +
 +Of particular interest for discrete optimization are properties of the original inequality system to define a lattice polytope, i.e., a polytope such that all of its vertices are integral (this can be tested by checking the property ''%%LATTICE%%''). One particularly interesting case occurs if the matrix defining the polytope is //totally unimodular// and the right hand side is integral.
 +
 +Using the polymake extension [[https://github.com/xammy/unimodularity-test/wiki/Polymake-Extension|Unimodularity]] by Matthias Walter, this can be checked as illustrated in the following examples.
 +
 +==== Example: Explicit Matrix ====
 +
 +In a first example, we directly create an integral matrix
 +
 +<code perl>
 +> $M=new Matrix<Integer>([[1,1,0,0],[1,0,1,0],[1,0,0,1]]);
 +</code>
 +The total unimodularity of this matrix can be checked as follows:
 +
 +<code perl>
 +> print totally_unimodular($M);
 +true
 +</code>
 +Thus, the given matrix is totally unimodular.
 +
 +==== Example: Matrix from Input ====
 +
 +In the second example, we reuse the file ''%%c3t%%'' from the example above. We read it into polymake:
 +
 +<code perl>
 +> $f=lp2poly('files/optimization/c3t.lp');
 +> $p = new Polytope<Rational>($f);
 +</code>
 +We now want to check whether the constraint matrix defined by the inequalities is totally unimodular (note that there are no equations in this example). Thus we first extract the inequality matrix without the first column (as an integer matrix) and then perform the test:
 +
 +<code perl>
 +> $A = new Matrix<Integer>($p->INEQUALITIES->minor(All, ~[0]));
 +> print totally_unimodular($A);
 +true
 +</code>
 +Thus, this matrix is totally unimodular as well.
 +
 +===== Total dual integrality =====
 +
 +Computations with respect to total dual integrality (TDI) can also be performed in polymake. Currently (August 2013), you need the perpetual beta version of polymake to access this functionality.
 +
 +The main functions are:
 +
 +  * The function ''%%totally_dual_integral%%'' takes an inequality system (as a matrix) and checks whether it is totally dual integral.
 +  * The function ''%%make_totally_dual_integral%%'' takes a polytope and returns a new polytope with inequalities that are TDI.
 +
 +Note that the input has to be full-dimensional in order to use these functions.
 +
 +To demonstrate the behavior of these functions, consider the 5-cycle example from above again:
 +
 +<code>
 +Maximize
 + obj: x#1 + x#2 + x#3 + x#4 + x#5
 +Subject to
 + edge_1: x#2 + x#1 <= 1
 + edge_2: x#3 + x#2 <= 1
 + edge_3: x#4 + x#3 <= 1
 + edge_4: x#5 + x#4 <= 1
 + edge_5: x#1 + x#5 <= 1
 +Bounds
 + 0 <= x#1 <= 1
 + 0 <= x#2 <= 1
 + 0 <= x#3 <= 1
 + 0 <= x#4 <= 1
 + 0 <= x#5 <= 1
 +General
 + x#1 x#2 x#3 x#4 x#5
 +End
 +</code>
 +Let us test whether the inequality system of this example is TDI. Thus, we first load the data as usual:
 +
 +<code perl>
 +> $f = lp2poly('files/optimization/stab.lp');
 +> $p = new Polytope<Rational>($f);
 +</code>
 +We now extract the corresponding inequality system and check it for TDIness:
 +
 +<code perl>
 +> $M = new Matrix<Rational>($p->INEQUALITIES);
 +> print totally_dual_integral($M);
 +false
 +</code>
 +The system is not TDI, which we expected from general theory, since we know that the polytope is not integral, but the system has integral coefficients. Consequently, let us construct a TDI-system for this polytope:
 +
 +<code perl>
 +> $q = make_totally_dual_integral($p);
 +> print_constraints($q);
 +Inequalities:
 +0: x5 >= 0
 +1: x4 >= 0
 +2: x3 >= 0
 +3: x2 >= 0
 +4: x1 >= 0
 +5: -x1 - x2 >= -1
 +6: -x1 - x5 >= -1
 +7: -x2 - x3 >= -1
 +8: -x3 - x4 >= -1
 +9: -x4 - x5 >= -1
 +10: -x1 - x2 - x3 - x4 - x5 >= -5/2
 +11: 0 >= -1
 +
 +</code>
 +As expected, the right hand side is non integral (otherwise, we know from general theory that the polytope would be integral as well). The result is now TDI:
 +
 +<code perl>
 +> $N = new Matrix<Rational>($q->INEQUALITIES);
 +> print totally_dual_integral($N);
 +true
 +</code>
 +//Note// that we need to take the inequalities instead of facets here, since facets are irredundant and thus might not be TDI, although the complete set of inequalities is TDI.
 +
 +===== Chvátal-Gomory Closure =====
 +
 +In the following we want to briefly show how closures of polytopes with respect to certain cutting-plane operators can be computed. We consider the two well-known cutting-plane operators here. The first one is the Chvátal-Gomory generator and the second one is the Lift-and-project operator as defined by Balas. For simplicity we will assume that the considered polytope is full-dimensional.
 +
 +==== Chvátal-Gomory Closure - Example 1 ====
 +
 +We first consider the polytope from the stable set problem from above:
 +
 +<code>
 +Maximize
 + obj: x#1 + x#2 + x#3 + x#4 + x#5
 +Subject to
 + edge_1: x#2 + x#1 <= 1
 + edge_2: x#3 + x#2 <= 1
 + edge_3: x#4 + x#3 <= 1
 + edge_4: x#5 + x#4 <= 1
 + edge_5: x#1 + x#5 <= 1
 +Bounds
 + 0 <= x#1 <= 1
 + 0 <= x#2 <= 1
 + 0 <= x#3 <= 1
 + 0 <= x#4 <= 1
 + 0 <= x#5 <= 1
 +General
 + x#1 x#2 x#3 x#4 x#5
 +End
 +</code>
 +As before we read in the file using ''%%lp2poly%%'':
 +
 +<code perl>
 +> $f = lp2poly('files/optimization/stab.lp');
 +> $p = new Polytope<Rational>($f);
 +</code>
 +The Chvátal-Gomory closure of a polytope can be computed with the function ''%%gc_closure%%''. The function takes a full-dimensional polytope and returns a new polytope. This contains the system of inequalities defining the closure in the property ''%%INEQUALITIES%%''. For our example, we obtain:
 +
 +<code perl>
 +> $g = gc_closure($p);
 +> print print_constraints($g);
 +Inequalities:
 +0: x5 >= 0
 +1: x4 >= 0
 +2: x3 >= 0
 +3: x2 >= 0
 +4: x1 >= 0
 +5: -x1 - x2 >= -1
 +6: -x1 - x5 >= -1
 +7: -x2 - x3 >= -1
 +8: -x3 - x4 >= -1
 +9: -x4 - x5 >= -1
 +10: -x1 - x2 - x3 - x4 - x5 >= -2
 +11: 0 >= -1
 +
 +</code>
 +Let us check whether the resulting polytope is integral:
 +
 +<code perl>
 +> print $g->LATTICE;
 +true
 +</code>
 +Thus, in this case, we have obtained the integer hull by one step of the Chvatal-Gomory-closure.
 +
 +==== Chvátal-Gomory Closure - Example 2 ====
 +
 +Let us now consider the classical example of a polytope with the vertices of simplex in d dimensions and the point 1/2 times (1, ..., 1). It can be shown that such a polytope has rank at least log(d) - 1, see [[http://www.box.net/shared/at1y8i3pq434bxt6m9xm|Pokutta, 2011]]]. In our example, we use d = 4:
 +
 +<code perl>
 +> $M = new Matrix<Rational>([[1,0,0,0,0],[1,1,0,0,0],[1,0,1,0,0],[1,0,0,1,0],[1,0,0,0,1],[1,1/2,1/2,1/2,1/2]]);
 +> $t = new Polytope<Rational>(POINTS => $M);
 +> $t1 = gc_closure($t);
 +> $t1->FACETS;
 +> print_constraints($t1);
 +Facets:
 +0: x4 >= 0
 +1: x3 >= 0
 +2: x2 >= 0
 +3: x1 >= 0
 +4: -x1 - x2 - x3 >= -1
 +5: -x1 - x2 - x4 >= -1
 +6: -x1 - x3 - x4 >= -1
 +7: -x2 - x3 - x4 >= -1
 +
 +> print $t1->LATTICE;
 +false
 +</code>
 +Thus, one round was not enough to produce an integral polytope. Indeed, the vertices are
 +
 +<code perl>
 +> $t1->VERTICES;
 +> print $t1->VERTICES;
 +1 1 0 0 0
 +1 0 0 0 0
 +1 1/3 1/3 1/3 1/3
 +1 0 1 0 0
 +1 0 0 1 0
 +1 0 0 0 1
 +</code>
 +However, one more round is enough:
 +
 +<code perl>
 +> $t2 = gc_closure($t1);
 +> $t2->FACETS;
 +> print_constraints($t2);
 +Facets:
 +0: x4 >= 0
 +1: x3 >= 0
 +2: x2 >= 0
 +3: x1 >= 0
 +4: -x1 - x2 - x3 - x4 >= -1
 +
 +> print $t2->LATTICE;
 +true
 +</code>
 +===== Lift-and-project closure =====
 +
 +The lift-and-project closure of a 0/1-polytope P can be generated as follows: for each coordinate compute the intersection of P with the pair of opposite cube faces and compute the convex hull. Then intersect the result with P. The following script defines a subroutine performing this operation - the code is somewhat complicated throught the fact that we need to check whether parts are empty.
 +
 +<code perl>
 +use application "polytope";
 +
 +sub lpclosure
 +{
 +    my $p = shift;
 +    my $d = $p->AMBIENT_DIM;
 +    my $q = new Polytope<Rational>($p);
 +    for (my $k = 0; $k < $d; $k = $k+1)
 +    {
 +        if ( $q->DIM == -1 )         # can stop as soon as $q is empty
 +        {
 +             return $q;
 +        }
 +    
 +        # create reversed opposite inequalities of 0/1-cube and corresponding polyhedra
 +        my $v1 = new Vector<Rational>(0 | -unit_vector($d, $k));
 +        my $v2 = new Vector<Rational>(-1 | unit_vector($d, $k));
 +    
 +        # create intersection of corresponding polyhedra with iterated polyhedron $q
 +        my $b1 = new Polytope<Rational>(INEQUALITIES => $v1 / $q->FACETS);
 +        my $b2 = new Polytope<Rational>(INEQUALITIES => $v2 / $q->FACETS);
 +    
 +        if ( ($b1->DIM > -1) && ($b2->DIM > -1) )
 +        {
 +            my $c = conv($b1, $b2);
 +            $q = intersection($q, $c);
 +        }
 +        elsif ( ($b1->DIM > -1) && ($b2->DIM == -1) )
 +        {
 +            $q = intersection($q, $b1);
 +        }
 +        elsif ( ($b1->DIM == -1) && ($b2->DIM > -1) )
 +        {
 +            $q = intersection($q, $b2);
 +        }
 +    }
 +    return $q;
 +}
 +> script("files/optimization/lpclosure.pl");
 +</code>
 +==== Lift-and-Project Closure - Example 1 ====
 +
 +For our well known stable set example, we get the following:
 +
 +<code perl>
 +> $q = lpclosure($p);
 +> $q->FACETS;
 +> print_constraints($q);
 +Facets:
 +0: -x1 - x2 >= -1
 +1: x3 >= 0
 +2: x2 >= 0
 +3: x4 >= 0
 +4: -x3 - x4 >= -1
 +5: -x4 - x5 >= -1
 +6: x1 >= 0
 +7: -x2 - x3 >= -1
 +8: -x1 - x2 - x3 - x4 - x5 >= -2
 +9: -x1 - x5 >= -1
 +10: x5 >= 0
 +
 +</code>
 +Thus, the lift-and-project closure in this case gives the integral hull (as we have seen above).
 +
 +==== Lift-and-Project Closure - Example 2 ====
 +
 +Let us now consider the same example as for CG-closures:
 +
 +<code perl>
 +> $M = new Matrix<Rational>([[1,0,0,0,0],[1,1,0,0,0],[1,0,1,0,0],[1,0,0,1,0],[1,0,0,0,1],[1,1/2,1/2,1/2,1/2]]);
 +> $p = new Polytope<Rational>(POINTS => $M);
 +> $q = lpclosure($p);
 +> $q->FACETS;
 +> print_constraints($q);
 +Facets:
 +0: x1 >= 0
 +1: x4 >= 0
 +2: -x1 - x2 - x3 - x4 >= -1
 +3: x2 >= 0
 +4: x3 >= 0
 +
 +</code>
 +Thus, we have obtained the integral hull in a single step of the lift-and-project closure as opposed to two steps in the CG-closure.
 +
 +==== Alternative Input: mps2poly ====
 +
 +Another input format that ''%%polymake%%'' can handle is the MPS-format used by soplex or gurobi. This format is more focused on the columns then on the rows of a linear program.
 +
 +<code>
 +Name          0/1-cube
 +Rows
 +  C0000000
 +  R0000000
 +COLUMNS
 +    x1        C0000000  1                        R0000000  -1
 +    x2        C0000000  1                        R0000000  -1
 +    x3        C0000000  1                        R0000000  -1
 +RHS
 +    B         R0000000  -2
 +BOUNDS
 + LO BND       x1        0
 + UP BND       x1        1
 + LO BND       x2        0
 + UP BND       x2        1
 + LO BND       x3        0
 + UP BND       x3        1
 +ENDDATA
 +</code>
 +This is again a description of the same 0/1-Cube like above.\\
 +If the example is contained in the file ''%%c3t.mps%%'', we can create again a polytope from it via:
 +
 +<code perl>
 +> $p = mps2poly('files/optimization/c3t.mps');
 +</code>
 +It gives use the same polytope with the same linear program as above.
 +
 +<code perl>
 +> print $p->VISUAL->DIRECTED_GRAPH;
 +</code>
 +<HTML>
 +<!--
 +polymake for knusper
 +Thu Mar  3 00:38:29 2022
 +p
 +-->
 +
 +
 +<html>
 +   <head>
 +      <meta charset=utf-8>
 +      <title>p</title>
 +      <style>
 +/*
 +// COMMON_CODE_BLOCK_BEGIN
 +*/
 +         html {overflow: scroll;}
 +         strong{font-size: 18px;}
 +         canvas { z-index: 8; }
 +         input[type='radio'] {margin-left:0;}
 +         input[type='checkbox'] {margin-right:7px; margin-left: 0px; padding-left:0px;}
 +         .group{padding-bottom: 15px;}
 +         .settings * {z-index: 11; }
 +         .settings{z-index: 10; font-family: Arial, Helvetica, sans-serif; margin-left: 30px; visibility: hidden; width: 14em; height: 96%; border: solid 1px silver; padding: 2px; overflow-y: scroll; box-sizing: border-box; background-color: white; position: absolute;}
 +         .indented{margin-left: 20px; margin-top: 10px; padding-bottom: 0px;} 
 +         .shownObjectsList{overflow: auto; max-width: 150px; max-height: 150px;}
 +         .showSettingsButton{visibility: visible; z-index: 12; position: absolute }
 +         .hideSettingsButton{visibility: hidden; z-index: 12; position: absolute; opacity: 0.5}
 +         button{margin-left: 0; margin-top: 10px}
 +         img{cursor: pointer;}
 +         .suboption{padding-top: 15px;}
 +         #model75625263093 { width: 100%; height: 100%; }
 +         .threejs_container { width: 100%; height: 75vh;}
 +         .settings{max-height: 74vh} 
 +         input[type=range] {
 +           -webkit-appearance: none;
 +           padding:0; 
 +           width:90%; 
 +           margin-left: auto;
 +           margin-right: auto;
 +           margin-top: 15px;
 +           margin-bottom: 15px;
 +           display: block;
 +         }
 +         input[type=range]:focus {
 +           outline: none;
 +         }
 +         input[type=range]::-webkit-slider-runnable-track {
 +           height: 4px;
 +           cursor: pointer;
 +           animate: 0.2s;
 +           box-shadow: 0px 0px 0px #000000;
 +           background: #E3E3E3;
 +           border-radius: 0px;
 +           border: 0px solid #000000;
 +         }
 +         input[type=range]::-webkit-slider-thumb {
 +           box-shadow: 1px 1px 2px #B8B8B8;
 +           border: 1px solid #ABABAB;
 +           height: 13px;
 +           width: 25px;
 +           border-radius: 20px;
 +           background: #E0E0E0;
 +           cursor: pointer;
 +           -webkit-appearance: none;
 +           margin-top: -5px;
 +         }
 +         input[type=range]:focus::-webkit-slider-runnable-track {
 +           background: #E3E3E3;
 +         }
 +         input[type=range]::-moz-range-track {
 +           height: 4px;
 +           cursor: pointer;
 +           animate: 0.2s;
 +           box-shadow: 0px 0px 0px #000000;
 +           background: #E3E3E3;
 +           border-radius: 0px;
 +           border: 0px solid #000000;
 +         }
 +         input[type=range]::-moz-range-thumb {
 +           box-shadow: 1px 1px 2px #B8B8B8;
 +           border: 1px solid #ABABAB;
 +           height: 13px;
 +           width: 25px;
 +           border-radius: 20px;
 +           background: #E0E0E0;
 +           cursor: pointer;
 +         }
 +         input[type=range]::-ms-track {
 +           height: 4px;
 +           cursor: pointer;
 +           animate: 0.2s;
 +           background: transparent;
 +           border-color: transparent;
 +           color: transparent;
 +         }
 +         input[type=range]::-ms-fill-lower {
 +           background: #E3E3E3;
 +           border: 0px solid #000000;
 +           border-radius: 0px;
 +           box-shadow: 0px 0px 0px #000000;
 +         }
 +         input[type=range]::-ms-fill-upper {
 +           background: #E3E3E3;
 +           border: 0px solid #000000;
 +           border-radius: 0px;
 +           box-shadow: 0px 0px 0px #000000;
 +         }
 +         input[type=range]::-ms-thumb {
 +           box-shadow: 1px 1px 2px #B8B8B8;
 +           border: 1px solid #ABABAB;
 +           height: 13px;
 +           width: 25px;
 +           border-radius: 20px;
 +           background: #E0E0E0;
 +           cursor: pointer;
 +         }
 +         input[type=range]:focus::-ms-fill-lower {
 +           background: #E3E3E3;
 +         }
 +         input[type=range]:focus::-ms-fill-upper {
 +           background: #E3E3E3;
 +         }
 +/*
 +// COMMON_CODE_BLOCK_END
 +*/
 + </style>
 +   </head>
 +<body>
 +   <div class='threejs_container'>
 + <div id='settings_2' class='settings'>
 + <div class=group id='explode_2'>
 + <strong>Explode</strong>
 + <input id='explodeRange_2' type='range' min='0.00001' max=6 step=0.01 value=0.00001>
 + <div class=indented><input id='explodeCheckbox_2' type='checkbox'>Automatic explosion</div>
 + <div class=suboption>Exploding speed</div>
 + <input id='explodingSpeedRange_2' type='range' min=0 max=0.5 step=0.001 value=0.05>
 + </div>
 +
 + <div class=group id='transparency_2' class='transparency'>
 + <strong>Transparency</strong>
 + <input id='transparencyRange_2' type='range' min=0 max=1 step=0.01 value=0>
 +            <div class=indented><input id='depthWriteCheckbox_2' type='checkbox'>depthWrite</div>
 + </div>
 +
 + <div class=group id='rotation_2'>
 + <strong>Rotation</strong>
 + <div class=indented>
 + <div><input type='checkbox' id='changeRotationX_2'> x-axis</div>
 + <div><input type='checkbox' id='changeRotationY_2'> y-axis</div>
 + <div><input type='checkbox' id='changeRotationZ_2'> z-axis</div>
 + <button id='resetButton_2'>Reset</button>
 + </div>
 +
 + <div class=suboption>Rotation speed</div>
 + <input id='rotationSpeedRange_2' type='range' min=0 max=5 step=0.01 value=2>
 + </div>
 +
 +
 + <div class=group id='display_2'>
 + <strong>Display</strong>
 + <div class=indented>
 + <div id='shownObjectTypesList_2' class='shownObjectsList'></div>
 + </div>
 + <div class=suboption>Objects</div>
 + <div class=indented>
 +    <div id='shownObjectsList_2' class='shownObjectsList'></div>
 + </div>
 + </div>
 +         
 +         <div class=group id='camera_2'>
 +            <strong>Camera</strong>
 +            <div class=indented>
 +               <form>
 +                  <select id="cameraType_2">
 +                     <option value='perspective' selected> Perspective<br></option>
 +                     <option value='orthographic' > Orthographic<br></option>
 +                  </select>
 +               </form>
 +            </div>
 +         </div>
 +
 + <div class=group id='svg_2'>
 + <strong>SVG</strong>
 + <div class=indented>
 + <form>
 + <input type="radio" name='screenshotMode' value='download' id='download_2' checked> Download<br>
 + <input type="radio" name='screenshotMode' value='tab' id='tab_2' > New tab<br>
 + </form>
 + <button id='takeScreenshot_2'>Screenshot</button>
 + </div>
 + </div>
 +
 + </div> <!-- end of settings -->
 + <img id='hideSettingsButton_2' class='hideSettingsButton' src='/kernelspecs/r118/polymake/close.svg' width=20px">
 + <img id='showSettingsButton_2' class='showSettingsButton' src='/kernelspecs/r118/polymake/menu.svg' width=20px">
 +<div id="model75625263093"></div>
 +</div>
 +   <script>
 +    requirejs.config({
 +      paths: {
 +        three: '/kernelspecs/r118/polymake/three',
 +        TrackballControls: '/kernelspecs/r118/polymake/TrackballControls',
 +        OrbitControls: '/kernelspecs/r118/polymake/OrbitControls',
 +        Projector: '/kernelspecs/r118/polymake/Projector',
 +        SVGRenderer: '/kernelspecs/r118/polymake/SVGRenderer',
 +        WEBGL: '/kernelspecs/r118/polymake/WebGL',
 +      },
 +      shim: {
 +        'three': { exports: 'THREE'},
 +        'SVGRenderer': { deps: [ 'three' ], exports: 'THREE.SVGRenderer' },
 +        'WEBGL': { deps: [ 'three' ], exports: 'THREE.WEBGL' },
 +        'Projector': { deps: [ 'three' ], exports: 'THREE.Projector' },
 +        'TrackballControls': { deps: [ 'three' ], exports: 'THREE.TrackballControls' },
 +        'OrbitControls': { deps: [ 'three' ], exports: 'THREE.OrbitControls' },
 +      }
 +    });
 +    
 +    require(['three'],function(THREE){
 +        window.THREE = THREE;
 +      require(['TrackballControls', 'OrbitControls', 'Projector', 'SVGRenderer', 'WEBGL'],
 +               function(TrackballControls, OrbitControls, Projector, SVGRenderer, WEBGL) {
 +    THREE.TrackballControls = TrackballControls;
 +    THREE.OrbitControls = OrbitControls;
 +    THREE.Projector = Projector;
 +    THREE.SVGRenderer = SVGRenderer;
 +    THREE.WEBGL = WEBGL;
 +
 +// COMMON_CODE_BLOCK_BEGIN
 +
 +const intervalLength = 25; // for automatic animations
 +const explodableModel = true; 
 +const modelContains = { points: false, pointlabels: false, lines: false, edgelabels: false, faces: false, arrowheads: false };
 +const foldables = [];
 +
 +var three = document.getElementById("model75625263093");
 +var scene = new THREE.Scene();
 +var renderer = new THREE.WebGLRenderer( { antialias: true } );
 +var svgRenderer = new THREE.SVGRenderer( { antialias: true } );
 +renderer.setPixelRatio( window.devicePixelRatio );
 +renderer.setClearColor(0xFFFFFF, 1);
 +svgRenderer.setClearColor(0xFFFFFF, 1);
 +three.appendChild(renderer.domElement);
 +
 +var frustumSize = 4;
 +var cameras = [new THREE.PerspectiveCamera(75, 1, 0.1, 1000), new THREE.OrthographicCamera()];
 +cameras.forEach(function(cam) {
 +    cam.position.set(0, 0, 5);
 +    cam.lookAt(0, 0, 0);  
 +    cam.up.set(0, 1, 0);         
 +});
 +var controls = [new THREE.TrackballControls(cameras[0], three), new THREE.OrbitControls(cameras[1], three)];
 +var camera, control;
 +
 +controls[0].zoomSpeed = 0.2;
 +controls[0].rotateSpeed = 4;
 +
 +
 +// class to allow move points together with labels and spheres
 +var PMPoint = function (x,y,z) {
 +   this.vector = new THREE.Vector3(x,y,z);
 +   this.sprite = null;
 +   this.sphere = null;
 +}
 +PMPoint.prototype.addLabel = function(labelsprite) {
 +   this.sprite = labelsprite;
 +   this.sprite.position.copy(this.vector);
 +}
 +PMPoint.prototype.addSphere = function(spheremesh) {
 +   this.sphere = spheremesh;
 +   this.sphere.position.copy(this.vector);
 +}
 +PMPoint.prototype.set = function(x,y,z) {
 +   this.vector.set(x,y,z);
 +   if (this.sprite) {
 +      this.sprite.position.copy(this.vector);
 +   }
 +   if (this.sphere) {
 +      this.sphere.position.copy(this.vector);
 +   }
 +}
 +PMPoint.prototype.radius = function() {
 +   if (this.sphere) {
 +      return this.sphere.geometry.parameters.radius;
 +   } else {
 +      return 0;
 +   }
 +};
 +// select the target node
 +var target = document.querySelector('#model75625263093');
 +
 +// create an observer instance
 +var observer = new MutationObserver(function(mutations) {
 +   mutations.forEach(function(mutation) {
 +      if (mutation.removedNodes && mutation.removedNodes.length > 0) {
 +         cancelAnimationFrame(renderId);
 +         observer.disconnect();
 +         console.log("cancelled frame "+renderId);
 +      }
 +   });
 +});
 +
 +// configuration of the observer:
 +var config = { childList: true, characterData: true }
 +
 +// pass in the target node, as well as the observer options
 +while (target) {
 +   if (target.className=="output") {
 +      observer.observe(target, config);
 +      break;
 +   }
 +   target = target.parentNode;
 +}
 +
 +// COMMON_CODE_BLOCK_END
 +
 +var obj0 = new THREE.Object3D();
 +obj0.name = "p";
 +obj0.userData.explodable = 1;
 +obj0.userData.points = [];
 +obj0.userData.points.push(new PMPoint(0, 0, 1));
 +obj0.userData.points.push(new PMPoint(1, 0, 0));
 +obj0.userData.points.push(new PMPoint(0, 0, 0));
 +obj0.userData.points.push(new PMPoint(0, 1, 0));
 +obj0.userData.points.push(new PMPoint(1, 1, 0));
 +obj0.userData.points.push(new PMPoint(1, 0, 1));
 +obj0.userData.points.push(new PMPoint(0, 1, 1));
 +
 +obj0.userData.pointradii = 0.02;
 +   <!-- Vertex style -->
 +obj0.userData.pointmaterial = new THREE.MeshBasicMaterial( { color: 0xFF0000, side: THREE.DoubleSide, transparent: false } );
 +obj0.userData.pointlabels = ["0", "1", "2", "3", "4", "5", "6"];
 +obj0.userData.edgeindices = [0, 2, 1, 2, 2, 3, 1, 4, 3, 4, 0, 5, 1, 5, 4, 5, 0, 6, 3, 6, 4, 6, 5, 6];
 +obj0.userData.facets = [[5, 4, 6], [6, 3, 2, 0], [1, 4, 5], [5, 0, 2, 1], [6, 4, 3], [3, 4, 1, 2], [5, 6, 0]];
 +   <!-- Facet style -->
 +obj0.userData.facetmaterial = new THREE.MeshBasicMaterial( { color: 0x77EC9E, depthFunc: THREE.LessDepth, depthWrite: false, opacity: 1, polygonOffset: true, polygonOffsetFactor: 1, polygonOffsetUnits: 0.5, side: THREE.DoubleSide, transparent: true } );
 +init_object(obj0);
 +scene.add(obj0);
 +
 +var obj1 = new THREE.Object3D();
 +obj1.name = "GRAPH directed with LP ";
 +obj1.userData.explodable = 1;
 +obj1.userData.points = [];
 +obj1.userData.points.push(new PMPoint(0, 0, 1));
 +obj1.userData.points.push(new PMPoint(1, 0, 0));
 +obj1.userData.points.push(new PMPoint(0, 0, 0));
 +obj1.userData.points.push(new PMPoint(0, 1, 0));
 +obj1.userData.points.push(new PMPoint(1, 1, 0));
 +obj1.userData.points.push(new PMPoint(1, 0, 1));
 +obj1.userData.points.push(new PMPoint(0, 1, 1));
 +
 +obj1.userData.edgeindices = [0, 5, 0, 6, 1, 4, 1, 5, 2, 0, 2, 1, 2, 3, 3, 4, 3, 6];
 +   <!-- Edge style -->
 +obj1.userData.edgematerial = new THREE.LineBasicMaterial( { color: 0x000000, linewidth: 1.5, transparent: false } );
 +obj1.userData.arrowstyle = 1;
 +init_object(obj1);
 +scene.add(obj1);
 +
 +// COMMON_CODE_BLOCK_BEGIN
 +function textSpriteMaterial(message, parameters) {
 +    if ( parameters === undefined ) parameters = {};
 +    var fontface = "Helvetica";
 +    var fontsize = parameters.hasOwnProperty("fontsize") ? parameters["fontsize"] : 15;
 +    fontsize = fontsize*10;
 +    var lines = message.split('\\n');
 +    var size = 512;
 +    for(var i = 0; i<lines.length; i++){
 +        var tmp = lines[i].length;
 +        while(tmp*fontsize > size){
 +           fontsize--;
 +        }
 +    }
 +    
 +    var canvas = document.createElement('canvas');
 +    canvas.width = size;
 +    canvas.height = size;
 +    var context = canvas.getContext('2d');
 +    context.fillStyle = "rgba(255, 255, 255, 0)";
 +    context.fill();
 +    context.font = fontsize + "px " + fontface;
 +    
 +    // text color
 +    context.fillStyle = "rgba(0, 0, 0, 1.0)";
 +     for(var i = 0; i<lines.length; i++){
 +        context.fillText(lines[i], size/2, size/2+i*fontsize);
 +     }
 +    
 +    // canvas contents will be used for a texture
 +    var texture = new THREE.Texture(canvas);
 +    texture.needsUpdate = true;
 +    
 +    var spriteMaterial = new THREE.SpriteMaterial({map: texture, depthTest: true, depthWrite: false, polygonOffset: true, polygonOffsetFactor: -1, polygonOffsetUnits: 1 });
 +    return spriteMaterial;
 +}
 +
 +
 +// ---------------------- INITIALIZING OBJECTS--------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +function init_object(obj) {
 +    if (obj.userData.hasOwnProperty("pointmaterial")) {
 +        init_points(obj);
 +        modelContains.points = true;
 +    }
 +    if (obj.userData.hasOwnProperty("pointlabels")) {
 +        init_pointlabels(obj);
 +        modelContains.pointlabels = true;
 +    }
 +    if (obj.userData.hasOwnProperty("edgematerial")) {
 +        init_lines(obj);
 +        modelContains.lines = true;
 +    }
 +    if (obj.userData.hasOwnProperty("edgelabels")) {
 +        init_edgelabels(obj);
 +        modelContains.edgelabels = true;
 +    }
 +    if (obj.userData.hasOwnProperty("arrowstyle")) {
 +        init_arrowheads(obj);
 +        modelContains.arrowheads = true;
 +    }
 +    if (obj.userData.hasOwnProperty("facetmaterial")) {
 +        init_faces(obj);
 +        modelContains.faces = true;
 +    }
 +}
 +
 +function init_points(obj) {
 +    var pointgroup = new THREE.Group();
 +    pointgroup.name = "points";
 +    var points = obj.userData.points;
 +    var radii = obj.userData.pointradii;
 +    var materials = obj.userData.pointmaterial;
 +    var geometry,material;
 +    if (!Array.isArray(radii)) {
 +        geometry = new THREE.SphereBufferGeometry(radii);  
 +    }
 +    if (!Array.isArray(materials)) {
 +        material = materials;
 +    }
 +    for (var i=0; i<points.length; i++) {
 +        var point = points[i];
 +        if (Array.isArray(radii)) {
 +            if (radii[i] == 0) {
 +                continue;
 +            }
 +            geometry = new THREE.SphereBufferGeometry(radii[i]);  
 +        } 
 +        if (Array.isArray(materials)) {
 +            material = materials[i];     
 +        } 
 +        var sphere = new THREE.Mesh(geometry, material);
 +        point.addSphere(sphere);
 +        pointgroup.add(sphere);
 +    }
 +    obj.add(pointgroup);
 +}
 +
 +function init_pointlabels(obj) {
 +    var points = obj.userData.points;
 +    var labels = obj.userData.pointlabels;
 +    var pointlabels = new THREE.Group();
 +    pointlabels.name = "pointlabels";
 +    if (Array.isArray(labels)) {
 +        for (var i=0; i<points.length; i++) {
 +            var point = points[i];
 +            var spriteMaterial = textSpriteMaterial( labels[i] );
 +         var sprite = new THREE.Sprite(spriteMaterial);
 +            point.addLabel(sprite);
 +            pointlabels.add(sprite);
 +        }
 +    } else {
 +        var spriteMaterial = textSpriteMaterial( labels );
 +        for (var i=0; i<points.length; i++) {
 +            var point = points[i];
 +         var sprite = new THREE.Sprite(spriteMaterial);
 +            point.addLabel(sprite);
 +            pointlabels.add(sprite);
 +        }
 +    }
 +    obj.add(pointlabels);
 +}
 +
 +function init_lines(obj) {
 +    var edgeindices = obj.userData.edgeindices;
 +    var points = obj.userData.points;
 +    var materials = obj.userData.edgematerial;
 +    var geometry = new THREE.BufferGeometry();
 +    var bufarr = new Float32Array( obj.userData.edgeindices.length * 3 );
 +    var bufattr = new THREE.Float32BufferAttribute( bufarr, 3 );
 +    var geometry = new THREE.BufferGeometry();
 +    geometry.setAttribute('position', bufattr);
 +    if (Array.isArray(materials)) {     
 +        for (var i=0; i<materials.length; i++) {
 +            geometry.addGroup(2*i,2,i);
 +        }
 +    }
 +    var lines = new THREE.LineSegments(geometry, materials);
 +    lines.name = "lines";
 +    obj.add(lines);
 +    updateEdgesPosition(obj);
 +}
 +
 +function init_edgelabels(obj) {
 +    var points = obj.userData.points;
 +    var edgeindices = obj.userData.edgeindices;
 +    var labels = obj.userData.edgelabels;
 +    var edgelabels = new THREE.Group();
 +    edgelabels.name = "edgelabels";
 +    if (Array.isArray(labels)) {
 +        for (var i=0; i<edgeindices.length/2; i++) {
 +            var spriteMaterial = textSpriteMaterial( labels[i] );
 +            var sprite = new THREE.Sprite(spriteMaterial);
 +            sprite.position.copy(new THREE.Vector3().addVectors(points[edgeindices[2*i]].vector,points[edgeindices[2*i+1]].vector).multiplyScalar(0.5));
 +            edgelabels.add(sprite);
 +        }
 +    } else {
 +        var spriteMaterial = textSpriteMaterial( labels );
 +        for (var i=0; i<edgeindices.length/2; i++) {
 +            var sprite = new THREE.Sprite(spriteMaterial);
 +            sprite.position.copy(new THREE.Vector3().addVectors(points[edgeindices[2*i]].vector,points[edgeindices[2*i+1]].vector).multiplyScalar(0.5));
 +            edgelabels.add(sprite);
 +        }
 +    }
 +    obj.add(edgelabels);
 +}
 +
 +function init_arrowheads(obj) {
 +    var arrowheads = new THREE.Group();
 +    arrowheads.name = "arrowheads";
 +    var arrowstyle = obj.userData.arrowstyle;
 +    var edgeindices = obj.userData.edgeindices;
 +    var edgematerials = obj.userData.edgematerial;
 +    var points = obj.userData.points;
 +    var material;
 +    if (!Array.isArray(edgematerials)) {
 +        material = new THREE.MeshBasicMaterial( {color: edgematerials.color} );
 +    }
 +
 +    for (var i=0; i<edgeindices.length; i=i+2) {
 +        var start = points[edgeindices[i]];
 +        var end = points[edgeindices[i+1]];
 +        var dist = start.vector.distanceTo( end.vector ) - start.radius() - end.radius();
 +        if (dist <= 0) {
 +            continue;
 +        }
 +        var dir = new THREE.Vector3().subVectors(end.vector,start.vector);
 +        dir.normalize();
 +        var axis = new THREE.Vector3().set(dir.z,0,-dir.x);
 +        axis.normalize();
 +        var radians = Math.acos( dir.y );
 +        var radius = dist/25;
 +        var height = dist/5;
 +        var geometry = new THREE.ConeBufferGeometry(radius,height);
 +        var position = new THREE.Vector3().addVectors(start.vector,dir.clone().multiplyScalar(start.radius()+dist-height/2));
 +        if (Array.isArray(edgematerials)) {
 +            material = new THREE.MeshBasicMaterial( {color: edgematerials[i].color} );
 +        }
 +        var cone = new THREE.Mesh( geometry, material );
 +        cone.quaternion.setFromAxisAngle(axis,radians);;
 +        cone.position.copy(position);;
 +        arrowheads.add(cone);
 +    }
 +    obj.add(arrowheads);
 +}
 +
 +function init_faces(obj) {
 +    var points = obj.userData.points;
 +    var facets = obj.userData.facets;
 +    obj.userData.triangleindices = [];
 +    for (var i=0; i<facets.length; i++) {
 +        facet = facets[i];
 +        for (var t=0; t<facet.length-2; t++) {
 +            obj.userData.triangleindices.push(facet[0],facet[t+1],facet[t+2]);  
 +        }
 +    }
 +    var bufarr = new Float32Array( obj.userData.triangleindices.length * 3 );
 +    var bufattr = new THREE.Float32BufferAttribute(bufarr,3);
 +    
 +    var materials = obj.userData.facetmaterial;
 +    var geometry = new THREE.BufferGeometry();
 +    var frontmaterials = [];
 +    var backmaterials = [];
 +    geometry.setAttribute('position',bufattr);
 +    if (Array.isArray(materials)) {
 +        var tricount = 0;
 +        var facet;
 +        for (var i=0; i<facets.length; i++) {
 +            facet = facets[i];
 +            geometry.addGroup(tricount,(facet.length-2)*3,i);
 +            tricount += (facet.length-2)*3;
 +        }
 +        for (var j=0; j<materials.length; j++) {
 +            var fmat = materials[j].clone()
 +            fmat.side = THREE.FrontSide;
 +            frontmaterials.push(fmat);
 +            var bmat = materials[j].clone()
 +            bmat.side = THREE.BackSide;
 +            backmaterials.push(bmat);
 +            obj.userData.facetmaterial = frontmaterials.concat(backmaterials);
 +        }
 +    } else if (materials instanceof THREE.Material) {
 +        frontmaterials = materials.clone()
 +        frontmaterials.side = THREE.FrontSide;
 +        backmaterials = materials.clone()
 +        backmaterials.side = THREE.BackSide;
 +        obj.userData.facetmaterial = [frontmaterials, backmaterials];
 +    }
 +    // duplicating the object with front and back should avoid transparency issues
 +    var backmesh = new THREE.Mesh(geometry, backmaterials);
 +    // meshname is used to show/hide objects
 +    backmesh.name = "backfaces";
 +    obj.add(backmesh);
 +    var frontmesh = new THREE.Mesh(geometry, frontmaterials);
 +    frontmesh.name = "frontfaces";
 +    obj.add(frontmesh);
 +    updateFacesPosition(obj);
 +}
 +// //INITIALIZING
 +
 +
 +function updateFacesPosition(obj) {
 +    var points = obj.userData.points;
 +    var indices = obj.userData.triangleindices;
 +    var faces = obj.getObjectByName("frontfaces");
 +    var ba = faces.geometry.getAttribute("position");
 +    for (var i=0; i<indices.length; i++) {
 +        ba.setXYZ(i, points[indices[i]].vector.x, points[indices[i]].vector.y ,points[indices[i]].vector.z); 
 +    }
 +    faces.geometry.attributes.position.needsUpdate = true;
 +    
 +}
 +
 +function updateEdgesPosition(obj) {
 +    var points = obj.userData.points;
 +    var indices = obj.userData.edgeindices;
 +    var lines = obj.getObjectByName("lines");
 +    var ba = lines.geometry.getAttribute("position"); 
 +    for (var i=0; i<indices.length; i++) {
 +        ba.setXYZ(i, points[indices[i]].vector.x, points[indices[i]].vector.y ,points[indices[i]].vector.z); 
 +    }
 +    lines.geometry.attributes.position.needsUpdate = true;
 +}
 +
 +function onWindowResize() {
 +    renderer.setSize( three.clientWidth, three.clientHeight );
 +    svgRenderer.setSize( three.clientWidth, three.clientHeight );
 +    updateCamera();
 +}
 +
 +function updateCamera() {
 +    var width = three.clientWidth;
 +    var height = three.clientHeight;
 +    var aspect = width / height;
 +    if (camera.type == "OrthographicCamera") {
 +        camera.left = frustumSize * aspect / - 2;
 +        camera.right = frustumSize * aspect / 2;
 +        camera.top = frustumSize / 2;
 +        camera.bottom = - frustumSize / 2;
 +    } else if (camera.type == "PerspectiveCamera") {
 +        camera.aspect = aspect;
 +    }
 +    camera.updateProjectionMatrix();
 +}
 +
 +function changeCamera(event) {
 +    var selindex = event.currentTarget.selectedIndex;
 +    camera = cameras[selindex];
 +    control = controls[selindex];
 +    control.enabled = true; 
 +    for (var i=0; i<controls.length; i++) {
 +        if (i!=selindex) {
 +            controls[i].enabled = false;
 +        }
 +    }
 +    updateCamera();
 +}
 +
 +var camtypenode = document.getElementById('cameraType_2');
 +camtypenode.onchange = changeCamera;
 +camtypenode.dispatchEvent(new Event('change'));
 +
 +onWindowResize();
 +window.addEventListener('resize', onWindowResize);
 +
 +
 +var xRotationEnabled = false;
 +var yRotationEnabled = false;
 +var zRotationEnabled = false;
 +var rotationSpeedFactor = 1;
 +var settingsShown = false;
 +var labelsShown = true;
 +var intervals = [];
 +var timeouts = [];
 +var explodingSpeed = 0.05;
 +var explodeScale = 0.000001;
 +var XMLS = new XMLSerializer();
 +var svgElement;
 +var renderId;
 +
 +var render = function () {
 +
 + renderId = requestAnimationFrame(render);
 +
 +// comment in for automatic explosion
 +// explode(updateFactor());
 +
 +    var phi = 0.02 * rotationSpeedFactor;
 +
 +    if (xRotationEnabled) {
 +        scene.rotation.x += phi;
 +    }
 +    if (yRotationEnabled) {
 +        scene.rotation.y += phi;
 +    }
 +    if (zRotationEnabled) {
 +        scene.rotation.z += phi;
 +    }
 +
 +    control.update();
 +    renderer.render(scene, camera);
 +};
 +
 +if ( THREE.WEBGL.isWebGLAvailable() ) {
 + render();
 +} else {
 + var warning = WEBGL.getWebGLErrorMessage();
 + three.appendChild( warning );
 +}
 +    
 +function changeTransparency() {
 +    var opacity = 1-Number(event.currentTarget.value);
 +    for (var i=0; i<scene.children.length; i++) {
 +        child = scene.children[i];
 +        if ( child.userData.hasOwnProperty("facetmaterial") ) {
 +            if (Array.isArray(child.userData.facetmaterial)) {
 +                for (var j=0; j<child.userData.facetmaterial.length; j++) {
 +                    child.userData.facetmaterial[j].opacity = opacity;
 +                }
 +            } else {
 +                child.userData.facetmaterial.opacity = opacity;
 +            }    
 +        }
 +    }
 +}
 +
 +function toggleDepthWrite(event) {
 +    depthwrite = event.currentTarget.checked;
 +    for (var i=0; i<scene.children.length; i++) {
 +        child = scene.children[i];
 +        if ( child.userData.hasOwnProperty("facetmaterial") ) {
 +            if (Array.isArray(child.userData.facetmaterial)) {
 +                for (var j=0; j<child.userData.facetmaterial.length; j++) {
 +                    child.userData.facetmaterial[j].depthWrite = depthwrite;
 +                }
 +            } else {
 +                child.userData.facetmaterial.depthWrite = depthWrite;
 +            }    
 +        }
 +    }
 +}
 +
 +function changeRotationX(event){
 +    xRotationEnabled = event.currentTarget.checked;
 +}
 +
 +function changeRotationY(event){
 +    yRotationEnabled = event.currentTarget.checked;
 +}
 +
 +function changeRotationZ(event){
 +    zRotationEnabled = event.currentTarget.checked;
 +}
 +
 +
 +function changeRotationSpeedFactor(event){
 +    rotationSpeedFactor = Number(event.currentTarget.value);
 +}
 +
 +function resetScene(){
 +    scene.rotation.set(0,0,0);
 +    camera.position.set(0,0,5);
 +    camera.up.set(0,1,0);
 +}
 +
 +function showSettings(event){
 +    document.getElementById('settings_2').style.visibility = 'visible';
 +    document.getElementById('showSettingsButton_2').style.visibility = 'hidden';
 +    document.getElementById('hideSettingsButton_2').style.visibility = 'visible';
 +    settingsShown = true;
 +}
 +
 +function hideSettings(event){
 +    document.getElementById('settings_2').style.visibility = 'hidden';
 +    document.getElementById('showSettingsButton_2').style.visibility = 'visible';
 +    document.getElementById('hideSettingsButton_2').style.visibility = 'hidden';
 +    settingsShown = false;
 +}
 +
 +
 +
 +var pos = 150* Math.PI;
 +
 +function updateFactor() {
 +    pos++;
 +    return Math.sin(.01*pos)+1;
 +}
 +
 +// ------------------------ FOLDING ------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +// rotate point p around axis defined by points p1 and p2 by given angle
 +function rotate(p, p1, p2, angle ){   
 +    angle = -angle;
 +    var x = p.x, y = p.y, z = p.z, 
 +    a = p1.x, b = p1.y, c = p1.z, 
 +    u = p2.x-p1.x, v = p2.y-p1.y, w = p2.z-p1.z;
 +    var result = [];
 +    var L = u*u + v*v + w*w;
 +    var sqrt = Math.sqrt;
 +    var cos = Math.cos;
 +    var sin = Math.sin;
 +
 +    result[0] = ((a*(v*v+w*w)-u*(b*v+c*w-u*x-v*y-w*z))*(1-cos(angle))+L*x*cos(angle)+sqrt(L)*(-c*v+b*w-w*y+v*z)*sin(angle))/L;
 +    result[1] = ((b*(u*u+w*w)-v*(a*u+c*w-u*x-v*y-w*z))*(1-cos(angle))+L*y*cos(angle)+sqrt(L)*(c*u-a*w+w*x-u*z)*sin(angle))/L;
 +    result[2] = ((c*(u*u+v*v)-w*(a*u+b*v-u*x-v*y-w*z))*(1-cos(angle))+L*z*cos(angle)+sqrt(L)*(-b*u+a*v-v*x+u*y)*sin(angle))/L;
 +
 +    return result;
 +}
 +
 +var fold = function(event){
 +    var obj = foldables[Number(event.currentTarget.name)];
 +    var foldvalue = Number(event.currentTarget.value);
 +    var scale = foldvalue - obj.userData.oldscale;
 +
 +    for (var j=0; j<obj.userData.axes.length; j++) {
 +        rotateVertices(obj, j, scale);
 +    }
 +    update(obj);
 +    obj.userData.oldscale += scale;
 +    lookAtBarycenter(obj);
 +}
 +
 +function lookAtBarycenter(obj){
 +    control.target = barycenter(obj);
 +}
 +
 +function barycenter(obj) {
 +    var center = new THREE.Vector3(0,0,0);
 +    var points = obj.userData.points;
 +    for (var i=0; i<points.length; i++){
 +        center.add(points[i].vector);
 +    }
 +    center.divideScalar(points.length);
 +    return center;
 +}
 +
 +function rotateVertices(obj, edge, scale) {
 +    var axes = obj.userData.axes;
 +    var subtrees = obj.userData.subtrees;
 +    var points = obj.userData.points;
 +    var angles = obj.userData.angles;
 +    if (edge < axes.length){
 +        for (var j=0; j<subtrees[edge].length; j++){
 +            var rotP = rotate(points[subtrees[edge][j]].vector, points[axes[edge][0]].vector,points[axes[edge][1]].vector, scale * (Math.PI - angles[edge]));
 +            points[subtrees[edge][j]].set(rotP[0],rotP[1],rotP[2]);
 +        }
 +    }
 +}
 +
 +function update(obj) {
 +   updateFacesPosition(obj);
 +   updateEdgesPosition(obj);
 +}
 +
 +if (foldables.length) {
 +    var settings = document.getElementById('settings_2');
 +    var foldDiv = document.createElement('div');
 +    foldDiv.id = 'fold_2';
 +    var title = document.createElement('strong');
 +    title.innerHTML = 'Fold';
 +    foldDiv.appendChild(title);
 +    foldDiv.className = 'group';
 +    for (var i=0; i<foldables.length; i++) {
 +        var range = document.createElement('input');
 +        range.type = 'range';
 +        range.min = 0;
 +        range.max = 1;
 +        range.value = 0;
 +        range.step = 0.001;
 +        range.name = String(i);
 +        range.oninput = fold;
 +        foldDiv.appendChild(range);
 +    }
 +    lookAtBarycenter(foldables[0]);
 +    settings.insertBefore(foldDiv,settings.childNodes[0]);
 +}
 +
 +    
 +// ---------------------- EXPLOSION ------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +if (explodableModel) {
 +    for (var i=0; i<scene.children.length; i++) {
 +        obj = scene.children[i];
 +        if ( obj.userData.explodable ) {
 +            computeCentroid(obj);
 +        }
 +    }
 +    document.getElementById('explodeRange_2').oninput = triggerExplode;
 +    document.getElementById('explodeCheckbox_2').onchange = triggerAutomaticExplode;
 +    document.getElementById('explodingSpeedRange_2').oninput = setExplodingSpeed;
 +    explode(0.000001);
 +}
 +
 +function computeCentroid(obj) {
 +    centroid = new THREE.Vector3();
 +    obj.userData.points.forEach(function(pmpoint) {
 +        centroid.add(pmpoint.vector);
 +    });
 +    centroid.divideScalar(obj.userData.points.length);
 +    obj.userData.centroid = centroid;
 +}
 +
 +function explode(factor) {
 +    for (var i=0; i<scene.children.length; i++) {
 +        var obj = scene.children[i];
 +        if (obj.userData.hasOwnProperty("centroid")) { 
 +            var c = obj.userData.centroid;
 +            obj.position.set(c.x*factor, c.y*factor, c.z*factor);
 +        }
 +    }
 +}
 +
 +function triggerExplode(event){
 +    explodeScale = Number(event.currentTarget.value);
 +    explode(explodeScale);
 +}
 +
 +function setExplodingSpeed(event){
 +    explodingSpeed = Number(event.currentTarget.value);
 +}
 +
 +function triggerAutomaticExplode(event){
 +    if (event.currentTarget.checked){
 +        startExploding();
 +    } else {
 +        clearIntervals();
 +    }
 +}
 +
 +function startExploding(){
 +    intervals.push(setInterval(explodingInterval, 25));
 +}
 +
 +
 +function explodingInterval(){
 +    explodeScale += explodingSpeed;
 +    if (explodeScale <= 6){ 
 +        explode(explodeScale);
 +    }
 +    else{
 +        explode(6);
 +        explodeScale = 6;
 +        clearIntervals();
 +        timeouts.push(setTimeout(startUnexploding, 3000));
 +    }
 +    document.getElementById('explodeRange_2').value = explodeScale;
 +}
 +
 +
 +function startUnexploding(){
 +    intervals.push(setInterval(unexplodingInterval, 25));
 +}
 +
 +function unexplodingInterval(){
 +    explodeScale -= explodingSpeed;
 +    if (explodeScale >= 0){
 +        explode(explodeScale);
 +    }
 +    else {
 +        explode(0);
 +        explodeScale = 0;
 +        clearIntervals();
 +        timeouts.push(setTimeout(startExploding, 3000));
 +    }
 +    document.getElementById('explodeRange_2').value = explodeScale;
 +}
 +
 +function clearIntervals(){
 +    intervals.forEach(function(interval){
 +        clearInterval(interval);
 +    });
 +    intervals = [];
 +    timeouts.forEach(function(timeout){
 +        clearTimeout(timeout);
 +    });
 +    timeouts = [];
 +}
 +
 +// ---------------------- DISPLAY --------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +const objectTypeInnerHTMLs = { points: "Points", pointlabels: "Point labels", lines: "Edges", edgelabels: "Edge labels", faces: "Faces", arrowheads: "Arrow heads" };
 +const objectTypeVisible = {};
 +Object.assign(objectTypeVisible,modelContains);
 +const sortedObjectTypeKeys = Object.keys(objectTypeInnerHTMLs).sort();
 +const shownObjectTypesList = document.getElementById('shownObjectTypesList_2');
 +
 +function setVisibility(bool,objname) {
 +    for (var i=0; i<scene.children.length; i++){
 +        var obj = scene.children[i].getObjectByName(objname);
 +        if (obj) {
 +            obj.visible = bool;
 +        }
 +    }
 +}
 +
 +function toggleObjectTypeVisibility(event){
 +    var name = event.currentTarget.name;
 +    var checked = event.currentTarget.checked;
 +    objectTypeVisible[name] = checked;
 +    if (name == "faces") {
 +        setVisibility(checked,"frontfaces");
 +        setVisibility(checked,"backfaces");
 +    } else {
 +        setVisibility(checked,name);
 +    }
 +}
 +
 +for (var i=0; i<sortedObjectTypeKeys.length; i++){
 +    var key = sortedObjectTypeKeys[i];
 +    if (modelContains[key]) {
 +        var objTypeNode = document.createElement('span');
 +        objTypeNode.innerHTML = objectTypeInnerHTMLs[key] + '<br>';
 +        var checkbox = document.createElement('input');
 +        checkbox.type = 'checkbox';
 +        checkbox.checked = true;
 +        checkbox.name = key;
 +        checkbox.onchange = toggleObjectTypeVisibility;
 +        shownObjectTypesList.appendChild(checkbox);
 +        shownObjectTypesList.appendChild(objTypeNode);
 +    }
 +}
 +
 +// ------------------------------------------------------
 +
 +function toggleObjectVisibility(event){
 +    var nr = Number(event.currentTarget.name);
 +    scene.children[nr].visible = event.currentTarget.checked;
 +}
 +
 +// append checkboxes for displaying or hiding objects
 +var shownObjectsList = document.getElementById('shownObjectsList_2');
 +for (var i=0; i<scene.children.length; i++){
 +    obj = scene.children[i];
 +    var objNode = document.createElement('span');
 +    objNode.innerHTML = obj.name + '<br>';
 +    var checkbox = document.createElement('input');
 +    checkbox.type = 'checkbox';
 +    checkbox.checked = true;
 +    checkbox.name = String(i);
 +    checkbox.onchange = toggleObjectVisibility;
 +    shownObjectsList.appendChild(checkbox);
 +    shownObjectsList.appendChild(objNode);
 +}
 +
 +// ---------------------- SVG ------------------------------------------------------
 +// ---------------------------------------------------------------------------------
 +
 +function takeSvgScreenshot() {
 +    if (objectTypeVisible["pointlabels"]) {
 +        setVisibility(false,"pointlabels");
 +    }
 +    if (objectTypeVisible["edgelabels"]) {
 +        setVisibility(false,"edgelabels");
 +    }
 +    svgRenderer.render(scene,camera);
 +    svgElement = XMLS.serializeToString(svgRenderer.domElement);
 +    
 +    if (objectTypeVisible["pointlabels"]) {
 +        setVisibility(true,"pointlabels");
 +    }
 +    if (objectTypeVisible["edgelabels"]) {
 +        setVisibility(true,"edgelabels");
 +    }
 +
 +    if (document.getElementById('tab_2').checked){
 +        //show in new tab
 +        var myWindow = window.open("","");
 +        myWindow.document.body.innerHTML = svgElement;
 +    } else{
 +        // download svg file 
 +        download("screenshot.svg", svgElement);
 +    }
 +}
 +
 +function download(filename, text) {
 +    var element = document.createElement('a');
 +    element.setAttribute('href', 'data:text/plain;charset=utf-8,' + encodeURIComponent(text));
 +    element.setAttribute('download', filename);
 +
 +    element.style.display = 'none';
 +    document.body.appendChild(element);
 +
 +    element.click();
 +
 +    document.body.removeChild(element);
 +}
 +
 +
 +document.getElementById('transparencyRange_2').oninput = changeTransparency;
 +document.getElementById('depthWriteCheckbox_2').onchange = toggleDepthWrite;
 +document.getElementById('changeRotationX_2').onchange = changeRotationX;
 +document.getElementById('changeRotationY_2').onchange = changeRotationY;
 +document.getElementById('changeRotationZ_2').onchange = changeRotationZ;
 +document.getElementById('resetButton_2').onclick = resetScene;
 +document.getElementById('rotationSpeedRange_2').oninput = changeRotationSpeedFactor;
 +document.getElementById('takeScreenshot_2').onclick = takeSvgScreenshot;
 +document.getElementById('showSettingsButton_2').onclick = showSettings;
 +document.getElementById('hideSettingsButton_2').onclick = hideSettings;
 +
 +
 +// ------------------ SHORTCUTS --------------------------------------------
 +// -------------------------------------------------------------------------
 +
 +/**
 + * http://www.openjs.com/scripts/events/keyboard_shortcuts/
 + * Version : 2.01.B
 + * By Binny V A
 + * License : BSD
 + */
 +shortcut = {
 + 'all_shortcuts':{},//All the shortcuts are stored in this array
 + 'add': function(shortcut_combination,callback,opt) {
 + //Provide a set of default options
 + var default_options = {
 + 'type':'keydown',
 + 'propagate':false,
 + 'disable_in_input':false,
 + 'target':document,
 + 'keycode':false
 + }
 + if(!opt) opt = default_options;
 + else {
 + for(var dfo in default_options) {
 + if(typeof opt[dfo] == 'undefined') opt[dfo] = default_options[dfo];
 + }
 + }
 +
 + var ele = opt.target;
 + if(typeof opt.target == 'string') ele = document.getElementById(opt.target);
 + var ths = this;
 + shortcut_combination = shortcut_combination.toLowerCase();
 +
 + //The function to be called at keypress
 + var func = function(e) {
 + e = e || window.event;
 +
 + if(opt['disable_in_input']) { //Don't enable shortcut keys in Input, Textarea fields
 + var element;
 + if(e.target) element=e.target;
 + else if(e.srcElement) element=e.srcElement;
 + if(element.nodeType==3) element=element.parentNode;
 +
 + if(element.tagName == 'INPUT' || element.tagName == 'TEXTAREA') return;
 + }
 +
 + //Find Which key is pressed
 + if (e.keyCode) code = e.keyCode;
 + else if (e.which) code = e.which;
 + var character = String.fromCharCode(code).toLowerCase();
 +
 + if(code == 188) character=","; //If the user presses , when the type is onkeydown
 + if(code == 190) character="."; //If the user presses , when the type is onkeydown
 +
 + var keys = shortcut_combination.split("+");
 + //Key Pressed - counts the number of valid keypresses - if it is same as the number of keys, the shortcut function is invoked
 + var kp = 0;
 +
 + //Work around for stupid Shift key bug created by using lowercase - as a result the shift+num combination was broken
 + var shift_nums = {
 + "`":"~",
 + "1":"!",
 + "2":"@",
 + "3":"#",
 + "4":"$",
 + "5":"%",
 + "6":"^",
 + "7":"&",
 + "8":"*",
 + "9":"(",
 + "0":")",
 + "-":"_",
 + "=":"+",
 + ";":":",
 + "'":"\"",
 + ",":"<",
 + ".":">",
 + "/":"?",
 + "\\":"|"
 + }
 + //Special Keys - and their codes
 + var special_keys = {
 + 'esc':27,
 + 'escape':27,
 + 'tab':9,
 + 'space':32,
 + 'return':13,
 + 'enter':13,
 + 'backspace':8,
 +
 + 'scrolllock':145,
 + 'scroll_lock':145,
 + 'scroll':145,
 + 'capslock':20,
 + 'caps_lock':20,
 + 'caps':20,
 + 'numlock':144,
 + 'num_lock':144,
 + 'num':144,
 +
 + 'pause':19,
 + 'break':19,
 +
 + 'insert':45,
 + 'home':36,
 + 'delete':46,
 + 'end':35,
 +
 + 'pageup':33,
 + 'page_up':33,
 + 'pu':33,
 +
 + 'pagedown':34,
 + 'page_down':34,
 + 'pd':34,
 +
 + 'left':37,
 + 'up':38,
 + 'right':39,
 + 'down':40,
 +
 + 'f1':112,
 + 'f2':113,
 + 'f3':114,
 + 'f4':115,
 + 'f5':116,
 + 'f6':117,
 + 'f7':118,
 + 'f8':119,
 + 'f9':120,
 + 'f10':121,
 + 'f11':122,
 + 'f12':123
 + }
 +
 + var modifiers = { 
 + shift: { wanted:false, pressed:false},
 + ctrl : { wanted:false, pressed:false},
 + alt  : { wanted:false, pressed:false},
 + meta : { wanted:false, pressed:false} //Meta is Mac specific
 + };
 +                        
 + if(e.ctrlKey) modifiers.ctrl.pressed = true;
 + if(e.shiftKey) modifiers.shift.pressed = true;
 + if(e.altKey) modifiers.alt.pressed = true;
 + if(e.metaKey)   modifiers.meta.pressed = true;
 +                        
 + for(var i=0; k=keys[i],i<keys.length; i++) {
 + //Modifiers
 + if(k == 'ctrl' || k == 'control') {
 + kp++;
 + modifiers.ctrl.wanted = true;
 +
 + } else if(k == 'shift') {
 + kp++;
 + modifiers.shift.wanted = true;
 +
 + } else if(k == 'alt') {
 + kp++;
 + modifiers.alt.wanted = true;
 + } else if(k == 'meta') {
 + kp++;
 + modifiers.meta.wanted = true;
 + } else if(k.length > 1) { //If it is a special key
 + if(special_keys[k] == code) kp++;
 +
 + } else if(opt['keycode']) {
 + if(opt['keycode'] == code) kp++;
 +
 + } else { //The special keys did not match
 + if(character == k) kp++;
 + else {
 + if(shift_nums[character] && e.shiftKey) { //Stupid Shift key bug created by using lowercase
 + character = shift_nums[character]; 
 + if(character == k) kp++;
 + }
 + }
 + }
 + }
 +
 + if(kp == keys.length && 
 + modifiers.ctrl.pressed == modifiers.ctrl.wanted &&
 + modifiers.shift.pressed == modifiers.shift.wanted &&
 + modifiers.alt.pressed == modifiers.alt.wanted &&
 + modifiers.meta.pressed == modifiers.meta.wanted) {
 + callback(e);
 +
 + if(!opt['propagate']) { //Stop the event
 + //e.cancelBubble is supported by IE - this will kill the bubbling process.
 + e.cancelBubble = true;
 + e.returnValue = false;
 +
 + //e.stopPropagation works in Firefox.
 + if (e.stopPropagation) {
 + e.stopPropagation();
 + e.preventDefault();
 + }
 + return false;
 + }
 + }
 + }
 + this.all_shortcuts[shortcut_combination] = {
 + 'callback':func, 
 + 'target':ele, 
 + 'event': opt['type']
 + };
 + //Attach the function with the event
 + if(ele.addEventListener) ele.addEventListener(opt['type'], func, false);
 + else if(ele.attachEvent) ele.attachEvent('on'+opt['type'], func);
 + else ele['on'+opt['type']] = func;
 + },
 +
 + //Remove the shortcut - just specify the shortcut and I will remove the binding
 + 'remove':function(shortcut_combination) {
 + shortcut_combination = shortcut_combination.toLowerCase();
 + var binding = this.all_shortcuts[shortcut_combination];
 + delete(this.all_shortcuts[shortcut_combination])
 + if(!binding) return;
 + var type = binding['event'];
 + var ele = binding['target'];
 + var callback = binding['callback'];
 +
 + if(ele.detachEvent) ele.detachEvent('on'+type, callback);
 + else if(ele.removeEventListener) ele.removeEventListener(type, callback, false);
 + else ele['on'+type] = false;
 + }
 +}
 +
 +shortcut.add("Alt+Left",function() {
 + var event = new Event('click');
 + if (settingsShown){
 + document.getElementById('hideSettingsButton_2').dispatchEvent(event);
 + } else {
 + document.getElementById('showSettingsButton_2').dispatchEvent(event);
 + }
 +});
 +
 +
 +// COMMON_CODE_BLOCK_END
 +
 +});});
 +      </script>
 +   </body>
 +</html>
 +</HTML>
 +==== Multiple LP's for one Polytope ====
 +
 +You can add more than one linear program to one polytope.
 +
 +<code perl>
 +> $p->add("LP", LINEAR_OBJECTIVE => [0,1,1,-1]);
 +</code>
 +The different linear programs are saved in an array. The first one can still be accessed like befor:
 +
 +<code perl>
 +> print $p->LP->LINEAR_OBJECTIVE;
 +0 1 1 1
 +</code>
 +The other ones can be accessed like the entries of an array.
 +
 +<code perl>
 +> print $p->LP->[1]->LINEAR_OBJECTIVE;
 +0 1 1 -1
 +</code>
 +You can find more about this [[user_guide:howto:scripting#multiple_subobjects|here]].
 +