/****************************************************************************
 *               parse.cpp
 *
 * This module implements a parser for the scene description files.
 *
 * from Persistence of Vision(tm) Ray Tracer version 3.6.
 * Copyright 1991-2003 Persistence of Vision Team
 * Copyright 2003-2004 Persistence of Vision Raytracer Pty. Ltd.
 *---------------------------------------------------------------------------
 * NOTICE: This source code file is provided so that users may experiment
 * with enhancements to POV-Ray and to port the software to platforms other
 * than those supported by the POV-Ray developers. There are strict rules
 * regarding how you are permitted to use this file. These rules are contained
 * in the distribution and derivative versions licenses which should have been
 * provided with this file.
 *
 * These licences may be found online, linked from the end-user license
 * agreement that is located at http://www.povray.org/povlegal.html
 *---------------------------------------------------------------------------
 * This program is based on the popular DKB raytracer version 2.12.
 * DKBTrace was originally written by David K. Buck.
 * DKBTrace Ver 2.0-2.12 were written by David K. Buck & Aaron A. Collins.
 *---------------------------------------------------------------------------
 *
 *===========================================================================
 * This file is part of MegaPOV, a modified and unofficial version of POV-Ray
 * For more information on MegaPOV visit our website:
 * http://megapov.inetart.net/
 *===========================================================================
 *
 * $RCSfile: parse.cpp,v $
 * $Revision: 1.47 $
 * $Author: chris $
 *
 *****************************************************************************/

#include <ctype.h>
#include <math.h>
#include <algorithm>

#ifdef __BORLANDC__
  #include <mem.h>
#endif

#include "frame.h"
#include "vector.h"
#include "parse.h"
#include "parstxtr.h"
#include "parsestr.h"
#include "atmosph.h"
#include "bezier.h"   
#include "blob.h"     
#include "boxes.h"
#include "bsphere.h"
#include "colour.h"
#include "cones.h"    
#include "csg.h"
#include "discs.h"
#include "express.h"  
#include "fractal.h"
#include "hfield.h"
#include "image.h"    
#include "interior.h"    
#include "isosurf.h" 
#include "lathe.h"    
#include "lightgrp.h"
#include "matrices.h"
#include "mesh.h"
#include "normal.h"
#include "objects.h"
#include "octree.h"
#include "pigment.h"
#include "planes.h"
#include "polygon.h"
#include "polysolv.h"
#include "poly.h"
#include "povray.h"   
#include "prism.h"    
#include "quadrics.h" 
#include "radiosit.h"      
#include "render.h"   
#include "sor.h"      
#include "spheres.h"  
#include "sphsweep.h"   /* Sphere sweep support */
#include "splines.h"
#include "super.h"
#include "targa.h"    
#include "texture.h"  
#include "tokenize.h" 
#include "torus.h"
#include "triangle.h" 
#include "truetype.h" 
#include "photons.h"
#include "fpmetric.h"
#include "pov_util.h"
#include "povmsgid.h"

#ifdef GLOW_PATCH
  #include "glow.h"
#endif

#ifdef MECHANICS_SIMULATION_PATCH    /* Christoph Hormann August 2002 */
  #include "function.h"
  #include "mechsim.h"
#endif

#ifdef CLOTHRAY_PATCH
#include "clothray.h"
#endif

#ifdef POST_PROCESS_PATCH
  #include "fnpovfpu.h"
  #include "postproc.h"
#endif

#ifdef MOTION_BLUR_PATCH
#include "mblur.h"
#endif

BEGIN_POV_NAMESPACE

/*****************************************************************************
* Local preprocessor defines
******************************************************************************/

/* Volume that is considered to be infinite. [DB 9/94] */

const DBL INFINITE_VOLUME = BOUND_HUGE;


/*****************************************************************************
* Local typedefs
******************************************************************************/



/*****************************************************************************
* Global variables
******************************************************************************/

extern int tempfunc; // GLOBAL VARIABLE

extern int backtraceFlag; // GLOBAL VARIABLE

extern PHOTON_OPTIONS photonOptions; // GLOBAL VARIABLE

#ifdef MECHANICS_SIMULATION_PATCH            /* Christoph Hormann August 2002 */
  extern MECHSIM_OPTIONS MechsimOptions; // GLOBAL VARIABLE
#endif
/*****************************************************************************
* External functions
******************************************************************************/



/*****************************************************************************
* Local variables
******************************************************************************/

short Not_In_Default; // GLOBAL VARIABLE
short Ok_To_Declare; // GLOBAL VARIABLE
short LValue_Ok; // GLOBAL VARIABLE

static TOKEN *Brace_Stack // GLOBAL VARIABLE
#ifdef FIX_MEMORY_LEAKS_PATCH
=NULL
#endif
;

#if !defined( STORE_TOKEN_POSITION_START_IN_FILE_PATCH ) && !defined( MOTION_BLUR_PATCH )
  static
#endif
int Brace_Index; // GLOBAL VARIABLE
static bool Destroying_Frame = false ; // GLOBAL VARIABLE

CAMERA *Default_Camera; // GLOBAL VARIABLE

#ifdef MOTION_BLUR_PATCH
static bool GlobalSettingsForMotionBlurParsed=false; // GLOBAL VARIABLE
#endif

/*****************************************************************************
* Static functions
******************************************************************************/

void Frame_Init(void);
void Parse_Coeffs(int order, DBL *Coeffs);

static OBJECT *Parse_Bicubic_Patch (void);
static OBJECT *Parse_Blob (void);
static OBJECT *Parse_Box (void);
static OBJECT *Parse_Cone (void);
static OBJECT *Parse_CSG (int CSG_Type);
static OBJECT *Parse_Light_Group (void);
static OBJECT *Parse_Cylinder (void);
static OBJECT *Parse_Disc (void);
static OBJECT *Parse_Julia_Fractal (void);
static OBJECT *Parse_HField (void);
static OBJECT *Parse_Lathe (void);
static OBJECT *Parse_Light_Source (void);
static OBJECT *Parse_Object_Id (void);
static OBJECT *Parse_Plane (void);
static OBJECT *Parse_Poly (int order);
static OBJECT *Parse_Polygon (void);
static OBJECT *Parse_Prism (void);
static OBJECT *Parse_Quadric (void);
static OBJECT *Parse_Smooth_Triangle (void);
static OBJECT *Parse_Sor (void);
static OBJECT *Parse_Sphere (void);
static OBJECT *Parse_Superellipsoid (void);
static OBJECT *Parse_Torus (void);
static OBJECT *Parse_Triangle (void);
static OBJECT *Parse_Mesh (void);
static OBJECT *Parse_Mesh2 (void);
static TEXTURE *Parse_Mesh_Texture (TEXTURE **t2, TEXTURE **t3);
static OBJECT *Parse_TrueType (void);
static void Parse_Blob_Element_Mods (BLOB_ELEMENT *Element);

#ifdef MOTION_BLUR_PATCH
  static OBJECT *Parse_Motion_Blur ();
#endif

#ifndef PIGMENT_CAMERA_VIEW_PATCH
  static void Parse_Camera (CAMERA **Camera_Ptr);
#endif

static bool Parse_Camera_Mods(CAMERA *New);
static void Parse_Frame (void);

void Link (OBJECT *New_Object,OBJECT **Object_List_Root);
void Link_To_Frame (OBJECT *Object);
void Post_Process (OBJECT *Object, OBJECT *Parent);

static void Parse_Global_Settings (void);
static void Global_Setting_Warn (void);

static void Set_CSG_Children_Flag (OBJECT*, unsigned int, unsigned int, unsigned int);
static void Set_CSG_Tree_Flag (OBJECT*, unsigned int, int);

/* NK layers - 1999 July 10 - for backwards compatiblity with layered textures */
static void Convert_Filter_To_Transmit(PIGMENT *Pigment);

static OBJECT *Parse_Isosurface (void);
static OBJECT *Parse_Parametric(void);

static OBJECT *Parse_Sphere_Sweep (void);       /* Sphere sweep support */
#ifndef REMOVE_NOT_USED_CODE_WARNINGS_PATCH
  static OBJECT *Parse_Mesh2 (void);
#endif
static int Parse_Three_UVCoords(UV_VECT UV1, UV_VECT UV2, UV_VECT UV3);

#ifdef MECHANICS_SIMULATION_PATCH            /* Christoph Hormann August 2002 */
  void Parse_Mech_Sim_Mass (int Group);
  void Parse_Mech_Sim_Connection (int Group);
  void Parse_Mech_Sim_Face (int Group);
  void Parse_Mech_Sim_Viscoelastic (void);
  void Parse_Mech_Sim (void);
#endif

#ifdef POST_PROCESS_PATCH
  void Parse_Post_Process (void);
#endif

#ifdef ALTERNATIVE_VECTOR_FUNCTION
  FUNCTION_PTR Parse_User_Function (void);
  VECTOR_FUNCTION *Parse_Vector_Function (int Input,int Output);
#endif

/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

/* Parse the file. */
void Parse()
{
	POVMSObject decobj;
	int i, imax, err;
	int l;

	Init_Random_Generators();

	Initialize_Tokenizer();
	Brace_Stack = (TOKEN *)POV_MALLOC(MAX_BRACES*sizeof (TOKEN), "brace stack");
	Brace_Index = 0;

	Default_Camera = Create_Camera();
	Default_Texture = Create_Texture();
	Default_Texture->Pigment = Create_Pigment();
	Default_Texture->Tnormal = NULL;
	Default_Texture->Finish  = Create_Finish();

	Not_In_Default = true;
	Ok_To_Declare = true;
	LValue_Ok = false;

#ifdef MOTION_BLUR_PATCH
  GlobalSettingsForMotionBlurParsed = false;
#endif

	Frame_Init ();

	Stage = STAGE_PARSING;

	l = 0;
	err = POVMSAttrList_Count(&opts.Declared_Variables, &l);
	if((err == 0) && (l > 0))
	{
		imax = l;
		for(i = 1; i <= imax; i++)
		{
			POVMSAttribute item;

			err = POVMSAttrList_GetNth(&opts.Declared_Variables, i, &decobj);
			if(err == 0)
				err = POVMSObject_Get(&decobj, &item, kPOVAttrib_Identifier);
			if(err == 0)
			{
				l = 0;
				err = POVMSAttr_Size(&item, &l);
				if(l > 0)
				{
					SYM_ENTRY *Temp_Entry;
					float floatval = 0.0;
					char *bufptr = (char *)POV_MALLOC(l, "declare");

					bufptr[0] = 0;
					err = POVMSUtil_GetFloat(&decobj, kPOVAttrib_Value, &floatval);
					if(err == 0)
						err = POVMSAttr_Get(&item, kPOVMSType_CString, bufptr, &l);
					if(err == 0)
					{
						Temp_Entry = Add_Symbol(true, bufptr, FLOAT_ID_TOKEN);
						Temp_Entry->Data = Create_Float();
						*((DBL *)(Temp_Entry->Data)) = floatval;

						POV_FREE(bufptr);
					}
				}
				(void)POVMSAttr_Delete(&item);
			}
		}
	}

	IncludeHeader(opts.Header_File_Name);

	Parse_Frame ();

	Stage = STAGE_CLEANUP_PARSE;

	Post_Media(Frame.Atmosphere);

	if (Frame.Objects == NULL)
	{
		Warning (0, "No objects in scene.");
	}

	Terminate_Tokenizer();
	Destroy_Textures(Default_Texture); 
	Destroy_Camera(Default_Camera); 

	POV_FREE (Brace_Stack);

	Destroy_Random_Generators();

	Default_Texture = NULL;
	Default_Camera = NULL;
	Brace_Stack = NULL;
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

/* Set up the fields in the frame to default values. */
void Frame_Init()
{
   Destroying_Frame = false ;
   Frame.Camera = Copy_Camera(Default_Camera);
   Frame.Number_Of_Light_Sources = 0;  
   Frame.Light_Sources = NULL;
   Frame.Light_Group_Lights = NULL;
   Frame.Objects = NULL;
   Frame.Atmosphere_IOR = 1.0;
   Frame.Atmosphere_Dispersion = 1.0;
   Frame.Antialias_Threshold = opts.Antialias_Threshold;

   #ifdef GLOW_PATCH
     Frame.Glows = NULL;
     Frame.NumOfGlows = 0;
   #endif

/* dmf -- the first is physically "more correct".  The second works better */
/*   Make_Colour (Frame.Irid_Wavelengths, 0.70, 0.52, 0.48); */
   Make_Colour (Frame.Irid_Wavelengths, 0.25, 0.18, 0.14);
   Make_Colour (Frame.Background_Colour, 0.0, 0.0, 0.0);
   Make_Colour (Frame.Ambient_Light, 1.0, 1.0, 1.0);

   /* Init atmospheric stuff. [DB 12/94] */

   Frame.Atmosphere = NULL;

   Frame.Fog = NULL;

   Frame.Rainbow = NULL;

   Frame.Skysphere = NULL;
  }



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Destroy_Frame()
{
    FOG *Fog, *Next_Fog;
    RAINBOW *Rainbow, *Next_Rainbow;

    // This is necessary as a user who hits CANCEL during any IO performed
    // by this routine (e.g. Destroy_Object(), which can complain about
    // isosurface max_gradient), will cause this routine to be entered again
    // before the relevent data member has been set to NULL (this is able
    // to happen since cancel will invoke a longjmp on most platforms).
    // This causes the currently-executing segment to be destroyed twice,
    // which is a Bad Thing(tm). [CJC 11/01]
    if (Destroying_Frame)
      return ;
    Destroying_Frame = true ;

    Destroy_Camera (Frame.Camera);
    Frame.Camera=NULL;

    /* Destroy fogs. [DB 12/94] */

    for (Fog = Frame.Fog; Fog != NULL;)
    {
      Next_Fog = Fog->Next;

      Destroy_Fog(Fog);

      Fog = Next_Fog;
    }
    
    Frame.Fog = NULL;

    /* Destroy rainbows. [DB 12/94] */

    for (Rainbow = Frame.Rainbow; Rainbow != NULL;)
    {
      Next_Rainbow = Rainbow->Next;

      Destroy_Rainbow(Rainbow);

      Rainbow = Next_Rainbow;
    }

    Frame.Rainbow = NULL;

    /* Destroy skysphere. [DB 12/94] */

    Destroy_Skysphere(Frame.Skysphere);

    Frame.Skysphere = NULL;

    /* Destroy atmosphere. [DB 1/95] */

    Destroy_Media(Frame.Atmosphere);

    #ifdef GLOW_PATCH
      if(Frame.NumOfGlows > 0)
      {
        Destroy_Glow_List(Frame.Glows);
        Frame.Glows = NULL;
        Frame.NumOfGlows = 0;
      }
    #endif

    Frame.Atmosphere = NULL;

    if (Frame.Objects != NULL) {
       Destroy_Object (Frame.Objects);
       Frame.Objects = NULL;
       Frame.Light_Sources = NULL;
    }

    if(Frame.Light_Group_Lights != NULL)
    {
      LIGHT_GROUP_LIGHT* Light_Group_Light;
      LIGHT_GROUP_LIGHT* Next;
      for (Light_Group_Light = Frame.Light_Group_Lights;
           Light_Group_Light != NULL; )
      {
        Next = Light_Group_Light->Next;
        POV_FREE(Light_Group_Light);
        Light_Group_Light = Next;
      }

      Frame.Light_Group_Lights = NULL;
    }
	#ifdef FIX_MEMORY_LEAKS_PATCH
		if (   Brace_Stack )
			POV_FREE(Brace_Stack);
	#endif   
    
}


/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Begin ()
{
  char *front;
   
  if(++Brace_Index >= MAX_BRACES)
  {
     Warning(0, "Too many nested '{' braces.");
     Brace_Index--;
  }

  Brace_Stack[Brace_Index]=Token.Token_Id;
  
  Get_Token();
  
  if(Token.Token_Id == LEFT_CURLY_TOKEN)
  {
     return;
  }
     
  front = Get_Token_String(Brace_Stack[Brace_Index]);
  Found_Instead_Error("Missing { after", front);
}
  

/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_End ()
{
   char *front;

   Get_Token();

   if(Token.Token_Id == RIGHT_CURLY_TOKEN)
   {
      if(--Brace_Index < 0)
      {
         Warning(0, "Possible '}' brace missmatch.");
         Brace_Index = 0;
      }
      return;
   }

   front = Get_Token_String (Brace_Stack[Brace_Index]);
   Found_Instead_Error("No matching } in", front);
}

/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Bicubic_Patch ()
{
   BICUBIC_PATCH *Object;
   int i, j;

   Parse_Begin ();

   if ( (Object = (BICUBIC_PATCH *)Parse_Object_Id()) != NULL)
   {
      return ((OBJECT *) Object);
   }
   
   Object = Create_Bicubic_Patch();

   EXPECT
     CASE_FLOAT
       Warning(150, "Should use keywords for bicubic parameters.");
       Object->Patch_Type = (int)Parse_Float();
       if (Object->Patch_Type == 2 ||
           Object->Patch_Type == 3)
       {
          Object->Flatness_Value = Parse_Float();
       }
       else
       {
          Object->Flatness_Value = 0.1;
       }
       Object->U_Steps = (int)Parse_Float();
       Object->V_Steps = (int)Parse_Float();
       EXIT
     END_CASE
       
     CASE (TYPE_TOKEN)
       Object->Patch_Type = (int)Parse_Float();
     END_CASE

     CASE (FLATNESS_TOKEN)
       Object->Flatness_Value = Parse_Float();
     END_CASE

     CASE (V_STEPS_TOKEN)
       Object->V_Steps = (int)Parse_Float();
     END_CASE

     CASE (U_STEPS_TOKEN)
       Object->U_Steps = (int)Parse_Float();
     END_CASE

     CASE (ACCURACY_TOKEN)
       Object->accuracy = Parse_Float();
     END_CASE
     
     CASE(UV_VECTORS_TOKEN)
       /* Store 4 ST coords for quadrilateral  */
       Parse_UV_Vect(Object->ST[0]);  Parse_Comma();
       Parse_UV_Vect(Object->ST[1]);  Parse_Comma();
       Parse_UV_Vect(Object->ST[2]);  Parse_Comma();
       Parse_UV_Vect(Object->ST[3]); 

       EXIT
     END_CASE

     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   if (Object->Patch_Type > 1)
   {
      Object->Patch_Type = 1;
      Warning(0, "Patch type no longer supported. Using type 1.");
   }

   if ((Object->Patch_Type < 0) || (Object->Patch_Type > MAX_PATCH_TYPE))
   {
     Error("Undefined bicubic patch type.");
   }

   Parse_Comma();

   for (i=0;i<4;i++)
   {
     for (j=0;j<4;j++)
     {
       Parse_Vector(Object->Control_Points[i][j]);
	   if(!((i == 3) && (j == 3)))
	     Parse_Comma();
     }
   }
   
   Precompute_Patch_Values(Object); /* interpolated mesh coords */

   Compute_Bicubic_Patch_BBox(Object);  

   Parse_Object_Mods ((OBJECT *)Object);

   return ((OBJECT *) Object);
}


/*****************************************************************************
*
* FUNCTION
*
*   Parse_Blob
*
* INPUT
*
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*
*   Dieter Bayer
*   
* DESCRIPTION
*
*   -
*
* CHANGES
*
*   Jul 1994 : Creation.
*
******************************************************************************/

static OBJECT *Parse_Blob()
{
  int npoints;
  DBL threshold;
  VECTOR Axis, Base, Apex;
  BLOB *Object;
  BLOB_LIST *blob_components, *blob_component;

  Parse_Begin();

  if ((Object = (BLOB *)Parse_Object_Id()) != NULL)
  {
    return ((OBJECT *) Object);
  }

  Object = Create_Blob();

  blob_components = NULL;

  npoints = 0;

  threshold = 1.0;

  EXPECT
    CASE (THRESHOLD_TOKEN)
      threshold = Parse_Float();
    END_CASE

    /*************************************************************************
     * Read sperical component (old syntax).
     *************************************************************************/

    CASE (COMPONENT_TOKEN)
      blob_component = Create_Blob_List_Element();

      blob_component->elem.Type = BLOB_SPHERE;

      blob_component->elem.c[2] = Parse_Float();

      Parse_Comma();

      blob_component->elem.rad2 = Parse_Float();

      Parse_Comma();

      blob_component->elem.rad2 = Sqr(blob_component->elem.rad2);

      Parse_Vector(blob_component->elem.O);

      /* Next component. */

      blob_component->next = blob_components;

      blob_components = blob_component;

      npoints++;
    END_CASE

    /*************************************************************************
     * Read sperical component (new syntax).
     *************************************************************************/

    CASE (SPHERE_TOKEN)
      blob_component = Create_Blob_List_Element();

      blob_component->elem.Type = BLOB_SPHERE;

      Parse_Begin();

      Parse_Vector(blob_component->elem.O);

      Parse_Comma();

      blob_component->elem.rad2 = Parse_Float();

      blob_component->elem.rad2 = Sqr(blob_component->elem.rad2);

      Parse_Comma();

      ALLOW(STRENGTH_TOKEN)

      blob_component->elem.c[2] = Parse_Float();

      Parse_Blob_Element_Mods(&blob_component->elem);

      /* Next component. */

      blob_component->next = blob_components;

      blob_components = blob_component;

      npoints++;
    END_CASE

    /*************************************************************************
     * Read cylindrical component.
     *************************************************************************/

    CASE (CYLINDER_TOKEN)
      blob_component = Create_Blob_List_Element();

      blob_component->elem.Type = BLOB_CYLINDER;

      blob_component->elem.Trans = Create_Transform();

      Parse_Begin();

      Parse_Vector(Base);

      Parse_Comma();

      Parse_Vector(Apex);

      Parse_Comma();

      blob_component->elem.rad2 = Parse_Float();

      blob_component->elem.rad2 = Sqr(blob_component->elem.rad2);

      Parse_Comma();

      ALLOW(STRENGTH_TOKEN)

      blob_component->elem.c[2] = Parse_Float();

      /* Calculate cylinder's coordinate system. */

      VSub(Axis, Apex, Base);

      VLength(blob_component->elem.len, Axis);

      if (blob_component->elem.len < EPSILON)
      {
        Error("Degenerate cylindrical component in blob.");
      }

      VInverseScaleEq(Axis, blob_component->elem.len);

      Compute_Coordinate_Transform(blob_component->elem.Trans, Base, Axis, 1.0, 1.0);

      Parse_Blob_Element_Mods(&blob_component->elem);

      /* Next component. */

      blob_component->next = blob_components;

      blob_components = blob_component;

      npoints++;
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  Create_Blob_Element_Texture_List(Object, blob_components, npoints);

  Parse_Object_Mods((OBJECT *)Object);

  /* The blob's texture has to be processed before Make_Blob() is called. */

  Post_Textures(Object->Texture);
  
  /* Finally, process the information */

  Make_Blob(Object, threshold, blob_components, npoints);

  return((OBJECT *)Object);
}



/*****************************************************************************
*
* FUNCTION
*
*   Parse_Blob_Element_Mods
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*
*   Dieter Bayer
*   
* DESCRIPTION
*
*   -
*
* CHANGES
*
*   Sep 1994 : Creation.
*
******************************************************************************/

static void Parse_Blob_Element_Mods(BLOB_ELEMENT *Element)
{
  VECTOR Local_Vector;
  MATRIX Local_Matrix;
  TRANSFORM Local_Trans;
  TEXTURE *Local_Texture;

  EXPECT
    CASE (TRANSLATE_TOKEN)
      Parse_Vector (Local_Vector);
      Translate_Blob_Element (Element, Local_Vector);
    END_CASE

    CASE (ROTATE_TOKEN)
      Parse_Vector (Local_Vector);
      Rotate_Blob_Element (Element, Local_Vector);
    END_CASE

    CASE (SCALE_TOKEN)
      Parse_Scale_Vector (Local_Vector);
      Scale_Blob_Element (Element, Local_Vector);
    END_CASE

    CASE (TRANSFORM_TOKEN)
      Transform_Blob_Element (Element, Parse_Transform(&Local_Trans));
    END_CASE

    CASE (MATRIX_TOKEN)
      Parse_Matrix (Local_Matrix);
      Compute_Matrix_Transform(&Local_Trans, Local_Matrix);
      Transform_Blob_Element (Element, &Local_Trans);
    END_CASE

    CASE (TEXTURE_TOKEN)
      Parse_Begin ();
      Local_Texture = Parse_Texture();
      Parse_End ();
      Link_Textures(&Element->Texture, Local_Texture);
    END_CASE

    CASE3 (PIGMENT_TOKEN, TNORMAL_TOKEN, FINISH_TOKEN)
      if (Element->Texture == NULL)
      {
        Element->Texture = Copy_Textures(Default_Texture);
      }
      else
      {
        if (Element->Texture->Type != PLAIN_PATTERN)
        {
          Link_Textures(&Element->Texture, Copy_Textures(Default_Texture));
        }
      }
      UNGET
      EXPECT
        CASE (PIGMENT_TOKEN)
          Parse_Begin ();
          Parse_Pigment(&Element->Texture->Pigment);
          Parse_End ();
        END_CASE

        CASE (TNORMAL_TOKEN)
          Parse_Begin ();
          Parse_Tnormal(&Element->Texture->Tnormal);
          Parse_End ();
        END_CASE

        CASE (FINISH_TOKEN)
          Parse_Finish(&Element->Texture->Finish);
        END_CASE

        OTHERWISE
          UNGET
          EXIT
        END_CASE
      END_EXPECT
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  Parse_End();

  /* Postprocess to make sure that HAS_FILTER will be set correctly. */

  Post_Textures(Element->Texture);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Box ()
{
   BOX *Object;
   DBL temp;

   Parse_Begin ();

   if ( (Object = (BOX *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);

   Object = Create_Box();

        Parse_Vector(Object->bounds[0]);     Parse_Comma();
        Parse_Vector(Object->bounds[1]);

    if (Object->bounds[0][X] > Object->bounds[1][X]) {
       temp = Object->bounds[0][X];
       Object->bounds[0][X] = Object->bounds[1][X];
       Object->bounds[1][X] = temp;
       }
    if (Object->bounds[0][Y] > Object->bounds[1][Y]) {
       temp = Object->bounds[0][Y];
       Object->bounds[0][Y] = Object->bounds[1][Y];
       Object->bounds[1][Y] = temp;
       }
    if (Object->bounds[0][Z] > Object->bounds[1][Z]) {
       temp = Object->bounds[0][Z];
       Object->bounds[0][Z] = Object->bounds[1][Z];
       Object->bounds[1][Z] = temp;
       }

   Compute_Box_BBox(Object);

   Parse_Object_Mods ((OBJECT *)Object);

   return ((OBJECT *) Object);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

#ifndef PIGMENT_CAMERA_VIEW_PATCH
static void Parse_Camera (CAMERA **Camera_Ptr)
#else
void Parse_Camera (CAMERA **Camera_Ptr, bool Extended_Output)
#endif
{
	int i;
	DBL Direction_Length = 1.0, Up_Length, Right_Length, Handedness;
	DBL k1, k2, k3;
	VECTOR tempv;
	MATRIX Local_Matrix;
	TRANSFORM Local_Trans;
	CAMERA *New;
	bool only_mods = false;

	Parse_Begin ();

	EXPECT
		CASE (CAMERA_ID_TOKEN)
			Destroy_Camera(*Camera_Ptr);
			*Camera_Ptr = Copy_Camera ((CAMERA *) Token.Data);
			if (opts.Language_Version >= 350)
				only_mods = true;
			EXIT
		END_CASE

		OTHERWISE
			UNGET
			EXIT
		END_CASE
	END_EXPECT

	New = *Camera_Ptr;

	if ((opts.Language_Version >= 350) && (only_mods == true))
	{
		TRANSFORM Backup_Trans;

		// keep a copy and clear it because this is a copy of a camera
		// and this will prevent that transforms are applied twice [trf]
		Backup_Trans = *New->Trans;
		Destroy_Transform(New->Trans);
		New->Trans = Create_Transform();

		EXPECT
			CASE (TRANSLATE_TOKEN)
				Parse_Vector (tempv);
	            Compute_Translation_Transform(&Local_Trans, tempv);
	            Compose_Transforms(New->Trans, &Local_Trans);
			END_CASE

			CASE (ROTATE_TOKEN)
				Parse_Vector (tempv);
	            Compute_Rotation_Transform(&Local_Trans, tempv);
	            Compose_Transforms(New->Trans, &Local_Trans);
			END_CASE

			CASE (SCALE_TOKEN)
				Parse_Scale_Vector(tempv);
				Compute_Scaling_Transform(&Local_Trans, tempv);
	            Compose_Transforms(New->Trans, &Local_Trans);
			END_CASE

			CASE (TRANSFORM_TOKEN)
				Parse_Transform(&Local_Trans);
	            Compose_Transforms(New->Trans, &Local_Trans);
			END_CASE

			CASE (MATRIX_TOKEN)
				Parse_Matrix(Local_Matrix);
				Compute_Matrix_Transform(&Local_Trans, Local_Matrix);
				Compose_Transforms(New->Trans, &Local_Trans);
			END_CASE

			OTHERWISE
				UNGET
				EXIT
			END_CASE
		END_EXPECT

		// apply camera transformations
		Transform_Camera(New, New->Trans);
		Compose_Transforms(&Backup_Trans, New->Trans);
	}
	else if (opts.Language_Version >= 350)
	{

		/*
		 * The camera statement in version 3.5 is a tiny bit more restrictive
		 * than in previous versions (Note: Backward compatibility is available
		 * with the version switch!).  It will always apply camera modifiers in
		 * the same order, regardless of the order in which they appeared in the
		 * camera statement.  The order is as follows:
		 *
		 * right
		 * direction
		 * angle (depends on right, changes direction-length)
		 * up
		 * sky
		 * location
		 * look_at (depends on location, right, direction, up, sky, changes right, up, direction)
		 * focal_point (depends on location)
		 *
		 * VERIFY: Is there a need to modify look_at to consider angle, right, up and/or direction??? [trf]
		 * VERIFY: Is there a need to modify angle to consider direction??? [trf]
		 */

		bool had_angle = false, had_up = false, had_right = false;
		VECTOR old_look_at, old_up, old_right, old_focal_point;
		DBL old_angle;

		Assign_Vector(old_look_at, New->Look_At);
		Make_Vector(New->Look_At, HUGE_VAL, HUGE_VAL, HUGE_VAL);
		Assign_Vector(old_up, New->Up);
		Make_Vector(New->Up, HUGE_VAL, HUGE_VAL, HUGE_VAL);
		Assign_Vector(old_right, New->Right);
		Make_Vector(New->Right, HUGE_VAL, HUGE_VAL, HUGE_VAL);
		Assign_Vector(old_focal_point, New->Focal_Point);
		Make_Vector(New->Focal_Point, HUGE_VAL, HUGE_VAL, HUGE_VAL);
		old_angle = New->Angle;
		New->Angle = HUGE_VAL;

		EXPECT
			CASE (PERSPECTIVE_TOKEN)
				New->Type = PERSPECTIVE_CAMERA;
			END_CASE

			CASE (ORTHOGRAPHIC_TOKEN)
				New->Type = ORTHOGRAPHIC_CAMERA;
			END_CASE

			CASE (FISHEYE_TOKEN)
				New->Type = FISHEYE_CAMERA;
			END_CASE

			CASE (ULTRA_WIDE_ANGLE_TOKEN)
				New->Type = ULTRA_WIDE_ANGLE_CAMERA;
			END_CASE

			CASE (OMNIMAX_TOKEN)
				New->Type = OMNIMAX_CAMERA;
			END_CASE

			CASE (PANORAMIC_TOKEN)
				New->Type = PANORAMIC_CAMERA;
			END_CASE

			CASE (SPHERICAL_TOKEN)
				New->Type = SPHERICAL_CAMERA;
			END_CASE

			CASE (CYLINDER_TOKEN)
				i = (int)Parse_Float();
				switch (i)
				{
					case 1: New->Type = CYL_1_CAMERA; break;
					case 2: New->Type = CYL_2_CAMERA; break;
					case 3: New->Type = CYL_3_CAMERA; break;
					case 4: New->Type = CYL_4_CAMERA; break;
					default: Error("Invalid cylinder camera type, valid are types 1 to 4."); break;
				}       
			END_CASE

			#ifdef USER_DEFINED_CAMERA_PATCH
			      CASE (USER_DEFINED_TOKEN)
			        parseUnofficialFeature(110);
			        New->Type = USER_DEFINED_CAMERA;
			      END_CASE
			#endif

 			OTHERWISE
				UNGET
				EXIT
			END_CASE
		END_EXPECT

		switch(New->Type)
		{
			case PERSPECTIVE_CAMERA:
				EXPECT
					CASE (ANGLE_TOKEN)
						New->Angle = Parse_Float();
						if (New->Angle < 0.0)
							Error("Negative viewing angle.");
					END_CASE

					CASE5(ORTHOGRAPHIC_TOKEN, FISHEYE_TOKEN, ULTRA_WIDE_ANGLE_TOKEN, OMNIMAX_TOKEN, PANORAMIC_TOKEN)
					CASE2(SPHERICAL_TOKEN, CYLINDER_TOKEN)
						Expectation_Error("perspective camera modifier");
					END_CASE

					OTHERWISE
						UNGET
						if(Parse_Camera_Mods(New) == false)
							EXIT
					END_CASE
				END_EXPECT
				break;
			case ORTHOGRAPHIC_CAMERA:
				EXPECT
					CASE (ANGLE_TOKEN)
						New->Angle = Allow_Float(0.0);
						if (New->Angle < 0.0)
							Error("Negative viewing angle.");
					END_CASE

					CASE5(PERSPECTIVE_TOKEN, FISHEYE_TOKEN, ULTRA_WIDE_ANGLE_TOKEN, OMNIMAX_TOKEN, PANORAMIC_TOKEN)
					CASE2(SPHERICAL_TOKEN, CYLINDER_TOKEN)
						Expectation_Error("orthographic camera modifier");
					END_CASE

					OTHERWISE
						UNGET
						if(Parse_Camera_Mods(New) == false)
							EXIT
					END_CASE
				END_EXPECT
				break;
			case FISHEYE_CAMERA:
				EXPECT
					CASE (ANGLE_TOKEN)
						New->Angle = Parse_Float();
						if (New->Angle < 0.0)
							Error("Negative viewing angle.");
					END_CASE

					CASE5(PERSPECTIVE_TOKEN, ORTHOGRAPHIC_TOKEN, ULTRA_WIDE_ANGLE_TOKEN, OMNIMAX_TOKEN, PANORAMIC_TOKEN)
					CASE2(SPHERICAL_TOKEN, CYLINDER_TOKEN)
						Expectation_Error("fisheye camera modifier");
					END_CASE

					OTHERWISE
						UNGET
						if(Parse_Camera_Mods(New) == false)
							EXIT
					END_CASE
				END_EXPECT
				break;
			case ULTRA_WIDE_ANGLE_CAMERA:
				EXPECT
					CASE (ANGLE_TOKEN)
						New->Angle = Parse_Float();
						if (New->Angle < 0.0)
							Error("Negative viewing angle.");
					END_CASE

					CASE5(PERSPECTIVE_TOKEN, ORTHOGRAPHIC_TOKEN, FISHEYE_TOKEN, OMNIMAX_TOKEN, PANORAMIC_TOKEN)
					CASE2(SPHERICAL_TOKEN, CYLINDER_TOKEN)
						Expectation_Error("ultra_wide_angle camera modifier");
					END_CASE

					OTHERWISE
						UNGET
						if(Parse_Camera_Mods(New) == false)
							EXIT
					END_CASE
				END_EXPECT
				break;
			case OMNIMAX_CAMERA:
				EXPECT
					CASE (ANGLE_TOKEN)
						New->Angle = Parse_Float();
						if (New->Angle < 0.0)
							Error("Negative viewing angle.");
					END_CASE

					CASE5(PERSPECTIVE_TOKEN, ORTHOGRAPHIC_TOKEN, FISHEYE_TOKEN, ULTRA_WIDE_ANGLE_TOKEN, PANORAMIC_TOKEN)
					CASE2(SPHERICAL_TOKEN, CYLINDER_TOKEN)
						Expectation_Error("omnimax camera modifier");
					END_CASE

					OTHERWISE
						UNGET
						if(Parse_Camera_Mods(New) == false)
							EXIT
					END_CASE
				END_EXPECT
				break;
			case PANORAMIC_CAMERA:
				EXPECT
					CASE (ANGLE_TOKEN)
						New->Angle = Parse_Float();
						if (New->Angle < 0.0)
							Error("Negative viewing angle.");
					END_CASE

					CASE5(PERSPECTIVE_TOKEN, ORTHOGRAPHIC_TOKEN, FISHEYE_TOKEN, ULTRA_WIDE_ANGLE_TOKEN, OMNIMAX_TOKEN)
					CASE2(SPHERICAL_TOKEN, CYLINDER_TOKEN)
						Expectation_Error("panoramic camera modifier");
					END_CASE

					OTHERWISE
						UNGET
						if(Parse_Camera_Mods(New) == false)
							EXIT
					END_CASE
				END_EXPECT
				break;
			case CYL_1_CAMERA:
			case CYL_2_CAMERA:
			case CYL_3_CAMERA:
			case CYL_4_CAMERA:
				EXPECT
					CASE (ANGLE_TOKEN)
						New->Angle = Parse_Float();
						if (New->Angle < 0.0)
							Error("Negative viewing angle.");
					END_CASE

					CASE6(PERSPECTIVE_TOKEN, ORTHOGRAPHIC_TOKEN, FISHEYE_TOKEN, ULTRA_WIDE_ANGLE_TOKEN, OMNIMAX_TOKEN, PANORAMIC_TOKEN)
					CASE(SPHERICAL_TOKEN)
						Expectation_Error("cylinder camera modifier");
					END_CASE

					OTHERWISE
						UNGET
						if(Parse_Camera_Mods(New) == false)
							EXIT
					END_CASE
				END_EXPECT
				break;
			case SPHERICAL_CAMERA:
				EXPECT
					CASE (ANGLE_TOKEN)
						New->H_Angle = Parse_Float();
						if (New->H_Angle < 0.0)
							Error("Negative horizontal angle not allowed.");
						Parse_Comma();
						New->V_Angle = Allow_Float(New->H_Angle * 0.5);
						if (New->V_Angle < 0.0)
							Error("Negative vertical angle not allowed.");
					END_CASE

					CASE6(PERSPECTIVE_TOKEN, ORTHOGRAPHIC_TOKEN, FISHEYE_TOKEN, ULTRA_WIDE_ANGLE_TOKEN, OMNIMAX_TOKEN, PANORAMIC_TOKEN)
					CASE(CYLINDER_TOKEN)
						Expectation_Error("spherical camera modifier");
					END_CASE

					OTHERWISE
						UNGET
						if(Parse_Camera_Mods(New) == false)
							EXIT
					END_CASE
				END_EXPECT
				break;

				#ifdef USER_DEFINED_CAMERA_PATCH
				      case USER_DEFINED_CAMERA:
				        parseUnofficialFeature(110);
				        GET(LOCATION_TOKEN)
				        Parse_Begin ();
				        New->f_Location = Parse_Vector_Function(2,3);
				        Parse_End ();
				        GET(DIRECTION_TOKEN)
				        Parse_Begin ();
				        New->f_Direction = Parse_Vector_Function(2,3);
				        Parse_End ();
				        break;
				#endif
		}

		// handle "up"
		if (New->Up[X] == HUGE_VAL)
		{
			Assign_Vector(New->Up, old_up); // restore default up
		}
		else
			had_up = true;

		// handle "right"
		if (New->Right[X] == HUGE_VAL)
		{
			Assign_Vector(New->Right, old_right); // restore default right
		}
		else
			had_right = true;

		// apply "angle"
		if (New->Angle != HUGE_VAL)
		{
			if ((New->Type == PERSPECTIVE_CAMERA) || (New->Type == ORTHOGRAPHIC_CAMERA))
			{
				if (New->Angle >= 180.0)
					Error("Viewing angle has to be smaller than 180 degrees.");

				if (New->Angle > 0.0)
				{
					VNormalize(New->Direction, New->Direction);
					VLength (Right_Length, New->Right);
					Direction_Length = Right_Length / tan(New->Angle * M_PI_360)/2.0;
					VScaleEq(New->Direction, Direction_Length);
				}
			}

			had_angle = true;
		}
		else
			New->Angle = old_angle; // restore default angle

		// apply "look_at"
		if (New->Look_At[X] != HUGE_VAL)
		{
			VLength (Direction_Length, New->Direction);
			VLength (Up_Length,        New->Up);
			VLength (Right_Length,     New->Right);
			VCross  (tempv,            New->Up, New->Direction);
			VDot    (Handedness,       tempv,   New->Right);

			Assign_Vector(New->Direction, New->Look_At);

			VSub          (New->Direction, New->Direction, New->Location);

			// Check for zero length direction vector.
			if (VSumSqr(New->Direction) < EPSILON)
				Error("Camera location and look_at point must be different.");

			VNormalize (New->Direction, New->Direction);

			// Save Right vector
			Assign_Vector (tempv, New->Right);

			VCross        (New->Right, New->Sky, New->Direction);

			// Avoid DOMAIN error (from Terry Kanakis)
			if((fabs(New->Right[X]) < EPSILON) &&
			   (fabs(New->Right[Y]) < EPSILON) &&
			   (fabs(New->Right[Z]) < EPSILON))
			{
				// Restore Right vector
				Assign_Vector (New->Right, tempv);
			}

			VNormalize (New->Right,     New->Right);
			VCross     (New->Up,        New->Direction, New->Right);
			VScale     (New->Direction, New->Direction, Direction_Length);

			if (Handedness > 0.0)
			{
				VScaleEq (New->Right, Right_Length);
			}
			else
			{
				VScaleEq (New->Right, -Right_Length);
			}

			VScaleEq(New->Up, Up_Length);
		}
		else
			Assign_Vector(New->Look_At, old_look_at); // restore default look_at

		// apply "orthographic"
		if (New->Type == ORTHOGRAPHIC_CAMERA)
		{
			// only if neither up nor right have been specified
			// or if angle has been specified regardless if up or right have been specified
			if (((had_up == false) && (had_right == false)) || (had_angle == true))
			{
				// resize right and up vector to get the same image 
				// area as we get with the perspective camera
				VSub(tempv, New->Look_At, New->Location);
				VLength(k1, tempv);
				VLength(k2, New->Direction);
				if ((k1 > EPSILON) && (k2 > EPSILON))
				{
					VScaleEq(New->Right, k1 / k2);
					VScaleEq(New->Up, k1 / k2);
				}
			}
		}

		// apply "focal_point"
		if (New->Focal_Point[X] != HUGE_VAL)
		{
			Assign_Vector(tempv, New->Focal_Point);
			VSubEq(tempv, New->Location);
			VLength (New->Focal_Distance, tempv);
		}
		else
			Assign_Vector(New->Focal_Point, old_focal_point); // restore default focal_point

		// apply camera transformations
		Transform_Camera(New, New->Trans);
	}
	else // old style syntax
	{
		EXPECT
			CASE (PERSPECTIVE_TOKEN)
				New->Type = PERSPECTIVE_CAMERA;
			END_CASE

			CASE (ORTHOGRAPHIC_TOKEN)
				New->Type = ORTHOGRAPHIC_CAMERA;
				// resize right and up vector to get the same image 
				// area as we get with the perspective camera 
				VSub(tempv, New->Look_At, New->Location);
				VLength(k1, tempv);
				VLength(k2, New->Direction);
				if ((k1 > EPSILON) && (k2 > EPSILON))
				{
					VScaleEq(New->Right, k1 / k2);
					VScaleEq(New->Up, k1 / k2);
				}
			END_CASE

			CASE (FISHEYE_TOKEN)
				New->Type = FISHEYE_CAMERA;
			END_CASE

			CASE (ULTRA_WIDE_ANGLE_TOKEN)
				New->Type = ULTRA_WIDE_ANGLE_CAMERA;
			END_CASE

			CASE (OMNIMAX_TOKEN)
				New->Type = OMNIMAX_CAMERA;
			END_CASE

			CASE (PANORAMIC_TOKEN)
				New->Type = PANORAMIC_CAMERA;
			END_CASE

			CASE (CYLINDER_TOKEN)
				i = (int)Parse_Float();
				switch (i)
				{
					case 1: New->Type = CYL_1_CAMERA; break;
					case 2: New->Type = CYL_2_CAMERA; break;
					case 3: New->Type = CYL_3_CAMERA; break;
					case 4: New->Type = CYL_4_CAMERA; break;
				}       
			END_CASE

			CASE (ANGLE_TOKEN)
				New->Angle = Parse_Float();

				if (New->Angle < 0.0)
					Error("Negative viewing angle.");

				if (New->Type == PERSPECTIVE_CAMERA)
				{
					if (New->Angle >= 180.0)
						Error("Viewing angle has to be smaller than 180 degrees.");

					VNormalize(New->Direction, New->Direction);
					VLength (Right_Length, New->Right);
					Direction_Length = Right_Length / tan(New->Angle * M_PI_360)/2.0;
					VScaleEq(New->Direction, Direction_Length);
				}
			END_CASE

			CASE (TNORMAL_TOKEN)
				Parse_Begin ();
				Parse_Tnormal(&(New->Tnormal));
				Parse_End ();
			END_CASE

			CASE (LOCATION_TOKEN)
				Parse_Vector(New->Location);
			END_CASE

			CASE (DIRECTION_TOKEN)
				Parse_Vector(New->Direction);
			END_CASE

			CASE (UP_TOKEN)
				Parse_Vector(New->Up);
			END_CASE

			CASE (RIGHT_TOKEN)
				Parse_Vector(New->Right);
			END_CASE

			CASE (SKY_TOKEN)
				Parse_Vector(New->Sky);
			END_CASE

			CASE (LOOK_AT_TOKEN)
				VLength (Direction_Length, New->Direction);
				VLength (Up_Length,        New->Up);
				VLength (Right_Length,     New->Right);
				VCross  (tempv,            New->Up, New->Direction);
				VDot    (Handedness,       tempv,   New->Right);

				Parse_Vector (New->Direction);
				Assign_Vector(New->Look_At, New->Direction);

				VSub          (New->Direction, New->Direction, New->Location);

				// Check for zero length direction vector.
				if (VSumSqr(New->Direction) < EPSILON)
					Error("Camera location and look_at point must be different.");

				VNormalize (New->Direction, New->Direction);

				// Save Right vector
				Assign_Vector (tempv, New->Right);

				VCross        (New->Right, New->Sky, New->Direction);

				// Avoid DOMAIN error (from Terry Kanakis)
				if((fabs(New->Right[X]) < EPSILON) &&
				   (fabs(New->Right[Y]) < EPSILON) &&
				   (fabs(New->Right[Z]) < EPSILON))
				{
					// Restore Right vector
					Assign_Vector (New->Right, tempv);
				}

				VNormalize (New->Right,     New->Right);
				VCross     (New->Up,        New->Direction, New->Right);
				VScale     (New->Direction, New->Direction, Direction_Length);

				if (Handedness > 0.0)
				{
					VScaleEq (New->Right, Right_Length);
				}
				else
				{
					VScaleEq (New->Right, -Right_Length);
				}

				VScaleEq(New->Up, Up_Length);
			END_CASE

			CASE (TRANSLATE_TOKEN)
				Parse_Vector (tempv);
				Translate_Camera (New, tempv);
			END_CASE

			CASE (ROTATE_TOKEN)
				Parse_Vector (tempv);
				Rotate_Camera (New, tempv);
			END_CASE

			CASE (SCALE_TOKEN)
				Parse_Scale_Vector (tempv);
				Scale_Camera (New, tempv);
			END_CASE

			CASE (TRANSFORM_TOKEN)
				Transform_Camera (New, Parse_Transform(&Local_Trans));
			END_CASE

			CASE (MATRIX_TOKEN)
				Parse_Matrix (Local_Matrix);
				Compute_Matrix_Transform(&Local_Trans, Local_Matrix);
				Transform_Camera (New, &Local_Trans);
			END_CASE

			CASE (BLUR_SAMPLES_TOKEN)
				New->Blur_Samples = Parse_Float();
				if (New->Blur_Samples <= 0)
					Error("Illegal number of focal blur samples.");
			END_CASE

			CASE (CONFIDENCE_TOKEN)
				k1 = Parse_Float();
				if ((k1 > 0.0) && (k1 < 1.0))
					New->Confidence = k1;
				else
					Warning(0, "Illegal confidence value. Default is used.");
			END_CASE

			CASE (VARIANCE_TOKEN)
				k1 = Parse_Float();
				if ((k1 >= 0.0) && (k1 <= 1.0))
					New->Variance = k1;
				else
					Warning(0, "Illegal variance value. Default is used.");
			END_CASE

			CASE (APERTURE_TOKEN)
				New->Aperture = Parse_Float();
			END_CASE

			CASE (FOCAL_POINT_TOKEN)
				Parse_Vector(New->Focal_Point);
				Assign_Vector(tempv, New->Focal_Point);
				VSubEq(tempv, New->Location);
				VLength (New->Focal_Distance, tempv);
			END_CASE

			OTHERWISE
				UNGET
				EXIT
			END_CASE
		END_EXPECT
	}

  #ifdef PIGMENT_CAMERA_VIEW_PATCH
    EXPECT
      CASE (OUTPUT_TOKEN)
        parseUnofficialFeature(110);
        if(Extended_Output)
        {
          New->Output_Type = (int)Parse_Float();
          if ( !( ( New->Output_Type >= COLOR_OUTPUT ) && (New->Output_Type < LAST_OUTPUT ) ) )
          {
            Error( "Output type out of range" );
          }
        }
        else
        {
          Error( "Output type not allowed here" );
        }
        EXIT
      END_CASE
      OTHERWISE
        UNGET
        EXIT
      END_CASE
    END_EXPECT
  #endif

	Parse_End ();

	// Make sure the focal distance hasn't been explicitly given
	if (New->Focal_Distance < 0.0)
		New->Focal_Distance = Direction_Length;
	if (New->Focal_Distance == 0.0)
		New->Focal_Distance = 1.0;

	// Print a warning message if vectors are not perpendicular. [DB 10/94]
	VDot(k1, New->Right, New->Up);
	VDot(k2, New->Right, New->Direction);
	VDot(k3, New->Up, New->Direction);

	if ((fabs(k1) > EPSILON) || (fabs(k2) > EPSILON) || (fabs(k3) > EPSILON))
	{
		Warning(0, "Camera vectors are not perpendicular.\n"
		           "Making look_at the last statement may help.");
	}
}


bool Parse_Camera_Mods(CAMERA *New)
{
	TRANSFORM Local_Trans;
	MATRIX Local_Matrix;
	VECTOR tempv;
	DBL k1;

	EXPECT_ONE
		CASE (TRANSLATE_TOKEN)
			Parse_Vector (tempv);
            Compute_Translation_Transform(&Local_Trans, tempv);
            Compose_Transforms(New->Trans, &Local_Trans);
		END_CASE

		CASE (ROTATE_TOKEN)
			Parse_Vector (tempv);
            Compute_Rotation_Transform(&Local_Trans, tempv);
            Compose_Transforms(New->Trans, &Local_Trans);
		END_CASE

		CASE (SCALE_TOKEN)
			Parse_Scale_Vector(tempv);
			Compute_Scaling_Transform(&Local_Trans, tempv);
            Compose_Transforms(New->Trans, &Local_Trans);
		END_CASE

		CASE (TRANSFORM_TOKEN)
			Parse_Transform(&Local_Trans);
            Compose_Transforms(New->Trans, &Local_Trans);
		END_CASE

		CASE (MATRIX_TOKEN)
			Parse_Matrix(Local_Matrix);
			Compute_Matrix_Transform(&Local_Trans, Local_Matrix);
			Compose_Transforms(New->Trans, &Local_Trans);
		END_CASE

		CASE (TNORMAL_TOKEN)
			Parse_Begin();
			Parse_Tnormal(&(New->Tnormal));
			Parse_End();
		END_CASE

		CASE (LOOK_AT_TOKEN)
			Parse_Vector(New->Look_At);
		END_CASE

		CASE (LOCATION_TOKEN)
			Parse_Vector(New->Location);
		END_CASE

		CASE (DIRECTION_TOKEN)
			Parse_Vector(New->Direction);
		END_CASE

		CASE (UP_TOKEN)
			Parse_Vector(New->Up);
		END_CASE

		CASE (RIGHT_TOKEN)
			Parse_Vector(New->Right);
		END_CASE

		CASE (SKY_TOKEN)
			Parse_Vector(New->Sky);
		END_CASE

	    CASE (BLUR_SAMPLES_TOKEN)
	      #ifdef AVOID_TYPE_CONVERSION_WARNINGS_PATCH
	        New->Blur_Samples = (int)Parse_Float();
	      #else
	        New->Blur_Samples = Parse_Float();
	      #endif
	      if (New->Blur_Samples <= 0)
	        Error("Illegal number of focal blur samples.");
	    END_CASE

		CASE (CONFIDENCE_TOKEN)
			k1 = Parse_Float();
			if ((k1 > 0.0) && (k1 < 1.0))
				New->Confidence = k1;
			else
				Warning(0, "Illegal confidence value. Default is used.");
		END_CASE

		CASE (VARIANCE_TOKEN)
			k1 = Parse_Float();
			if ((k1 >= 0.0) && (k1 <= 1.0))
				New->Variance = k1;
			else
				Warning(0, "Illegal variance value. Default is used.");
		END_CASE

		CASE (APERTURE_TOKEN)
			New->Aperture = Parse_Float();
		END_CASE

		CASE (FOCAL_POINT_TOKEN)
			Parse_Vector(New->Focal_Point);
		END_CASE

		OTHERWISE
			UNGET
			return false;
		END_CASE
	END_EXPECT

	return true;
}


/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_CSG (int CSG_Type)
{
   CSG *Object;
   OBJECT *Local;
   int Object_Count = 0;
   int Light_Source_Union = true;

   Parse_Begin ();

   if ( (Object = (CSG *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);

   if (CSG_Type & CSG_UNION_TYPE)
     Object = Create_CSG_Union ();
   else
     if (CSG_Type & CSG_MERGE_TYPE)
       Object = Create_CSG_Merge ();
     else
       Object = Create_CSG_Intersection ();

   Object->Children = NULL;

   while ((Local = Parse_Object ()) != NULL)
     {
      if ((CSG_Type & CSG_INTERSECTION_TYPE) && (Local->Type & PATCH_OBJECT))
      { // Do not issue this warning in case it's a mesh with Inside_Vector
        if (Local->Methods != &Mesh_Methods || !((MESH*)Local)->has_inside_vector)
          Warning(0, "Patch objects not allowed in intersection.");
      }
      Object_Count++;
      if ((CSG_Type & CSG_DIFFERENCE_TYPE) && (Object_Count > 1))
        Invert_Object (Local);
      Object->Type |=  (Local->Type & CHILDREN_FLAGS);
      if (!(Local->Type & LIGHT_SOURCE_OBJECT))
         Light_Source_Union = false;
      Local->Type |= IS_CHILD_OBJECT;
      Link(Local, &Object->Children);
   }

   if (Light_Source_Union)
   {
     Object->Type |= LT_SRC_UNION_OBJECT;
   }
   
   if ((Object_Count < 2) && (opts.Language_Version >= 150))
     Warning(150, "Should have at least 2 objects in csg.");

   Compute_CSG_BBox((OBJECT *)Object);

   Parse_Object_Mods ((OBJECT *)Object);

   if (CSG_Type & CSG_DIFFERENCE_TYPE)
   {
     Object->Type |= CSG_DIFFERENCE_OBJECT;
   }

   return ((OBJECT *) Object);
  }




/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Cone ()
{
   CONE *Object;

   Parse_Begin ();

   if ( (Object = (CONE *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);

   Object = Create_Cone();

        Parse_Vector(Object->apex);  Parse_Comma ();
        Object->apex_radius = Parse_Float();  Parse_Comma ();

        Parse_Vector(Object->base);  Parse_Comma ();
   Object->base_radius = Parse_Float();

   EXPECT
     CASE(OPEN_TOKEN)
       Clear_Flag(Object, CLOSED_FLAG);
       EXIT
     END_CASE
     
     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   /* Compute run-time values for the cone */
   Compute_Cone_Data((OBJECT *)Object);

   Compute_Cone_BBox(Object);  

   Parse_Object_Mods ((OBJECT *)Object);

   return ((OBJECT *) Object);
}


/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Cylinder ()
{
   CONE *Object;

   Parse_Begin ();

   if ( (Object = (CONE *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);

   Object = Create_Cylinder();

        Parse_Vector(Object->apex);  Parse_Comma ();
        Parse_Vector(Object->base);  Parse_Comma ();
   Object->apex_radius = Parse_Float();
   Object->base_radius = Object->apex_radius;

   EXPECT
     CASE(OPEN_TOKEN)
       Clear_Flag(Object, CLOSED_FLAG);
       EXIT
     END_CASE

     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   Compute_Cylinder_Data((OBJECT *)Object);

   Compute_Cone_BBox(Object);

   Parse_Object_Mods ((OBJECT *)Object);

   return ((OBJECT *) Object);
  }


/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Disc ()
{
   DISC *Object;
   DBL tmpf;

   Parse_Begin ();

   if ( (Object = (DISC *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);

   Object = Create_Disc();

        Parse_Vector(Object->center); Parse_Comma ();
        Parse_Vector(Object->normal); Parse_Comma ();
   VNormalize(Object->normal, Object->normal);

   tmpf = Parse_Float(); Parse_Comma ();
   Object->oradius2 = tmpf * tmpf;

   EXPECT
     CASE_FLOAT
       tmpf = Parse_Float();
       Object->iradius2 = tmpf * tmpf;
     END_CASE

     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   /* Calculate info needed for ray-disc intersections */
   VDot(tmpf, Object->center, Object->normal);
   Object->d = -tmpf;

   Compute_Disc(Object);  

   Parse_Object_Mods ((OBJECT *)Object);

   return ((OBJECT *) Object);
  }



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_HField ()
{
  VECTOR Local_Vector;
  DBL Temp_Water_Level;
  HFIELD *Object;
  IMAGE *Image;

  Parse_Begin ();

  if ( (Object = (HFIELD *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);

   Object = Create_HField();

   Image = Parse_Image (HF_FILE);
   Image->Use_Colour_Flag = false;

   Make_Vector(Object->bounding_corner1, 0.0, 0.0, 0.0);

   if (Image->File_Type == POT_FILE)
   {
     Object->bounding_corner2[X] = Image->width/2.0 - 1.0;
   }
   else
   {
     Object->bounding_corner2[X] = Image->width - 1.0;
   }

   Object->bounding_corner2[Y] = 65536.0;
   Object->bounding_corner2[Z] = Image->height - 1.0;

   Make_Vector(Local_Vector,
     1.0 / (Object->bounding_corner2[X]),
     1.0 / (Object->bounding_corner2[Y]),
     1.0 / (Object->bounding_corner2[Z]));

   Compute_Scaling_Transform(Object->Trans, Local_Vector);

   EXPECT
     CASE (WATER_LEVEL_TOKEN)
       Temp_Water_Level = Parse_Float();
       if (opts.Language_Version < 200)
         Temp_Water_Level /=256.0;
       ((HFIELD *) Object)->bounding_corner1[Y] = 65536.0 * Temp_Water_Level;
     END_CASE

     CASE (SMOOTH_TOKEN)
       Set_Flag(Object, SMOOTHED_FLAG);
     END_CASE

     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   Parse_Object_Mods ((OBJECT *)Object);

   Compute_HField(Object, Image);

   Compute_HField_BBox(Object);

   Destroy_Image (Image);

   return ((OBJECT *) Object);
  }



/*****************************************************************************
*
* FUNCTION  	Parse_Isosurface
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR			  R. Suzuki
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Isosurface()
{
	ISOSURFACE *Object;
	DBL temp;

	Parse_Begin();

	if ((Object = (ISOSURFACE *)Parse_Object_Id()) != NULL)
		return ((OBJECT *)Object);
      
	Object = Create_IsoSurface();

	Get_Token();
	if(Token.Token_Id != FUNCTION_TOKEN)
		Parse_Error(FUNCTION_TOKEN);

	Object->Function = Parse_Function();

	EXPECT
		CASE(CONTAINED_BY_TOKEN)
			Parse_Begin();
			{ 
				int Exit_Flag2 = false; 

				while (!Exit_Flag2) 
				{
					Get_Token();  
					switch(Token.Token_Id) 
					{
						CASE(BOX_TOKEN)
							Object->container_shape = 0;

							Parse_Begin();

							Parse_Vector(Object->container.box.corner1);
							Parse_Comma();
							Parse_Vector(Object->container.box.corner2);

							Parse_End();

							if (Object->container.box.corner1[X] > Object->container.box.corner2[X])
							{
								temp = Object->container.box.corner1[X];
								Object->container.box.corner1[X] = Object->container.box.corner2[X];
								Object->container.box.corner2[X] = temp;
							}
							if (Object->container.box.corner1[Y] > Object->container.box.corner2[Y])
							{
								temp = Object->container.box.corner1[Y];
								Object->container.box.corner1[Y] = Object->container.box.corner2[Y];
								Object->container.box.corner2[Y] = temp;
							}
							if (Object->container.box.corner1[Z] > Object->container.box.corner2[Z])
							{
								temp = Object->container.box.corner1[Z];
								Object->container.box.corner1[Z] = Object->container.box.corner2[Z];
								Object->container.box.corner2[Z] = temp;
							}

							if (Object->Trans != NULL)
								Compute_IsoSurface_BBox(Object);

							Exit_Flag2 = true;
						END_CASE

						CASE(SPHERE_TOKEN)
							Object->container_shape = 1; 

							Parse_Begin();

							Parse_Vector(Object->container.sphere.center);
							Parse_Comma();
							Object->container.sphere.radius = Parse_Float();

							Parse_End();

							Make_BBox(Object->BBox,
							          Object->container.sphere.center[X] - Object->container.sphere.radius,
							          Object->container.sphere.center[Y] - Object->container.sphere.radius,
							          Object->container.sphere.center[Z] - Object->container.sphere.radius,
							          2.0 * Object->container.sphere.radius,
							          2.0 * Object->container.sphere.radius,
							          2.0 * Object->container.sphere.radius);

							if (Object->Trans != NULL)
								Compute_IsoSurface_BBox(Object);

							Exit_Flag2 = true;
						END_CASE

						OTHERWISE
							UNGET
							Exit_Flag2 = true;
		 			}
		 		}
		 	}
			Parse_End();
		END_CASE

		CASE(THRESHOLD_TOKEN)
			Object->threshold = Parse_Float();
		END_CASE

		CASE(ACCURACY_TOKEN)
			Object->accuracy = Parse_Float();
		END_CASE

		CASE(MAX_GRADIENT_TOKEN)
			Object->max_gradient = Parse_Float();
		END_CASE

  #ifdef MESSAGE_PATCH
    CASE(MESSAGE_TOKEN)
      parseUnofficialFeature(100);
      Object->mginfo->message=(unsigned char)Parse_Float(); //on off
    END_CASE

  #endif

		CASE(MAX_TRACE_TOKEN)
			Object->max_trace = (short)Parse_Float();
		END_CASE

		CASE(EVALUATE_TOKEN)
			Object->eval = true;
			Object->eval_param[0] = Parse_Float();
			Parse_Comma();
			Object->eval_param[1] = Parse_Float();
			Parse_Comma();
			Object->eval_param[2] = Parse_Float();
		END_CASE

		CASE(OPEN_TOKEN)
			Object->closed = false;
		END_CASE

		CASE(ALL_INTERSECTIONS_TOKEN)
			Object->max_trace = ISOSURFACE_MAXTRACE;
		END_CASE

		OTHERWISE
			UNGET
			EXIT
		END_CASE
	END_EXPECT

	if (Object->accuracy <= 0.0)
	{
		Warning(0, "Isosurface 'accuracy' is not positive. Using 0.001 (default).");
		Object->accuracy = 0.001;
	}
	if (Object->max_gradient <= 0.0)
	{
		Warning(0, "Isosurface 'max_gradient' is not positive. Using 1.1 (default).");
		Object->max_gradient = 1.1;
	}
	if (Object->max_trace > ISOSURFACE_MAXTRACE)
	{
		Warning(0, "Isosurface 'max_trace' exceeds maximum of %d. Using maximum.", (int)ISOSURFACE_MAXTRACE);
		Object->max_trace = ISOSURFACE_MAXTRACE;
	}
	if (Object->max_trace < 1)
	{
		Warning(0, "Isosurface 'max_trace' is not positive. Using 1 (default).");
		Object->max_trace = 1;
	}

	Parse_Object_Mods ((OBJECT *)Object);

	return ((OBJECT *)Object);
}


/*****************************************************************************
*
* FUNCTION
*
*   Parse_Julia_Fractal
*
* INPUT None
*   
* OUTPUT Fractal Objecstructure filledt
*   
* RETURNS 
*
*   OBJECT * -
*   
* AUTHOR
*
*   Pascal Massimino
*   
* DESCRIPTION
*
*   -
*
* CHANGES
*
*   Dec 1994 : Adopted to version 3.0. [DB]
*   Sept 1995 : Total rewrite for new syntax [TW]
*
******************************************************************************/

static OBJECT *Parse_Julia_Fractal ()
{
  FRACTAL *Object;
  DBL P;

  Parse_Begin();

  if ( (Object = (FRACTAL *)Parse_Object_Id()) != NULL)
    return((OBJECT *)Object);

  Object = Create_Fractal();

  Parse_Vector4D(Object->Julia_Parm); 

  EXPECT

    CASE(MAX_ITERATION_TOKEN)
      Object->n = (int)floor(Parse_Float()); 

      if (Object->n <= 0)
      {
        Object->n = 1;
      }
    END_CASE

    CASE(SLICE_TOKEN)
      Parse_Vector4D(Object->Slice);
      Parse_Comma();
      Object->SliceDist = Parse_Float(); 

      /* normalize slice vector */
      V4D_Dot(P,Object->Slice, Object->Slice);
      if (fabs(P) < EPSILON)
      {
        Error("Slice vector is zero.");
      }
      if (fabs(Object->Slice[T]) < EPSILON)
      {
        Error("Slice t component is zero.");
      }
      P = sqrt(P);
      V4D_InverseScaleEq(Object->Slice, P);      

    END_CASE

    CASE(PRECISION_TOKEN)
      P = Parse_Float(); 
      if ( P < 1.0 )
      {
        P = 1.0;
      }
      Object->Precision = 1.0 / P;
    END_CASE
      
    CASE(FLOAT_FUNCT_TOKEN)
      switch(Token.Function_Id)
      {
        case EXP_TOKEN:
          Object->Sub_Type = EXP_STYPE;
          break;
        case LN_TOKEN:
          Object->Sub_Type = LN_STYPE;
          break;
        case SIN_TOKEN:
          Object->Sub_Type = SIN_STYPE;
          break;
        case ASIN_TOKEN:
          Object->Sub_Type = ASIN_STYPE;
          break;
        case COS_TOKEN:
          Object->Sub_Type = COS_STYPE;
          break;
        case ACOS_TOKEN:
          Object->Sub_Type = ACOS_STYPE;
          break;
        case TAN_TOKEN:
          Object->Sub_Type = TAN_STYPE;
          break;
        case ATAN_TOKEN:
          Object->Sub_Type = ATAN_STYPE;
          break;
        case COSH_TOKEN:
          Object->Sub_Type = COSH_STYPE;
          break;
        case SINH_TOKEN:
          Object->Sub_Type = SINH_STYPE;
          break;
        case TANH_TOKEN:
          Object->Sub_Type = TANH_STYPE;
          break;
        case ATANH_TOKEN:
          Object->Sub_Type = ATANH_STYPE;
          break;
        case ACOSH_TOKEN:
          Object->Sub_Type = ACOSH_STYPE;
          break;
        case ASINH_TOKEN:
          Object->Sub_Type = ASINH_STYPE;
          break;
        default: Expectation_Error ("fractal keyword");
      }    
    END_CASE

    /* if any of the next become supported by the expression parser,
     * then their handling would need to move above to the FLOAT_FUNCT_TOKEN
     * case above.
     */

    CASE(SQR_TOKEN)
      Object->Sub_Type = SQR_STYPE;
    END_CASE

    CASE(PWR_TOKEN)
      Object->Sub_Type = PWR_STYPE;
      Parse_Float_Param2(&Object->exponent.x,&Object->exponent.y);
    END_CASE

    CASE(CUBE_TOKEN)
      Object->Sub_Type = CUBE_STYPE;
    END_CASE

    CASE(RECIPROCAL_TOKEN)
      Object->Sub_Type = RECIPROCAL_STYPE;
    END_CASE

    CASE(HYPERCOMPLEX_TOKEN)
      Object->Algebra = HYPERCOMPLEX_TYPE;
    END_CASE

    CASE(QUATERNION_TOKEN)
      Object->Algebra = QUATERNION_TYPE;
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE

  END_EXPECT

  Parse_Object_Mods((OBJECT *)Object);

  SetUp_Fractal(Object);

  return((OBJECT *)Object);
}
/*****************************************************************************
*
* FUNCTION
*
*   Parse_Lathe
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*
*   OBJECT * -
*   
* AUTHOR
*
*   Dieter Bayer
*   
* DESCRIPTION
*
*   Read a lathe primitive.
*
* CHANGES
*
*   Jun 1994 : Creation.
*
******************************************************************************/

static OBJECT *Parse_Lathe()
{
  int i;
  LATHE *Object;
  UV_VECT *Points;

  Parse_Begin();

  if ((Object = (LATHE *)Parse_Object_Id()) != NULL)
  {
    return((OBJECT *)Object);
  }

  Object = Create_Lathe();

  /* Determine kind of spline used and aspect ratio. */

  EXPECT
    CASE(LINEAR_SPLINE_TOKEN)
      Object->Spline_Type = LINEAR_SPLINE;
    END_CASE

    CASE(QUADRATIC_SPLINE_TOKEN)
      Object->Spline_Type = QUADRATIC_SPLINE;
    END_CASE

    CASE(CUBIC_SPLINE_TOKEN)
      Object->Spline_Type = CUBIC_SPLINE;
    END_CASE

    CASE(BEZIER_SPLINE_TOKEN)
      Object->Spline_Type = BEZIER_SPLINE;
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  /* Get number of points. */

  Object->Number = (int)Parse_Float();

  switch (Object->Spline_Type)
  {
    case LINEAR_SPLINE :

      if (Object->Number < 2)
      {
        Error("Lathe with linear splines must have at least two points.");
      }

      break;

    case QUADRATIC_SPLINE :

      if (Object->Number < 3)
      {
        Error("Lathe with quadratic splines must have at least three points.");
      }

      break;

    case CUBIC_SPLINE :

      if (Object->Number < 4)
      {
        Error("Prism with cubic splines must have at least four points.");
      }

      break;

    case BEZIER_SPLINE :

      if ((Object->Number & 3) != 0)
      {
        Error("Lathe with Bezier splines must have four points per segment.");
      }

      break;
  }

  /* Get temporary points describing the rotated curve. */

  Points = (UV_VECT *)POV_MALLOC(Object->Number*sizeof(UV_VECT), "temporary lathe points");

  /* Read points (x : radius; y : height; z : not used). */

  for (i = 0; i < Object->Number; i++)
  {
    Parse_Comma();

    Parse_UV_Vect(Points[i]);

    if ((i > 0) && (i < Object->Number - 1) && (Points[i][X] < 0.0))
    {
      Error("Incorrect point in lathe.");
    }
  }

  /* Compute spline segments. */

  Compute_Lathe(Object, Points);

  /* Compute bounding box. */

  Compute_Lathe_BBox(Object);

  /* Parse object's modifiers. */

  Parse_Object_Mods((OBJECT *)Object);

  /* Destroy temporary points. */

  POV_FREE (Points);

  return((OBJECT *) Object);
}



/*****************************************************************************
*
* FUNCTION
*
*   Parse_Light_Group
*
* INPUT
*
*   -
*
* OUTPUT
*
* RETURNS
*
*   Light group object
*   
* AUTHOR
*
*   Thorsten Froehlich [trf]
*   
* DESCRIPTION
*
*   Parse light_group object
*
* CHANGES
*
*   Jun 2000 : Creation.
*
******************************************************************************/

static OBJECT *Parse_Light_Group()
{
	CSG *Object;
	OBJECT *Local;
	VECTOR Local_Vector;
	MATRIX Local_Matrix;
	TRANSFORM Local_Trans;

	Parse_Begin();

	Object = Create_CSG_Union();

	Object->Type |= LIGHT_GROUP_OBJECT;
	Object->Children = NULL;
	Set_Flag(Object, NO_GLOBAL_LIGHTS_FLAG);

	while((Local = Parse_Object ()) != NULL)
	{
		// prevent light sources from being added to Frame.Light_Sources
 		if((Local->Type & LIGHT_SOURCE_OBJECT) == LIGHT_SOURCE_OBJECT)
			Local->Type |= LIGHT_GROUP_LIGHT_OBJECT;
		Local->Type |= IS_CHILD_OBJECT;
		Link(Local, &Object->Children);
	}

	Promote_Local_Lights(Object); // in lightgrp.cpp [trf]

	Compute_CSG_BBox((OBJECT *)Object);

	// Note: We cannot use Parse_Object_Mods here because 
	// it would allow all kinds of modifiers. However,
	// changing it to not allow those would slow it down,
	// so the bits of code needed are just duplicated
	// here. [trf]
	// Parse_Object_Mods((OBJECT *)Object);

	EXPECT
		CASE (TRANSLATE_TOKEN)
			Parse_Vector (Local_Vector);
			Compute_Translation_Transform(&Local_Trans, Local_Vector);
			Translate_Object ((OBJECT *)Object, Local_Vector, &Local_Trans);
		END_CASE

		CASE (ROTATE_TOKEN)
			Parse_Vector (Local_Vector);
			Compute_Rotation_Transform(&Local_Trans, Local_Vector);
			Rotate_Object ((OBJECT *)Object, Local_Vector, &Local_Trans);
		END_CASE

		CASE (SCALE_TOKEN)
			Parse_Scale_Vector (Local_Vector);
			Compute_Scaling_Transform(&Local_Trans, Local_Vector);
			Scale_Object ((OBJECT *)Object, Local_Vector, &Local_Trans);
		END_CASE

		CASE (TRANSFORM_TOKEN)
			Transform_Object ((OBJECT *)Object, Parse_Transform(&Local_Trans));
		END_CASE

		CASE (MATRIX_TOKEN)
			Parse_Matrix (Local_Matrix);
			Compute_Matrix_Transform(&Local_Trans, Local_Matrix);
			Transform_Object ((OBJECT *)Object, &Local_Trans);
		END_CASE

		CASE (GLOBAL_LIGHTS_TOKEN)
			Bool_Flag (Object, NO_GLOBAL_LIGHTS_FLAG, !(Allow_Float(1.0) > 0.5));
		END_CASE

		CASE(PHOTONS_TOKEN)
			Parse_Begin();
			EXPECT
				CASE(TARGET_TOKEN)
					Object->Ph_Density = Allow_Float(1.0);
					if (Object->Ph_Density > 0)
					{
						Set_Flag(Object,PH_TARGET_FLAG);
						CheckPassThru((OBJECT *)Object, PH_TARGET_FLAG);
					}
					else
					{
						Clear_Flag(Object, PH_TARGET_FLAG);
					}
				END_CASE

				CASE(REFRACTION_TOKEN)
					if((int)Parse_Float())
					{ 
						Set_Flag(Object, PH_RFR_ON_FLAG);
						Clear_Flag(Object, PH_RFR_OFF_FLAG);
						CheckPassThru((OBJECT *)Object, PH_RFR_ON_FLAG);
					}
					else
					{ 
						Clear_Flag(Object, PH_RFR_ON_FLAG);
						Set_Flag(Object, PH_RFR_OFF_FLAG);
					}
				END_CASE

				CASE(REFLECTION_TOKEN)
					if((int)Parse_Float())
					{ 
						Set_Flag(Object, PH_RFL_ON_FLAG); 
						Clear_Flag(Object, PH_RFL_OFF_FLAG); 
					}
					else
					{ 
						Clear_Flag(Object, PH_RFL_ON_FLAG); 
						Set_Flag(Object, PH_RFL_OFF_FLAG); 
					}
				END_CASE

				CASE(PASS_THROUGH_TOKEN)
					if((int)Allow_Float(1.0))
					{
						Set_Flag(Object, PH_PASSTHRU_FLAG);
						CheckPassThru((OBJECT *)Object, PH_PASSTHRU_FLAG);
					}
					else
					{
						Clear_Flag(Object, PH_PASSTHRU_FLAG);
					}
				END_CASE

				CASE(COLLECT_TOKEN)
					Bool_Flag (Object, PH_IGNORE_PHOTONS_FLAG, !(Allow_Float(1.0) > 0.0));
				END_CASE

				OTHERWISE
					UNGET
					EXIT
				END_CASE
			END_EXPECT
			Parse_End();
		END_CASE
		OTHERWISE
			UNGET
			EXIT
		END_CASE
	END_EXPECT

	Set_CSG_Children_Flag((OBJECT *)Object, Test_Flag((OBJECT *)Object, NO_GLOBAL_LIGHTS_FLAG),
	                      NO_GLOBAL_LIGHTS_FLAG, NO_GLOBAL_LIGHTS_SET_FLAG);

	Parse_End();

	return ((OBJECT *)Object);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Light_Source ()
{
   DBL Len;
   VECTOR Local_Vector;
   MATRIX Local_Matrix;
   TRANSFORM Local_Trans;
   LIGHT_SOURCE *Object;
   /* NK phmap */
   BLEND_MAP *Map;
   SNGL Value;
   int i;
   /* NK ---- */

   Parse_Begin ();

   if ( (Object = (LIGHT_SOURCE *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);
      
   Object = Create_Light_Source ();

   Parse_Vector(Object->Center);

   Parse_Comma();

   Parse_Colour (Object->Colour);

   EXPECT
     /* NK phmap */
     CASE (COLOUR_MAP_TOKEN)
       Destroy_Blend_Map(Object->blend_map);
       Map = Object->blend_map = Parse_Colour_Map ();

       Make_Colour(Object->Colour, 0, 0, 0);
       for (i = 0; i < Map->Number_Of_Entries; i++)
       {
         Value = Map->Blend_Map_Entries[i].value;
         CRGBAddScaledEq(Object->Colour,Value,Map->Blend_Map_Entries[i].Vals.Colour);
       }
     END_CASE

     CASE(PHOTONS_TOKEN)
       Parse_Begin();
       EXPECT
#ifdef GLOBAL_PHOTONS
         CASE(GLOBAL_TOKEN)
           Object->Ph_Density = Allow_Float(1.0);
           if (Object->Ph_Density > 0)
           {
             Set_Flag(Object, PH_TARGET_FLAG);
             /*CheckPassThru(Object, PH_TARGET_FLAG);*/
           }
           else
           {
             Clear_Flag(Object, PH_TARGET_FLAG);
           }
         END_CASE
#endif

         CASE(REFRACTION_TOKEN)
           if((int)Parse_Float())
           { 
             Set_Flag(Object, PH_RFR_ON_FLAG); 
             Clear_Flag(Object, PH_RFR_OFF_FLAG); 
           }
           else
           { 
             Clear_Flag(Object, PH_RFR_ON_FLAG); 
             Set_Flag(Object, PH_RFR_OFF_FLAG); 
           }
         END_CASE

         CASE(REFLECTION_TOKEN)
           if((int)Parse_Float())
           { 
             Set_Flag(Object, PH_RFL_ON_FLAG); 
             Clear_Flag(Object, PH_RFL_OFF_FLAG); 
           }
           else
           {
             Clear_Flag(Object, PH_RFL_ON_FLAG); 
             Set_Flag(Object, PH_RFL_OFF_FLAG); 
           }
         END_CASE

         CASE (AREA_LIGHT_TOKEN)
           Object->Photon_Area_Light = true;
         END_CASE

         OTHERWISE
           UNGET
           EXIT
         END_CASE
       END_EXPECT
       Parse_End();
     END_CASE

     CASE (LOOKS_LIKE_TOKEN)
       if (Object->Children != NULL)
         Error("Only one looks_like allowed per light_source.");
       Parse_Begin ();
       Object->Type &= ~(int)PATCH_OBJECT;
       if ((Object->Children = Parse_Object ()) == NULL)
         Expectation_Error ("object");
       Compute_Translation_Transform(&Local_Trans, Object->Center);
       Translate_Object (Object->Children, Object->Center, &Local_Trans);
       Parse_Object_Mods (Object->Children);
       Set_Flag(Object->Children, NO_SHADOW_FLAG);
       Set_Flag(Object, NO_SHADOW_FLAG);
       Object->Type |= (Object->Children->Type & CHILDREN_FLAGS);
       Set_Flag(Object, PH_PASSTHRU_FLAG);
     END_CASE

     CASE (PROJECTED_THROUGH_TOKEN)
       if (Object->Projected_Through_Object != NULL)
         Error("Only one projected through allowed per light_source.");
       Parse_Begin ();
       Object->Type &= ~(int)PATCH_OBJECT;
       if ((Object->Projected_Through_Object = Parse_Object ()) == NULL)
         Expectation_Error ("object");
       Parse_Object_Mods (Object->Projected_Through_Object);
       Set_Flag(Object, NO_SHADOW_FLAG);
       Set_Flag(Object, PH_PASSTHRU_FLAG);
     END_CASE

     CASE (FILL_LIGHT_TOKEN)
       Object->Light_Type = FILL_LIGHT_SOURCE;
     END_CASE

     CASE (PARALLEL_TOKEN)
       Object->Parallel= true;
     END_CASE

     CASE (SPOTLIGHT_TOKEN)
       Object->Light_Type = SPOT_SOURCE;
       Object->Radius = cos(30 * M_PI_180);
       Object->Falloff = cos(45 * M_PI_180);
       Object->Coeff = 0;
     END_CASE

     CASE (CYLINDER_TOKEN)
       Object->Light_Type = CYLINDER_SOURCE;
       Object->Radius = 0.75;
       Object->Falloff = 1;
       Object->Coeff = 0;
       Object->Parallel = true;
     END_CASE

     CASE (POINT_AT_TOKEN)
       if ((Object->Light_Type == SPOT_SOURCE) || (Object->Light_Type == CYLINDER_SOURCE)
           || Object->Parallel)
       {
         Parse_Vector(Object->Points_At);
       }
       else
       {
         Not_With ("point_at","standard light source");
       }
     END_CASE

     CASE (TIGHTNESS_TOKEN)
       if ((Object->Light_Type == SPOT_SOURCE) || (Object->Light_Type == CYLINDER_SOURCE))
         Object->Coeff = Parse_Float();
       else
         Not_With ("tightness","standard light source");
     END_CASE

     CASE (RADIUS_TOKEN)
       if ((Object->Light_Type == SPOT_SOURCE) || (Object->Light_Type == CYLINDER_SOURCE))
       {
         Object->Radius = Parse_Float();
         if (Object->Light_Type == SPOT_SOURCE)
         {
           Object->Radius  = cos(Object->Radius * M_PI_180);
         }
       }
       else
         Not_With ("radius","standard light source");
     END_CASE

     CASE (FALLOFF_TOKEN)
       if ((Object->Light_Type == SPOT_SOURCE) || (Object->Light_Type == CYLINDER_SOURCE))
       {
         Object->Falloff = Parse_Float();
         if (Object->Light_Type == SPOT_SOURCE)
         {
           Object->Falloff = cos(Object->Falloff * M_PI_180);
         }
       }
       else
         Not_With ("falloff","standard light source");
     END_CASE

     CASE (FADE_DISTANCE_TOKEN)
       Object->Fade_Distance = Parse_Float();
     END_CASE

     CASE (FADE_POWER_TOKEN)
       Object->Fade_Power = Parse_Float();
     END_CASE

     CASE (AREA_LIGHT_TOKEN)
       Object->Area_Light = true;
       Parse_Vector (Object->Axis1); Parse_Comma ();
       Parse_Vector (Object->Axis2); Parse_Comma ();
       Object->Area_Size1 = (int)Parse_Float(); Parse_Comma ();
       Object->Area_Size2 = (int)Parse_Float();
       Object->Light_Grid = Create_Light_Grid (Object->Area_Size1, Object->Area_Size2);
     END_CASE

     CASE (JITTER_TOKEN)
       Object->Jitter = true;
     END_CASE

     /* Orient area lights to the point [ENB 9/97] */
     CASE (ORIENT_TOKEN)
       Object->Orient = true;
       if (!(Object->Area_Light))
       {
         Warning(0,"Orient only affects area_light");
       }
     END_CASE

     /* Circular area lights [ENB 9/97] */
     CASE (CIRCULAR_TOKEN)
       Object->Circular = true;
       if (!(Object->Area_Light))
       {
         Warning(0,"Circular only affects area_light");
       }
     END_CASE

     CASE (ADAPTIVE_TOKEN)
       Object->Adaptive_Level = (int)Parse_Float();
     END_CASE

#ifdef AREA_LIGHT_REFLECTION_PATCH
     CASE (MAX_TRACE_LEVEL_TOKEN)
			 parseUnofficialFeature(120);
       Object->Area_Light_Max_Trace_Level = (int)Parse_Float();
       if (!(Object->Area_Light))
       {
         Warning(0,"max_trace_level only affects area_light");
       }
     END_CASE
#endif

     CASE (MEDIA_ATTENUATION_TOKEN)
       Object->Media_Attenuation = Allow_Float(1.0) > 0.0;
     END_CASE

     CASE (MEDIA_INTERACTION_TOKEN)
       Object->Media_Interaction = Allow_Float(1.0) > 0.0;
     END_CASE

     CASE (TRANSLATE_TOKEN)
       Parse_Vector (Local_Vector);
       Compute_Translation_Transform(&Local_Trans, Local_Vector);
       Translate_Object ((OBJECT *)Object, Local_Vector, &Local_Trans);
     END_CASE

     CASE (ROTATE_TOKEN)
       Parse_Vector (Local_Vector);
       Compute_Rotation_Transform(&Local_Trans, Local_Vector);
       Rotate_Object ((OBJECT *)Object, Local_Vector, &Local_Trans);
     END_CASE

     CASE (SCALE_TOKEN)
       Parse_Scale_Vector (Local_Vector);
       Compute_Scaling_Transform(&Local_Trans, Local_Vector);
       Scale_Object ((OBJECT *)Object, Local_Vector, &Local_Trans);
     END_CASE

     CASE (TRANSFORM_TOKEN)
       Transform_Object ((OBJECT *)Object, Parse_Transform(&Local_Trans));
     END_CASE

     CASE (MATRIX_TOKEN)
       Parse_Matrix (Local_Matrix);
       Compute_Matrix_Transform(&Local_Trans, Local_Matrix);
       Transform_Object ((OBJECT *)Object, &Local_Trans);
     END_CASE

      #ifdef GLOW_PATCH
        CASE(GLOW_TOKEN)
        {
          GLOW * NewGlow = NULL;
          parseUnofficialFeature(100);
          NewGlow = Create_Glow();
          Assign_Colour(NewGlow->Colour, Object->Colour);
          Assign_Vector(NewGlow->Center, Object->Center);
          Parse_Glow(NewGlow);
          Add_Glow(NewGlow);/*Add glow to Frame*/
          /*Add glow to list of glow pointers in object*/
          if(Object->glowList == NULL)
          {
            Object->glowList =(GLOW**) POV_MALLOC(sizeof(GLOW_PTR), "glow pointer array");
          }
          else
          {
            GLOW_PTR * newGlowList =(GLOW**) POV_MALLOC((Object->numOfGlows+2)*sizeof(GLOW_PTR), "glow pointer array");
            POV_MEMCPY(newGlowList, Object->glowList, (Object->numOfGlows+1)*sizeof(GLOW_PTR));
            POV_FREE(Object->glowList);
            Object->glowList = newGlowList;
          }
          Object->glowList[Object->numOfGlows] = NewGlow;
          Object->numOfGlows++;
        }
        END_CASE
      #endif

     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   Parse_End ();


   VSub(Object->Direction, Object->Points_At, Object->Center);

   VLength(Len, Object->Direction);

   if (Len > EPSILON)
   {
     VInverseScaleEq(Object->Direction, Len);
   }

   /* Make sure that circular light sources are larger than 1 by x [ENB 9/97] */
   if (Object->Circular)
   {
      if ((Object->Area_Size1 <= 1) || (Object->Area_Size2 <= 1))
      {
         Error("Circular area light must have more than 1 point per axis");
      }
   }

   return ((OBJECT *) Object);
}



/*****************************************************************************
*
* FUNCTION
*
*   Parse_Mesh
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*
*   OBJECT
*   
* AUTHOR
*
*   Dieter Bayer
*   
* DESCRIPTION
*
*   Read a triangle mesh.
*
* CHANGES
*
*   Feb 1995 : Creation.
*
******************************************************************************/

static OBJECT *Parse_Mesh()
{
  /* NK 1998 - added all sorts of uv variables*/
  int i;
  int number_of_normals, number_of_textures, number_of_triangles, number_of_vertices, number_of_uvcoords;
  int max_normals, max_textures, max_triangles, max_vertices, max_uvcoords;
  DBL l1, l2, l3;
  VECTOR D1, D2, P1, P2, P3, N1, N2, N3, N;
  UV_VECT UV1, UV2, UV3;
  SNGL_VECT *Normals, *Vertices;
  TEXTURE **Textures;
  UV_VECT *UVCoords;
  MESH *Object;
  MESH_TRIANGLE *Triangles;
  int fully_textured=true;
  /* NK 1998 */
  VECTOR Inside_Vect;
  TEXTURE *t2, *t3;
  bool foundZeroNormal=false;

  Make_Vector(Inside_Vect, 0, 0, 0);

  Parse_Begin();

  if ((Object = (MESH *)Parse_Object_Id()) != NULL)
  {
    return((OBJECT *)Object);
  }

  /* Create object. */

  Object = Create_Mesh();

  /* Allocate temporary normals, textures, triangles and vertices. */

  max_normals = 256;

  max_vertices = 256;

  max_textures = 16;

  max_triangles = 256;

  Normals = (SNGL_VECT *)POV_MALLOC(max_normals*sizeof(SNGL_VECT), "temporary triangle mesh data");

  Textures = (TEXTURE **)POV_MALLOC(max_textures*sizeof(TEXTURE *), "temporary triangle mesh data");

  Triangles = (MESH_TRIANGLE *)POV_MALLOC(max_triangles*sizeof(MESH_TRIANGLE), "temporary triangle mesh data");

  Vertices = (SNGL_VECT *)POV_MALLOC(max_vertices*sizeof(SNGL_VECT), "temporary triangle mesh data");

  /* Read raw triangle file. */

  number_of_normals = 0;

  number_of_textures = 0;

  number_of_triangles = 0;

  number_of_vertices = 0;

  /* NK 1998 */
  max_uvcoords = 256;
  UVCoords = (UV_VECT *)POV_MALLOC(max_uvcoords*sizeof(UV_VECT), "temporary triangle mesh data");
  number_of_uvcoords = 0;
  /* NK ---- */

  /* Create hash tables. */

  Create_Mesh_Hash_Tables();

  EXPECT
    CASE(TRIANGLE_TOKEN)
      Parse_Begin();

      Parse_Vector(P1);  Parse_Comma();
      Parse_Vector(P2);  Parse_Comma();
      Parse_Vector(P3);

      if (!Mesh_Degenerate(P1, P2, P3))
      {
        if (number_of_triangles >= max_triangles)
        {
          if (max_triangles >= INT_MAX/2)
          {
            Error("Too many triangles in triangle mesh.");
          }

          max_triangles *= 2;

          Triangles = (MESH_TRIANGLE *)POV_REALLOC(Triangles, max_triangles*sizeof(MESH_TRIANGLE), "triangle triangle mesh data");
        }

        /* Init triangle. */

        Init_Mesh_Triangle(&Triangles[number_of_triangles]);

        Triangles[number_of_triangles].P1 = Mesh_Hash_Vertex(&number_of_vertices, &max_vertices, &Vertices, P1);
        Triangles[number_of_triangles].P2 = Mesh_Hash_Vertex(&number_of_vertices, &max_vertices, &Vertices, P2);
        Triangles[number_of_triangles].P3 = Mesh_Hash_Vertex(&number_of_vertices, &max_vertices, &Vertices, P3);

        /* NK 1998 */
        Parse_Three_UVCoords(UV1,UV2,UV3);
        Triangles[number_of_triangles].UV1 = Mesh_Hash_UV(&number_of_uvcoords, &max_uvcoords, &UVCoords,UV1);
        Triangles[number_of_triangles].UV2 = Mesh_Hash_UV(&number_of_uvcoords, &max_uvcoords, &UVCoords,UV2);
        Triangles[number_of_triangles].UV3 = Mesh_Hash_UV(&number_of_uvcoords, &max_uvcoords, &UVCoords,UV3);
        /* NK ---- */

        /* NK */
        /* read possibly three instead of only one texture */
        /* read these before compute!!! */
        t2 = t3 = NULL;
        Triangles[number_of_triangles].Texture = Mesh_Hash_Texture(&number_of_textures, &max_textures, &Textures, Parse_Mesh_Texture(&t2,&t3));
        if (t2) Triangles[number_of_triangles].Texture2 = Mesh_Hash_Texture(&number_of_textures, &max_textures, &Textures, t2);
        if (t3) Triangles[number_of_triangles].Texture3 = Mesh_Hash_Texture(&number_of_textures, &max_textures, &Textures, t3);
        if (t2 || t3) Triangles[number_of_triangles].ThreeTex = true;

        Compute_Mesh_Triangle(&Triangles[number_of_triangles], false, P1, P2, P3, N);

        Triangles[number_of_triangles].Normal_Ind = Mesh_Hash_Normal(&number_of_normals, &max_normals, &Normals, N);


        if (Triangles[number_of_triangles].Texture < 0)
        {
          fully_textured = false;
        }

        number_of_triangles++;
      }
      /* NK degenerate fix */
      else
      {
        /* parse the uv and texture info - even though we'll just throw it
           away.  why?  if not we get a parse error - we should just ignore the
           degenerate triangle */
        t2=t3=NULL;
        Parse_Three_UVCoords(UV1,UV2,UV3);
        Parse_Mesh_Texture(&t2,&t3);
      }

      Parse_End();
    END_CASE

    CASE(SMOOTH_TRIANGLE_TOKEN)
      Parse_Begin();

      Parse_Vector(P1);  Parse_Comma();
      Parse_Vector(N1);  Parse_Comma();
      if(fabs(N1[X])<EPSILON && fabs(N1[Y])<EPSILON && fabs(N1[Z])<EPSILON)
      {
        N1[X] = 1.0;  // make it nonzero
        if(!foundZeroNormal)
          Warning(0,"Normal vector in mesh cannot be zero - changing it to <1,0,0>.");
        foundZeroNormal = true;
      }

      Parse_Vector(P2);  Parse_Comma();
      Parse_Vector(N2);  Parse_Comma();
      if(fabs(N2[X])<EPSILON && fabs(N2[Y])<EPSILON && fabs(N2[Z])<EPSILON)
      {
        N2[X] = 1.0;  // make it nonzero
        if(!foundZeroNormal)
          Warning(0,"Normal vector in mesh cannot be zero - changing it to <1,0,0>.");
        foundZeroNormal = true;
      }

      Parse_Vector(P3);  Parse_Comma();
      Parse_Vector(N3);
      if(fabs(N3[X])<EPSILON && fabs(N3[Y])<EPSILON && fabs(N3[Z])<EPSILON)
      {
        N3[X] = 1.0;  // make it nonzero
        if(!foundZeroNormal)
          Warning(0,"Normal vector in mesh cannot be zero - changing it to <1,0,0>.");
        foundZeroNormal = true;
      }

      VLength(l1, N1);
      VLength(l2, N2);
      VLength(l3, N3);

      if ((l1 != 0.0) && (l2 != 0.0) && (l3 != 0.0) && (!Mesh_Degenerate(P1, P2, P3)))
      {
        if (number_of_triangles >= max_triangles)
        {
          if (max_triangles >= INT_MAX/2)
          {
            Error("Too many triangles in triangle mesh.");
          }

          max_triangles *= 2;

          Triangles = (MESH_TRIANGLE *)POV_REALLOC(Triangles, max_triangles*sizeof(MESH_TRIANGLE), "triangle triangle mesh data");
        }

        VInverseScaleEq(N1, l1);
        VInverseScaleEq(N2, l2);
        VInverseScaleEq(N3, l3);

        /* Init triangle. */

        Init_Mesh_Triangle(&Triangles[number_of_triangles]);

        Triangles[number_of_triangles].P1 = Mesh_Hash_Vertex(&number_of_vertices, &max_vertices, &Vertices, P1);
        Triangles[number_of_triangles].P2 = Mesh_Hash_Vertex(&number_of_vertices, &max_vertices, &Vertices, P2);
        Triangles[number_of_triangles].P3 = Mesh_Hash_Vertex(&number_of_vertices, &max_vertices, &Vertices, P3);

        /* Check for equal normals. */

        VSub(D1, N1, N2);
        VSub(D2, N1, N3);

        VDot(l1, D1, D1);
        VDot(l2, D2, D2);

        /* NK 1998 */
        Parse_Three_UVCoords(UV1,UV2,UV3);
        Triangles[number_of_triangles].UV1 = Mesh_Hash_UV(&number_of_uvcoords, &max_uvcoords, &UVCoords,UV1);
        Triangles[number_of_triangles].UV2 = Mesh_Hash_UV(&number_of_uvcoords, &max_uvcoords, &UVCoords,UV2);
        Triangles[number_of_triangles].UV3 = Mesh_Hash_UV(&number_of_uvcoords, &max_uvcoords, &UVCoords,UV3);

        /* read possibly three instead of only one texture */
        /* read these before compute!!! */
        t2 = t3 = NULL;
        Triangles[number_of_triangles].Texture = Mesh_Hash_Texture(&number_of_textures, &max_textures, &Textures, Parse_Mesh_Texture(&t2,&t3));
        if (t2) Triangles[number_of_triangles].Texture2 = Mesh_Hash_Texture(&number_of_textures, &max_textures, &Textures, t2);
        if (t3) Triangles[number_of_triangles].Texture3 = Mesh_Hash_Texture(&number_of_textures, &max_textures, &Textures, t3);
        if (t2 || t3) Triangles[number_of_triangles].ThreeTex = true;

        if ((fabs(l1) > EPSILON) || (fabs(l2) > EPSILON))
        {
          /* Smooth triangle. */

          Triangles[number_of_triangles].N1 = Mesh_Hash_Normal(&number_of_normals, &max_normals, &Normals, N1);
          Triangles[number_of_triangles].N2 = Mesh_Hash_Normal(&number_of_normals, &max_normals, &Normals, N2);
          Triangles[number_of_triangles].N3 = Mesh_Hash_Normal(&number_of_normals, &max_normals, &Normals, N3);

          Compute_Mesh_Triangle(&Triangles[number_of_triangles], true, P1, P2, P3, N);
        }
        else
        {
          /* Flat triangle. */

          Compute_Mesh_Triangle(&Triangles[number_of_triangles], false, P1, P2, P3, N);
        }

        Triangles[number_of_triangles].Normal_Ind = Mesh_Hash_Normal(&number_of_normals, &max_normals, &Normals, N);

        if (Triangles[number_of_triangles].Texture < 0)
        {
          fully_textured = false;
        }

        number_of_triangles++;
      }
      /* NK degenerate fix */
      else
      {
        /* parse the uv and texture info - even though we'll just throw it
           away.  why?  if not we get a parse error - we should just ignore the
           degenerate triangle */
        t2=t3=NULL;
        Parse_Three_UVCoords(UV1,UV2,UV3);
        Parse_Mesh_Texture(&t2,&t3);
      }

      Parse_End();
    END_CASE

    /* NK 1998 */
    CASE(INSIDE_VECTOR_TOKEN)
      Parse_Vector(Inside_Vect);

    END_CASE
    /* NK ---- */

    #ifdef TRANSPARENT_INVERTED_NORMALS_PATCH
      CASE(INVERTED_NORMALS_TOKEN)
         parseUnofficialFeature(100);
        Object->inverted_normals = ( Allow_Float( (DBL) true ) != (DBL) false );
      END_CASE
    #endif

  OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  /* Destroy hash tables. */

  Destroy_Mesh_Hash_Tables();

  /* If there are no triangles something went wrong. */

  if (number_of_triangles == 0)
  {
    Error("No triangles in triangle mesh.");
  }

  /* Init triangle mesh data. */

  Object->Data = (MESH_DATA *)POV_MALLOC(sizeof(MESH_DATA), "triangle mesh data");


  Object->Data->References = 1;

  Object->Data->Tree = NULL;
  /* NK 1998 */
 
  if( (fabs(Inside_Vect[X]) < EPSILON) &&  (fabs(Inside_Vect[Y]) < EPSILON) &&  (fabs(Inside_Vect[Z]) < EPSILON))
    Object->has_inside_vector=false;
  else
  {
    VNormalize(Object->Data->Inside_Vect, Inside_Vect);
    Object->has_inside_vector=true;
  }
  
  Object->Data->Normals   = NULL;

  /* [LSK] Removed "Data->" */
  Object->Textures  = NULL;
  
  Object->Data->Triangles = NULL;
  Object->Data->Vertices  = NULL;

  /* Allocate memory for normals, textures, triangles and vertices. */

  Object->Number_Of_Textures = number_of_textures;

  Object->Data->Number_Of_Normals = number_of_normals;

  Object->Data->Number_Of_Triangles = number_of_triangles;

  Object->Data->Number_Of_Vertices = number_of_vertices;

  Object->Data->Normals = (SNGL_VECT *)POV_MALLOC(number_of_normals*sizeof(SNGL_VECT), "triangle mesh data");

  if (number_of_textures)
  {
    Set_Flag(Object, MULTITEXTURE_FLAG);

 	/* [LSK] Removed "Data->" */
   Object->Textures = (TEXTURE **)POV_MALLOC(number_of_textures*sizeof(TEXTURE *), "triangle mesh data");
  }

  Object->Data->Triangles = (MESH_TRIANGLE *)POV_MALLOC(number_of_triangles*sizeof(MESH_TRIANGLE), "triangle mesh data");

  Object->Data->Vertices = (SNGL_VECT *)POV_MALLOC(number_of_vertices*sizeof(SNGL_VECT), "triangle mesh data");

  /* Copy normals, textures, triangles and vertices into mesh. */

  for (i = 0; i < number_of_normals; i++)
  {
    Assign_Vector(Object->Data->Normals[i], Normals[i]);
  }

  for (i = 0; i < number_of_textures; i++)
  {
	/* [LSK] Removed "Data->" */
    Object->Textures[i] = Copy_Textures(Textures[i]);
    Post_Textures(Object->Textures[i]);

    /* now free the texture, in order to decrement the reference count */
    Destroy_Textures(Textures[i]);
  }
  
  if (fully_textured)
  {
    Object->Type |= TEXTURED_OBJECT;
  }

  for (i = 0; i < number_of_triangles; i++)
  {
    Object->Data->Triangles[i] = Triangles[i];
  }

  for (i = 0; i < number_of_vertices; i++)
  {
    Assign_Vector(Object->Data->Vertices[i], Vertices[i]);
  }

  /* NK 1998 */
  /* do the four steps above, but for UV coordinates*/
  Object->Data->UVCoords  = NULL;
  Object->Data->Number_Of_UVCoords = number_of_uvcoords;
  Object->Data->UVCoords = (UV_VECT *)POV_MALLOC(number_of_uvcoords*sizeof(UV_VECT), "triangle mesh data");
  for (i = 0; i < number_of_uvcoords; i++)
  {
    Assign_UV_Vect(Object->Data->UVCoords[i], UVCoords[i]);
  }
  POV_FREE(UVCoords);
  /* NK ---- */

  /* Free temporary memory. */

  POV_FREE(Normals);
  POV_FREE(Textures);
  POV_FREE(Triangles);
  POV_FREE(Vertices);

/*
  Render_Info("Mesh: %ld bytes: %ld vertices, %ld normals, %ld textures, %ld triangles\n",
    Object->Data->Number_Of_Normals*sizeof(SNGL_VECT)+
    Object->Number_Of_Textures*sizeof(TEXTURE *)+
    Object->Data->Number_Of_Triangles*sizeof(MESH_TRIANGLE)+
    Object->Data->Number_Of_Vertices*sizeof(SNGL_VECT),
    Object->Data->Number_Of_Vertices,
    Object->Data->Number_Of_Normals,
    Object->Number_Of_Textures,
    Object->Data->Number_Of_Triangles);
*/

  /* Create bounding box. */

  Compute_Mesh_BBox(Object);

  /* Parse object modifiers. */

  Parse_Object_Mods((OBJECT *)Object);

  /* Create bounding box tree. */

  Build_Mesh_BBox_Tree(Object);

  return((OBJECT *)Object);
}

/*****************************************************************************
*
* FUNCTION
*
*   Parse_Mesh2
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
*   OBJECT
*
* AUTHOR
*
*   Nathan Kopp
*
* DESCRIPTION
*
*   Read a triangle mesh - syntax version 2.
*
* CHANGES
*
*   Feb 1998 : Creation.
*
******************************************************************************/
static OBJECT *Parse_Mesh2()
{
  int i;
  int number_of_normals, number_of_textures, number_of_triangles, number_of_vertices, number_of_uvcoords;
  int number_of_normal_indices;
  int a,b,c;
  int n1, n2, n3;
  int found_normal_indices = 0;
  int found_uv_indices = 0;
  bool fully_textured = true;
  bool foundZeroNormal = false;

  DBL l1, l2;
  VECTOR D1, D2, P1, P2, P3, N1, N;
  VECTOR Inside_Vect;

  UV_VECT UV1;
  SNGL_VECT *Normals = NULL;
  SNGL_VECT *Vertices = NULL;
  TEXTURE **Textures = NULL;
  UV_VECT *UVCoords = NULL;
  MESH *Object;
  MESH_TRIANGLE *Triangles;

  Make_Vector(Inside_Vect, 0, 0, 0);

  Parse_Begin();

  if ((Object = (MESH *)Parse_Object_Id()) != NULL)
    return((OBJECT *)Object);

  /* Create object. */
  Object = Create_Mesh();

  /* normals, uvcoords, and textures are optional */
  number_of_vertices = 0;
  number_of_uvcoords = 0;
  number_of_textures = 0;
  number_of_normals = 0;
  number_of_normal_indices = 0;


  /* -----------------  Get the Normals & UV Vectors & Textures ------------ */
  EXPECT
    /* -------------------  Get the Vertices ------------------- */
    CASE(VERTEX_VECTORS_TOKEN)
      Parse_Begin();

      number_of_vertices = (int)Parse_Float(); Parse_Comma();

      if (number_of_vertices<=0)
        Error("No vertices in triangle mesh.");

      /* allocate memory for vertices */
      Vertices = (SNGL_VECT *)POV_MALLOC(number_of_vertices*sizeof(SNGL_VECT), "triangle mesh data");

      for(i=0; i<number_of_vertices; i++)
      {
        Parse_Vector(P1); Parse_Comma();
        Assign_Vector(Vertices[i], P1);
      }
      Parse_End();
    END_CASE

    CASE(NORMAL_VECTORS_TOKEN)
      Parse_Begin();
      number_of_normals = (int)Parse_Float(); Parse_Comma();

      if (number_of_normals>0)
      {
        Normals = (SNGL_VECT *)POV_MALLOC(number_of_normals*sizeof(SNGL_VECT), "triangle mesh data");

        /* leave space in the array for the raw triangle normals */
        for(i=0; i<number_of_normals; i++)
        {
          Parse_Vector(N1); Parse_Comma();
          if(fabs(N1[X])<EPSILON && fabs(N1[Y])<EPSILON && fabs(N1[Z])<EPSILON)
          {
            N1[X] = 1.0;  // make it nonzero
            if(!foundZeroNormal)
              Warning(0,"Normal vector in mesh2 cannot be zero - changing it to <1,0,0>.");
            foundZeroNormal = true;
          }
          VNormalizeEq(N1);
          Assign_Vector(Normals[i], N1);
        }
      }

      Parse_End();
    END_CASE

    CASE(UV_VECTORS_TOKEN)
      Parse_Begin();
      number_of_uvcoords = (int)Parse_Float(); Parse_Comma();

      if (number_of_uvcoords>0)
      {
        UVCoords = (UV_VECT *)POV_MALLOC(number_of_uvcoords*sizeof(UV_VECT), "triangle mesh data");

        for(i=0; i<number_of_uvcoords; i++)
        {
          Parse_UV_Vect(UV1); Parse_Comma();
          Assign_UV_Vect(UVCoords[i], UV1);
        }
      }

      Parse_End();
    END_CASE

    /*OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  EXPECT*/
    CASE(TEXTURE_LIST_TOKEN)
      Parse_Begin();

      number_of_textures = (int)Parse_Float();  Parse_Comma();

      if (number_of_textures>0)
      {
        Textures = (TEXTURE **)POV_MALLOC(number_of_textures*sizeof(TEXTURE *), "triangle mesh data");

        for(i=0; i<number_of_textures; i++)
        {
          /*
          GET(TEXTURE_ID_TOKEN)
          Textures[i] = Copy_Texture_Pointer((TEXTURE *)Token.Data);
          */
          GET(TEXTURE_TOKEN);
          Parse_Begin();
          Textures[i] = Parse_Texture();
          Post_Textures(Textures[i]);
          Parse_End();
          Parse_Comma();
        }
      }

      Parse_End();
      EXIT
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE

  END_EXPECT

  if (number_of_vertices == 0)
    Error("Vertex vectors not found in mesh2");

  /* first make sure we at least have one UV coordinate */
  if (number_of_uvcoords == 0)
  {
    number_of_uvcoords = 1;
    UVCoords = (UV_VECT *)POV_MALLOC(number_of_uvcoords*sizeof(UV_VECT), "triangle mesh data");
    UVCoords[0][U] = 0;
    UVCoords[0][V] = 0;
  }

  /* -------------------  Get the Faces ------------------- */
  GET(FACE_INDICES_TOKEN)
      Parse_Begin();

  /* number faces is mandatory, so we ask how many there are */
  #ifdef AVOID_TYPE_CONVERSION_WARNINGS_PATCH
    number_of_triangles = (int)Parse_Float(); Parse_Comma();
  #else
    number_of_triangles = Parse_Float(); Parse_Comma();
  #endif

  if (number_of_triangles == 0)
  {
    Error("No triangles in triangle mesh.");
  }

  /* allocate memory for triangles */
  Triangles = (MESH_TRIANGLE *)POV_MALLOC(number_of_triangles*sizeof(MESH_TRIANGLE), "triangle mesh data");

  /* start reading triangles */

      for(i=0; i<number_of_triangles; i++)
      {
        /* read in the indices vector */
        Parse_Vector(P1); Parse_Comma();

        /* convert the vector to integers */
        a = (int)P1[X];
        b = (int)P1[Y];
        c = (int)P1[Z];

        /* a--;b--;c--; use this to start external stuff at 1 */
        if ( a<0 || b<0 || c<0 ||
             a>=number_of_vertices || b>=number_of_vertices ||
             c>=number_of_vertices)
        {
          Error("Mesh face index out of range.");
        }

        /* Init triangle. */
        Init_Mesh_Triangle(&Triangles[i]);

        /* assign the vertices */
        Triangles[i].P1 = a;
        Triangles[i].P2 = b;
        Triangles[i].P3 = c;

        /* look for a texture index */
        EXPECT
          CASE_FLOAT
            #ifdef AVOID_TYPE_CONVERSION_WARNINGS_PATCH
              Triangles[i].Texture = (int)Parse_Float(); Parse_Comma();
            #else
     	       Triangles[i].Texture = Parse_Float(); Parse_Comma();
            #endif
            if (Triangles[i].Texture >= number_of_textures ||
                Triangles[i].Texture < 0)
              Error("Texture index out of range in mesh2.");
            EXIT
          END_CASE

          OTHERWISE
            Triangles[i].Texture = -1;
			fully_textured = false;
            EXIT
            UNGET
          END_CASE
        END_EXPECT
        /* look for a texture index */
        EXPECT
          CASE_FLOAT
            #ifdef AVOID_TYPE_CONVERSION_WARNINGS_PATCH
              Triangles[i].Texture2 = (int)Parse_Float(); Parse_Comma();
            #else
              Triangles[i].Texture2 = Parse_Float(); Parse_Comma();
            #endif
            if (Triangles[i].Texture2 >= number_of_textures ||
                Triangles[i].Texture2 < 0)
              Error("Texture index out of range in mesh2.");
            Triangles[i].ThreeTex = true;
            EXIT
          END_CASE
          OTHERWISE
            Triangles[i].Texture2 = -1;
            EXIT
            UNGET
          END_CASE
        END_EXPECT
        /* look for a texture index */
        EXPECT
          CASE_FLOAT
            #ifdef AVOID_TYPE_CONVERSION_WARNINGS_PATCH
              Triangles[i].Texture3 = (int)Parse_Float(); Parse_Comma();
            #else
              Triangles[i].Texture3 = Parse_Float(); Parse_Comma();
            #endif
            if (Triangles[i].Texture3 >= number_of_textures ||
                Triangles[i].Texture3 < 0)
              Error("Texture index out of range in mesh2.");
            Triangles[i].ThreeTex = true;
            EXIT
          END_CASE
          OTHERWISE
            Triangles[i].Texture3 = -1;
            EXIT
            UNGET
          END_CASE
        END_EXPECT

      }

      Parse_End();

  /* now we get the uv_indices & normal_indices in either order */

  EXPECT
    CASE(UV_INDICES_TOKEN)
      if (found_uv_indices)
      {
        Error("Only one uv_indices section is allowed in mesh2");
      }
      found_uv_indices = 1;
      Parse_Begin();

      if (Parse_Float() != number_of_triangles)
        Error("Number of uv indices must equal number of faces.");
      Parse_Comma();

      for (i=0; i<number_of_triangles; i++)
      {
        /* read in the indices vector */
        Parse_Vector(P1); Parse_Comma();

        /* convert the vector to integers */
        a = (int)P1[X];
        b = (int)P1[Y];
        c = (int)P1[Z];

        /* a--;b--;c--; use this to start external stuff at 1 */
        if ( a<0 || b<0 || c<0 ||
             a>=number_of_uvcoords || b>=number_of_uvcoords ||
             c>=number_of_uvcoords)
        {
          Error("Mesh UV index out of range.");
        }

        /* assign the uv coordinate */
        Triangles[i].UV1 = a;
        Triangles[i].UV2 = b;
        Triangles[i].UV3 = c;
      }
      Parse_End();
      /*EXIT*/
    END_CASE

  /*
    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  EXPECT
  */
    CASE(NORMAL_INDICES_TOKEN)
      if (found_normal_indices)
      {
        Error("Only one normal_indices section is allowed in mesh2");
      }
      found_normal_indices = 1;
      Parse_Begin();

      /*
      Change - if fewer normals than triangles, then no problem - the
      rest will be flat triangles.

      if (Parse_Float() != number_of_triangles)
        Error("Number of normal indices must equal number of faces.");
      */
      #ifdef AVOID_TYPE_CONVERSION_WARNINGS_PATCH
        number_of_normal_indices = (int)Parse_Float();
      #else
        number_of_normal_indices = Parse_Float();
      #endif
      if (number_of_normal_indices > number_of_triangles)
        Error("Number of normal indices cannot be more than the number of faces.");

      Parse_Comma();

      for (i=0; i<number_of_normal_indices; i++)
      {
        /* read in the indices vector */
        Parse_Vector(P1); Parse_Comma();

        /* convert the vector to integers */
        a = (int)P1[X];
        b = (int)P1[Y];
        c = (int)P1[Z];

        /* a--;b--;c--; use this to start external stuff at 1 */
        if ( a<0 || b<0 ||
             c<0 ||
             a>=number_of_normals || b>=number_of_normals ||
             c>=number_of_normals)
        {
          Error("Mesh normal index out of range.");
        }

        /* assign the uv coordinate */
        Triangles[i].N1 = a;
        Triangles[i].N2 = b;
        Triangles[i].N3 = c;
      }
      Parse_End();
      /*EXIT*/
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  /* ----------------------------------------------------- */
  /* ----------------------------------------------------- */

  EXPECT
    CASE(INSIDE_VECTOR_TOKEN)
      Parse_Vector(Inside_Vect);
    END_CASE

    #ifdef TRANSPARENT_INVERTED_NORMALS_PATCH
      CASE(INVERTED_NORMALS_TOKEN)
         parseUnofficialFeature(100);
        Object->inverted_normals = ( Allow_Float( (DBL) true ) != (DBL) false );
      END_CASE
    #endif

    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  if (fully_textured)
  {
    Object->Type |= TEXTURED_OBJECT;
  }

  if (!found_uv_indices)
  {
    if (number_of_uvcoords==number_of_vertices)
    {
      for (i=0; i<number_of_triangles; i++)
      {
        Triangles[i].UV1 = Triangles[i].P1;
        Triangles[i].UV2 = Triangles[i].P2;
        Triangles[i].UV3 = Triangles[i].P3;
      }
    }
    else if (number_of_uvcoords==1)
    {
      for (i=0; i<number_of_triangles; i++)
      {
        Triangles[i].UV1 = 0;
        Triangles[i].UV2 = 0;
        Triangles[i].UV3 = 0;
      }
    }
    else
    {
      Error("Missing uv_indicies section in mesh2.");
    }
  }

  if (!found_normal_indices)
  {
    if (number_of_normals==number_of_vertices)
    {
      /* If number of normals matches number of vertices, then assume
         that the normal_indices are the same as the triangle indices
         (left out for file size reasons).
         So, we pretend that we read in some normal_indices
      */
      number_of_normal_indices = number_of_triangles;

      for (i=0; i<number_of_triangles; i++)
      {
        Triangles[i].N1 = Triangles[i].P1;
        Triangles[i].N2 = Triangles[i].P2;
        Triangles[i].N3 = Triangles[i].P3;
      }
    }
    else if (number_of_normals)
    {
      Error("Missing normal_indicies section in mesh2.");
    }
  }

  /* ---------------- Compute Triangle Normals ---------------- */

  /* reallocate the normals stuff */
  if (!number_of_normals)
    Normals = (SNGL_VECT *)POV_MALLOC(number_of_triangles*sizeof(SNGL_VECT), "triangle mesh data");
  else
    Normals = (SNGL_VECT *)POV_REALLOC(Normals, (number_of_normals+number_of_triangles)*sizeof(SNGL_VECT), "triangle mesh data");

  for (i=0; i<number_of_triangles; i++)
  {
    a = Triangles[i].P1;
    b = Triangles[i].P2;
    c = Triangles[i].P3;
    n1 = Triangles[i].N1;
    n2 = Triangles[i].N2;
    n3 = Triangles[i].N3;

    Assign_Vector(P1, Vertices[a]);
    Assign_Vector(P2, Vertices[b]);
    Assign_Vector(P3, Vertices[c]);

    Triangles[i].Smooth = false;

    /* compute the normal (check for smoothness) */
    /* if number_of_normal_indices > 0, then the first triangles
       are smooth and the rest are flat */
    if (i<number_of_normal_indices)
    {
      /* Check for equal normals. */
      VSub(D1, Normals[n1], Normals[n2]);
      VSub(D2, Normals[n1], Normals[n3]);

      VDot(l1, D1, D1);
      VDot(l2, D2, D2);

      if ((fabs(l1) > EPSILON) || (fabs(l2) > EPSILON))
      {
        /* Smooth triangle. */
        Compute_Mesh_Triangle(&Triangles[i], true, P1, P2, P3, N);
        Triangles[i].Smooth = true;
      }
      else
      {
        /* Flat triangle. */
        Compute_Mesh_Triangle(&Triangles[i], false, P1, P2, P3, N);
      }
    }
    else
    {
      /* Flat triangle. */
      Compute_Mesh_Triangle(&Triangles[i], false, P1, P2, P3, N);
    }

    /* assign the triangle normal that we just computed */
    Triangles[i].Normal_Ind = i+number_of_normals;
    Assign_Vector(Normals[i+number_of_normals], N);
  }

  /* now remember how many normals we really have */
  number_of_normals += number_of_triangles;

  /* ----------------------------------------------------- */

  /* Init triangle mesh data. */
  Object->Data = (MESH_DATA *)POV_MALLOC(sizeof(MESH_DATA), "triangle mesh data");
  Object->Data->References = 1;
  Object->Data->Tree = NULL;
  /* NK 1998 */
  /*YS* 31/12/1999 */
 
  if( (fabs(Inside_Vect[X]) < EPSILON) &&  (fabs(Inside_Vect[Y]) < EPSILON) &&  (fabs(Inside_Vect[Z]) < EPSILON))
    Object->has_inside_vector=false;
  else
  {
    VNormalize(Object->Data->Inside_Vect, Inside_Vect);
    Object->has_inside_vector=true;
  }
  /*YS*/

  /* copy pointers to normals, triangles, textures, and vertices. */
  Object->Data->Normals   = Normals;
  Object->Data->Triangles = Triangles;
  Object->Data->Vertices  = Vertices;
  Object->Data->UVCoords  = UVCoords;
  /* [LSK] Removed "Data->" */
  Object->Textures  = Textures;

  /* copy number of for normals, textures, triangles and vertices. */
  Object->Data->Number_Of_Normals = number_of_normals;
  Object->Data->Number_Of_Triangles = number_of_triangles;
  Object->Data->Number_Of_Vertices = number_of_vertices;
  Object->Data->Number_Of_UVCoords  = number_of_uvcoords;
  Object->Number_Of_Textures = number_of_textures;

  if (number_of_textures)
  {
    Set_Flag(Object, MULTITEXTURE_FLAG);
  }

  /* Create bounding box. */
  Compute_Mesh_BBox(Object);

  /* Parse object modifiers. */
  Parse_Object_Mods((OBJECT *)Object);

  /* Create bounding box tree. */
  Build_Mesh_BBox_Tree(Object);

/*
  Render_Info("Mesh2: %ld bytes: %ld vertices, %ld normals, %ld textures, %ld triangles\n",
    Object->Data->Number_Of_Normals*sizeof(SNGL_VECT)+
    Object->Number_Of_Textures*sizeof(TEXTURE *)+
    Object->Data->Number_Of_Triangles*sizeof(MESH_TRIANGLE)+
    Object->Data->Number_Of_Vertices*sizeof(SNGL_VECT),
    Object->Data->Number_Of_Vertices,
    Object->Data->Number_Of_Normals,
    Object->Number_Of_Textures,
    Object->Data->Number_Of_Triangles);
*/

  return((OBJECT *)Object);
}


/*****************************************************************************
*
* FUNCTION
*
*   Parse_Mesh_Texture
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*
*   OBJECT
*   
* AUTHOR
*
*   Dieter Bayer
*   
* DESCRIPTION
*
*   Read an individual triangle mesh texture.
*
* CHANGES
*
*   Feb 1995 : Creation.
*
******************************************************************************/

static TEXTURE *Parse_Mesh_Texture (TEXTURE **t2, TEXTURE **t3)
{
  TEXTURE *Texture;

  Texture = NULL;

  EXPECT
    CASE(TEXTURE_TOKEN)
      Parse_Begin();

      GET(TEXTURE_ID_TOKEN);

      Texture = (TEXTURE *)Token.Data;

      Parse_End();
    END_CASE

    /* NK */
    CASE(TEXTURE_LIST_TOKEN)
      Parse_Begin();
      
      GET(TEXTURE_ID_TOKEN);
      Texture = (TEXTURE *)Token.Data;

      Parse_Comma();

      GET(TEXTURE_ID_TOKEN);
      *t2 = (TEXTURE *)Token.Data;

      Parse_Comma();

      GET(TEXTURE_ID_TOKEN);
      *t3 = (TEXTURE *)Token.Data;

      Parse_End();
      EXIT
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  return(Texture);
}



/*****************************************************************************
*
* FUNCTION
*
*   Parse_Parametric
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Parametric(void)
{
	PARAMETRIC  *Object;
	DBL         temp;
	char        PrecompFlag = 0;
	int         PrecompDepth = 1;
	UV_VECT     tempUV;

	Parse_Begin();

	if((Object = (PARAMETRIC *)Parse_Object_Id()) != NULL)
		return ((OBJECT *)Object);

	Object = Create_Parametric();

	EXPECT
		CASE(FUNCTION_TOKEN)
			Object->Function[0]= Parse_Function();
			EXIT
		END_CASE
		OTHERWISE
			Object->Function[0]= Parse_FunctionContent();
			EXIT
		END_CASE
	END_EXPECT

	Parse_Comma();

	EXPECT
		CASE(FUNCTION_TOKEN)
			Object->Function[1]= Parse_Function();
			EXIT
		END_CASE
		OTHERWISE
			Object->Function[1]= Parse_FunctionContent();
			EXIT
		END_CASE
	END_EXPECT

	Parse_Comma();

	EXPECT
		CASE(FUNCTION_TOKEN)
			Object->Function[2]= Parse_Function();
			EXIT
		END_CASE
		OTHERWISE
			Object->Function[2]= Parse_FunctionContent();
			EXIT
		END_CASE
	END_EXPECT

	Parse_UV_Vect(tempUV);
	Object->umin = tempUV[U];
	Object->vmin = tempUV[V];
	Parse_Comma();

	Parse_UV_Vect(tempUV);
	Object->umax = tempUV[U];
	Object->vmax = tempUV[V];

	if(Object->umin>Object->umax)
	{
		temp = Object->umin;
		Object->umin = Object->umax;
		Object->umax = temp;
	}
	if(Object->vmin>Object->vmax)
	{
		temp = Object->vmin;
		Object->vmin = Object->vmax;
		Object->vmax = temp;
	}

	EXPECT
		CASE(ACCURACY_TOKEN)
			Object->accuracy= Parse_Float();
		END_CASE

    CASE(MAX_GRADIENT_TOKEN)
      Object->max_gradient = Parse_Float();
    END_CASE

		CASE(PRECOMPUTE_TOKEN)
			PrecompDepth= Parse_Float();
			Parse_Comma();

			EXPECT
				CASE(VECTOR_FUNCT_TOKEN)
					if(Token.Function_Id != X_TOKEN)
					{
						UNGET
					}
					else
						PrecompFlag |= OK_X;
					EXIT
				END_CASE

				OTHERWISE
					UNGET
					EXIT
				END_CASE
			END_EXPECT

			Parse_Comma();

			EXPECT
				CASE(VECTOR_FUNCT_TOKEN)
					if(Token.Function_Id != Y_TOKEN)
					{
						UNGET
					}
					else
						PrecompFlag |= OK_Y;
					EXIT
				END_CASE

				OTHERWISE
					UNGET
					EXIT
				END_CASE
			END_EXPECT

			Parse_Comma();

			EXPECT
				CASE(VECTOR_FUNCT_TOKEN)
					if(Token.Function_Id != Z_TOKEN)
					{
						UNGET
					}
					else
						PrecompFlag |= OK_Z;
					EXIT
				END_CASE

				OTHERWISE
					UNGET
					EXIT
				END_CASE
			END_EXPECT
		END_CASE

		CASE(CONTAINED_BY_TOKEN)
			Parse_Begin();
			{ 
				int Exit_Flag2 = false; 

				while (!Exit_Flag2) 
				{
					Get_Token();  
					switch(Token.Token_Id) 
					{
						CASE(BOX_TOKEN)
							Object->container_shape = 0;

							Parse_Begin();

							Parse_Vector((Object->container.box.corner1));
							Parse_Comma();
							Parse_Vector((Object->container.box.corner2));

							Parse_End();

							if (Object->container.box.corner1[X] > Object->container.box.corner2[X])
							{
								temp = Object->container.box.corner1[X];
								Object->container.box.corner1[X] = Object->container.box.corner2[X];
								Object->container.box.corner2[X] = temp;
							}
							if (Object->container.box.corner1[Y] > Object->container.box.corner2[Y])
							{
								temp = Object->container.box.corner1[Y];
								Object->container.box.corner1[Y] = Object->container.box.corner2[Y];
								Object->container.box.corner2[Y] = temp;
							}
							if (Object->container.box.corner1[Z] > Object->container.box.corner2[Z])
							{
								temp = Object->container.box.corner1[Z];
								Object->container.box.corner1[Z] = Object->container.box.corner2[Z];
								Object->container.box.corner2[Z] = temp;
							}

							if (Object->Trans != NULL)
								Compute_Parametric_BBox(Object);

							Exit_Flag2 = true;
						END_CASE

						CASE(SPHERE_TOKEN)
							Object->container_shape = 1; 

							Parse_Begin();

							Parse_Vector(Object->container.sphere.center);
							Parse_Comma();
							Object->container.sphere.radius = Parse_Float();

							Parse_End();

							Make_BBox(Object->BBox,
							          Object->container.sphere.center[X] - Object->container.sphere.radius,
							          Object->container.sphere.center[Y] - Object->container.sphere.radius,
							          Object->container.sphere.center[Z] - Object->container.sphere.radius,
							          2.0 * Object->container.sphere.radius,
							          2.0 * Object->container.sphere.radius,
							          2.0 * Object->container.sphere.radius);

							if (Object->Trans != NULL)
								Compute_Parametric_BBox(Object);

							Exit_Flag2 = true;
						END_CASE

						OTHERWISE
							UNGET
							Exit_Flag2 = true;
			 		}
			 	}
			}
			Parse_End();
		END_CASE

		OTHERWISE
			UNGET
			EXIT
		END_CASE
	END_EXPECT

	Parse_Object_Mods((OBJECT *)Object);

	if(PrecompFlag != 0)
	{
		Object->PData = Precompute_Parametric_Values(Object, PrecompFlag, PrecompDepth);
	}

	return ((OBJECT *)Object);
}


/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Plane ()
{
   DBL len;
   PLANE *Object;

   Parse_Begin ();

   if ( (Object = (PLANE *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);

   Object = Create_Plane();

   Parse_Vector(Object->Normal_Vector);   Parse_Comma();
   VLength(len, Object->Normal_Vector);
   if (len < EPSILON)
   {
     Error("Degenerate plane normal.");
   }
   VInverseScaleEq(Object->Normal_Vector, len);
   Object->Distance = -Parse_Float();

   Compute_Plane_BBox(Object);

   Parse_Object_Mods ((OBJECT *)Object);

   return ((OBJECT *) Object);
  }



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Poly (int order)
{
   POLY *Object;

   Parse_Begin ();

   if ( (Object = (POLY *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);

   if (order == 0)
     {
      order = (int)Parse_Float();      Parse_Comma();
      if (order < 2 || order > MAX_ORDER)
        Error("Order of poly is out of range.");
     }

   Object = Create_Poly(order);

   Parse_Coeffs(Object->Order, &(Object->Coeffs[0]));

   Compute_Poly_BBox(Object);

   Parse_Object_Mods ((OBJECT *)Object);

   return ((OBJECT *) Object);
}


/*****************************************************************************
*
* FUNCTION
*
*   Parse_Polygon
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
*   OBJECT * -
*
* AUTHOR
*
*   Dieter Bayer
*   
* DESCRIPTION
*
*   -
*
* CHANGES
*
*   May 1994 : Creation.
*
*   Oct 1994 : Modified to use new polygon data structure. [DB]
*
******************************************************************************/

static OBJECT *Parse_Polygon()
{
  int i, closed = false;
  int Number;
  POLYGON *Object;
  VECTOR *Points;
  VECTOR P;

  Parse_Begin();

  if ((Object = (POLYGON *)Parse_Object_Id()) != NULL)
  {
    return((OBJECT *) Object);
  }

  Object = Create_Polygon();

  Number = (int)Parse_Float();

  if (Number < 3)
  {
    Error("Polygon needs at least three points.");
  }

  Points = (VECTOR *)POV_MALLOC((Number+1)*sizeof(VECTOR), "temporary polygon points");

  for (i = 0; i < Number; i++)
  {
    Parse_Comma();

    Parse_Vector(Points[i]);
  }

  /* Check for closed polygons. */

  Assign_Vector(P, Points[0]);

  for (i = 1; i < Number; i++)
  {
    closed = false;

    if ((fabs(P[X] - Points[i][X]) < EPSILON) &&
        (fabs(P[Y] - Points[i][Y]) < EPSILON) &&
        (fabs(P[Z] - Points[i][Z]) < EPSILON))
    {
      i++;

      if (i < Number)
      {
        Assign_Vector(P, Points[i]);
      }

      closed = true;
    }
  }

  if (!closed)
  {
    Warning(0, "Polygon not closed. Closing it.");

    Assign_Vector(Points[Number], P);

    Number++;
  }

  Compute_Polygon(Object, Number, Points);

  POV_FREE (Points);

  Parse_Object_Mods ((OBJECT *)Object);

  return((OBJECT *) Object);
}



/*****************************************************************************
*
* FUNCTION
*
*   Parse_Prism
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
*   OBJECT * -
*
* AUTHOR
*
*   Dieter Bayer
*
* DESCRIPTION
*
*   -
*
* CHANGES
*
*   May 1994 : Creation.
*
******************************************************************************/

static OBJECT *Parse_Prism()
{
  int i, closed = false;
  DBL h;
  int loopStart = 0;

  PRISM *Object;
  UV_VECT *Points;
  UV_VECT P;

  Parse_Begin();

  if ((Object = (PRISM *)Parse_Object_Id()) != NULL)
  {
    return((OBJECT *) Object);
  }

  Object = Create_Prism();

  /* 
   * Determine kind of spline used (linear, quadratic, cubic) 
   * and type of sweeping (linear, conic).
   */

  EXPECT
    CASE(LINEAR_SPLINE_TOKEN)
      Object->Spline_Type = LINEAR_SPLINE;
    END_CASE

    CASE(QUADRATIC_SPLINE_TOKEN)
      Object->Spline_Type = QUADRATIC_SPLINE;
    END_CASE

    CASE(CUBIC_SPLINE_TOKEN)
      Object->Spline_Type = CUBIC_SPLINE;
    END_CASE

    CASE(BEZIER_SPLINE_TOKEN)
      Object->Spline_Type = BEZIER_SPLINE;
    END_CASE

    CASE(LINEAR_SWEEP_TOKEN)
      Object->Sweep_Type = LINEAR_SWEEP;
    END_CASE

    CASE(CONIC_SWEEP_TOKEN)
      Object->Sweep_Type = CONIC_SWEEP;
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  /* Read prism heights. */

  Object->Height1 = Parse_Float(); Parse_Comma();
  Object->Height2 = Parse_Float(); Parse_Comma();

  if (Object->Height1 > Object->Height2)
  {
    h = Object->Height1;
    Object->Height1 = Object->Height2;
    Object->Height2 = h;
  }

  /* Get number of points = number of segments. */

  Object->Number = (int)Parse_Float();

  switch (Object->Spline_Type)
  {
    case LINEAR_SPLINE :

      if (Object->Number < 3)
      {
        Error("Prism with linear splines must have at least three points.");
      }

      break;

    case QUADRATIC_SPLINE :

      if (Object->Number < 5)
      {
        Error("Prism with quadratic splines must have at least five points.");
      }

      break;

    case CUBIC_SPLINE :

      if (Object->Number < 6)
      {
        Error("Prism with cubic splines must have at least six points.");
      }

      break;

    case BEZIER_SPLINE :

      if ((Object->Number & 3) != 0)
      {
        Error("Prism with Bezier splines must have four points per segment.");
      }

      break;
  }

  /* Allocate Object->Number points for the prism. */

  Points = (UV_VECT *)POV_MALLOC((Object->Number+1) * sizeof(UV_VECT), "temporary prism points");

  /* Read points (x, y : coordinate of 2d point; z : not used). */

  for (i = 0; i < Object->Number; i++)
  {
    Parse_Comma();

    Parse_UV_Vect(Points[i]);
  }

  /* Closed or not closed that's the question. */

  EXPECT
    CASE(OPEN_TOKEN)
      Clear_Flag(Object, CLOSED_FLAG);
      EXIT
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  /* Check for closed prism. */

  if ((Object->Spline_Type == LINEAR_SPLINE) ||
      (Object->Spline_Type == QUADRATIC_SPLINE) ||
      (Object->Spline_Type == CUBIC_SPLINE))
  {
    switch (Object->Spline_Type)
    {
      case LINEAR_SPLINE :

        i = 1;

        Assign_UV_Vect(P, Points[0]);

        break;

      case QUADRATIC_SPLINE :
      case CUBIC_SPLINE :

        i = 2;

        Assign_UV_Vect(P, Points[1]);

        break;
    }

    for ( ; i < Object->Number; i++)
    {
      closed = false;

      if ((fabs(P[X] - Points[i][X]) < EPSILON) &&
          (fabs(P[Y] - Points[i][Y]) < EPSILON))
      {
        switch (Object->Spline_Type)
        {
          case LINEAR_SPLINE :

            i++;

            if (i < Object->Number)
            {
              Assign_UV_Vect(P, Points[i]);
            }

            break;

          case QUADRATIC_SPLINE :

            i += 2;

            if (i < Object->Number)
            {
              Assign_UV_Vect(P, Points[i]);
            }

            break;

          case CUBIC_SPLINE :

            i += 3;

            if (i < Object->Number)
            {
              Assign_UV_Vect(P, Points[i]);
            }

            break;
        }

        closed = true;
      }
    }
  }
  else
  {
    closed = true;

    loopStart = 0;
    
	for (i = 4; i < Object->Number; i += 4)
    {
      if ((fabs(Points[i][X] - Points[i-1][X]) > EPSILON) ||
          (fabs(Points[i][Y] - Points[i-1][Y]) > EPSILON))
      {
		//. this is a different point.  Check if we have a loop.
		if ((fabs(Points[i-1][X] - Points[loopStart][X]) > EPSILON) ||
			(fabs(Points[i-1][Y] - Points[loopStart][Y]) > EPSILON))
		{
	  	  closed = false;
		  break;
		}

		loopStart = i;
      }
    }
	if ((fabs(Points[Object->Number-1][X] - Points[loopStart][X]) > EPSILON) ||
		(fabs(Points[Object->Number-1][Y] - Points[loopStart][Y]) > EPSILON))
	{
	  closed = false;
	}
  }

  if (!closed)
  {
    if (Object->Spline_Type == LINEAR_SPLINE)
    {
      Assign_UV_Vect(Points[Object->Number], P);

      Object->Number++;

      Warning(0, "Linear prism not closed. Closing it.");
    }
    else
    {
      Set_Flag(Object, DEGENERATE_FLAG);

      Warning(0, "Prism not closed. Ignoring it.");
    }
  }

  /* Compute spline segments. */

  Compute_Prism(Object, Points);

  /* Compute bounding box. */

  Compute_Prism_BBox(Object);

  /* Parse object's modifiers. */

  Parse_Object_Mods((OBJECT *)Object);

  /* Destroy temporary points. */

  POV_FREE (Points);

  return((OBJECT *) Object);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Quadric ()
{
   VECTOR Min, Max;
   QUADRIC *Object;

   Parse_Begin ();

   if ( (Object = (QUADRIC *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);

   Object = Create_Quadric();

        Parse_Vector(Object->Square_Terms);     Parse_Comma();
        Parse_Vector(Object->Mixed_Terms);      Parse_Comma();
        Parse_Vector(Object->Terms);            Parse_Comma();
   Object->Constant = Parse_Float();

   Make_Vector(Min, -BOUND_HUGE, -BOUND_HUGE, -BOUND_HUGE);
   Make_Vector(Max,  BOUND_HUGE,  BOUND_HUGE,  BOUND_HUGE);

   Compute_Quadric_BBox(Object, Min, Max);  

   Parse_Object_Mods ((OBJECT *)Object);

   return ((OBJECT *) Object);
  }



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Smooth_Triangle ()
{
   SMOOTH_TRIANGLE *Object;
   short degen;
   DBL vlen;

   degen=false;

   Parse_Begin ();

   if ( (Object = (SMOOTH_TRIANGLE *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);

   Object = Create_Smooth_Triangle();

   Parse_Vector (Object->P1);    Parse_Comma();
   Parse_Vector (Object->N1);    Parse_Comma();

   VLength(vlen,Object->N1);

   if (vlen == 0.0)
     degen=true;
   else
     VNormalize (Object->N1, Object->N1);

   Parse_Vector (Object->P2);    Parse_Comma();
   Parse_Vector (Object->N2);    Parse_Comma();

   VLength(vlen,Object->N2);

   if (vlen == 0.0)
     degen=true;
   else
     VNormalize (Object->N2, Object->N2);

   Parse_Vector (Object->P3);    Parse_Comma();
   Parse_Vector (Object->N3);

   VLength(vlen,Object->N3);

   if (vlen == 0.0)
     degen=true;
   else
     VNormalize (Object->N3, Object->N3);

   if (!degen)
   {
     degen=!Compute_Triangle ((TRIANGLE *) Object,true);
   }

   if (degen)
   {
     Warning(0, "Degenerate triangle. Please remove.");
   }

   Compute_Triangle_BBox((TRIANGLE *)Object);

   Parse_Object_Mods ((OBJECT *)Object);

   return ((OBJECT *) Object);
 }



/*****************************************************************************
*
* FUNCTION
*
*   Parse_Sor
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*
*   OBJECT * -
*
* AUTHOR
*
*   Dieter Bayer
*   
* DESCRIPTION
*
*   Read a surface of revolution primitive.
*
* CHANGES
*
*   May 1994 : Creation.
*
******************************************************************************/

static OBJECT *Parse_Sor()
{
  int i;
  SOR *Object;
  UV_VECT *Points;

  Parse_Begin();

  if ((Object = (SOR *)Parse_Object_Id()) != NULL)
  {
    return((OBJECT *)Object);
  }

  Object = Create_Sor();

  /* Get number of points. */

  Object->Number = (int)Parse_Float();

  if (Object->Number <4)
  {
    Error("Surface of revolution must have at least four points.");
  }

  /* Get temporary points describing the rotated curve. */

  Points = (UV_VECT *)POV_MALLOC(Object->Number*sizeof(UV_VECT), "temporary surface of revolution points");

  /* Read points (x : radius; y : height; z : not used). */

  for (i = 0; i < Object->Number; i++)
  {
    Parse_Comma();

    Parse_UV_Vect(Points[i]);

    if ((Points[i][X] < 0.0) ||
        ((i > 1 ) && (i < Object->Number - 1) && (Points[i][Y] <= Points[i-1][Y])))
    {
      Error("Incorrect point in surface of revolution.");
    }
  }

  /* Closed or not closed that's the question. */

  EXPECT
    CASE(OPEN_TOKEN)
      Clear_Flag(Object, CLOSED_FLAG);
      EXIT
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  /* There are Number-3 segments! */

  Object->Number -= 3;

  /* Compute spline segments. */

  Compute_Sor(Object, Points);

  /* Compute bounding box. */

  Compute_Sor_BBox(Object);

  /* Parse object's modifiers. */

  Parse_Object_Mods((OBJECT *)Object);

  /* Destroy temporary points. */

  POV_FREE (Points);

  return ((OBJECT *) Object);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Sphere()
{
  SPHERE *Object;

  Parse_Begin();

  if ((Object = (SPHERE *)Parse_Object_Id()) != NULL)
  {
    return ((OBJECT *) Object);
  }

  Object = Create_Sphere();

  Parse_Vector(Object->Center);

  Parse_Comma();

  Object->Radius = Parse_Float();

  Compute_Sphere_BBox(Object);  

  Parse_Object_Mods((OBJECT *)Object);

  return((OBJECT *)Object);
}



/*****************************************************************************
*
* FUNCTION
*
*       Parse_Sphere_Sweep
*
* INPUT
*
*   -
*
* OUTPUT
*
*   -
*
* RETURNS
*
*   Object
*
* AUTHOR
*
*       Jochen Lippert
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Sphere_Sweep()
{
    SPHERE_SWEEP    *Object;
    int                             i;
    
    Parse_Begin();
    
    if ((Object = (SPHERE_SWEEP *)Parse_Object_Id()) != NULL)
    {
        return ((OBJECT *) Object);
    }
    
    Object = Create_Sphere_Sweep();
    
    /* Get type of interpolation */
    EXPECT
        CASE(LINEAR_SPLINE_TOKEN)
            Object->Interpolation = LINEAR_SPHERE_SWEEP;
            EXIT
        END_CASE
        CASE(CUBIC_SPLINE_TOKEN)
            Object->Interpolation = CATMULL_ROM_SPLINE_SPHERE_SWEEP;
            EXIT
        END_CASE
        CASE(B_SPLINE_TOKEN)
            Object->Interpolation = B_SPLINE_SPHERE_SWEEP;
            EXIT
        END_CASE
        OTHERWISE
            UNGET
            EXIT
        END_CASE
    END_EXPECT
    
    if (Object->Interpolation == -1)
    {
        Error("Invalid type of interpolation.");
    }
    
    Parse_Comma();
    
    /* Get number of modeling spheres */
    Object->Num_Modeling_Spheres = (int)Parse_Float();
    
    Object->Modeling_Sphere =
        (SPHSWEEP_SPH *)POV_MALLOC(Object->Num_Modeling_Spheres * sizeof(SPHSWEEP_SPH),
        "sphere sweep modeling spheres");
    
    for (i = 0; i < Object->Num_Modeling_Spheres; i++)
    {
        Parse_Comma();
        Parse_Vector(Object->Modeling_Sphere[i].Center);
        Parse_Comma();
        Object->Modeling_Sphere[i].Radius = Parse_Float();
    }
    
    EXPECT
        CASE(TOLERANCE_TOKEN)
            Object->Depth_Tolerance = Parse_Float();
            EXIT
        END_CASE
        OTHERWISE
            UNGET
            EXIT
        END_CASE
    END_EXPECT
    
    Compute_Sphere_Sweep(Object);
    
    Compute_Sphere_Sweep_BBox(Object);
    
    Parse_Object_Mods((OBJECT *)Object);
    
    return ((OBJECT *)Object);
}



/*****************************************************************************
*
* FUNCTION
*
*   Parse_Superellipsoid
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*
*   OBJECT * -
*   
* AUTHOR
*
*   Dieter Bayer
*   
* DESCRIPTION
*
*   Read a superellipsoid primitive.
*
* CHANGES
*
*   Oct 1994 : Creation.
*
******************************************************************************/

static OBJECT *Parse_Superellipsoid()
{
  UV_VECT V1;
  SUPERELLIPSOID *Object;

  Parse_Begin();

  if ((Object = (SUPERELLIPSOID *)Parse_Object_Id()) != NULL)
  {
    return((OBJECT *)Object);
  }

  Object = Create_Superellipsoid();

  Parse_UV_Vect(V1);

  /* The x component is e, the y component is n. */

  Object->Power[X] = 2.0  / V1[X];
  Object->Power[Y] = V1[X] / V1[Y];
  Object->Power[Z] = 2.0  / V1[Y];

  /* Compute bounding box. */

  Compute_Superellipsoid_BBox(Object);

  /* Parse object's modifiers. */

  Parse_Object_Mods((OBJECT *)Object);

  return((OBJECT *) Object);
}


/*****************************************************************************
*
* FUNCTION
*
*   Parse_Torus
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*
*   OBJECT
*
* AUTHOR
*
*   Dieter Bayer
*   
* DESCRIPTION
*
*   -
*
* CHANGES
*
*   Jul 1994 : Creation.
*
******************************************************************************/

static OBJECT *Parse_Torus()
{
  TORUS *Object;

  Parse_Begin();

  if ((Object = (TORUS *)Parse_Object_Id()) != NULL)
  {
    return((OBJECT *)Object);
  }

  Object = Create_Torus();

  /* Read in the two radii. */

  Object->R = Parse_Float(); /* Big radius */

  Parse_Comma();

  Object->r = Parse_Float(); /* Little radius */

  Compute_Torus_BBox(Object);

  Parse_Object_Mods ((OBJECT *)Object);

  return ((OBJECT *) Object);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Triangle()
{
  TRIANGLE *Object;

  Parse_Begin();

  if ((Object = (TRIANGLE *)Parse_Object_Id()) != NULL)
  {
    return((OBJECT *) Object);
  }

  Object = Create_Triangle();

  Parse_Vector(Object->P1);    Parse_Comma();
  Parse_Vector(Object->P2);    Parse_Comma();
  Parse_Vector(Object->P3);

  /* Note that Compute_Triangle also computes the bounding box. */

  if (!Compute_Triangle(Object, false))
  {
    Warning(0, "Degenerate triangle. Please remove.");
  }

  Parse_Object_Mods((OBJECT *)Object);

  return((OBJECT *)Object);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_TrueType ()
{
   OBJECT *Object;
   char *filename;
   UCS2 *text_string;
   DBL depth;
   VECTOR offset;
   TRANSFORM Local_Trans;

  #ifdef FONT_ALIGNMENT_PATCH
    int h_position=H_ALIGN_LEFT;
    int v_position=V_ALIGN_NONE;
  #endif
   Parse_Begin ();
   
   GET(TTF_TOKEN);

   if ( (Object = (OBJECT *)Parse_Object_Id()) != NULL)
      return ((OBJECT *) Object);
      
   /*** Object = Create_TTF(); */

   /* Parse the TrueType font file name */
   filename = Parse_C_String(true);
   Parse_Comma();

   /* Parse the text string to be rendered */
   text_string = Parse_String();
   Parse_Comma();

   /* Get the extrusion depth */
   depth = Parse_Float(); Parse_Comma ();

   /* Get the offset vector */
   Parse_Vector(offset);

   /* Process all this good info */
   Object = (OBJECT *)Create_CSG_Union ();


    #ifdef FONT_ALIGNMENT_PATCH
      EXPECT
        CASE (H_ALIGN_LEFT_TOKEN)
          parseUnofficialFeature(100);
          h_position =H_ALIGN_LEFT;
        END_CASE

        CASE(H_ALIGN_CENTER_TOKEN)
          parseUnofficialFeature(100);
          h_position = H_ALIGN_CENTER;
        END_CASE

        CASE(H_ALIGN_RIGHT_TOKEN)
          parseUnofficialFeature(100);
          h_position = H_ALIGN_RRIGHT;
        END_CASE

        CASE (V_ALIGN_TOP_TOKEN)
          parseUnofficialFeature(100);
          v_position =V_ALIGN_TOP;
        END_CASE

        CASE(V_ALIGN_CENTER_TOKEN)
          parseUnofficialFeature(100);
          v_position = V_ALIGN_CENTER;
        END_CASE

        CASE(V_ALIGN_BOTTOM_TOKEN)
          parseUnofficialFeature(100);
          v_position = V_ALIGN_BOTTOM;
        END_CASE


        OTHERWISE
          UNGET
          EXIT
        END_CASE
      END_EXPECT

      ProcessNewTTF((OBJECT *)Object, filename, text_string, depth, offset, h_position, v_position );
      #else
     ProcessNewTTF((OBJECT *)Object, filename, text_string, depth, offset);
      #endif

   /* Free up the filename and text string memory */
   POV_FREE (filename);
   POV_FREE (text_string);

   /**** Compute_TTF_BBox(Object); */
   Compute_CSG_BBox((OBJECT *)Object);

/* This tiny rotation should fix cracks in text that lies along an axis */
   Make_Vector(offset, 0.001, 0.001, 0.001);
   Compute_Rotation_Transform(&Local_Trans, offset);
   Rotate_Object ((OBJECT *)Object, offset, &Local_Trans);

   /* Get any rotate/translate or texturing stuff */
   Parse_Object_Mods ((OBJECT *)Object);

   return ((OBJECT *) Object);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

OBJECT *Parse_Object ()
{
   OBJECT *Object = NULL;

   EXPECT

     CASE (ISOSURFACE_TOKEN)
       Object = Parse_Isosurface ();
       EXIT
     END_CASE

     CASE (PARAMETRIC_TOKEN)
       Object = Parse_Parametric ();
       EXIT
     END_CASE

     CASE (JULIA_FRACTAL_TOKEN)
       Object = Parse_Julia_Fractal ();
       EXIT
     END_CASE

     CASE (SPHERE_TOKEN)
       Object = Parse_Sphere ();
       EXIT
     END_CASE

     CASE (SPHERE_SWEEP_TOKEN)
       Object = Parse_Sphere_Sweep ();
       EXIT
     END_CASE

     CASE (PLANE_TOKEN)
       Object = Parse_Plane ();
       EXIT
     END_CASE

     CASE (CONE_TOKEN)
       Object = Parse_Cone ();
       EXIT
     END_CASE

     CASE (CYLINDER_TOKEN)
       Object = Parse_Cylinder ();
       EXIT
     END_CASE

     CASE (DISC_TOKEN)
       Object = Parse_Disc ();
       EXIT
     END_CASE

     CASE (QUADRIC_TOKEN)
       Object = Parse_Quadric ();
       EXIT
     END_CASE

     CASE (CUBIC_TOKEN)
       Object = Parse_Poly (3);
       EXIT
     END_CASE

     CASE (QUARTIC_TOKEN)
       Object = Parse_Poly (4);
       EXIT
     END_CASE

     CASE (POLY_TOKEN)
       Object = Parse_Poly (0);
       EXIT
     END_CASE

     CASE (TORUS_TOKEN)
       Object = Parse_Torus ();
       EXIT
     END_CASE

     /* Parse lathe primitive. [DB 8/94] */

     CASE (LATHE_TOKEN)
       Object = Parse_Lathe();
       EXIT
     END_CASE

     /* Parse polygon primitive. [DB 8/94] */

     CASE (POLYGON_TOKEN)
       Object = Parse_Polygon();
       EXIT
     END_CASE

     /* Parse prism primitive. [DB 8/94] */

     CASE (PRISM_TOKEN)
       Object = Parse_Prism();
       EXIT
     END_CASE

     /* Parse surface of revolution primitive. [DB 8/94] */

     CASE (SOR_TOKEN)
       Object = Parse_Sor();
       EXIT
     END_CASE

     /* Parse superellipsoid primitive. [DB 11/94] */

     CASE (SUPERELLIPSOID_TOKEN)
       Object = Parse_Superellipsoid();
       EXIT
     END_CASE

     /* Parse triangle mesh primitive. [DB 2/95] */

     CASE (MESH_TOKEN)
       Object = Parse_Mesh();
       EXIT
     END_CASE

     /* NK 1998 Parse triangle mesh primitive - syntax version 2. */
     CASE (MESH2_TOKEN)
       Object = Parse_Mesh2();
       EXIT
     END_CASE
     /* NK ---- */

     CASE (TEXT_TOKEN)
       Object = Parse_TrueType ();
       EXIT
     END_CASE

     CASE (OBJECT_ID_TOKEN)
       Object = Copy_Object((OBJECT *) Token.Data);
       EXIT
     END_CASE

     CASE (UNION_TOKEN)
       Object = Parse_CSG (CSG_UNION_TYPE);
       EXIT
     END_CASE

     CASE (LIGHT_GROUP_TOKEN)
       Object = Parse_Light_Group ();
       EXIT
     END_CASE

     CASE (COMPOSITE_TOKEN)
       Warning(150, "Use union instead of composite.");
       Object = Parse_CSG (CSG_UNION_TYPE);
       EXIT
     END_CASE

     CASE (MERGE_TOKEN)
       Object = Parse_CSG (CSG_MERGE_TYPE);
       EXIT
     END_CASE

     CASE (INTERSECTION_TOKEN)
       Object = Parse_CSG (CSG_INTERSECTION_TYPE);
       EXIT
     END_CASE

     CASE (DIFFERENCE_TOKEN)
       Object = Parse_CSG (CSG_DIFFERENCE_TYPE+CSG_INTERSECTION_TYPE);
       EXIT
     END_CASE

     CASE (BICUBIC_PATCH_TOKEN)
       Object = Parse_Bicubic_Patch ();
       EXIT
     END_CASE

     CASE (TRIANGLE_TOKEN)
       Object = Parse_Triangle ();
       EXIT
     END_CASE

     CASE (SMOOTH_TRIANGLE_TOKEN)
       Object = Parse_Smooth_Triangle ();
       EXIT
     END_CASE

     CASE (HEIGHT_FIELD_TOKEN)
       Object = Parse_HField ();
       EXIT
     END_CASE

     CASE (BOX_TOKEN)
       Object = Parse_Box ();
       EXIT
     END_CASE

     CASE (BLOB_TOKEN)
       Object = Parse_Blob ();
       EXIT
     END_CASE

     CASE (LIGHT_SOURCE_TOKEN)
       Object = Parse_Light_Source ();
       EXIT
     END_CASE

#ifdef MOTION_BLUR_PATCH
     CASE (MOTION_BLUR_TOKEN)
     parseUnofficialFeature(110);
       Object = Parse_Motion_Blur();
       EXIT
     END_CASE
#endif

     CASE (OBJECT_TOKEN)
       Parse_Begin ();
       Object = Parse_Object ();
       if (!Object)
         Expectation_Error ("object");
       Parse_Object_Mods ((OBJECT *)Object);
       EXIT
     END_CASE

     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   return ((OBJECT *) Object);
  }
  


/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Default ()
{
   TEXTURE *Local_Texture;
   PIGMENT *Local_Pigment;
   TNORMAL *Local_Tnormal;
   FINISH  *Local_Finish;

   Not_In_Default = false;
   Parse_Begin();

   EXPECT
     CASE (TEXTURE_TOKEN)
       Local_Texture = Default_Texture;
       Parse_Begin ();
       Default_Texture = Parse_Texture();
       Parse_End ();
       if (Default_Texture->Type != PLAIN_PATTERN)
         Error("Default texture cannot be material map or tiles.");
       if (Default_Texture->Next != NULL)
         Error("Default texture cannot be layered.");
       Destroy_Textures(Local_Texture);
     END_CASE

     CASE (PIGMENT_TOKEN)
       Local_Pigment = Copy_Pigment((Default_Texture->Pigment));
       Parse_Begin ();
       Parse_Pigment (&Local_Pigment);
       Parse_End ();
       Destroy_Pigment(Default_Texture->Pigment);
       Default_Texture->Pigment = Local_Pigment;
     END_CASE

     CASE (TNORMAL_TOKEN)
       Local_Tnormal = Copy_Tnormal((Default_Texture->Tnormal));
       Parse_Begin ();
       Parse_Tnormal (&Local_Tnormal);
       Parse_End ();
       Destroy_Tnormal(Default_Texture->Tnormal);
       Default_Texture->Tnormal = Local_Tnormal;
     END_CASE

     CASE (FINISH_TOKEN)
       Local_Finish = Copy_Finish((Default_Texture->Finish));
       Parse_Finish (&Local_Finish);
       Destroy_Finish(Default_Texture->Finish);
       Default_Texture->Finish = Local_Finish;
     END_CASE

     CASE (CAMERA_TOKEN)
       Parse_Camera (&Default_Camera);
     END_CASE

     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   Parse_End();

   Not_In_Default = true;
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static void Parse_Frame ()
{
   OBJECT *Object;
   RAINBOW  *Local_Rainbow;
   FOG  *Local_Fog;
   SKYSPHERE  *Local_Skysphere;
   int i;
   bool had_camera = false;

   EXPECT
     CASE (RAINBOW_TOKEN)
       Local_Rainbow = Parse_Rainbow();
       Local_Rainbow->Next = Frame.Rainbow;
       Frame.Rainbow = Local_Rainbow;
     END_CASE

     CASE (SKYSPHERE_TOKEN)
       if (opts.Options & OUTPUT_ALPHA)
       {
         Warning(0,"Alpha channel output is enabled, so the skysphere will only be visible in reflections.");
       }
       Local_Skysphere = Parse_Skysphere();
       if (Frame.Skysphere != NULL)
       {
         Warning(0, "Only one sky-sphere allowed (last one will be used).");
         Destroy_Skysphere(Frame.Skysphere);
       }
       Frame.Skysphere = Local_Skysphere;
       for (i=0; i<Local_Skysphere->Count; i++)
       {
         Post_Pigment(Local_Skysphere->Pigments[i]);
       }
     END_CASE

     CASE (FOG_TOKEN)
       Local_Fog = Parse_Fog();
       Local_Fog->Next = Frame.Fog;
       Frame.Fog = Local_Fog;
     END_CASE

     CASE (MEDIA_TOKEN)
       Parse_Media(&Frame.Atmosphere);
     END_CASE

     CASE (BACKGROUND_TOKEN)
       if (opts.Options & OUTPUT_ALPHA)
       {
         Warning(0,"Alpha channel output is enabled, so background color will only be visible in reflections.");
       }
       Parse_Begin();
       Parse_Colour (Frame.Background_Colour);
       Parse_End();
     END_CASE

     CASE (CAMERA_TOKEN)
       if (opts.Language_Version >= 350)
       {
         if (had_camera == true)
           Warning(0, "More than one camera in scene. Ignoring previous camera(s).");
         had_camera = true;

         Destroy_Camera(Frame.Camera);
         Frame.Camera = Copy_Camera(Default_Camera);
       }

       Parse_Camera(&Frame.Camera);
     END_CASE

      #ifdef CLOTHRAY_PATCH
        CASE (SIMCLOTH_TOKEN)
          parseUnofficialFeature(100);
          Parse_Cloth_Sys();
        END_CASE
      #endif

     CASE (DECLARE_TOKEN)
       UNGET
       Warning(295,"Should have '#' before 'declare'.");
       Parse_Directive (false);
     END_CASE

     CASE (INCLUDE_TOKEN)
       UNGET
       Warning(295,"Should have '#' before 'include'.");
       Parse_Directive (false);
     END_CASE

     CASE (FLOAT_FUNCT_TOKEN)
       switch(Token.Function_Id)
       {
          case VERSION_TOKEN:
            UNGET
            Parse_Directive (false);
            UNGET
            break;
            
          default:
            UNGET
            Expectation_Error ("object or directive");
            break;
       }
     END_CASE

     CASE (MAX_TRACE_LEVEL_TOKEN)
       if (opts.Language_Version >= 350)
       {
         PossibleError("'max_trace_level' should be in global_settings block.\n"
                       "Future versions may not support 'max_trace_level' outside global_settings.");
       }
       Global_Setting_Warn();
       Max_Trace_Level = (int)Parse_Float();
       Max_Trace_Level = max(1, Max_Trace_Level);
       Had_Max_Trace_Level = true;
       if(Max_Trace_Level > MAX_TRACE_LEVEL_LIMIT)
       {
         Warning(0, "Maximum max_trace_level is %d but %d was specified.\n"
                    "Going to use max_trace_level %d.",
                    MAX_TRACE_LEVEL_LIMIT, Max_Trace_Level, MAX_TRACE_LEVEL_LIMIT);
         Max_Trace_Level = MAX_TRACE_LEVEL_LIMIT;
       }
     END_CASE

     CASE (MAX_INTERSECTIONS)
       if (opts.Language_Version >= 350)
       {
         PossibleError("'max_intersections' should be in global_settings block.\n"
                       "Future versions may not support 'max_intersections' outside global_settings.");
       }
       Global_Setting_Warn();
       Max_Intersections = (int)Parse_Float();
       Max_Intersections = max((unsigned int)2, Max_Intersections);
     END_CASE

     CASE (DEFAULT_TOKEN)
       Parse_Default();
     END_CASE

     CASE (END_OF_FILE_TOKEN)
       EXIT
     END_CASE

     CASE (GLOBAL_SETTINGS_TOKEN)
       Parse_Global_Settings();
     END_CASE

    #ifdef GLOW_PATCH
        // YS sept 18 2000 added this here too to allow
        // for a scene with only glow
        CASE(GLOW_TOKEN)
        parseUnofficialFeature(100);
        {
          GLOW * NewGlow = Create_Glow();
          Parse_Glow(NewGlow);
          Add_Glow(NewGlow);
        }
        END_CASE
    #endif

     OTHERWISE
       UNGET
       Object = Parse_Object();
       if (Object == NULL)
         Expectation_Error ("object or directive");
       Post_Process (Object, NULL);
       Link_To_Frame (Object);
     END_CASE
   END_EXPECT
  }



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static void Parse_Global_Settings()
{
   Parse_Begin();
   EXPECT
#ifdef MOTION_BLUR_PATCH
     CASE (MOTION_BLUR_TOKEN)
        parseUnofficialFeature(110);

       opts.motionBlurCount = (int)Parse_Float();
       Parse_Comma();
       opts.blurClockDelta = Parse_Float();

       if(opts.motionBlurCount<=1)
         Error("motion blur count must be greater than 1.\n");
       if(opts.blurClockDelta<=0)
         Error("blur clock delta must be greater than zero.\n");
		GlobalSettingsForMotionBlurParsed=true;
     END_CASE
#endif
     CASE (IRID_WAVELENGTH_TOKEN)
       Parse_Colour (Frame.Irid_Wavelengths);
     END_CASE
     CASE (CHARSET_TOKEN)
       EXPECT
         CASE (ASCII_TOKEN)
           opts.String_Encoding = 0; // ASCII
           EXIT
         END_CASE
         CASE (UTF8_TOKEN)
           opts.String_Encoding = 1; // UTF8
           EXIT
         END_CASE
         CASE (SYS_TOKEN)
           opts.String_Encoding = 2; // System Specific
           EXIT
         END_CASE
         OTHERWISE
           Expectation_Error ("charset type");
         END_CASE
       END_EXPECT
     END_CASE
     CASE (ASSUMED_GAMMA_TOKEN)
     {
       DBL AssumedGamma;
       AssumedGamma = Parse_Float ();

       if (fabs(AssumedGamma - opts.DisplayGamma) < 0.1)
       {
         opts.GammaFactor = 1.0;
         opts.Options &= ~GAMMA_CORRECT; /* turn off gamma correction */
       }
       else
       {
         opts.GammaFactor = AssumedGamma/opts.DisplayGamma;
         opts.Options |= GAMMA_CORRECT; /* turn on gamma correction */
       }
     }
     END_CASE

     CASE (MAX_TRACE_LEVEL_TOKEN)
     {
       int Trace_Level = (int)Parse_Float();
       Max_Trace_Level = max(1, Trace_Level);
       Had_Max_Trace_Level = true;
       if(Max_Trace_Level > MAX_TRACE_LEVEL_LIMIT)
       {
         Warning(0, "Maximum max_trace_level is %d but %d was specified.\n"
                    "Going to use max_trace_level %d.",
                    MAX_TRACE_LEVEL_LIMIT, Max_Trace_Level, MAX_TRACE_LEVEL_LIMIT);
         Max_Trace_Level = MAX_TRACE_LEVEL_LIMIT;
       }
     }
     END_CASE

     CASE (ADC_BAILOUT_TOKEN)
       ADC_Bailout = Parse_Float ();
     END_CASE

     CASE (NUMBER_OF_WAVES_TOKEN)
       Number_Of_Waves = (int) Parse_Float ();
       if(Number_Of_Waves <=0)
       {
         Warning(0, "Illegal Value: Number_Of_Waves = 0. Changed to 1.");
         Number_Of_Waves = 1;
       }
     END_CASE

     CASE (MAX_INTERSECTIONS)
       Max_Intersections = (unsigned int) Parse_Float ();
       Max_Intersections = max((unsigned int)2, Max_Intersections);
     END_CASE

     CASE (NOISE_GENERATOR_TOKEN)
       opts.Noise_Generator = (int) Parse_Float();
       if (opts.Noise_Generator < 1 || opts.Noise_Generator > 3)
         Error ("Value for noise_generator in global_settings must be 1, 2, or 3.");
     END_CASE

     CASE (AMBIENT_LIGHT_TOKEN)
       Parse_Colour (Frame.Ambient_Light);
     END_CASE

     /* NK phmap */
     CASE (PHOTONS_TOKEN)
       Parse_Begin();
       EXPECT
         CASE(RADIUS_TOKEN)
           photonOptions.photonMap.minGatherRad = Allow_Float(-1.0);
           Parse_Comma();
           photonOptions.photonMap.minGatherRadMult = Allow_Float(1.0);
           Parse_Comma();
           photonOptions.mediaPhotonMap.minGatherRadMult = Allow_Float(-1.0);
           Parse_Comma();
           photonOptions.mediaPhotonMap.minGatherRadMult = Allow_Float(1.0);
         END_CASE

         CASE(SPACING_TOKEN)
           photonOptions.surfaceSeparation = Parse_Float();
         END_CASE

#ifdef GLOBAL_PHOTONS
         CASE(GLOBAL_TOKEN)
           photonOptions.globalCount = (int)Parse_Float();
         END_CASE
#endif
			 
         CASE (EXPAND_THRESHOLDS_TOKEN)
           photonOptions.expandTolerance = Parse_Float(); Parse_Comma();
           #ifdef AVOID_TYPE_CONVERSION_WARNINGS_PATCH
             photonOptions.minExpandCount = (int)Parse_Float();
           #else
             photonOptions.minExpandCount = Parse_Float();
           #endif
           if (photonOptions.expandTolerance<0.0)
           {
             Warning(100,"The first parameter of expand_thresholds must be greater than or equal to 0.\nSetting it to 0 now.");
             photonOptions.expandTolerance = 0.0;
           }
           if (photonOptions.minExpandCount<0)
           {
             Warning(100,"The second parameter of expand_thresholds must be greater than or equal to 0.\nSetting it to 0 now.");
             photonOptions.minExpandCount = 0;
           }
         END_CASE

         CASE (GATHER_TOKEN)
           photonOptions.minGatherCount = (int)Parse_Float();
           Parse_Comma();
           photonOptions.maxGatherCount = (int)Parse_Float();
         END_CASE

         CASE (JITTER_TOKEN)
           photonOptions.jitter = Parse_Float();
         END_CASE

         CASE (COUNT_TOKEN)
           photonOptions.surfaceCount = (int)Parse_Float();
         END_CASE

         CASE (AUTOSTOP_TOKEN)
           photonOptions.autoStopPercent = Parse_Float();
         END_CASE

         CASE (ADC_BAILOUT_TOKEN)
           photonOptions.ADC_Bailout = Parse_Float ();
         END_CASE

         CASE (MAX_TRACE_LEVEL_TOKEN)
           #ifdef AVOID_TYPE_CONVERSION_WARNINGS_PATCH
             photonOptions.Max_Trace_Level = (int)Parse_Float();
           #else
             photonOptions.Max_Trace_Level = Parse_Float();
           #endif
         END_CASE

         CASE(LOAD_FILE_TOKEN)
           if(photonOptions.fileName)
           {
             if(photonOptions.loadFile)
               Warning(100,"Filename already given, using new name");
             else
               Warning(100,"Cannot both load and save photon map. Now switching to load mode.");
             POV_FREE(photonOptions.fileName);
           }
           photonOptions.fileName = Parse_C_String(true);
           photonOptions.loadFile = true;
         END_CASE

         CASE(SAVE_FILE_TOKEN)
           if(photonOptions.fileName)
           {
             if(!photonOptions.loadFile)
               Warning(100,"Filename already given, using new name");
             else
               Warning(100,"Cannot both load and save photon map. Now switching to save mode.");
             POV_FREE(photonOptions.fileName);
           }
           photonOptions.fileName = Parse_C_String(true);
           photonOptions.loadFile = false;
         END_CASE

         CASE(MEDIA_TOKEN)
           photonOptions.maxMediaSteps = (int)Parse_Float(); Parse_Comma();
           if (photonOptions.maxMediaSteps<0)
             Error("max media steps must be non-negative.");

           photonOptions.mediaSpacingFactor = Allow_Float(1.0);
           if (photonOptions.mediaSpacingFactor <= 0.0)
             Error("media spacing factor must be greater than zero.");
         END_CASE

         OTHERWISE
           UNGET
           EXIT
         END_CASE
       END_EXPECT

       /* max_gather_count = 0  means no photon maps */
       if (photonOptions.maxGatherCount > 0)
         photonOptions.photonsEnabled = true;
       else
         photonOptions.photonsEnabled = false;

       if(photonOptions.photonsEnabled)
       {
         /* check for range errors */
         if (photonOptions.minGatherCount < 0)
           Error("min_gather_count cannot be negative.");
         if (photonOptions.maxGatherCount < 0)
           Error("max_gather_count cannot be negative.");
         if (photonOptions.minGatherCount > photonOptions.maxGatherCount)
           Error("min_gather_count must be less than max_gather_count.");
       }

       Parse_End();

       if((photonOptions.photonsEnabled == false) && (opts.Quality < 9))
       {
         photonOptions.photonsEnabled = false;
         Warning(0, "A photons{}-block has been found but photons remain disabled because\n"
                    "the output quality is set to 8 or less.");
       }
     END_CASE
 

     CASE (RADIOSITY_TOKEN)
       Experimental_Flag |= EF_RADIOS;

       /* enable radiosity only if the user includes a "radiosity" token */
       opts.Radiosity_Enabled = (opts.Quality >= 9);
       if((opts.Radiosity_Enabled == false) || (opts.Language_Version < 350))
       {
         Warning(0, "In POV-Ray 3.5 and later a radiosity{}-block will automatically\n"
                      "turn on radiosity if the output quality is set to 9 or higher.\n"
                      "Read the documentation to find out more about radiosity changes!");
       }

       Parse_Begin();
       EXPECT
         CASE(LOAD_FILE_TOKEN)
           if(opts.Radiosity_Load_File_Name)
           {
             Warning(100,"Filename already given, using new name");
             POV_FREE(opts.Radiosity_Load_File_Name);
           }
           opts.Radiosity_Load_File_Name = Parse_C_String(true);
         END_CASE

         CASE(SAVE_FILE_TOKEN)
           if(opts.Radiosity_Save_File_Name)
           {
             Warning(100,"Filename already given, using new name");
             POV_FREE(opts.Radiosity_Save_File_Name);
           }
           opts.Radiosity_Save_File_Name = Parse_C_String(true);
         END_CASE

         CASE(ALWAYS_SAMPLE_TOKEN)
           opts.Radiosity_Add_On_Final_Trace = (int)Parse_Float();
         END_CASE

         CASE (PRETRACE_START_TOKEN)
           opts.radPretraceStart = Parse_Float();
         END_CASE

         CASE (PRETRACE_END_TOKEN)
           opts.radPretraceEnd = Parse_Float();
         END_CASE

#ifdef NEW_RADIOSITY_PRETRACE_PATCH /* Christoph Hormann March 2005 */

         CASE (ADAPTIVE_TOKEN)
			     parseUnofficialFeature(120);
           opts.radPretraceAdaptive = (int)Allow_Float(0.0);
         END_CASE

#endif

         CASE (BRIGHTNESS_TOKEN)
           if ((opts.Radiosity_Brightness = Parse_Float()) <= 0.0)
           {
              Error("Radiosity brightness must be a positive number.");
           }
         END_CASE

         CASE (COUNT_TOKEN)
           if (( opts.Radiosity_Count = (int)Parse_Float()) <= 0)
           {
             Error("Radiosity count must be a positive number.");
           }
#ifdef RAD_CUSTOM_SAMPLES_PATCH    /* Christoph Hormann January 2003 */
           if ((opts.Radiosity_Count > 1600) && (opts.Rad_Samples_Type == 0))
           {
             Error("Radiosity count by default can not be more than 1600,\nif you want to use higher counts you need to use custom\nray directions.");
             opts.Radiosity_Count = 1600;
           }
           if ((opts.Radiosity_Count > opts.Rad_Samples_Count) && (opts.Rad_Samples_Type == -1))
           {
             Error("Radiosity count specified is higher than the number\nof sample ray directions.");
             opts.Radiosity_Count = opts.Rad_Samples_Count;
           }
#else
           if ( opts.Radiosity_Count > 1600)
           {
             Error("Radiosity count can not be more than 1600.");
             opts.Radiosity_Count = 1600;
           }
#endif
         END_CASE

         CASE (DISTANCE_MAXIMUM_TOKEN)
           if (( opts.Radiosity_Dist_Max = Parse_Float()) < 0.0)
           {
             Error("Radiosity distance maximum must be a positive number.");
           }
         END_CASE

         CASE (ERROR_BOUND_TOKEN)

#ifdef RADIOSITY_ADAPTIVE_RADIUS_PATCH

					 EXPECT
					 	 /* old syntax */
				 		 CASE_EXPRESS
           	 	 if (( opts.Radiosity_Error_Bound = Parse_Float()) <= 0.0)
           	 	 {
           	 	   Error("Radiosity error bound must be a positive number.");
           	 	 }
							 EXIT
						 END_CASE

						 /* new syntax */
						 CASE(LEFT_CURLY_TOKEN)
               parseUnofficialFeature(120);
							 UNGET
							 Parse_Begin();

           	 	 if (( opts.Radiosity_Error_Bound = Parse_Float()) <= 0.0)
           	 	 {
           	 	   Error("Radiosity error bound must be a positive number.");
           	 	 }

							 EXPECT
								 CASE (ADAPTIVE_TOKEN)
								 	 if (( opts.Radiosity_Error_Bound_Multiplier = Allow_Float(1.5)) <= 1.0)
								 	 {
								 		 Error("Radiosity adaptive error bound multiplier must be > 1.");
								 	 }
							 	 	 Parse_Comma();
							 	 	 opts.Radiosity_Error_Bound_Max = Allow_Float(10.0);

								 END_CASE

								 OTHERWISE
								 	 UNGET
								 	 EXIT
								 END_CASE

							 END_EXPECT

							 Parse_End();
							 EXIT
						 END_CASE

						 OTHERWISE
							 UNGET
							 EXIT
						 END_CASE

					 END_EXPECT

#else
           if (( opts.Radiosity_Error_Bound = Parse_Float()) <= 0.0)
           {
             Error("Radiosity error bound must be a positive number.");
           }
#endif

         END_CASE

         CASE (GRAY_THRESHOLD_TOKEN)
           opts.Radiosity_Gray = Parse_Float();
           if (( opts.Radiosity_Gray < 0.0) || ( opts.Radiosity_Gray > 1.0))
           {
              Error("Radiosity gray threshold must be from 0.0 to 1.0.");
           }
         END_CASE

         CASE (LOW_ERROR_FACTOR_TOKEN)
           if (( opts.Radiosity_Low_Error_Factor = Parse_Float()) <= 0.0)
           {
             Error("Radiosity low error factor must be a positive number.");
           }
         END_CASE

         CASE (MINIMUM_REUSE_TOKEN)
           if (( opts.Radiosity_Min_Reuse = Parse_Float()) < 0.0)
           {
              Error("Radiosity min reuse can not be a negative number.");
           }
         END_CASE

         CASE (NEAREST_COUNT_TOKEN)
           opts.Radiosity_Nearest_Count = (int)Parse_Float();
           if (( opts.Radiosity_Nearest_Count < 1) ||
               ( opts.Radiosity_Nearest_Count > MAX_NEAREST_COUNT))
           {
              Error("Radiosity nearest count must be a value from 1 to %d.", (int)MAX_NEAREST_COUNT);
           }
         END_CASE

         CASE (RECURSION_LIMIT_TOKEN)
           opts.Radiosity_Recursion_Limit = (int)Parse_Float();
           if ((opts.Radiosity_Recursion_Limit < 1) || (opts.Radiosity_Recursion_Limit > 20))
           {
              Error("Radiosity recursion limit must be in the range 1 to 20.");
           }
         END_CASE

         /* NK radiosity */
         CASE (MAX_SAMPLE_TOKEN)
           opts.Maximum_Sample_Brightness = Parse_Float();
         END_CASE

         CASE (ADC_BAILOUT_TOKEN)
           if (( opts.Radiosity_ADC_Bailout = Parse_Float()) <= 0)
           {
              Error("ADC Bailout must be a positive number.");
           }
         END_CASE

         CASE (TNORMAL_TOKEN)
           opts.Radiosity_Use_Normal = (int)Parse_Float();
         END_CASE

         CASE (MEDIA_TOKEN)
           opts.Radiosity_Use_Media = (int)Parse_Float();
         END_CASE
#ifdef MOTION_BLUR_PATCH
         CASE (MOTION_BLUR_TOKEN)
         parseUnofficialFeature(110);
           opts.Radiosity_Use_Motion_Blur = (int)Parse_Float();
         END_CASE
#endif

#ifdef RAD_CUSTOM_SAMPLES_PATCH    /* Christoph Hormann January 2003 */
         CASE (SAMPLES_TOKEN)
           parseUnofficialFeature(110);

           EXPECT
             /* internal type */
             CASE_EXPRESS
               opts.Rad_Samples_Type = (int)Parse_Float();
               if (opts.Rad_Samples_Type<0)
                 Error("radiosity samples type must be >=0.");

               EXIT
             END_CASE

             /* custom sample directions */
             CASE(LEFT_CURLY_TOKEN)
               UNGET
               Parse_Begin();

               opts.Rad_Samples_Type = -1;
               opts.Rad_Samples_Count = (int)Parse_Float();
               if (opts.Rad_Samples_Count<1)
                 Error("Number of radiosity sample directions must be >0.");

               opts.Rad_Samples = (VECTOR *)POV_MALLOC(opts.Rad_Samples_Count*sizeof (VECTOR), "radiosity sample directions");

               int i;

               for (i=0; i<opts.Rad_Samples_Count; i++)
               {
                 Parse_Comma();
                 Parse_Vector(opts.Rad_Samples[i]);
                 VNormalizeEq(opts.Rad_Samples[i]);
               }

               EXPECT
                 CASE(WEIGHT_TOKEN)
                   opts.Rad_Samples_Weight = (int)Parse_Float();
                 END_CASE
                 OTHERWISE
                   UNGET
                   EXIT
                 END_CASE
               END_EXPECT

               Parse_End();

             END_CASE
             OTHERWISE
               UNGET
               EXIT
             END_CASE
           END_EXPECT

         END_CASE /* Samples */

#endif

#ifdef RAD_SAMPLES_RANDOMIZE_PATCH    /* Christoph Hormann January 2003 */
         CASE (RANDOMIZE_TOKEN)
           parseUnofficialFeature(110);
           opts.Rad_Samples_Random = (int)Parse_Float();
         END_CASE
#endif

#ifdef RADIOSITY_VISUALIZE_PATCH
         CASE (SHOW_SAMPLES_TOKEN)
           parseUnofficialFeature(120);
           opts.Radiosity_Show_Sample_Points_Radius = Parse_Float();
					 if (opts.Radiosity_Show_Sample_Points_Radius > 0.0)
					 {
						 Parse_Comma();
						 Parse_Colour (opts.Radiosity_Show_Sample_Points);
					 }
         END_CASE
         CASE (SHOW_LOW_COUNT_TOKEN)
           parseUnofficialFeature(120);
				   opts.Radiosity_Show_Low_Count = true;
				   Parse_Colour (opts.Radiosity_Low_Count_Color);
					 Parse_Comma();
				   Parse_Colour (opts.Radiosity_Gather_Final_Color);
         END_CASE
#endif

         OTHERWISE
           UNGET
           EXIT
         END_CASE
       END_EXPECT
       Parse_End();
     END_CASE

      #ifdef MECHANICS_SIMULATION_PATCH            /* Christoph Hormann August 2002 */
        CASE (MECHSIM_TOKEN)
          parseUnofficialFeature(100);
          Parse_Mech_Sim();
        END_CASE
      #endif

      #ifdef POST_PROCESS_PATCH
        CASE (POST_PROCESS_TOKEN)
          parseUnofficialFeature(110);
          Parse_Post_Process();
        END_CASE
      #endif

      #ifdef FILM_EXPOSURE_SIMULATION_PATCH
       CASE (EXPOSURE_TOKEN)
         parseUnofficialFeature(100);
         opts.Exposure = Parse_Float ();
       END_CASE
       CASE (EXPOSURE_GAIN_TOKEN)
         parseUnofficialFeature(100);
         opts.Exposure_Gain = Parse_Float ();
        END_CASE
      #endif

#ifdef CUSTOM_TONE_MAPPING_PATCH    /* Christoph Hormann September 2004 */
		 CASE (TONE_MAPPING_TOKEN)
			 parseUnofficialFeature(120);
			 Parse_Begin();

			 Get_Token();
			 if(Token.Token_Id != FUNCTION_TOKEN)
				 Parse_Error(FUNCTION_TOKEN);

			 opts.Mapping_Function = Parse_Function();

			 EXPECT

				 CASE (INVERSE_TOKEN)

					 EXPECT

						 CASE (FUNCTION_TOKEN)
				 	 		 opts.Mapping_Function_Inverse = Parse_Function();
							 EXIT
						 END_CASE

						 OTHERWISE
							 UNGET
							 opts.Mapping_Function_Invert = (int)Parse_Float();
							 EXIT
						 END_CASE

					 END_EXPECT

				 END_CASE

				 OTHERWISE
					 UNGET
					 EXIT
				 END_CASE

			 END_EXPECT

       opts.Use_Mapping_Function = true;

			 Parse_End();

       END_CASE
#endif

     CASE (HF_GRAY_16_TOKEN)
       if (Allow_Float(1.0)>EPSILON)     
       {
         opts.Options |= HF_GRAY_16;
         opts.PaletteOption = GREY;        /* Force gray scale preview */
         opts.Output_File_Type |= HF_FTYPE;
       }
     END_CASE

     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT
   Parse_End();
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Object_Mods (OBJECT *Object)
{
   DBL V1, V2;
   VECTOR Min, Max;  
   VECTOR Local_Vector;
   MATRIX Local_Matrix;
   TRANSFORM Local_Trans;
   BBOX BBox;
   OBJECT *Sib;
   TEXTURE *Local_Texture;
   TEXTURE *Local_Int_Texture;
   MATERIAL Local_Material;
   OBJECT *Temp1_Object;
   OBJECT *Temp2_Object;
   COLOUR Local_Colour;

   EXPECT
     CASE(UV_MAPPING_TOKEN)
       /* if no texture than allow uv_mapping
          otherwise, warn user */
       if (Object->Texture == NULL)
       {
         Set_Flag(Object, UV_FLAG);
       }
       else
       {
         Error ("uv_mapping must be specified before texture.");
       }
     END_CASE

     CASE(SPLIT_UNION_TOKEN)
       if (Object->Methods != &CSG_Union_Methods)
       {
         Error("split_union found in non-union object.\n");
       }

       ((CSG*)Object)->do_split = (int)Parse_Float();
     END_CASE

     CASE(PHOTONS_TOKEN)
       Parse_Begin();
       EXPECT
         CASE(TARGET_TOKEN)
           Object->Ph_Density = Allow_Float(1.0);
           if (Object->Ph_Density > 0)
           {
             Set_Flag(Object,PH_TARGET_FLAG);
             CheckPassThru(Object, PH_TARGET_FLAG);
           }
           else
           {
             Clear_Flag(Object, PH_TARGET_FLAG);
           }
         END_CASE

         CASE(REFRACTION_TOKEN)
           if((int)Parse_Float())
           { 
             Set_Flag(Object, PH_RFR_ON_FLAG);
             Clear_Flag(Object, PH_RFR_OFF_FLAG);
             CheckPassThru(Object, PH_RFR_ON_FLAG);
           }
           else
           { 
             Clear_Flag(Object, PH_RFR_ON_FLAG);
             Set_Flag(Object, PH_RFR_OFF_FLAG);
           }
         END_CASE

         CASE(REFLECTION_TOKEN)
           if((int)Parse_Float())
           { 
             Set_Flag(Object, PH_RFL_ON_FLAG); 
             Clear_Flag(Object, PH_RFL_OFF_FLAG); 
           }
           else
           { 
             Clear_Flag(Object, PH_RFL_ON_FLAG); 
             Set_Flag(Object, PH_RFL_OFF_FLAG); 
           }
         END_CASE

         CASE(PASS_THROUGH_TOKEN)
           if((int)Allow_Float(1.0))
           {
             Set_Flag(Object, PH_PASSTHRU_FLAG);
             CheckPassThru(Object, PH_PASSTHRU_FLAG);
           }
           else
           {
             Clear_Flag(Object, PH_PASSTHRU_FLAG);
           }
         END_CASE

         CASE(COLLECT_TOKEN)
           Bool_Flag (Object, PH_IGNORE_PHOTONS_FLAG, !(Allow_Float(1.0) > 0.0));
         END_CASE

         OTHERWISE
           UNGET
           EXIT
         END_CASE
       END_EXPECT
       Parse_End();

     END_CASE

     CASE(CUTAWAY_TEXTURES_TOKEN)
       if (Object->Methods != &CSG_Intersection_Methods)
       {
         Error("cutaway_textures can only be used with intersection and difference.");
       }

       Set_Flag(Object, MULTITEXTURE_FLAG);
     END_CASE

     CASE_COLOUR
       Parse_Colour (Local_Colour);
       if (opts.Language_Version < 150)
       {
         if (Object->Texture != NULL)
         {
           if (Object->Texture->Type == PLAIN_PATTERN)
           {
             if (opts.Quality_Flags & Q_QUICKC)
             {
              Assign_Colour(Object->Texture->Pigment->Colour,Local_Colour);
              END_CASE
             }
           }
         }
       }
       Warning(0, "Quick color belongs in texture. Color ignored.");
     END_CASE

     CASE (TRANSLATE_TOKEN)
       Parse_Vector (Local_Vector);
       Compute_Translation_Transform(&Local_Trans, Local_Vector);
       Translate_Object (Object, Local_Vector, &Local_Trans);
     END_CASE

     CASE (ROTATE_TOKEN)
       Parse_Vector (Local_Vector);
       Compute_Rotation_Transform(&Local_Trans, Local_Vector);
       Rotate_Object (Object, Local_Vector, &Local_Trans);
     END_CASE

     CASE (SCALE_TOKEN)
       Parse_Scale_Vector (Local_Vector);
       Compute_Scaling_Transform(&Local_Trans, Local_Vector);
       Scale_Object (Object, Local_Vector, &Local_Trans);
     END_CASE

     CASE (TRANSFORM_TOKEN)
       Transform_Object(Object, Parse_Transform(&Local_Trans));
     END_CASE

     CASE (MATRIX_TOKEN)
       Parse_Matrix (Local_Matrix);
       Compute_Matrix_Transform(&Local_Trans, Local_Matrix);
       Transform_Object (Object, &Local_Trans);
     END_CASE

     CASE (BOUNDED_BY_TOKEN)
       Parse_Begin ();
       if (Object->Bound != NULL)
         if (Object->Clip == Object->Bound)
           Error ("Cannot add bounds after linking bounds and clips.");

       EXPECT
         CASE (CLIPPED_BY_TOKEN)
           if (Object->Bound != NULL)
             Error ("Cannot link clips with previous bounds.");
           Object->Bound = Object->Clip;
           EXIT
         END_CASE

         OTHERWISE
           UNGET
           Temp1_Object = Temp2_Object = Parse_Bound_Clip ();
           while (Temp2_Object->Sibling != NULL)
             Temp2_Object = Temp2_Object->Sibling;
           Temp2_Object->Sibling = Object->Bound;
           Object->Bound = Temp1_Object;
           EXIT
         END_CASE
       END_EXPECT

       Parse_End ();
     END_CASE

     CASE (CLIPPED_BY_TOKEN)
       Parse_Begin ();
       if (Object->Clip != NULL)
         if (Object->Clip == Object->Bound)
           Error ("Cannot add clips after linking bounds and clips.");

       EXPECT
         CASE (BOUNDED_BY_TOKEN)
           if (Object->Clip != NULL)
             Error ("Cannot link bounds with previous clips.");
           Object->Clip = Object->Bound;
           EXIT
         END_CASE

         OTHERWISE
           UNGET
           Temp1_Object = Temp2_Object = Parse_Bound_Clip ();
           while (Temp2_Object->Sibling != NULL)
             Temp2_Object = Temp2_Object->Sibling;
           Temp2_Object->Sibling = Object->Clip;
           Object->Clip = Temp1_Object;

           /* Compute quadric bounding box before transformations. [DB 8/94] */

           if (Object->Methods == &Quadric_Methods)
           {
             Make_Vector(Min, -BOUND_HUGE, -BOUND_HUGE, -BOUND_HUGE);
             Make_Vector(Max,  BOUND_HUGE,  BOUND_HUGE,  BOUND_HUGE);

             Compute_Quadric_BBox((QUADRIC *)Object, Min, Max);
           }
           EXIT
         END_CASE
       END_EXPECT

       Parse_End ();
     END_CASE

     CASE (TEXTURE_TOKEN)
       Object->Type |= TEXTURED_OBJECT;
       Parse_Begin ();
       Local_Texture = Parse_Texture ();
       Parse_End ();
       Link_Textures(&(Object->Texture), Local_Texture);
     END_CASE

     CASE (INTERIOR_TEXTURE_TOKEN)
       Object->Type |= TEXTURED_OBJECT;
       Parse_Begin ();
       Local_Int_Texture = Parse_Texture ();
       Parse_End ();
       Link_Textures(&(Object->Interior_Texture), Local_Int_Texture);
     END_CASE

     CASE (INTERIOR_TOKEN)
       Parse_Interior((INTERIOR **)(&Object->Interior));
     END_CASE

     CASE (MATERIAL_TOKEN)
       Local_Material.Texture  = Object->Texture;
       Local_Material.Interior_Texture  = Object->Interior_Texture;
       Local_Material.Interior = Object->Interior;
       Parse_Material(&Local_Material);
       Object->Texture  = Local_Material.Texture;
       if ( Object->Texture ) 
       {
         Object->Type |= TEXTURED_OBJECT;
       }
       Object->Interior_Texture  = Local_Material.Interior_Texture;
       Object->Interior = Local_Material.Interior;
     END_CASE

     CASE3 (PIGMENT_TOKEN, TNORMAL_TOKEN, FINISH_TOKEN)
       Object->Type |= TEXTURED_OBJECT;
       if (Object->Texture == NULL)
         Object->Texture = Copy_Textures(Default_Texture);
       else
         if (Object->Texture->Type != PLAIN_PATTERN)
           Link_Textures(&(Object->Texture), Copy_Textures(Default_Texture));
       UNGET
       EXPECT
         CASE (PIGMENT_TOKEN)
           Parse_Begin ();
           Parse_Pigment ( &(Object->Texture->Pigment) );
           Parse_End ();
         END_CASE

         CASE (TNORMAL_TOKEN)
           Parse_Begin ();
           Parse_Tnormal ( &(Object->Texture->Tnormal) );
           Parse_End ();
         END_CASE

         CASE (FINISH_TOKEN)
           Parse_Finish ( &(Object->Texture->Finish) );
         END_CASE

         OTHERWISE
           UNGET
           EXIT
         END_CASE
       END_EXPECT
     END_CASE

     CASE (INVERSE_TOKEN)
       if (Object->Type & PATCH_OBJECT)
         Warning (0, "Cannot invert a patch object.");
       Invert_Object (Object);
     END_CASE

     CASE (STURM_TOKEN)
       if (!(Object->Type & STURM_OK_OBJECT))
         Not_With ("sturm","this object");
       Bool_Flag (Object, STURM_FLAG, (Allow_Float(1.0) > 0.0));
     END_CASE

     /* Object-Ray Options
        Do not intersect with camera rays [ENB 9/97] */
     CASE (NO_IMAGE_TOKEN)
       Bool_Flag (Object, NO_IMAGE_FLAG, (Allow_Float(1.0) > 0.0));
     END_CASE

     /* Object-Ray Options
        Do not intersect with reflection rays [ENB 9/97] */
     CASE (NO_REFLECTION_TOKEN)
       Bool_Flag (Object, NO_REFLECTION_FLAG, (Allow_Float(1.0) > 0.0));
     END_CASE

     CASE (NO_SHADOW_TOKEN)
       Set_Flag(Object, NO_SHADOW_FLAG);
     END_CASE

#ifdef NO_RADIOSITY_PATCH /* Christoph Hormann January 2004 */
     CASE (NO_RADIOSITY_TOKEN)
       parseUnofficialFeature(110);
       Bool_Flag (Object, NO_RADIOSITY_FLAG, (Allow_Float(1.0) > 0.0));
     END_CASE
#endif

     CASE (LIGHT_SOURCE_TOKEN)
       Error("Light source must be defined using new syntax.");
     END_CASE

     CASE(HIERARCHY_TOKEN)
       if (!(Object->Type & HIERARCHY_OK_OBJECT))
         Not_With ("hierarchy", "this object");
       Bool_Flag (Object, HIERARCHY_FLAG, (Allow_Float(1.0) > 0.0));
     END_CASE

     CASE(HOLLOW_TOKEN)
       Bool_Flag (Object, HOLLOW_FLAG, (Allow_Float(1.0) > 0.0));
       Set_Flag (Object, HOLLOW_SET_FLAG);
       if ((Object->Methods == &CSG_Intersection_Methods) ||
           (Object->Methods == &CSG_Merge_Methods) ||
           (Object->Methods == &CSG_Union_Methods))
       {
         Set_CSG_Children_Flag(Object, Test_Flag(Object, HOLLOW_FLAG), HOLLOW_FLAG, HOLLOW_SET_FLAG);
       }
     END_CASE

     /* NK 1998 double_illuminate */
     CASE(DOUBLE_ILLUMINATE_TOKEN)
       Bool_Flag (Object, DOUBLE_ILLUMINATE_FLAG, (Allow_Float(1.0) > 0.0));
       if ((Object->Methods == &CSG_Intersection_Methods) ||
           (Object->Methods == &CSG_Merge_Methods) ||
           (Object->Methods == &CSG_Union_Methods))
       {
         Set_CSG_Tree_Flag(Object, DOUBLE_ILLUMINATE_FLAG,Test_Flag(Object, DOUBLE_ILLUMINATE_FLAG));
       }     
     END_CASE
	 /* NK ---- */

     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   /*
    * Assign bounding objects' bounding box to object
    * if object's bounding box is larger. [DB 9/94]
    */

   if (Object->Bound != NULL)
   {
     /* Get bounding objects bounding box. */

     Make_Vector(Min, -BOUND_HUGE, -BOUND_HUGE, -BOUND_HUGE);
     Make_Vector(Max,  BOUND_HUGE,  BOUND_HUGE,  BOUND_HUGE);

     for (Sib = Object->Bound; Sib != NULL; Sib = Sib->Sibling)
     {
       if (!Test_Flag(Sib, INVERTED_FLAG))
       {
         Min[X] = max(Min[X], (DBL)(Sib->BBox.Lower_Left[X]));
         Min[Y] = max(Min[Y], (DBL)(Sib->BBox.Lower_Left[Y]));
         Min[Z] = max(Min[Z], (DBL)(Sib->BBox.Lower_Left[Z]));
         Max[X] = min(Max[X], (DBL)(Sib->BBox.Lower_Left[X] + Sib->BBox.Lengths[X]));
         Max[Y] = min(Max[Y], (DBL)(Sib->BBox.Lower_Left[Y] + Sib->BBox.Lengths[Y]));
         Max[Z] = min(Max[Z], (DBL)(Sib->BBox.Lower_Left[Z] + Sib->BBox.Lengths[Z]));
       }
     }

     Make_BBox_from_min_max(BBox, Min, Max);

     /* Get bounding boxes' volumes. */

     BOUNDS_VOLUME(V1, BBox);
     BOUNDS_VOLUME(V2, Object->BBox);

     if (V1 < V2)
     {
       Object->BBox = BBox;
     }
   }

   /*
    * Assign clipping objects' bounding box to object
    * if object's bounding box is larger. [DB 9/94]
    */

   if (Object->Clip != NULL)
   {
     /* Get clipping objects bounding box. */

     Make_Vector(Min, -BOUND_HUGE, -BOUND_HUGE, -BOUND_HUGE);
     Make_Vector(Max,  BOUND_HUGE,  BOUND_HUGE,  BOUND_HUGE);

     for (Sib = Object->Clip; Sib != NULL; Sib = Sib->Sibling)
     {
       if (!Test_Flag(Sib, INVERTED_FLAG))
       {
         Min[X] = max(Min[X], (DBL)(Sib->BBox.Lower_Left[X]));
         Min[Y] = max(Min[Y], (DBL)(Sib->BBox.Lower_Left[Y]));
         Min[Z] = max(Min[Z], (DBL)(Sib->BBox.Lower_Left[Z]));
         Max[X] = min(Max[X], (DBL)(Sib->BBox.Lower_Left[X] + Sib->BBox.Lengths[X]));
         Max[Y] = min(Max[Y], (DBL)(Sib->BBox.Lower_Left[Y] + Sib->BBox.Lengths[Y]));
         Max[Z] = min(Max[Z], (DBL)(Sib->BBox.Lower_Left[Z] + Sib->BBox.Lengths[Z]));
       }
     }

     Make_BBox_from_min_max(BBox, Min, Max);

     /* Get bounding boxes' volumes. */

     BOUNDS_VOLUME(V1, BBox);
     BOUNDS_VOLUME(V2, Object->BBox);

     if (V1 < V2)
     {
       Object->BBox = BBox;
     }
   }

   if((Object->Texture ==NULL)&&(Object->Interior_Texture != NULL))
       Error("Interior texture requires an exterior texture.");

   Parse_End ();
  }



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Matrix(MATRIX Matrix)
{
  int i, j;

  EXPECT
    CASE (LEFT_ANGLE_TOKEN)
      Matrix[0][0] = Parse_Float();
      for (i = 0; i < 4; i++)
      {
        for (j = !i ? 1 : 0; j < 3; j++)
        {
          Parse_Comma();

          Matrix[i][j] = Parse_Float();
        }

        Matrix[i][3] = (i != 3 ? 0.0 : 1.0);
      }
      GET (RIGHT_ANGLE_TOKEN);

      /* Check to see that we aren't scaling any dimension by zero */
      for (i = 0; i < 3; i++)
      {
        if (fabs(Matrix[0][i]) < EPSILON && fabs(Matrix[1][i]) < EPSILON &&
            fabs(Matrix[2][i]) < EPSILON)
        {
          Warning(0,"Illegal matrix column: Scale by 0.0. Changed to 1.0.");
          Matrix[i][i] = 1.0;
        }
      }
      EXIT
    END_CASE

    OTHERWISE
      Parse_Error (LEFT_ANGLE_TOKEN);
    END_CASE
  END_EXPECT
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

TRANSFORM *Parse_Transform(TRANSFORM *Trans)
{
  Get_Token();
  if(Token.Token_Id == TRANSFORM_ID_TOKEN)
  {
    /* using old "transform TRANS_IDENT" syntax */
    if(Trans == NULL)
       Trans=Create_Transform();
    else
    {
       MIdentity (Trans->matrix);
       MIdentity (Trans->inverse);
    }
    Compose_Transforms(Trans, (TRANSFORM *)Token.Data);
  }
  else
  {
    /* using new "transform {TRANS}" syntax */
    Unget_Token();
    Trans = Parse_Transform_Block(Trans);
  }
  return Trans;
}
/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

TRANSFORM *Parse_Transform_Block(TRANSFORM *New)
{
   MATRIX Local_Matrix;
   TRANSFORM Local_Trans;
   VECTOR Local_Vector;
   bool isInverse = false;

#ifdef NORMAL_TRANSFORM_PATCH            /* Christoph Hormann January 2003 */
   bool isNormal = false;
#endif

   Parse_Begin();
   if(New == NULL)
      New = Create_Transform();
   else
   {
      MIdentity (New->matrix);
      MIdentity (New->inverse);
   }

   EXPECT
     CASE(INVERSE_TOKEN)
       isInverse = true;
     END_CASE

     CASE(TRANSFORM_ID_TOKEN)
       Compose_Transforms(New, (TRANSFORM *)Token.Data);
     END_CASE

     CASE (TRANSFORM_TOKEN)
       Compose_Transforms(New, Parse_Transform(&Local_Trans));
     END_CASE

     CASE (TRANSLATE_TOKEN)
       Parse_Vector(Local_Vector);
       Compute_Translation_Transform(&Local_Trans, Local_Vector);
       Compose_Transforms(New, &Local_Trans);
     END_CASE

     CASE (ROTATE_TOKEN)
       Parse_Vector(Local_Vector);
       Compute_Rotation_Transform(&Local_Trans, Local_Vector);
       Compose_Transforms(New, &Local_Trans);
     END_CASE

     CASE (SCALE_TOKEN)
       Parse_Scale_Vector(Local_Vector);
       Compute_Scaling_Transform(&Local_Trans, Local_Vector);
       Compose_Transforms(New, &Local_Trans);
     END_CASE

     CASE (MATRIX_TOKEN)
       Parse_Matrix(Local_Matrix);
       Compute_Matrix_Transform(&Local_Trans, Local_Matrix);
       Compose_Transforms(New, &Local_Trans);
     END_CASE

#ifdef NORMAL_TRANSFORM_PATCH            /* Christoph Hormann January 2003 */
     CASE (TNORMAL_TOKEN)
			 parseUnofficialFeature(110);
       isNormal = true;
     END_CASE
#endif

     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   Parse_End();

   if(isInverse == true)
   {
       MInvers(New->matrix, New->matrix);
       MInvers(New->inverse, New->inverse);
   }

#ifdef NORMAL_TRANSFORM_PATCH            /* Christoph Hormann January 2003 */
   if(isNormal == true)
   {
     /*
       For transforming normal vectors we need the transpose of the inverse of the
       original matrix.
     */

     POV_MEMCPY(New->matrix, New->inverse, sizeof(MATRIX));
     MTranspose (New->matrix);
     MInvers(New->inverse, New->matrix);

   }
#endif

   return (New);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

OBJECT *Parse_Bound_Clip ()
{
   VECTOR Local_Vector;
   MATRIX Local_Matrix;
   TRANSFORM Local_Trans;
   OBJECT *First, *Current, *Prev;

   First = Prev = NULL;

   while ((Current = Parse_Object ()) != NULL)
     {
      if (Current->Type & (TEXTURED_OBJECT+PATCH_OBJECT))
#ifdef INSIDE_FUNCTION_FOR_MESHES_PATCH
        if (Current->Methods != &Mesh_Methods || !((MESH*)Current)->has_inside_vector)
#endif
        Error ("Illegal texture or patch in clip, bound or object pattern.");
      if (First == NULL)
        First = Current;
      if (Prev != NULL)
        Prev->Sibling = Current;
      Prev = Current;
     }

   EXPECT
     CASE (TRANSLATE_TOKEN)
       Parse_Vector (Local_Vector);
       Compute_Translation_Transform(&Local_Trans, Local_Vector);
       for (Current = First; Current != NULL; Current = Current->Sibling)
       {
         Translate_Object (Current, Local_Vector, &Local_Trans);
       }
     END_CASE

     CASE (ROTATE_TOKEN)
       Parse_Vector (Local_Vector);
       Compute_Rotation_Transform(&Local_Trans, Local_Vector);
       for (Current = First; Current != NULL; Current = Current->Sibling)
       {
         Rotate_Object (Current, Local_Vector, &Local_Trans);
       }
     END_CASE

     CASE (SCALE_TOKEN)
       Parse_Scale_Vector (Local_Vector);
       Compute_Scaling_Transform(&Local_Trans, Local_Vector);
       for (Current = First; Current != NULL; Current = Current->Sibling)
       {
         Scale_Object (Current, Local_Vector, &Local_Trans);
       }
     END_CASE

     CASE (TRANSFORM_TOKEN)
       Parse_Transform(&Local_Trans);

       for (Current = First; Current != NULL; Current = Current->Sibling)
       {
         Transform_Object (Current, &Local_Trans);
       }
     END_CASE

     CASE (MATRIX_TOKEN)
       Parse_Matrix (Local_Matrix);
       Compute_Matrix_Transform(&Local_Trans, Local_Matrix);
       for (Current = First; Current != NULL; Current = Current->Sibling)
       {
         Transform_Object (Current, &Local_Trans);
       }
     END_CASE

     OTHERWISE
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   if (First==NULL)
   {
      Expectation_Error("object");
   }

   return (First);
  }


/*****************************************************************************
*
* FUNCTION
*
*   Parse_Three_UVCoords
*
* INPUT
*
* OUTPUT
*
*   UV1..UV3 are the uv
*
* RETURNS
*
*   1 for successful read, 0 if UV_VECTORS_TOKEN not found
*
* AUTHOR
*
*   Nathan Kopp
*
* DESCRIPTION
*
*   Look for UV_VECTORS_TOKEN and then read in three UV coordinates
*
******************************************************************************/

static int Parse_Three_UVCoords(UV_VECT UV1, UV_VECT UV2, UV_VECT UV3)
{
  int Return_Value;

  EXPECT
    CASE(UV_VECTORS_TOKEN)
      Parse_UV_Vect(UV1);  Parse_Comma();
      Parse_UV_Vect(UV2);  Parse_Comma();
      Parse_UV_Vect(UV3);

      Return_Value = 1;
      EXIT
    END_CASE

    OTHERWISE
      UV1[0] = UV1[1] = 0.0;
      UV2[0] = UV2[1] = 0.0;
      UV3[0] = UV3[1] = 0.0;
      Return_Value = 0;
      UNGET
      EXIT
    END_CASE

  END_EXPECT

  return(Return_Value);
}


/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Comma (void)
{
   Get_Token();
   if (Token.Token_Id != COMMA_TOKEN)
   {
      UNGET;
   }
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Semi_Colon (bool force_semicolon)
{
   Get_Token();
   if (Token.Token_Id != SEMI_COLON_TOKEN)
   {
      UNGET;
      if ((opts.Language_Version >= 350) && (force_semicolon == true))
      {
         Error("All #declares of float, vector, and color require semi-colon ';' at end if the\n"
               "language version is set to 3.5 or higher.\n"
               "Either add the semi-colon or set the language version to 3.1 or lower.");
      }
      else if (opts.Language_Version >= 310)
      {
         PossibleError("All #version and #declares of float, vector, and color require semi-colon ';' at end.");
      }
   }
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Coeffs(int order, DBL *Coeffs)
{
   int i;

   EXPECT
     CASE (LEFT_ANGLE_TOKEN)
       Coeffs[0] = Parse_Float();
       for (i = 1; i < term_counts(order); i++)
         {
          Parse_Comma();
          Coeffs[i] = Parse_Float();
         }
       GET (RIGHT_ANGLE_TOKEN);
       EXIT
     END_CASE

     OTHERWISE
       Parse_Error (LEFT_ANGLE_TOKEN);
     END_CASE
   END_EXPECT
  }


/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static OBJECT *Parse_Object_Id ()
{
   OBJECT *Object;

   EXPECT
     CASE (OBJECT_ID_TOKEN)
       Warn_State(OBJECT_ID_TOKEN, OBJECT_TOKEN);
       Object = Copy_Object((OBJECT *) Token.Data);
       Parse_Object_Mods (Object);
       EXIT
     END_CASE

     OTHERWISE
       Object = NULL;
       UNGET
       EXIT
     END_CASE
   END_EXPECT

   return (Object);
}


/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/
#ifdef SET_PATCH
void Parse_Declare(bool is_local, bool after_hash,bool Modify_Only)
#else
void Parse_Declare(bool is_local, bool after_hash)
#endif

{
  int Previous=-1;
  int Local_Index
                  #ifdef AVOID_MIGHT_BE_USED_UNINITIALIZED_WARNINGS_PATCH
                  = 0
                  #endif
                  ,Local_Flag;
  SYM_ENTRY *Temp_Entry = NULL;
  bool allow_redefine = true;

  Ok_To_Declare = false;

  if ((opts.Language_Version >= 350) && (after_hash == false))
  {
    PossibleError("'declare' should be changed to '#declare'.\n"
                  "Future versions may not support 'declare' and may require '#declare'.");
  }
#ifdef SET_PATCH
  Local_Flag = (Token.Token_Id==LOCAL_TOKEN );
  if(Local_Flag)
  {
     Local_Index=Table_Index;
  }
  else
  {
     if(Token.Token_Id==SET_TOKEN || Modify_Only == true)
     {
      /*Search for most recent table with an identifier with this name
      Parts copied from Read_Symbol*/
      int Idx;
      int Found = false;

      Get_Token();
      /* Search tables from newest to oldest */
      for(Idx=Table_Index; Idx > 0; Idx--)
      {
          if(Find_Symbol(Idx, Token.Token_String) != NULL)
          {
            Local_Index = Idx;
            Found = true;
            break;
          }
      }
      if(Found != true)
      {
        Error("#set cannot assign to uninitialized identifier %s.\n", Token.Token_String);
      }
      Unget_Token();
     }
     else
     {
      Local_Index=1;
     }
  }
 #else
  Local_Flag=(Token.Token_Id==LOCAL_TOKEN);
  if (Local_Flag)
  {
     Local_Index=Table_Index;
  }
  else
  {
     Local_Index=1;
  }

#endif

  LValue_Ok = true;

  EXPECT
    CASE (IDENTIFIER_TOKEN)
      allow_redefine  = !Token.is_array_elem;
      Temp_Entry = Add_Symbol (Local_Index,Token.Token_String,IDENTIFIER_TOKEN);
      Token.NumberPtr = &(Temp_Entry->Token_Number);
      Token.DataPtr   = &(Temp_Entry->Data);
      Previous        = Token.Token_Id;
      EXIT
    END_CASE

    CASE2 (FUNCT_ID_TOKEN, VECTFUNCT_ID_TOKEN)
      if((!Token.is_array_elem) || (*(Token.DataPtr) != NULL))
        Error("Redeclaring functions is not allowed - #undef the function first!");
      // fall through

    // These are also used in Parse_Directive UNDEF_TOKEN section and Parse_Macro,
    // and all three functions should accept exactly the same identifiers! [trf]
    CASE4 (TNORMAL_ID_TOKEN, FINISH_ID_TOKEN, TEXTURE_ID_TOKEN, OBJECT_ID_TOKEN)
    CASE4 (COLOUR_MAP_ID_TOKEN, TRANSFORM_ID_TOKEN, CAMERA_ID_TOKEN, PIGMENT_ID_TOKEN)
    CASE4 (SLOPE_MAP_ID_TOKEN, NORMAL_MAP_ID_TOKEN, TEXTURE_MAP_ID_TOKEN, COLOUR_ID_TOKEN)
    CASE4 (PIGMENT_MAP_ID_TOKEN, MEDIA_ID_TOKEN, STRING_ID_TOKEN, INTERIOR_ID_TOKEN)
    CASE4 (DENSITY_MAP_ID_TOKEN, ARRAY_ID_TOKEN, DENSITY_ID_TOKEN, UV_ID_TOKEN)
    CASE4 (VECTOR_4D_ID_TOKEN, RAINBOW_ID_TOKEN, FOG_ID_TOKEN, SKYSPHERE_ID_TOKEN)
    CASE2 (MATERIAL_ID_TOKEN, SPLINE_ID_TOKEN )
      allow_redefine  = !Token.is_array_elem;
      if (Local_Flag && (Token.Table_Index != Table_Index))
      {
        Temp_Entry = Add_Symbol (Local_Index,Token.Token_String,IDENTIFIER_TOKEN);
        Token.NumberPtr = &(Temp_Entry->Token_Number);
        Token.DataPtr   = &(Temp_Entry->Data);
        Previous        = IDENTIFIER_TOKEN;
      }
      else
      {
        Previous        = Token.Token_Id;
      }
      EXIT
    END_CASE

    CASE (EMPTY_ARRAY_TOKEN)
      allow_redefine  = !Token.is_array_elem;
      Previous = Token.Token_Id;
      EXIT
    END_CASE

    CASE2 (VECTOR_FUNCT_TOKEN, FLOAT_FUNCT_TOKEN)
      allow_redefine  = !Token.is_array_elem;
      switch(Token.Function_Id)
        {
         case VECTOR_ID_TOKEN:
         case FLOAT_ID_TOKEN:
           if (Local_Flag && (Token.Table_Index != Table_Index))
           {
              Temp_Entry = Add_Symbol (Local_Index,Token.Token_String,IDENTIFIER_TOKEN);
              Token.NumberPtr = &(Temp_Entry->Token_Number);
              Token.DataPtr   = &(Temp_Entry->Data);
           }
           Previous           = Token.Function_Id;
           break;

         default:
           Parse_Error(IDENTIFIER_TOKEN);
           break;
        }
      EXIT
    END_CASE

    OTHERWISE
      allow_redefine  = !Token.is_array_elem;
      Parse_Error(IDENTIFIER_TOKEN);
    END_CASE
  END_EXPECT

  LValue_Ok = false;

  GET (EQUALS_TOKEN);
  Ok_To_Declare = true;  
  if (!Parse_RValue (Previous, Token.NumberPtr, Token.DataPtr, Temp_Entry, false, true, is_local, allow_redefine, MAX_NUMBER_OF_TABLES))
  {
    Expectation_Error("RValue to declare");
  }
  if ( after_hash ) 
  {
    Ok_To_Declare = false;
    ALLOW( SEMI_COLON_TOKEN );
    Ok_To_Declare = true;
  }

}

int Parse_RValue (int Previous, int *NumberPtr, void **DataPtr, SYM_ENTRY *sym, bool ParFlag, bool SemiFlag, bool is_local, bool allow_redefine, int old_table_index)
{
  EXPRESS Local_Express;
  COLOUR *Local_Colour;
  PIGMENT *Local_Pigment;
  TNORMAL *Local_Tnormal;
  FINISH *Local_Finish;
  TEXTURE *Local_Texture, *Temp_Texture;
  TRANSFORM *Local_Trans;
  OBJECT *Local_Object;
  CAMERA *Local_Camera;
  IMEDIA *Local_Media;
  PIGMENT *Local_Density;
  INTERIOR *Local_Interior;
  MATERIAL *Local_Material;
  void *Temp_Data;
  POV_PARAM *New_Par;
  int Found=true;
  int Temp_Count=30000;
  int Old_Ok=Ok_To_Declare;
  int Terms;
  bool function_identifier ;
  bool callable_identifier ;
  bool had_callable_identifier ;
  
  EXPECT
    CASE4 (TNORMAL_ID_TOKEN, FINISH_ID_TOKEN, TEXTURE_ID_TOKEN, OBJECT_ID_TOKEN)
    CASE4 (COLOUR_MAP_ID_TOKEN, TRANSFORM_ID_TOKEN, CAMERA_ID_TOKEN, PIGMENT_ID_TOKEN)
    CASE4 (SLOPE_MAP_ID_TOKEN,NORMAL_MAP_ID_TOKEN,TEXTURE_MAP_ID_TOKEN,ARRAY_ID_TOKEN)
    CASE4 (PIGMENT_MAP_ID_TOKEN, MEDIA_ID_TOKEN,INTERIOR_ID_TOKEN,DENSITY_ID_TOKEN)
    CASE4 (DENSITY_MAP_ID_TOKEN, RAINBOW_ID_TOKEN, FOG_ID_TOKEN, SKYSPHERE_ID_TOKEN)
    CASE2 (MATERIAL_ID_TOKEN, STRING_ID_TOKEN) 
      if ((ParFlag) && (Token.Table_Index <= old_table_index))
      {
        // pass by reference
        New_Par            = (POV_PARAM *)POV_MALLOC(sizeof(POV_PARAM),"parameter");
        New_Par->NumberPtr = Token.NumberPtr;
        New_Par->DataPtr   = Token.DataPtr;
        New_Par->Table_Index = Token.Table_Index;
        *NumberPtr = PARAMETER_ID_TOKEN;
        *DataPtr   = (void *)New_Par;
      }
      else
      {
        // pass by value
        Temp_Data  = (void *) Copy_Identifier((void *)*Token.DataPtr,*Token.NumberPtr);
        *NumberPtr = *Token.NumberPtr;
        Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
        *DataPtr   = Temp_Data;
      }
      EXIT
    END_CASE

    CASE (IDENTIFIER_TOKEN)
      if (ParFlag)
      {
         Error("Cannot pass uninitialized identifier as macro parameter.\nInitialize identifier first.");
      }
      else
      {
         Error("Cannot assign uninitialized identifier.");
      }
      EXIT
    END_CASE

    CASE_COLOUR
      if((Token.Token_Id != COLOUR_ID_TOKEN) || (opts.Language_Version < 350))
      {
         Local_Colour  = Create_Colour();
         Ok_To_Declare = false;
         Parse_Colour (*Local_Colour);
         if (SemiFlag)
         {
            Parse_Semi_Colon(true);
         }
         Ok_To_Declare = true;
         *NumberPtr    = COLOUR_ID_TOKEN;
         Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
         *DataPtr      = (void *) Local_Colour;
         EXIT
         END_CASE
      }
      // intentional to allow color dot expressions as macro parameters if #version is 3.5 or higher [trf]

    CASE_VECTOR
      // It seems very few people understand what is going on here, so let me try to
      // explain it. All comments below are mine and they are based on how I think it
      // works and understand it. As I didn't write most of the code I cannot really
      // tell for sure, so if anybody finds incorrect comments please let me know!
      // BTW, when saying #declare it always implies "or #local" :-)

      // determine the type of the first identifier
      function_identifier = (Token.Token_Id==FUNCT_ID_TOKEN) || (Token.Token_Id==VECTFUNCT_ID_TOKEN);
      callable_identifier = (Token.Token_Id==FUNCT_ID_TOKEN) || (Token.Token_Id==VECTFUNCT_ID_TOKEN) || (Token.Token_Id==SPLINE_ID_TOKEN);

      // don't allow #declares from here
      Ok_To_Declare = false;

      // if what follows could be a function/spline call or
      // is a macro parameter taking a float, vector or ids
      // of a float, vector or color then count the tokens
      // found between now and the time when the function
      // Parse_Unknown_Vector returns
      if (callable_identifier || (ParFlag && 
          (((Token.Token_Id==FLOAT_FUNCT_TOKEN) && (Token.Function_Id==FLOAT_ID_TOKEN)) ||
           ((Token.Token_Id==VECTOR_FUNCT_TOKEN) &&  (Token.Function_Id==VECTOR_ID_TOKEN)) ||
           (Token.Token_Id==VECTOR_4D_ID_TOKEN) || (Token.Token_Id==UV_ID_TOKEN) || (Token.Token_Id==COLOUR_ID_TOKEN))))
      {
         Temp_Count = token_count;
      }

      // assume no callable identifier (that is a function or spline identifier) has been found
      had_callable_identifier = false;

      // parse the expression and determine if it was a callable identifier
      Terms = Parse_Unknown_Vector (Local_Express, true, &had_callable_identifier);

      // if in a #declare force a semicolon at the end
      if (SemiFlag)
         Parse_Semi_Colon(true);

      // get the number of tokens found
      Temp_Count -= token_count;

#ifdef MACRO_PARAMETER_BUGFIX_PATCH
      // no tokens have been found or a fucntion call had no parameters in parenthesis
      if (!((Temp_Count==-1) || (Temp_Count==2500)) && had_callable_identifier)
         Error("Identifier expected, incomplete function call or spline call found instead.");

      // only one identifier token has been found so pass it by reference
      if (((Temp_Count==-1) || (Temp_Count==2500)) && (Token.Table_Index <= old_table_index))
#else
      // no tokens have been found or a fucntion call had no parameters in parenthesis
      if (!((Temp_Count==-1) || (Temp_Count==1000)) && had_callable_identifier)
         Error("Identifier expected, incomplete function call or spline call found instead.");

      // only one identifier token has been found so pass it by reference
      if (((Temp_Count==-1) || (Temp_Count==1000)) && (Token.Table_Index <= old_table_index))
#endif
      {
         // It is important that functions are passed by value and not by reference! [trf]
         if(!(ParFlag) || (ParFlag && function_identifier))
         {
           // pass by value
           Temp_Data  = (void *) Copy_Identifier((void *)*Token.DataPtr,*Token.NumberPtr);
           *NumberPtr = *Token.NumberPtr;
           Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
           *DataPtr   = Temp_Data;
         }
         else
         {
           // pass by reference
           New_Par            = (POV_PARAM *)POV_MALLOC(sizeof(POV_PARAM),"parameter");
           New_Par->NumberPtr = Token.NumberPtr;
           New_Par->DataPtr   = Token.DataPtr;
           New_Par->Table_Index = Token.Table_Index;

           *NumberPtr = PARAMETER_ID_TOKEN;
           *DataPtr   = (void *)New_Par;
         }
      }
      else // an expression has been found, so create a new identifier
      {
         switch(Terms)
         {
           case 1:
            *NumberPtr = FLOAT_ID_TOKEN;
            Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
            *DataPtr   = (void *) Create_Float();
            *((DBL *)*DataPtr)  = Local_Express[X];
            break;
            
           case 2:
            *NumberPtr = UV_ID_TOKEN;
            Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
            *DataPtr   = (void *) Create_UV_Vect();
            Assign_UV_Vect((DBL *)*DataPtr, Local_Express);
            break;
            
           case 3:
            *NumberPtr = VECTOR_ID_TOKEN;
            Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
            *DataPtr   = (void *) Create_Vector();
            Assign_Vector((DBL *)*DataPtr, Local_Express);
            break;
            
           case 4:
            *NumberPtr = VECTOR_4D_ID_TOKEN;
            Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
            *DataPtr   = (void *) Create_Vector_4D();
            Assign_Vector_4D((DBL *)*DataPtr, Local_Express);
            break;
            
           case 5:
            *NumberPtr    = COLOUR_ID_TOKEN;
            Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
            *DataPtr      = (void *) Create_Colour();
            Assign_Colour_Express((COLC*)(*DataPtr), Local_Express);
            break;
         }
      } 

      // allow #declares again
      Ok_To_Declare = true;
      EXIT
    END_CASE

    CASE (PIGMENT_TOKEN)
      Local_Pigment = Copy_Pigment((Default_Texture->Pigment));
      Parse_Begin ();
      Parse_Pigment (&Local_Pigment);
      Parse_End ();
      *NumberPtr = PIGMENT_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = (void *)Local_Pigment;
      EXIT
    END_CASE

    CASE (TNORMAL_TOKEN)
      Local_Tnormal = Copy_Tnormal((Default_Texture->Tnormal));
      Parse_Begin ();
      Parse_Tnormal (&Local_Tnormal);
      Parse_End ();
      *NumberPtr = TNORMAL_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = (void *) Local_Tnormal;
      EXIT
    END_CASE

    CASE (FINISH_TOKEN)
      Local_Finish = Copy_Finish((Default_Texture->Finish));
      Parse_Finish (&Local_Finish);
      *NumberPtr = FINISH_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = (void *) Local_Finish;
      EXIT
    END_CASE

    CASE (CAMERA_TOKEN)
      Local_Camera = Copy_Camera(Default_Camera);
      Parse_Camera (&Local_Camera);
      *NumberPtr = CAMERA_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = (void *) Local_Camera;
      EXIT
    END_CASE

    CASE (TEXTURE_TOKEN)
      Parse_Begin ();
      Local_Texture = Parse_Texture ();
      Parse_End ();
      Temp_Texture=NULL;
      Link_Textures(&Temp_Texture, Local_Texture);
      Ok_To_Declare = false;
      EXPECT
        CASE (TEXTURE_TOKEN)
          Parse_Begin ();
          Local_Texture = Parse_Texture ();
          Parse_End ();
          Link_Textures(&Temp_Texture, Local_Texture);
        END_CASE

        OTHERWISE
          UNGET
          EXIT
        END_CASE
      END_EXPECT

      *NumberPtr    = TEXTURE_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr      = (void *)Temp_Texture;
      Ok_To_Declare = true;
      EXIT
    END_CASE

    CASE (COLOUR_MAP_TOKEN)
      Temp_Data=(void *) Parse_Colour_Map ();
      *NumberPtr = COLOUR_MAP_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (PIGMENT_MAP_TOKEN)
      Temp_Data  = (void *) Parse_Blend_Map (PIGMENT_TYPE,NO_PATTERN);
      *NumberPtr = PIGMENT_MAP_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (SPLINE_TOKEN)
      Experimental_Flag |= EF_SPLINE;
      Parse_Begin();
      Temp_Data=(char *) Parse_Spline();
      Parse_End();
      *NumberPtr = SPLINE_ID_TOKEN;
      Test_Redefine( Previous, NumberPtr, *DataPtr , allow_redefine);
      *DataPtr = (void *)Temp_Data;
      EXIT
    END_CASE

    CASE (DENSITY_MAP_TOKEN)
      Temp_Data  = (void *) Parse_Blend_Map (DENSITY_TYPE,NO_PATTERN);
      *NumberPtr = DENSITY_MAP_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (SLOPE_MAP_TOKEN)
      Temp_Data  = (void *) Parse_Blend_Map (SLOPE_TYPE,NO_PATTERN);
      *NumberPtr = SLOPE_MAP_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (TEXTURE_MAP_TOKEN)
      Temp_Data  = (void *) Parse_Blend_Map (TEXTURE_TYPE,NO_PATTERN);
      *NumberPtr = TEXTURE_MAP_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (NORMAL_MAP_TOKEN)
      Temp_Data  = (void *) Parse_Blend_Map (NORMAL_TYPE,NO_PATTERN);
      *NumberPtr = NORMAL_MAP_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (RAINBOW_TOKEN)
      Temp_Data  = (void *) Parse_Rainbow();
      *NumberPtr = RAINBOW_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (FOG_TOKEN)
      Temp_Data  = (void *) Parse_Fog();
      *NumberPtr = FOG_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (MEDIA_TOKEN)
      Local_Media = NULL;
      Parse_Media(&Local_Media);
      Temp_Data  = (void *)Local_Media;
      *NumberPtr = MEDIA_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (DENSITY_TOKEN)
      Local_Density = NULL;
      Parse_Begin ();
      Parse_Media_Density_Pattern (&Local_Density);
      Parse_End ();
      *NumberPtr = DENSITY_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = (void *)Local_Density;
      EXIT
    END_CASE

    CASE (INTERIOR_TOKEN)
      Local_Interior = NULL;
      Parse_Interior(&Local_Interior);
      Temp_Data  = (void *)Local_Interior;
      *NumberPtr = INTERIOR_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (MATERIAL_TOKEN)
      Local_Material = Create_Material();
      Parse_Material(Local_Material);
      Temp_Data  = (void *)Local_Material;
      *NumberPtr = MATERIAL_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (SKYSPHERE_TOKEN)
      Temp_Data  = (void *) Parse_Skysphere();
      *NumberPtr = SKYSPHERE_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (FUNCTION_TOKEN)
      // Do NOT allow to redefine functions! [trf]
      //   #declare foo = function(x) { x }
      //   #declare foo = function(x) { foo(x) } // Error!
      // Reason: Code like this would be unreadable but possible. Is it 
      // a recursive function or not? - It is not recursive because the
      // foo in the second line refers to the first function, which is
      // not logical. Further, recursion is not supported in POV-Ray 3.5
      // anyway. However, allowing such code now would cause problems
      // implementing recursive functions after POV-Ray 3.5!
      if(sym != NULL)
         Temp_Data  = (void *)Parse_DeclareFunction(NumberPtr, sym->Token_Name, is_local);
      else
         Temp_Data  = (void *)Parse_DeclareFunction(NumberPtr, NULL, is_local);
      Test_Redefine(Previous, NumberPtr, *DataPtr, false);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (TRANSFORM_TOKEN)
      Local_Trans = Parse_Transform ();
      *NumberPtr  = TRANSFORM_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr    = (void *) Local_Trans;
      EXIT
    END_CASE

    #ifdef NTH_OUTPUT_FILE_NAME_PATCH
      CASE (OUTPUT_FILENAME_TOKEN)
      parseUnofficialFeature(110);
    #endif
    #ifdef TIME_DATE_PATCH
      CASE(GET_DATE_TOKEN) /*YS get date and time string*/
      parseUnofficialFeature(100);
    #endif
    CASE5 (STRING_LITERAL_TOKEN,CHR_TOKEN,SUBSTR_TOKEN,STR_TOKEN,VSTR_TOKEN)
    CASE3 (CONCAT_TOKEN,STRUPR_TOKEN,STRLWR_TOKEN)
      UNGET
      Temp_Data  = Parse_String();
      *NumberPtr = STRING_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    CASE (ARRAY_TOKEN)
      Temp_Data  = (void *) Parse_Array_Declare();
      *NumberPtr = ARRAY_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr   = Temp_Data;
      EXIT
    END_CASE

    OTHERWISE
      UNGET
      Local_Object = Parse_Object ();
      Found=(Local_Object!=NULL);
      *NumberPtr   = OBJECT_ID_TOKEN;
      Test_Redefine(Previous,NumberPtr,*DataPtr, allow_redefine);
      *DataPtr     = (void *) Local_Object;
      EXIT
    END_CASE

  END_EXPECT
  
  Ok_To_Declare=Old_Ok;
  return(Found);
}

void Destroy_Ident_Data(void *Data, int Type)
{
  int i;
  POV_ARRAY *a;
  DATA_FILE *Temp_File;

  if(Data == NULL)
     return;

  switch(Type)
  {
     case COLOUR_ID_TOKEN:
       Destroy_Colour((COLOUR *)Data);
       break;
     case VECTOR_ID_TOKEN:
       Destroy_Vector((VECTOR *)Data);
       break;
     case UV_ID_TOKEN:
       Destroy_UV_Vect((UV_VECT *)Data);
       break;
     case VECTOR_4D_ID_TOKEN:
       Destroy_Vector_4D((VECTOR_4D *)Data);
       break;
     case FLOAT_ID_TOKEN:
       Destroy_Float((DBL *)Data);
       break;
     case PIGMENT_ID_TOKEN:
     case DENSITY_ID_TOKEN:
       Destroy_Pigment((PIGMENT *)Data);
       break;
     case TNORMAL_ID_TOKEN:
       Destroy_Tnormal((TNORMAL *)Data);
       break;
     case FINISH_ID_TOKEN:
       Destroy_Finish(Data);
       break;
     case MEDIA_ID_TOKEN:
       Destroy_Media((IMEDIA *)Data);
       break;
     case INTERIOR_ID_TOKEN:
       Destroy_Interior((INTERIOR *)Data);
       break;
     case MATERIAL_ID_TOKEN:
       Destroy_Material((MATERIAL *)Data);
       break;
     case TEXTURE_ID_TOKEN:
       Destroy_Textures((TEXTURE *)Data);
       break;
     case OBJECT_ID_TOKEN:
       Destroy_Object((OBJECT *)Data);
       break;
     case COLOUR_MAP_ID_TOKEN:
     case PIGMENT_MAP_ID_TOKEN:
     case SLOPE_MAP_ID_TOKEN:
     case TEXTURE_MAP_ID_TOKEN:
     case NORMAL_MAP_ID_TOKEN:
     case DENSITY_MAP_ID_TOKEN:
       Destroy_Blend_Map((BLEND_MAP *)Data);
       break;
     case TRANSFORM_ID_TOKEN:
       Destroy_Transform((TRANSFORM *)Data);
       break;
     case CAMERA_ID_TOKEN:
       Destroy_Camera((CAMERA *)Data);
       break;
     case RAINBOW_ID_TOKEN:
       Destroy_Rainbow((RAINBOW *)Data);
       break;
     case FOG_ID_TOKEN:
       Destroy_Fog((FOG *)Data);
       break;
     case SKYSPHERE_ID_TOKEN:
       Destroy_Skysphere((SKYSPHERE *)Data);
       break;
     case MACRO_ID_TOKEN:
     case TEMPORARY_MACRO_ID_TOKEN:
       Destroy_Macro((POV_MACRO *)Data);
       break;
     case STRING_ID_TOKEN:
         POV_FREE(Data);
       break;
     case ARRAY_ID_TOKEN:
       a = (POV_ARRAY *)Data;
       for(i=0; i<a->Total; i++)
       {
         Destroy_Ident_Data(a->DataPtrs[i], a->Type);
       }
       if(a->DataPtrs != NULL)
         POV_FREE(a->DataPtrs);
       POV_FREE(a);
       break;
     case PARAMETER_ID_TOKEN:
       POV_FREE(Data);
       break;
     case FILE_ID_TOKEN:
       Temp_File = (DATA_FILE *)Data;
       if(Temp_File->In_File != NULL)
         delete Temp_File->In_File;
       if(Temp_File->Out_File != NULL)
         delete Temp_File->Out_File;
       POV_FREE(Data);
       break;
     case FUNCT_ID_TOKEN:
     case VECTFUNCT_ID_TOKEN:
       Destroy_Function((FUNCTION_PTR)Data);
       break;
     case SPLINE_ID_TOKEN:
       Destroy_Spline((SPLINE *)Data);
       break;
     default:
       Error("Do not know how to free memory for identifier type %d", Type);
   }
}





/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Link(OBJECT *New_Object, OBJECT **Object_List_Root)
{
	New_Object->Sibling = *Object_List_Root;
	*Object_List_Root = New_Object;
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Link_Textures (TEXTURE **Old_Textures, TEXTURE *New_Textures)
{
   TEXTURE *Layer;
   
   if (New_Textures == NULL)
     return;

   if ((*Old_Textures) != NULL)
   {
      if ((*Old_Textures)->Type != PLAIN_PATTERN) 
      {
         Error("Cannot layer over a patterned texture.");
      }
   }
   for (Layer = New_Textures ;
        Layer->Next != NULL ;
        Layer = (TEXTURE *)Layer->Next)
  {
    /* NK layers - 1999 June 10 - for backwards compatiblity with layered textures */
    if(opts.Language_Version<=310)
      Convert_Filter_To_Transmit(Layer->Pigment);
  }

  /* NK layers - 1999 Nov 16 - for backwards compatiblity with layered textures */
  if ((opts.Language_Version<=310) && (*Old_Textures!=NULL))
    Convert_Filter_To_Transmit(Layer->Pigment);

   Layer->Next = (TPATTERN *)*Old_Textures;
   *Old_Textures = New_Textures;

   if ((New_Textures->Type != PLAIN_PATTERN) && (New_Textures->Next != NULL))
   {
      Error("Cannot layer a patterned texture over another.");
   }
  }



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

char *Get_Token_String (TOKEN Token_Id)
{
  register int i;

  for (i = 0 ; i < LAST_TOKEN ; i++)
     if (Reserved_Words[i].Token_Number == Token_Id)
        return (Reserved_Words[i].Token_Name);
  return ("");
}

/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
*   Return a list of keywords seperated with a '\n'. The last keyword does not
*   have a LF after it. The caller is responsible for freeing the memory. This
*   function is intended to be used by GUI implementations that need a keyword
*   list for syntax-coloring edit windows.
*
* CHANGES
*
******************************************************************************/

char *Get_Reserved_Words (
                           #ifdef CUSTOM_KEYWORDS_IN_WINPOV_PATCH
                             bool use_tokens,
                           #endif
                           const char *additional_words)
{
  int length = 0 ;
  int i ;

  #ifdef CUSTOM_KEYWORDS_IN_WINPOV_PATCH
  if (use_tokens)
  #endif
  for (i = 0; i < LAST_TOKEN; i++)
  {
    if (!isalpha (Reserved_Words [i].Token_Name [0]))
      continue ;
    if (strchr (Reserved_Words [i].Token_Name, ' ') != NULL)
      continue ;
    length += strlen (Reserved_Words[i].Token_Name) + 1 ;
  }

  length += strlen (additional_words) ;

  char *result = (char *) POV_MALLOC (++length, "Keyword List") ;
  strcpy (result, additional_words) ;
  char *s = result + strlen (additional_words) ;

  #ifdef CUSTOM_KEYWORDS_IN_WINPOV_PATCH
  if (use_tokens)
  #endif
  for (i = 0 ; i < LAST_TOKEN ; i++)
  {
    if (!isalpha (Reserved_Words [i].Token_Name [0]))
      continue ;
    if (strchr (Reserved_Words [i].Token_Name, ' ') != NULL)
      continue ;
    s += sprintf (s, "%s\n", Reserved_Words[i].Token_Name) ;
  }
  *--s = '\0' ;

  return (result) ;
}

/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
*    Changed macro ALLOW_REDEFINE to parameter which defaults to true [trf]
*
******************************************************************************/

void Test_Redefine(TOKEN Previous, TOKEN *NumberPtr, void *Data, bool allow_redefine)
{
	if ((Previous == IDENTIFIER_TOKEN) || (Previous == EMPTY_ARRAY_TOKEN))
	{
		return;
	}
	/* NK 1998 - allow user to redefine all identifiers! */
	if( allow_redefine)
	{
		Destroy_Ident_Data(Data, Previous);
	}
	else
	{
		if (Previous == *NumberPtr)
		{
			Destroy_Ident_Data(Data,*NumberPtr);
		}
		else
		{
			char *oldt, *newt;

			oldt = Get_Token_String (Previous);
			newt = Get_Token_String (*NumberPtr);
			*NumberPtr = Previous;

			Error ("Attempted to redefine %s as %s.", oldt, newt);
		}
	}
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Error(TOKEN Token_Id)
{
   Expectation_Error(Get_Token_String(Token_Id));
}
  
/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Found_Instead_Error(const char *exstr, const char *extokstr)
{
   char *found;
  
   Stage = STAGE_FOUND_INSTEAD;

   switch(Token.Token_Id)
   {
      case IDENTIFIER_TOKEN:
         Error("%s '%s', undeclared identifier '%s' found instead", exstr, extokstr, Token.Token_String);
         break;
      case VECTOR_FUNCT_TOKEN:
         found = Get_Token_String(Token.Function_Id);
         Error("%s '%s', vector function '%s' found instead", exstr, extokstr, found);
         break;
      case FLOAT_FUNCT_TOKEN:
         found = Get_Token_String(Token.Function_Id);
         Error("%s '%s', float function '%s' found instead", exstr, extokstr, found);
         break;
      case COLOUR_KEY_TOKEN:
         found = Get_Token_String(Token.Function_Id);
         Error("%s '%s', color keyword '%s' found instead", exstr, extokstr, found);
         break;
      default:
         found = Get_Token_String(Token.Token_Id);
         Error("%s '%s', %s found instead", exstr, extokstr, found);
   }
}

  
/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Warn_State(TOKEN Token_Id, TOKEN Type)
{
	char *str;

	str = (char *)POV_MALLOC(160, "global setting warning string");

	strcpy(str, "Found '");
	strcat(str, Get_Token_String (Token_Id));
	strcat(str, "' that should be in '");
	strcat(str, Get_Token_String (Type));
	strcat(str, "' statement.");
	if(150 >= opts.Language_Version)
		Warning(0, str);
	POV_FREE(str);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void MAError (const char *str, long size)
{
  Error ("Out of memory.  Cannot allocate %ld bytes for %s.", size, str);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Post_Process (OBJECT *Object,OBJECT  *Parent)
{
  DBL Volume;
  OBJECT *Sib;
  FINISH *Finish;

  if (Object == NULL)
  {
    return;
  }

  if (Object->Type & LT_SRC_UNION_OBJECT) 
  {
    for (Sib = ((CSG *)Object)->Children; Sib != NULL; Sib = Sib->Sibling)
    {
      Post_Process(Sib, Object);
    }
    return;
  }

  /* Promote texture etc. from parent to children. */

  if (Parent != NULL)
  {
    if (Object->Texture == NULL)
    {
      Object->Texture = Copy_Texture_Pointer(Parent->Texture);
      /* NK 1998 copy uv_mapping flag if and only if we copy the texture */
      if (Test_Flag(Parent, UV_FLAG))
        Set_Flag(Object, UV_FLAG);
      /* NK ---- */
    }
    if (Object->Interior_Texture == NULL)
    {
      Object->Interior_Texture = Copy_Texture_Pointer(Parent->Interior_Texture);
      if(Test_Flag(Parent, UV_FLAG))
        Set_Flag(Object, UV_FLAG);
    }
    if (Object->Interior == NULL)
    {
      Object->Interior = Copy_Interior_Pointer(Parent->Interior);
    }

    if (Test_Flag(Parent, NO_REFLECTION_FLAG))
    {
      Set_Flag(Object, NO_REFLECTION_FLAG);
    }
    if (Test_Flag(Parent, NO_IMAGE_FLAG))
    {
      Set_Flag(Object, NO_IMAGE_FLAG);
    }
    if (Test_Flag(Parent, NO_SHADOW_FLAG))
    {
      Set_Flag(Object, NO_SHADOW_FLAG);
    }
#ifdef NO_RADIOSITY_PATCH /* Christoph Hormann January 2004 */
    if (Test_Flag(Parent, NO_RADIOSITY_FLAG))
    {
      Set_Flag(Object, NO_RADIOSITY_FLAG);
    }
#endif
    if (Test_Flag(Parent, MULTITEXTURE_FLAG))
    {
      Set_Flag(Object, MULTITEXTURE_FLAG);
    }

    /* NK phmap */
    /* promote photon mapping flags to child */
    if (Test_Flag(Parent, PH_TARGET_FLAG))
    {
      Set_Flag(Object, PH_TARGET_FLAG);
      Object->Ph_Density = Parent->Ph_Density;
      CheckPassThru(Object, PH_TARGET_FLAG);
    }

    if(Test_Flag(Parent, PH_PASSTHRU_FLAG))
    {
      Set_Flag(Object, PH_PASSTHRU_FLAG);
      CheckPassThru(Object, PH_PASSTHRU_FLAG);
    }

    if (Test_Flag(Parent, PH_RFL_ON_FLAG))
    {
      Set_Flag(Object, PH_RFL_ON_FLAG);
      Clear_Flag(Object, PH_RFL_OFF_FLAG);
    }
    else if (Test_Flag(Parent, PH_RFL_OFF_FLAG))
    {
      Set_Flag(Object, PH_RFL_OFF_FLAG);
      Clear_Flag(Object, PH_RFL_ON_FLAG);
    }

    if (Test_Flag(Parent, PH_RFR_ON_FLAG))
    {
      Set_Flag(Object, PH_RFR_ON_FLAG);
      Clear_Flag(Object, PH_RFR_OFF_FLAG);
      CheckPassThru(Object, PH_RFR_ON_FLAG);
    }
    else if (Test_Flag(Parent, PH_RFR_OFF_FLAG))
    {
      Set_Flag(Object, PH_RFR_OFF_FLAG);
      Clear_Flag(Object, PH_RFR_ON_FLAG);
    }

    if(Test_Flag(Parent, PH_IGNORE_PHOTONS_FLAG))
    {
      Set_Flag(Object, PH_IGNORE_PHOTONS_FLAG);
    }

#ifdef MOTION_BLUR_PATCH
    if(Parent->TimeStamp)
    {
      Object->TimeStamp = Parent->TimeStamp;
    }
#endif
  }

  if (Object->Interior != NULL)
  {
     Post_Media(Object->Interior->IMedia);
  }

  if ((Object->Texture == NULL) &&
      !(Object->Type & TEXTURED_OBJECT) &&
      !(Object->Type & LIGHT_SOURCE_OBJECT))
  {
    if (Parent)
    {
      if(Parent->Methods != &CSG_Intersection_Methods ||
         !Test_Flag(Parent, MULTITEXTURE_FLAG))
      {
        Object->Texture = Copy_Textures(Default_Texture);
      }
    }
    else
    Object->Texture = Copy_Textures(Default_Texture);
  }

  if(!(Object->Type & LIGHT_GROUP_OBJECT) &&
     !(Object->Type & LIGHT_GROUP_LIGHT_OBJECT))
  {
    Post_Textures(Object->Texture);  /*moved cey 6/97*/

    if (Object->Interior_Texture)
    {
      Post_Textures(Object->Interior_Texture);
    }
  }

  if(Object->Type & LIGHT_SOURCE_OBJECT)
  {
    DBL len1,len2;
    LIGHT_SOURCE *Light = (LIGHT_SOURCE *)Object;


    /* check some properties of the orient light sources */
    if (Light->Orient)
    {
      if(!Light->Circular)
      {
        Light->Circular = true;
        Warning(0,"Orient can only be used with circular area lights. This area light is now circular.");
      }

      VLength(len1,Light->Axis1);
      VLength(len2,Light->Axis2);

      if(fabs(len1-len2)>EPSILON)
      {
        Warning(0, "When using orient, the two axes of the area light must be of equal length.\nOnly the length of the first axis will be used.");
        /*
        the equalization is actually done in the lighting code, since only the length of
        Axis1 will be used
        */
      }

      if(Light->Area_Size1 != Light->Area_Size2)
      {
        Warning(0, "When using orient, the two sample sizes for the area light should be equal.");
      }
    }

    /* Make sure that circular light sources are larger than 1 by x [ENB 9/97] */
    if (Light->Circular)
    {
      if ((Light->Area_Size1 <= 1) || (Light->Area_Size2 <= 1))
      {
        Warning(0, "Circular area lights must have more than one sample along each axis.");
        Light->Circular = false;
      }
    }
  }

  if (Object->Type & LIGHT_SOURCE_OBJECT)
  {
    // post-process the light source
    if (((LIGHT_SOURCE *)Object)->Projected_Through_Object != NULL) 
    {
      if (((LIGHT_SOURCE *)Object)->Projected_Through_Object->Interior != NULL)
      {
         Destroy_Interior(((LIGHT_SOURCE *)Object)->Projected_Through_Object->Interior);
         ((LIGHT_SOURCE *)Object)->Projected_Through_Object->Interior=NULL;
         Warning(0,"Projected through objects can not have interior, interior removed.");
      }
      if (((LIGHT_SOURCE *)Object)->Projected_Through_Object->Texture != NULL) 
      {
         Destroy_Textures(((LIGHT_SOURCE *)Object)->Projected_Through_Object->Texture);
         ((LIGHT_SOURCE *)Object)->Projected_Through_Object->Texture = NULL;
         Warning(0,"Projected through objects can not have texture, texture removed.");
      }
    }

    // only global light sources are in Frame.Light_Sources list [trf]
    if(!(Object->Type & LIGHT_GROUP_LIGHT_OBJECT))
    {
      // add this light to the frame's list of global light sources
      ((LIGHT_SOURCE *)Object)->Next_Light_Source = Frame.Light_Sources;

      Frame.Light_Sources = (LIGHT_SOURCE *)Object;

      Frame.Number_Of_Light_Sources++;
    }
    else
    {
      // Put it into the frame's list of light-source lights

      // NOTE - We use LIGHT_GROUP_LIGHT to create this linked list, because we can't
      //        use the Next_Light_Source pointer that is part of the light source.
      //        We can't use that pointer because it is used for light group purposes,
      //        and here we just want a list of all lights that are in any light group.

      LIGHT_GROUP_LIGHT *Node = (LIGHT_GROUP_LIGHT*)POV_MALLOC(sizeof(LIGHT_GROUP_LIGHT),"light_group_light");
      
      Node->Light = (LIGHT_SOURCE *)Object;

      // link node to list
      Node->Next = Frame.Light_Group_Lights;
      Frame.Light_Group_Lights = Node;
    }
  }
  else
  {
    // post-process the object

    /* If there is no interior create one. */

    if (Object->Interior == NULL)
    {
      Object->Interior = Create_Interior();
    }

    /* Promote hollow flag to interior. */

    Object->Interior->hollow = (Test_Flag(Object, HOLLOW_FLAG) != false);

    /* Promote finish's IOR to interior IOR. */

    if (Object->Texture != NULL)
    {
      if (Object->Texture->Type == PLAIN_PATTERN)
      {
        if ((Finish = Object->Texture->Finish) != NULL)
        {
          if (Finish->Temp_IOR >= 0.0)
          {
            Object->Interior->IOR = Finish->Temp_IOR;
            Object->Interior->Dispersion = Finish->Temp_Dispersion;
          }
          if (Finish->Temp_Caustics >= 0.0)
          {
            Object->Interior->Caustics = Finish->Temp_Caustics;
          }

          Object->Interior->Old_Refract = Finish->Temp_Refract;
        }
      }
    }

    /* If there is no IOR specified use the atmopshere ior. */

    if (Object->Interior->IOR == 0.0)
    {
      Object->Interior->IOR = Frame.Atmosphere_IOR;
      Object->Interior->Dispersion = Frame.Atmosphere_Dispersion;
    }
  }

  if (Object->Type & IS_COMPOUND_OBJECT)
  {
    for (Sib = ((CSG *)Object)->Children; Sib != NULL; Sib = Sib->Sibling)
    {
      Post_Process(Sib, Object);
    }
  }
  /* Test wether the object is finite or infinite. [DB 9/94] */

  BOUNDS_VOLUME(Volume, Object->BBox);

  if (Volume > INFINITE_VOLUME)
  {
    Set_Flag(Object, INFINITE_FLAG);
  }

  /* Test wether the object is opaque or not. [DB 8/94] */

  if ((Object->Methods != &Blob_Methods) &&
      (Object->Methods != &Mesh_Methods) &&
      (Test_Opacity(Object->Texture)) && 
      ((Object->Interior_Texture==NULL) || 
        Test_Opacity(Object->Interior_Texture)))
  {
    Set_Flag(Object, OPAQUE_FLAG);
  }
  else
  {
    /* Objects with multiple textures have to be handled separately. */

    if (Object->Methods == &Blob_Methods)
    {
      Test_Blob_Opacity((BLOB *)Object);
    }

    if (Object->Methods == &Mesh_Methods)
    {
      Test_Mesh_Opacity((MESH *)Object);
    }
  }
}

/*****************************************************************************
*
* FUNCTION
*
*   Link_To_Frame
*
* INPUT
*
*   Object - Pointer to object
*   
* OUTPUT
*
*   Object
*   
* RETURNS
*   
* AUTHOR
*
*   POV-Ray Team
*   
* DESCRIPTION
*
*   -
*
* CHANGES
*
*   Sep 1994 : Added optional splitting of bounded unions if children are
*              finite. Added removing of unnecessary bounding. [DB]
*
******************************************************************************/

void Link_To_Frame(OBJECT *Object)
{
  int finite;
  DBL Volume;
  OBJECT *This_Sib, *Next_Sib;

  if (Object == NULL)           /* patches a memory addressing error jdm mar/95 */
    return;

  /* Remove bounding object if object is cheap to intersect. [DB 8/94]  */

  if ((opts.Options & REMOVE_BOUNDS) && (Object->Bound != NULL))
  {
    if ((Object->Methods != &CSG_Union_Methods)        &&
        (Object->Methods != &CSG_Intersection_Methods) &&
        (Object->Methods != &CSG_Merge_Methods)        &&
        (Object->Methods != &Poly_Methods)             &&
        (Object->Methods != &TTF_Methods)			   &&
		(Object->Methods != &Quadric_Methods ||
			((QUADRIC *)Object)->Automatic_Bounds))
    {
      /* Destroy only, if bounding object is not used as clipping object. */

      if (Object->Bound != Object->Clip)
      {
        Destroy_Object(Object->Bound);
      }

      Object->Bound = NULL;

      Warning(0, "Unnecessary bounding object removed.");
    }
  }

  /*
   * [CJC 8/01]
   *
   * if all children of a union have the no_shadow flag, then the union should
   * have it as well.
   */
  if (Object->Methods == &CSG_Union_Methods)
  {
    for (This_Sib = ((CSG * ) Object)->Children ; This_Sib != NULL ; This_Sib = This_Sib->Sibling)
      if (This_Sib->Methods != &Light_Source_Methods && !Test_Flag (This_Sib, NO_SHADOW_FLAG))
        break ;
    if (This_Sib == NULL)
      Set_Flag (Object, NO_SHADOW_FLAG) ;
  }

  /*
   * Link the object to the frame if it's not a CSG union object,
   * if it's clipped or if bounding slabs aren't used.
   */

  if ((Object->Methods != &CSG_Union_Methods) || (Object->Clip != NULL) || (!opts.Use_Slabs))
  {
    Link(Object, &(Frame.Objects));

    return;
  }

  /*
   * [DB 8/94]
   *
   * The object is a CSG union object. It will be split if all siblings are
   * finite, i.e. the volume of the bounding box doesn't exceed a threshold.
   */

  /* NK phmap - added code so union is not split up if it is
                flagged for hi-density photon mapping...
            maybe we SHOULD split it anyways... do speed tests later */
  if(!((CSG *)Object)->do_split)
  {
    Link(Object, &(Frame.Objects));
    return;
  }

  if (Object->Bound != NULL)
  {
    /* Test if all siblings are finite. */

    finite = true;

    for (This_Sib = ((CSG *)Object)->Children; This_Sib != NULL; This_Sib = This_Sib->Sibling)
    {
      BOUNDS_VOLUME(Volume, This_Sib->BBox);

      if (Volume > BOUND_HUGE)
      {
        finite = false;

        break;
      }
    }

    /*
     * If the union has infinite children or splitting is not used and
     * the object is not a light group link the union to the frame.
     */

    if (((!finite) || !(opts.Options & SPLIT_UNION)) && ((Object->Type & LIGHT_GROUP_OBJECT) != LIGHT_GROUP_OBJECT))
    {
      if (finite)
      {
        Warning(0, "CSG union unnecessarily bounded.");
      }

      Link(Object, &(Frame.Objects));

      return;
    }

    Warning(0, "Bounded CSG union split.");
  }

  /* Link all siblings of a union to the frame. */

  for (This_Sib = ((CSG *)Object)->Children; This_Sib != NULL; This_Sib = Next_Sib)
  {
    /* Link_To_Frame() changes Sibling so save it */

    Next_Sib = This_Sib->Sibling;

    /* Sibling is no longer inside a CSG object. */

    This_Sib->Type &= ~IS_CHILD_OBJECT;

    Link_To_Frame (This_Sib);
  }

/*
  Object->Texture = NULL;
*/

  Object->Sibling = NULL;

  ((CSG *)Object)->Children = NULL;

  Destroy_Object (Object);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Only_In(const  char *s1, const char *s2)
{
  Error("Keyword '%s' can only be used in a %s statement.",s1,s2);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Not_With(const char *s1, const char *s2)
{
  Error("Keyword '%s' cannot be used with %s.",s1,s2);
}

void Warn_Compat(int f, const char *syn)
{
  char isNotText[] = "is not";
  char mayNotText[] = "may not be";
  char *text;

  if (f)
  {
    text = isNotText;
  }
  else
  {
    text = mayNotText;
  }
    
  Warning(0,"%s\n"
            "  Use of this syntax %s backwards compatable with earlier versions of POV-Ray.\n"
            "  The #version directive or +MV switch will not help.",
            syn, text);
}

/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
*   Mar 1996 : Add line number info to warning message  [AED]
*
******************************************************************************/

static void Global_Setting_Warn()
{
	char *str;

	str = (char *)POV_MALLOC(strlen(Token.Token_String) + 80, "global setting warning string");

	if (opts.Language_Version >= 300)
	{
		strcpy(str, "'");
		strcat(str, Token.Token_String);
		strcat(str, "' should be in 'global_settings{...}' statement.");
		PossibleError(str);
	}

	POV_FREE(str);
}

/*****************************************************************************
*
* FUNCTION
*
*  Set_CSG_Children_Hollow
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static void Set_CSG_Children_Flag(OBJECT *Object, unsigned int f, unsigned int  flag, unsigned int  set_flag)
{
  OBJECT *Sib;

  for (Sib = ((CSG *)Object)->Children; Sib != NULL; Sib = Sib->Sibling)
  {
    if (!Test_Flag(Sib, set_flag))
    {
      if ((Sib->Methods == &CSG_Intersection_Methods) ||
          (Sib->Methods == &CSG_Merge_Methods) ||
          (Sib->Methods == &CSG_Union_Methods))
      {
        Set_CSG_Children_Flag(Sib, f, flag, set_flag);
      }
      else
      {
        Sib->Flags = (Sib->Flags & (~flag)) | f;
      }
    }
  }
}



/*****************************************************************************
*
* FUNCTION
*
*  Set_CSG_Tree_Flag
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

static void Set_CSG_Tree_Flag(OBJECT *Object, unsigned int f, int val)
{
  OBJECT *Sib;

  for (Sib = ((CSG *)Object)->Children; Sib != NULL; Sib = Sib->Sibling)
  {
    if ((Sib->Methods == &CSG_Intersection_Methods) ||
        (Sib->Methods == &CSG_Merge_Methods) ||
        (Sib->Methods == &CSG_Union_Methods))
    {
      Set_CSG_Tree_Flag(Sib, f, val);
    }
    Bool_Flag (Sib, f, val);
  }
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void *Copy_Identifier (void *Data, int Type)
{
  int i;
  POV_ARRAY *a, *na;
  VECTOR *vp;
  DBL *dp;
  UV_VECT *uvp;
  VECTOR_4D *v4p;
  int len;
  void *New=NULL;
  
  if (Data==NULL)
  {
     return(NULL);
  }

  switch (Type)
  {
     case COLOUR_ID_TOKEN:
       New = (void *)Copy_Colour(*(COLOUR *)Data);
       break;
     case VECTOR_ID_TOKEN:
       vp = Create_Vector();
       Assign_Vector((*vp),(*((VECTOR *)Data)));
       New=vp;
       break;
     case UV_ID_TOKEN:
       uvp = Create_UV_Vect();
       Assign_UV_Vect((*uvp),(*((UV_VECT *)Data)));
       New=uvp;
       break;
     case VECTOR_4D_ID_TOKEN:
       v4p = Create_Vector_4D();
       Assign_Vector_4D((*v4p),(*((VECTOR_4D *)Data)));
       New=v4p;
       break;
     case FLOAT_ID_TOKEN:
       dp = Create_Float();
       *dp = *((DBL *)Data);
       New = dp;
       break;
     case PIGMENT_ID_TOKEN:
     case DENSITY_ID_TOKEN:
       New = (void *)Copy_Pigment((PIGMENT *)Data);
       break;
     case TNORMAL_ID_TOKEN:
       New = (void *)Copy_Tnormal((TNORMAL *)Data);
       break;
     case FINISH_ID_TOKEN:
       New = (void *)Copy_Finish((FINISH *)Data);
       break;
     case MEDIA_ID_TOKEN:
       New = (void *)Copy_Media((IMEDIA *)Data);
       break;
     case INTERIOR_ID_TOKEN:
       New = (void *)Copy_Interior((INTERIOR *)Data);
       break;
     case MATERIAL_ID_TOKEN:
       New = (void *)Copy_Material((MATERIAL *)Data);
       break;
     case TEXTURE_ID_TOKEN:
       New = (void *)Copy_Textures((TEXTURE *)Data);
       break;
     case OBJECT_ID_TOKEN:
       New = (void *)Copy_Object((OBJECT *)Data);
       break;
     case COLOUR_MAP_ID_TOKEN:
     case PIGMENT_MAP_ID_TOKEN:
     case SLOPE_MAP_ID_TOKEN:
     case TEXTURE_MAP_ID_TOKEN:
     case NORMAL_MAP_ID_TOKEN:
     case DENSITY_MAP_ID_TOKEN:
       New = (void *)Copy_Blend_Map((BLEND_MAP *)Data);
       break;
     case TRANSFORM_ID_TOKEN:
       New = (void *)Copy_Transform((TRANSFORM *)Data);
       break;
     case CAMERA_ID_TOKEN:
       New = (void *)Copy_Camera((CAMERA *)Data);
       break;
     case RAINBOW_ID_TOKEN:
       New = (void *)Copy_Rainbow((RAINBOW *)Data);
       break;
     case FOG_ID_TOKEN:
       New = (void *)Copy_Fog((FOG *)Data);
       break;
     case SKYSPHERE_ID_TOKEN:
       New = (void *)Copy_Skysphere((SKYSPHERE *)Data);
       break;
     case STRING_ID_TOKEN:
       //New = (void *)POV_STRDUP((char *)Data);
       len = UCS2_strlen((UCS2 *)(Data)) + 1;
       New = (UCS2 *)POV_MALLOC(len * sizeof(UCS2), "UCS2 String");
       POV_MEMMOVE((void *)New, (void *)(Data), len * sizeof(UCS2));
       break;
     case ARRAY_ID_TOKEN:
       a=(POV_ARRAY *)Data;
       na=(POV_ARRAY *)POV_MALLOC(sizeof(POV_ARRAY),"array");
       *na=*a;
       na->DataPtrs = (void **)POV_MALLOC(sizeof(void *)*(a->Total),"array");
       for (i=0; i<a->Total; i++)
       {
         na->DataPtrs[i] = (void *)Copy_Identifier (a->DataPtrs[i],a->Type);
       }
       New = (void *)na;
       break;
     case FUNCT_ID_TOKEN:
     case VECTFUNCT_ID_TOKEN:
       New = (void *)Copy_Function((FUNCTION_PTR )Data);
       break;
     case SPLINE_ID_TOKEN:
       New = (void *)Copy_Spline((SPLINE *)Data);
       break;
     default:
       Error("Cannot copy identifier");
   }
   return(New);
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

/* NK layers - 1999 June 10 - for backwards compatiblity with layered textures */
void Convert_Filter_To_Transmit(PIGMENT *Pigment)
{
  int i;
  BLEND_MAP *Map;

  if (Pigment==NULL) return;

  switch (Pigment->Type)
  {
    case PLAIN_PATTERN:
      Pigment->Colour[pTRANSM]+=Pigment->Colour[pFILTER];
      Pigment->Colour[pFILTER]=0;
      break;

    default:
      if (Pigment->Blend_Map != NULL)
      {
        Map = Pigment->Blend_Map;
        /* go through blend map */
        if ((Map->Type == PIGMENT_TYPE) || (Map->Type == DENSITY_TYPE))
        {
           for (i = 0; i < Map->Number_Of_Entries; i++)
           {
             Convert_Filter_To_Transmit(Map->Blend_Map_Entries[i].Vals.Pigment);
           }
        }
        else
        {
           for (i = 0; i < Map->Number_Of_Entries; i++)
           {
             Map->Blend_Map_Entries[i].Vals.Colour[pTRANSM]+=Map->Blend_Map_Entries[i].Vals.Colour[pFILTER];
             Map->Blend_Map_Entries[i].Vals.Colour[pFILTER]=0;
           }
        }

      }

      break;
  }
}



/*****************************************************************************
*
* FUNCTION
*
* INPUT
*   
* OUTPUT
*   
* RETURNS
*   
* AUTHOR
*   
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Expectation_Error(const char *s)
{
  Found_Instead_Error("Expected", s);
}


//patches from here on
#ifdef GLOW_PATCH
/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Glow(GLOW * Glow)
{
  MATRIX Matrix;
  int Warped = false;
  TRANSFORM * Trans = Create_Transform();
  TRANSFORM * Temp_Trans = Create_Transform();
  TRANSFORM * Trans2;
  VECTOR Local_Vector;
  TPATTERN * Warps = (TPATTERN *)POV_MALLOC(sizeof(TPATTERN), "TPATTERN");
  Init_TPat_Fields(Warps);
  Warps->Type = BOZO_PATTERN;/*Set it to a pattern type to allow turbulence*/

  Parse_Begin();

  EXPECT
    CASE(SIZE_TOKEN)
      Glow->Size = Parse_Float();
    END_CASE

    CASE(RADIUS_TOKEN)
      Glow->Cutoff_Radius = Parse_Float();
    END_CASE

    CASE(FADE_POWER_TOKEN)
      Glow->fade_pow = Parse_Float();
    END_CASE

    CASE(TYPE_TOKEN)
      Glow->Glow_Type = (int)Parse_Float();
      if(Glow->Glow_Type < 0 || Glow->Glow_Type > 3)
        Error("\nUnknown glow type.\n");
    END_CASE

    CASE(LOCATION_TOKEN)
      Parse_Vector(Glow->Center);
    END_CASE

    CASE(COLOUR_TOKEN)
      Parse_Colour(Glow->Colour);
    END_CASE

    CASE (WARP_TOKEN)
      Warped = true;
      Parse_Warp(&(Warps->Warps));
    END_CASE

    CASE (TRANSLATE_TOKEN)
      Parse_Vector(Local_Vector);
      Compute_Translation_Transform(Temp_Trans, Local_Vector);
      Compose_Transforms(Trans, Temp_Trans);
    END_CASE

    CASE (ROTATE_TOKEN)
      Parse_Vector(Local_Vector);
      Compute_Rotation_Transform(Temp_Trans, Local_Vector);
      Compose_Transforms(Trans, Temp_Trans);
    END_CASE

    CASE (SCALE_TOKEN)
      Parse_Scale_Vector(Local_Vector);
      Compute_Scaling_Transform(Trans, Local_Vector);
      Compose_Transforms(Trans, Temp_Trans);
    END_CASE

    CASE (MATRIX_TOKEN)
      Parse_Matrix(Matrix);
      Compute_Matrix_Transform(Temp_Trans, Matrix);
      Compose_Transforms(Trans, Temp_Trans);
    END_CASE

    CASE (TRANSFORM_TOKEN)
      Trans2 = Parse_Transform();
      Compose_Transforms(Trans, Trans2);
      Destroy_Transform(Trans2);
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  Parse_End();
  if(Warped == true)
  {
    Glow->Warps = Warps;
  }
  else
  {
    Destroy_TPat_Fields(Warps);
    POV_FREE(Warps);
  }
  Transform_Glow(Glow, Trans);
  Destroy_Transform(Trans);
  Destroy_Transform(Temp_Trans);
}
#endif


#ifdef MECHANICS_SIMULATION_PATCH            /* Christoph Hormann August 2002 */

/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Mech_Sim_Mass (int Group)
{

  Parse_Begin();

  if ( MechsimOptions.Data.mass_count >= MechsimOptions.Data.max_masses )  // allocate new space if necessary
  {
    if (MechsimOptions.Data.max_masses >= INT_MAX/2)
    {
      Error("Too many masses in simulation data");
    }
    MechsimOptions.Data.max_masses *= 2;
    MechsimOptions.Data.Masses = (MECHSIM_MASS *)POV_REALLOC(MechsimOptions.Data.Masses, MechsimOptions.Data.max_masses*sizeof (MECHSIM_MASS), "mechanics simulation mass data block");
  }

  Parse_Vector(MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Position); Parse_Comma();  // Position parameter
  Parse_Vector(MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Velocity); Parse_Comma();  // Velocity parameter
  MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Radius = Parse_Float();                 // Radius parameter

  MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Mass = 0.0;                             // Mass default value
  MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Flag = 0;                               // Flag default value
  MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Attach = -1;                            // Attach default value
  MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Force = -1;                             // Custom force default value
  MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Group = Group;                          // Group as in parameter
  MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Index = MechsimOptions.Data.mass_count; // The number of the mass

  EXPECT

    CASE (MASS_TOKEN)                                                                                // Mass parameter (optional)
      MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Mass = Parse_Float();
    END_CASE

    CASE (DENSITY_TOKEN)                                                                             // Density parameter (optional)
      MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Mass =
          (4.0/3.0)*M_PI*MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Radius*
          MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Radius*
          MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Radius*Parse_Float();
    END_CASE

    CASE (FIXED_TOKEN)                                                                               // 'fixed' parameter (optional)
      if (Allow_Float(1.0) > 0.0)
        MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Flag |= MECHSIM_FLAG_FIXED;
    END_CASE

    CASE (ATTACH_TOKEN)                                                                              // 'attach' parameter (optional)
      MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Attach = (int)Parse_Float();
      if ((MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Attach > MechsimOptions.Attachments.count-1) ||
          (MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Attach < 0))
      {
        Error("Attachment index out of range (0..%d)", MechsimOptions.Attachments.count-1);
      }
    END_CASE

    CASE (FORCE_TOKEN)                                                                              // 'force' parameter (optional)
      MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Force = (int)Parse_Float();
      if ((MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Force > MechsimOptions.Forces.count-1) ||
          (MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Force < 0))
      {
        Error("Custom force index out of range (0..%d)", MechsimOptions.Forces.count-1);
      }
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE

  END_EXPECT

	if ((MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Flag & MECHSIM_FLAG_FIXED) && 
			(MechsimOptions.Data.Masses[MechsimOptions.Data.mass_count].Attach >= 0))
		Error("Mass cannot be both attached and fixed.");

  MechsimOptions.Data.mass_count++;

  Parse_End();
}

/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Mech_Sim_Connection (int Group)
{
  int Idx1, Idx2;

  Parse_Begin();

  if ( MechsimOptions.Data.connect_count >= MechsimOptions.Data.max_connects )  // allocate new space if necessary
  {
    if (MechsimOptions.Data.max_connects >= INT_MAX/2)
    {
      Error("Too many connections in simulation data");
    }
    MechsimOptions.Data.max_connects *= 2;
    MechsimOptions.Data.Connects = (MECHSIM_CONNECT *)POV_REALLOC(MechsimOptions.Data.Connects, MechsimOptions.Data.max_connects*sizeof (MECHSIM_CONNECT), "mechanics simulation connection data block");
  }

  Idx1 = (int)Parse_Float(); Parse_Comma();
  Idx2 = (int)Parse_Float();

  if ((Idx1>=MechsimOptions.Data.mass_count) || (Idx2>=MechsimOptions.Data.mass_count))
  {
    Error("Connection mass index out of range (%d)", MechsimOptions.Data.mass_count);
  }

  MechsimOptions.Data.Connects[MechsimOptions.Data.connect_count].Idx1 = Idx1;                      // First mass index parameter
  MechsimOptions.Data.Connects[MechsimOptions.Data.connect_count].Idx2 = Idx2;                      // Second mass index parameter

  //MechsimOptions.Data.Connects[MechsimOptions.Data.connect_count].Length = 0.0;

  VDist(MechsimOptions.Data.Connects[MechsimOptions.Data.connect_count].Length,                     // Default length: current distance
        MechsimOptions.Data.Masses[Idx1].Position,
        MechsimOptions.Data.Masses[Idx2].Position);
  MechsimOptions.Data.Connects[MechsimOptions.Data.connect_count].Stiffness = 0.0;                  // Default Stiffness
  MechsimOptions.Data.Connects[MechsimOptions.Data.connect_count].Damping = 0.0;                    // Default Damping
  MechsimOptions.Data.Connects[MechsimOptions.Data.connect_count].Group = Group;                    // Group as in parameter
  MechsimOptions.Data.Connects[MechsimOptions.Data.connect_count].Index = MechsimOptions.Data.connect_count; // The number of the connection

  EXPECT

    CASE (LENGTH_TOKEN)                                                                             // Length parameter (optional)
      MechsimOptions.Data.Connects[MechsimOptions.Data.connect_count].Length = Parse_Float();
    END_CASE

    CASE (STIFFNESS_TOKEN)                                                                          // Stiffness parameter (optional)
      MechsimOptions.Data.Connects[MechsimOptions.Data.connect_count].Stiffness = Parse_Float();
    END_CASE

    CASE (DAMPING_TOKEN)                                                                            // Damping parameter (optional)
      MechsimOptions.Data.Connects[MechsimOptions.Data.connect_count].Damping = Parse_Float();
    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE

  END_EXPECT


  MechsimOptions.Data.connect_count++;

  Parse_End();
}

/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Mech_Sim_Face (int Group)
{
  int Idx1, Idx2, Idx3;

  Parse_Begin();

  if ( MechsimOptions.Data.faces_count >= MechsimOptions.Data.max_faces )  // allocate new space if necessary
  {
    if (MechsimOptions.Data.max_faces >= INT_MAX/2)
    {
      Error("Too many faces in simulation data");
    }
    MechsimOptions.Data.max_faces *= 2;
    MechsimOptions.Data.Faces = (MECHSIM_FACE *)POV_REALLOC(MechsimOptions.Data.Faces, MechsimOptions.Data.max_faces*sizeof (MECHSIM_FACE), "mechanics simulation face data block");
  }

  Idx1 = (int)Parse_Float(); Parse_Comma();
  Idx2 = (int)Parse_Float(); Parse_Comma();
  Idx3 = (int)Parse_Float();

  if ((Idx1>=MechsimOptions.Data.mass_count) || (Idx2>=MechsimOptions.Data.mass_count) || (Idx3>=MechsimOptions.Data.mass_count) ||
      (Idx1<0) || (Idx2<0) || (Idx3<0))
  {
    Error("Face mass index out of range (%d)", MechsimOptions.Data.mass_count);
  }

  MechsimOptions.Data.Faces[MechsimOptions.Data.faces_count].Idx1 = Idx1;                      // First mass index parameter
  MechsimOptions.Data.Faces[MechsimOptions.Data.faces_count].Idx2 = Idx2;                      // Second mass index parameter
  MechsimOptions.Data.Faces[MechsimOptions.Data.faces_count].Idx3 = Idx3;                      // Third mass index parameter

  MechsimOptions.Data.Faces[MechsimOptions.Data.faces_count].Group = Group;                    // Group as in parameter
  MechsimOptions.Data.Faces[MechsimOptions.Data.faces_count].Index = MechsimOptions.Data.faces_count; // The number of the face
  MechsimOptions.Data.faces_count++;

  Parse_End();
}

/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Mech_Sim_Viscoelastic ()
{
  int Idx1, Idx2, i;

  Parse_Begin();

  if ( MechsimOptions.Data.ve_count >= MechsimOptions.Data.max_viscoelastics )  // allocate new space if necessary
  {
    if (MechsimOptions.Data.max_viscoelastics >= INT_MAX/2)
    {
      Error("Too many viscoelastics in simulation data");
    }
    MechsimOptions.Data.max_viscoelastics *= 2;
    MechsimOptions.Data.Viscoelastics = (MECHSIM_VISCOELASTIC *)POV_REALLOC(MechsimOptions.Data.Viscoelastics, 
																																						MechsimOptions.Data.max_viscoelastics*sizeof (MECHSIM_VISCOELASTIC), 
																																						"mechanics simulation viscoelastic data block");
  }

  Idx1 = (int)Parse_Float(); Parse_Comma();
  Idx2 = (int)Parse_Float();

  if ((Idx1>=MechsimOptions.Data.mass_count) || (Idx2>=MechsimOptions.Data.mass_count))
  {
    Error("Connection mass index out of range (%d)", MechsimOptions.Data.mass_count);
  }

  MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Idx1 = Idx1;                      // First mass index parameter
  MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Idx2 = Idx2;                      // Second mass index parameter

  //MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Length = 0.0;

  VDist(MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Length,                     // Default length: current distance
        MechsimOptions.Data.Masses[Idx1].Position,
        MechsimOptions.Data.Masses[Idx2].Position);
  MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Stiffness0 = 0.0;                 // Default Stiffness0
  MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Accuracy = 1;                     // Default Accuracy
  MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Index = MechsimOptions.Data.ve_count; // The number of the viscoelastic

  MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].element_count = 0;

  EXPECT

    CASE (LENGTH_TOKEN)                                                                             // Length parameter (optional)
      MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Length = Parse_Float();
    END_CASE

    CASE (STIFFNESS_TOKEN)                                                                          // Stiffness0 parameter (optional)
      MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Stiffness0 = Parse_Float();
    END_CASE

    CASE (ACCURACY_TOKEN)                                                                          // Accuracy parameter (optional)
      MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Accuracy = Parse_Float();

		  if (MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Accuracy < 1)
			{
				Error("Viscoelastic element accuracy has to be >= 1 (%d)", MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Accuracy);
			}
    END_CASE

    CASE (ELEMENT_TOKEN)                                                                            // Element

			Parse_Begin();

			i = MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].element_count;

			if (i == 0)
			{  // first element
				MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Element = 
											(MECHSIM_VE_ELEMENT *)POV_MALLOC(sizeof (MECHSIM_VE_ELEMENT), "mechanics simulation viscoelastic data");
			}
		  else
			{
				MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Element =
					(MECHSIM_VE_ELEMENT *)POV_REALLOC(MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Element, 
																						(i+1)*sizeof (MECHSIM_VE_ELEMENT), "mechanics simulation viscoelastic data");

			}

      MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Element[i].Stiffness = Parse_Float();
			Parse_Comma();
      MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Element[i].Damping = Parse_Float();
			Parse_Comma();
      MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].Element[i].Y = Allow_Float(0.0);

			MechsimOptions.Data.Viscoelastics[MechsimOptions.Data.ve_count].element_count++;

			Parse_End();

    END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE

  END_EXPECT


  MechsimOptions.Data.ve_count++;

  Parse_End();
}

/*****************************************************************************
*
* FUNCTION
*
* INPUT
*
* OUTPUT
*
* RETURNS
*
* AUTHOR
*
* DESCRIPTION
*
* CHANGES
*
******************************************************************************/

void Parse_Mech_Sim ()
{
  IStream *fdi;
  DBL Time;
	VECTOR Fixed;
  MECHSIM_ENVIRONMENT *Env;
  int Current_Group = 0;
	bool Time_Step_Given = false;
	bool Step_Count_Given = false;
	bool Time_Given = false;

  Parse_Begin();

  /* Allocate initial amount of data fields */
  MechsimOptions.Data.max_masses = 20;
  MechsimOptions.Data.max_connects = 20;
  MechsimOptions.Data.max_faces = 20;
  MechsimOptions.Data.max_viscoelastics = 20;

  MechsimOptions.Data.Masses = (MECHSIM_MASS *)POV_MALLOC(MechsimOptions.Data.max_masses*sizeof (MECHSIM_MASS), 
																													"mechanics simulation mass data block");
  MechsimOptions.Data.Connects = (MECHSIM_CONNECT *)POV_MALLOC(MechsimOptions.Data.max_connects*sizeof (MECHSIM_CONNECT), 
																															 "mechanics simulation connection data block");
  MechsimOptions.Data.Faces = (MECHSIM_FACE *)POV_MALLOC(MechsimOptions.Data.max_faces*sizeof (MECHSIM_FACE), 
																												 "mechanics simulation connection data block");
  MechsimOptions.Data.Viscoelastics = (MECHSIM_VISCOELASTIC *)POV_MALLOC(MechsimOptions.Data.max_viscoelastics*sizeof (MECHSIM_VISCOELASTIC), 
																																				 "mechanics simulation viscoelastic data block");

  MechsimOptions.Enabled = true;

  EXPECT
    CASE (BOUNDING_TOKEN)         // ----- Bounding method (integer) -----
      MechsimOptions.Bounding = (unsigned int)Parse_Float();
    END_CASE

    CASE (METHOD_TOKEN)         // ----- Method (integer) -----
      MechsimOptions.Method = (unsigned int)Parse_Float();
    END_CASE

    CASE (FIXED_TOKEN)         // ----- Fixed (integer vector) -----
      Parse_Vector(Fixed);
		  MechsimOptions.Fixed[0] = (Fixed[X] > EPSILON);
		  MechsimOptions.Fixed[1] = (Fixed[Y] > EPSILON);
		  MechsimOptions.Fixed[2] = (Fixed[Z] > EPSILON);
      if (MechsimOptions.Fixed[0] && MechsimOptions.Fixed[1] && MechsimOptions.Fixed[2])
            Warning(0, "MechSim: Simulation is fixed in all directions\n       results may be not as expected.");
    END_CASE

    CASE (GRAVITY_TOKEN)        // ----- Gravity parameter (vector) -----
      Parse_Vector(MechsimOptions.gravity);
    END_CASE

    CASE (TIME_STEP_TOKEN)      // ----- Time Step length (float) -----
      MechsimOptions.Time_Step = Parse_Float();
		  Time_Step_Given=true;
		  if (Step_Count_Given)
				MechsimOptions.End_Time = MechsimOptions.Start_Time + MechsimOptions.Step_Count*MechsimOptions.Time_Step;
      MechsimOptions.Timing_Given = true;
    END_CASE

    CASE (STEP_COUNT_TOKEN)     // ----- Step Count (integer) -----
      MechsimOptions.Step_Count = (int)Parse_Float();
		  Step_Count_Given = true;
		  if (Time_Step_Given)
				MechsimOptions.End_Time = MechsimOptions.Start_Time + MechsimOptions.Step_Count*MechsimOptions.Time_Step;
      MechsimOptions.Timing_Given = true;
    END_CASE

    CASE (TIME_TOKEN)           // ----- Total Time (float) -----
      Time = Parse_Float();
		  MechsimOptions.End_Time = MechsimOptions.Start_Time+Time;
		  if (Step_Count_Given)
				MechsimOptions.Time_Step = ((MechsimOptions.End_Time-MechsimOptions.Start_Time)/MechsimOptions.Step_Count);
			else if (Time_Step_Given)
				MechsimOptions.Step_Count = (int)((MechsimOptions.End_Time-MechsimOptions.Start_Time)/MechsimOptions.Time_Step);
				
      MechsimOptions.Timing_Given = true;
    END_CASE

    CASE (START_TIME_TOKEN)     // ----- Start Time (float) -----
      MechsimOptions.Start_Time = Parse_Float();
		  MechsimOptions.Time = MechsimOptions.Start_Time;
      MechsimOptions.End_Time += MechsimOptions.Start_Time;
      MechsimOptions.Timing_Given = true;
    END_CASE

    CASE (END_TIME_TOKEN)       // ----- End Time (float) -----
      MechsimOptions.End_Time = Parse_Float();
      MechsimOptions.Timing_Given = true;
    END_CASE

    CASE (ACCURACY_TOKEN)     // ----- Adaptive time stepping accuracy (float) -----
      MechsimOptions.Accuracy = Parse_Float();
		  Parse_Comma();
			MechsimOptions.Time_Step_Min = Allow_Float(1.0e-6);
    END_CASE

    CASE (INTERACTION_TOKEN)    // ----- Interaction -----
      Parse_Begin();

      EXPECT

        CASE (FUNCTION_TOKEN)       // Interaction Function (parameters: x, y, z, Distance)
          if (MechsimOptions.Interactions.max_count==0)
          {
            MechsimOptions.Interactions.max_count = 10;
            MechsimOptions.Interactions.Function =
              (FUNCTION_PTR *)POV_MALLOC(MechsimOptions.Interactions.max_count*sizeof (FUNCTION_PTR),
              "mechanics simulation interaction functions");
          }
          else if ( MechsimOptions.Interactions.count >= MechsimOptions.Interactions.max_count )  // allocate new space if necessary
          {
            if (MechsimOptions.Interactions.max_count >= INT_MAX/2)
            {
              Error("Too many interaction functions");
            }
            MechsimOptions.Interactions.max_count *= 2;
            MechsimOptions.Interactions.Function =
              (FUNCTION_PTR *)POV_REALLOC(MechsimOptions.Interactions.Function,
              MechsimOptions.Interactions.max_count*sizeof (FUNCTION_PTR),
              "mechanics simulation interaction functions");
          }

          MechsimOptions.Interactions.Function[MechsimOptions.Interactions.count] = Parse_ParameterFunction();

          MechsimOptions.Interactions.count++;

        END_CASE

        OTHERWISE
          UNGET
          EXIT
        END_CASE

      END_EXPECT

      Parse_End();
    END_CASE

    CASE (FIELD_TOKEN)    // ----- Fields -----
      Parse_Begin();

      EXPECT

        CASE (FUNCTION_TOKEN)       // Field Function (parameters: x, y, z, return: color)
          if (MechsimOptions.Fields.max_count==0)
          {
            MechsimOptions.Fields.max_count = 10;
            MechsimOptions.Fields.Function =
              (FUNCTION_PTR *)POV_MALLOC(MechsimOptions.Fields.max_count*sizeof (FUNCTION_PTR),
              "mechanics simulation field functions");
          }
          else if ( MechsimOptions.Fields.count >= MechsimOptions.Fields.max_count )  // allocate new space if necessary
          {
            if (MechsimOptions.Fields.max_count >= INT_MAX/2)
            {
              Error("Too many field functions");
            }
            MechsimOptions.Fields.max_count *= 2;
            MechsimOptions.Fields.Function =
              (FUNCTION_PTR *)POV_REALLOC(MechsimOptions.Fields.Function,
              MechsimOptions.Fields.max_count*sizeof (FUNCTION_PTR),
              "mechanics simulation fields functions");
          }

          MechsimOptions.Fields.Function[MechsimOptions.Fields.count] = Parse_PigmentFunction();

          MechsimOptions.Fields.count++;

        END_CASE

        OTHERWISE
          UNGET
          EXIT
        END_CASE

      END_EXPECT

      Parse_End();
    END_CASE

    CASE (FORCE_TOKEN)    // ----- Forces -----
      Parse_Begin();

      EXPECT

        CASE (FUNCTION_TOKEN)       // Force Function (parameters: x, y, z, return: color)
          if (MechsimOptions.Forces.max_count==0)
          {
            MechsimOptions.Forces.max_count = 10;
            MechsimOptions.Forces.Function =
              (FUNCTION_PTR *)POV_MALLOC(MechsimOptions.Forces.max_count*sizeof (FUNCTION_PTR),
              "mechanics simulation force functions");
          }
          else if ( MechsimOptions.Forces.count >= MechsimOptions.Forces.max_count )  // allocate new space if necessary
          {
            if (MechsimOptions.Forces.max_count >= INT_MAX/2)
            {
              Error("Too many force functions");
            }
            MechsimOptions.Forces.max_count *= 2;
            MechsimOptions.Forces.Function =
              (FUNCTION_PTR *)POV_REALLOC(MechsimOptions.Forces.Function,
              MechsimOptions.Forces.max_count*sizeof (FUNCTION_PTR),
              "mechanics simulation force functions");
          }

          MechsimOptions.Forces.Function[MechsimOptions.Forces.count] = Parse_PigmentFunction();

          MechsimOptions.Forces.count++;

        END_CASE

        OTHERWISE
          UNGET
          EXIT
        END_CASE

      END_EXPECT

      Parse_End();
    END_CASE

    CASE (ATTACH_TOKEN)    // ----- Attachment -----

			/* get space (function and fixed) */
			if (MechsimOptions.Attachments.max_count==0)
			{
				MechsimOptions.Attachments.max_count = 10;
				MechsimOptions.Attachments.Function =
					(FUNCTION_PTR *)POV_MALLOC(MechsimOptions.Attachments.max_count*sizeof (FUNCTION_PTR),
																		 "mechanics simulation attachment functions");
				MechsimOptions.Attachments.Fixed =
					(INT_VECT *)POV_MALLOC(MechsimOptions.Attachments.max_count*sizeof (INT_VECT),
																 "mechanics simulation attachment vector");
			}
			else if ( MechsimOptions.Attachments.count >= MechsimOptions.Attachments.max_count )  // allocate new space if necessary
			{
				if (MechsimOptions.Attachments.max_count >= INT_MAX/2)
				{
					Error("Too many attachment functions");
				}
				MechsimOptions.Attachments.max_count *= 2;
				MechsimOptions.Attachments.Function =
					(FUNCTION_PTR *)POV_REALLOC(MechsimOptions.Attachments.Function,
																			MechsimOptions.Attachments.max_count*sizeof (FUNCTION_PTR),
																			"mechanics simulation attachment functions");
				MechsimOptions.Attachments.Fixed =
					(INT_VECT *)POV_REALLOC(MechsimOptions.Attachments.Fixed,
																	MechsimOptions.Attachments.max_count*sizeof (INT_VECT),
																	"mechanics simulation attachment vector");
			}

		  /* default: all directions attached */
		  MechsimOptions.Attachments.Fixed[MechsimOptions.Attachments.count][0] = 1;
		  MechsimOptions.Attachments.Fixed[MechsimOptions.Attachments.count][1] = 1;
		  MechsimOptions.Attachments.Fixed[MechsimOptions.Attachments.count][2] = 1;

      EXPECT
			
        CASE_EXPRESS

				  Parse_Vector(Fixed);
		      MechsimOptions.Attachments.Fixed[MechsimOptions.Attachments.count][0] = (Fixed[X] > EPSILON);
		      MechsimOptions.Attachments.Fixed[MechsimOptions.Attachments.count][1] = (Fixed[Y] > EPSILON);
		      MechsimOptions.Attachments.Fixed[MechsimOptions.Attachments.count][2] = (Fixed[Z] > EPSILON);
					if (!MechsimOptions.Attachments.Fixed[MechsimOptions.Attachments.count][0] && 
							!MechsimOptions.Attachments.Fixed[MechsimOptions.Attachments.count][1] && 
							!MechsimOptions.Attachments.Fixed[MechsimOptions.Attachments.count][2])
            Warning(0, "MechSim: Attachment %d is not fixed in any direction.", MechsimOptions.Attachments.count);
					

        END_CASE

        CASE(LEFT_CURLY_TOKEN)
          UNGET
          Parse_Begin();

		      EXPECT

						CASE (FUNCTION_TOKEN)       // Attachment Function (spline function) (parameters: time, return: 3d vector)

              MechsimOptions.Attachments.Function[MechsimOptions.Attachments.count] = Parse_SplineFunction();

							MechsimOptions.Attachments.count++;

							EXIT

						END_CASE

						OTHERWISE
							Expectation_Error ("function");
						END_CASE

					END_EXPECT

					Parse_End();
					EXIT
        END_CASE

				OTHERWISE
					UNGET
          EXIT
				END_CASE

			END_EXPECT

		END_CASE // -- end attach --

    CASE (ENVIRONMENT_TOKEN)    // ----- Environment -----
      Parse_Begin();

      if (MechsimOptions.Env_max_count==0)
			{
				MechsimOptions.Env_max_count = 10;
				MechsimOptions.Environment = 
					(MECHSIM_ENVIRONMENT *)POV_MALLOC(MechsimOptions.Env_max_count*sizeof (MECHSIM_ENVIRONMENT),
																						"mechanics simulation environment data");
			}
			else if ( MechsimOptions.Env_count >= MechsimOptions.Env_max_count )  // allocate new space if necessary
			{
				if (MechsimOptions.Env_max_count >= INT_MAX/2)
				{
					Error("Too many environments");
				}
				MechsimOptions.Env_max_count *= 2;
				MechsimOptions.Environment = 
					(MECHSIM_ENVIRONMENT *)POV_REALLOC(MechsimOptions.Environment,
																						 MechsimOptions.Env_max_count*sizeof (MECHSIM_ENVIRONMENT),
																						 "mechanics simulation environment data");
			}

			Env = &MechsimOptions.Environment[MechsimOptions.Env_count];

      Env->Function = NULL;   // default values
      Env->Object = NULL;
      Env->Stiffness = 0.0;
      Env->Damping = 0.0;
      Env->Friction = 0.0;
      Env->Friction_Excess = 1.0;
      Env->Method = MECHSIM_ENV_METHOD_FORCE;
      Env->Next = NULL;

      EXPECT

        CASE (FUNCTION_TOKEN)       // Environment Function (for collision detection)
          /* ParseParameterFunction was added in function.cpp because
             we want to be able to add an additional time parameter
             just like with Parse_DeclareFunction */
          Env->Function = Parse_ParameterFunction();
          if (Env->Function == NULL) Expectation_Error ("function");
        END_CASE

        CASE (OBJECT_TOKEN)         // alternatively: Environment object
          Env->Object = Parse_Object();
          if(Env->Object == NULL) Expectation_Error ("object");

          if ((MechsimOptions.Method != MECHSIM_METHOD_EULER) && (MechsimOptions.Method != MECHSIM_METHOD_GRADIENT))
            Warning(0, "MechSim: Using object based environments with higher order\n       integration methods is not recommended.");

        END_CASE

        CASE (STIFFNESS_TOKEN)      // Environment Stiffness
          Env->Stiffness = Parse_Float();
        END_CASE

        CASE (DAMPING_TOKEN)        // Environment Damping
          Env->Damping = Parse_Float();
        END_CASE

        CASE (FRICTION_TOKEN)       // Environment Friction
          Env->Friction = Parse_Float();
          Parse_Comma();
          Env->Friction_Excess = Allow_Float(1.0);
        END_CASE

        CASE (METHOD_TOKEN)        // Environment calculation method
          Env->Method = (unsigned int)Parse_Float();
        END_CASE

        OTHERWISE
          UNGET
          EXIT
        END_CASE

      END_EXPECT

      MechsimOptions.Env_count++;

      Parse_End();

    END_CASE // -- end environment --

    CASE (COLLISION_TOKEN)    // ----- Collision -----
      Parse_Begin();

      // collision parameter, see 'mechsim.h' for the relevant constants
      MechsimOptions.Calc_Collision = (int)Allow_Float(1.0);

      if ((MechsimOptions.Calc_Collision < 0) || (MechsimOptions.Calc_Collision > MECHSIM_COLLISION_GROUP))
      {
        Warning(0, "Unknown mechanics simulation collision flag.");
      }

      Parse_Comma();
      // optional face collision parameter
      MechsimOptions.Calc_Collision_Faces = (int)Allow_Float(-1.0);

      // default: same as MechsimOptions.Calc_Collision
      if (MechsimOptions.Calc_Collision_Faces < 0) MechsimOptions.Calc_Collision_Faces = MechsimOptions.Calc_Collision;

      if ((MechsimOptions.Calc_Collision_Faces > MECHSIM_COLLISION_GROUP))
      {
        Warning(0, "Unknown mechanics simulation collision flag.");
      }

      Parse_Comma();
      // optional connection collision parameter
      MechsimOptions.Calc_Collision_Connections = (int)Allow_Float(-1.0);

      // default: No collision tests
      if (MechsimOptions.Calc_Collision_Connections < 0) MechsimOptions.Calc_Collision_Connections = MECHSIM_COLLISION_NONE;

      if ((MechsimOptions.Calc_Collision_Connections > MECHSIM_COLLISION_GROUP))
      {
        Warning(0, "Unknown mechanics simulation collision flag.");
      }

      EXPECT

        CASE (STIFFNESS_TOKEN)      // Collision Stiffness
          MechsimOptions.Collision_Stiffness = Parse_Float();
        END_CASE

        CASE (DAMPING_TOKEN)        // Collision Damping
          MechsimOptions.Collision_Damping = Parse_Float();
        END_CASE

        CASE (FRICTION_TOKEN)
          MechsimOptions.Collision_Friction = Parse_Float();
          Parse_Comma();
          MechsimOptions.Collision_Friction_Excess = Allow_Float(1.0);
        END_CASE

        OTHERWISE
          UNGET
          EXIT
        END_CASE

      END_EXPECT
      Parse_End();
    END_CASE // -- end collision --

    CASE (TOPOLOGY_TOKEN)      // ----- Topology -----
      Parse_Begin();
      EXPECT

        CASE(GROUP_TOKEN)          // topology group

          Get_Token();
          if(Token.Token_Id == LEFT_PAREN_TOKEN)   // explicit group index given
          {
            Current_Group = (int)Parse_Float();
            Get_Token();
            if(Token.Token_Id != RIGHT_PAREN_TOKEN) Expectation_Error(")");
          }
          else                                     // otherwise count up group index
          {
            Unget_Token();
            Current_Group++;
          }

          Parse_Begin();
          EXPECT

            CASE(MASS_TOKEN)           // point mass definition
              Parse_Mech_Sim_Mass(Current_Group);
            END_CASE

            CASE(CONNECTION_TOKEN)     // connection definition
              Parse_Mech_Sim_Connection(Current_Group);
            END_CASE

            CASE(FACE_TOKEN)           // face definition
              Parse_Mech_Sim_Face(Current_Group);
            END_CASE

            CASE(VISCOELASTIC_TOKEN)   // viscoelastic definition
							Parse_Mech_Sim_Viscoelastic();
						END_CASE

            CASE(LOAD_FILE_TOKEN)        // load topology data from separate file
              MechsimOptions.Load_File_Name = Parse_C_String(true);

              fdi = New_Checked_IStream(MechsimOptions.Load_File_Name, POV_File_Unknown);
              if ( fdi != NULL)
              {
                Mechsim_read_file(fdi, Current_Group);
                delete fdi;
              }
              POV_FREE(MechsimOptions.Load_File_Name);
              MechsimOptions.Load_File_Name = NULL;

            END_CASE
            /*
            CASE(SAVE_FILE_TOKEN)        // save topology data after simulation
              MechsimOptions.Save_File_Name = Parse_C_String(true);
              MechsimOptions.Save_File = true;
            END_CASE
            */
            OTHERWISE
              UNGET
              EXIT
            END_CASE

          END_EXPECT
          Parse_End();

        END_CASE // -- end group --

        // everything without group gets group index 0

        CASE(MASS_TOKEN)           // point mass definition
          Parse_Mech_Sim_Mass(0);
        END_CASE

        CASE(CONNECTION_TOKEN)     // connection definition
          Parse_Mech_Sim_Connection(0);
        END_CASE

        CASE(FACE_TOKEN)           // face definition
          Parse_Mech_Sim_Face(0);
        END_CASE

        CASE(VISCOELASTIC_TOKEN)           // viscoelastic definition
          Parse_Mech_Sim_Viscoelastic();
        END_CASE

        CASE(LOAD_FILE_TOKEN)        // load topology data from separate file
          MechsimOptions.Load_File_Name = Parse_C_String(true);

          fdi = New_Checked_IStream(MechsimOptions.Load_File_Name, POV_File_Unknown);
          if ( fdi != NULL)
          {
            Mechsim_read_file(fdi, 0);
            delete fdi;
          }
          POV_FREE(MechsimOptions.Load_File_Name);
          MechsimOptions.Load_File_Name = NULL;

        END_CASE

        CASE(SAVE_FILE_TOKEN)        // save topology data after simulation
          MechsimOptions.Save_File_Name = Parse_C_String(true);
          Parse_Comma();
          MechsimOptions.Save_File_Type = (int)Allow_Float(2.0);
          if ((MechsimOptions.Save_File_Type < 2) || (MechsimOptions.Save_File_Type > 3))
          {
            Error("file type %d not supported\n", MechsimOptions.Save_File_Type);
          }
					/* Save_File_Type only supported for backwards compatibility - it is ignored now */
          MechsimOptions.Save_File_Type = 4;
          MechsimOptions.Save_File = true;
        END_CASE

        OTHERWISE
          UNGET
          EXIT
        END_CASE

      END_EXPECT

      Parse_End();

    END_CASE // -- end topology --


    OTHERWISE
      UNGET
      EXIT
    END_CASE

  END_EXPECT
  Parse_End();

  Mechsim_Simulate();

}

#endif

#ifdef ALTERNATIVE_VECTOR_FUNCTION

VECTOR_FUNCTION *Parse_Vector_Function (int Input,int Output)
{
  VECTOR_FUNCTION *New = Create_Vector_Function();

  EXPECT
    CASE(PIGMENT_TOKEN)
      New->IsPigment = true;
      EXIT
    END_CASE
    OTHERWISE
      UNGET
      New->IsPigment = false;
      EXIT
    END_CASE
  END_EXPECT

  New->InputDimension = Input;
  New->OutputDimension = Output;
  if(New->IsPigment)
  {
    Parse_Begin ();
    New->Data.Pigment = Create_Pigment();
    Parse_Pigment(&(New->Data.Pigment));
    Post_Pigment(New->Data.Pigment);
    Parse_End ();
  }
  else
    for(int i=0; i<Output; i++)
      New->Data.Function[i] = Parse_User_Function();

  return New;
}

FUNCTION_PTR Parse_User_Function()
{
  FUNCTION_PTR Function;

  EXPECT
    CASE(FUNCTION_TOKEN)
      Function = Parse_Function();
      EXIT
    END_CASE
    OTHERWISE
      Function = Parse_FunctionContent();
      EXIT
    END_CASE
  END_EXPECT

  Parse_Comma();

  return Function;
}

#endif

#ifdef POST_PROCESS_PATCH

void Parse_Post_Process()
{
  POST_PROCESS **Post_Process_Cursor = &First_Post_Process;
  POST_PROCESS *Post_Process=Create_Post_Process();
  int i;

  Post_Process->prev = NULL;

  i = 1;
  while(*Post_Process_Cursor)
  {
    i++;
    Post_Process->prev = *Post_Process_Cursor;
    Post_Process_Cursor = &((*Post_Process_Cursor)->next);
  }

  Parse_Begin();

  Post_Process->Function[PP_COLOUR_RED]   = Parse_User_Function();
  Post_Process->Function[PP_COLOUR_GREEN] = Parse_User_Function();
  Post_Process->Function[PP_COLOUR_BLUE]  = Parse_User_Function();
  Post_Process->Function[PP_COLOUR_ALPHA] = Parse_User_Function();
  Post_Process->Order = i;

  EXPECT
    CASE(SAVE_FILE_TOKEN)
      Post_Process->Save_File_Name = Parse_C_String(true);
      EXIT
    END_CASE
    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

  Parse_End();

  *Post_Process_Cursor = Post_Process;
}
#endif

#ifdef MOTION_BLUR_PATCH
struct InputFileData
{
	POV_BASE_NAMESPACE::ITextStream *In_File;
	bool R_Flag;
	#ifdef UNOFFICIAL_BLOCKING_PATCH
	  int unofficialVersion;
	#endif
};
extern InputFileData *Input_File;
  extern int CS_Index;
  static int beginCS_Index;
static OBJECT *Parse_Motion_Blur ()
{
  MOTION_BLUR *Object;
  OBJECT *Local;
  DBL saveClock,clockDelta;
  int i;
	if ( GlobalSettingsForMotionBlurParsed==false)
		Error("No global_settings for motion_blur found.\nMotion_blur in global_settings must preceed all motion blur objects.\n");

  int Type=0; //standard type ==shake

  Object = Create_Motion_Blur ();
  Parse_Begin();

  EXPECT
     CASE (TYPE_TOKEN)
       Type= (int)Parse_Float();
        if ( Type != 0 && Type != 1)
          Error("For type only 0 (default) and 1 allowed.");
     END_CASE

    OTHERWISE
      UNGET
      EXIT
    END_CASE
  END_EXPECT

 POV_BASE_NAMESPACE::ITextStream::FilePos File_Pos;
  File_Pos =  Token.File_Pos;
Token.TokenFileChangedAfterSaveFilePos=false;

  Object->Children = NULL;
  saveClock = opts.FrameSeq.Clock_Value;

  if (opts.motionBlurCount>1)
  {
    switch (Type)
    {
    case 1:
       clockDelta = (opts.blurClockDelta)/(opts.motionBlurCount-1);
      break;
      default:  //default megapov 0.x behavior
        opts.FrameSeq.Clock_Value -= opts.blurClockDelta/2.0;
        clockDelta = (opts.blurClockDelta)/(opts.motionBlurCount-1);
      break;
    }
  }
  else
    clockDelta = 0;

  beginCS_Index = CS_Index;
  for(i=1; i<=opts.motionBlurCount; i++)
  {
	  if ( Token.TokenFileChangedAfterSaveFilePos)
    	Error("File error! Possibly a reference to an include file in a motion_blur{} block.");

	  Input_File->In_File->seekg(File_Pos);
	  Token.Unget_Token=false;
	  if (Token.Token_Id==RIGHT_CURLY_TOKEN)
    	++Brace_Index;

    Local = Parse_Object ();

    if(Local)
    {
      Object->Type |= (Local->Type & CHILDREN_FLAGS);

      if ((Local->Type & MOTION_BLUR_OBJECT))
      {
        Error("motion_blur objects cannot be nested\n");
      }
      Local->Type |= IS_CHILD_OBJECT;
      Local->TimeStamp = i;
      Link(Local,  &Object->Children);
      /* now parse the mods for motion_blur */
      Parse_Object_Mods(Local);
    }
    opts.FrameSeq.Clock_Value += clockDelta;
    if (CS_Index != beginCS_Index)
    {
      Error("close brace for motion_blur cannot be enclosed in an #if, #case, or #while block.\n");
    }
  }

  if (!Object->Children)
    Error("No object inside motion_blur.");

  Object->Type |= MOTION_BLUR_OBJECT;

  opts.FrameSeq.Clock_Value = saveClock;

  Compute_Motion_Blur_BBox((OBJECT *)Object);

  return ((OBJECT *) Object);
}


#endif

END_POV_NAMESPACE
