/**********************************************************************
 * $Id$
 *
 * PostGIS - Spatial Types for PostgreSQL
 * http://postgis.refractions.net
 * Copyright 2001-2006 Refractions Research Inc.
 *
 * This is free software; you can redistribute and/or modify it under
 * the terms of the GNU General Public Licence. See the COPYING file.
 * 
 **********************************************************************/

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include "liblwgeom.h"

LWMCURVE *
lwmcurve_deserialize(uchar *srl)
{
        LWMCURVE *result;
        LWGEOM_INSPECTED *insp;
        int stype;
        int type = lwgeom_getType(srl[0]);
        int i;

        if(type != MULTICURVETYPE)
        {
                lwerror("lwmcurve_deserialize called on NON multicurve: %d", type);
                return NULL;
        }

        insp = lwgeom_inspect(srl);

        result = lwalloc(sizeof(LWMCURVE));
        result->type = insp->type;
        result->SRID = insp->SRID;
        result->ngeoms = insp->ngeometries;
        result->geoms = lwalloc(sizeof(LWCURVE *)*insp->ngeometries);

        if(lwgeom_hasBBOX(srl[0]))
        {
                result->bbox = lwalloc(sizeof(BOX2DFLOAT4));
                memcpy(result->bbox, srl+1, sizeof(BOX2DFLOAT4));
        }
        else result->bbox = NULL;

        for(i = 0; i < insp->ngeometries; i++)
        {
                stype = lwgeom_getType(insp->sub_geoms[i][0]);
                if(stype == CURVETYPE)
                {
                        result->geoms[i] = (LWGEOM *)lwcurve_deserialize(insp->sub_geoms[i]);
                }
                else if(stype == LINETYPE)
                {
                        result->geoms[i] = (LWGEOM *)lwline_deserialize(insp->sub_geoms[i]);
                }
                else
                {
                        lwerror("Only Circular and Line strings are currenly permitted in a MultiCurve.");
                        free(result);
                        free(insp);
                        return NULL;
                }
                        
                if(TYPE_NDIMS(result->geoms[i]->type) != TYPE_NDIMS(result->type))
                {
                        lwerror("Mixed diminsions (multicurve: %d, curve %d:%d)",
                                TYPE_NDIMS(result->type), i,
                                TYPE_NDIMS(result->geoms[i]->type));
                        free(result);
                        free(insp);
                        return NULL;
                }
        }
        return result;
}

/*
 * Add 'what' to this multicurve at position 'where'.
 * where=0 == prepend
 * where=-1 == append
 * Returns a MULTICURVE or a COLLECTION
 */
LWGEOM *
lwmcurve_add(const LWMCURVE *to, uint32 where, const LWGEOM *what)
{
        LWCOLLECTION *col;
        LWGEOM **geoms;
        int newtype;
        uint32 i;

        if(where == -1) where = to->ngeoms;
        else if(where < -1 || where > to->ngeoms)
        {
                lwerror("lwmcurve_add: add position out of range %d..%d",
                        -1, to->ngeoms);
                return NULL;
        }

        /* dimensions compatibility are checked by caller */

        /* Construct geoms array */
        geoms = lwalloc(sizeof(LWGEOM *)*(to->ngeoms+1));
        for(i = 0; i < where; i++)
        {
                geoms[i] = lwgeom_clone((LWGEOM *)to->geoms[i]);
        }
        geoms[where] = lwgeom_clone(what);
        for(i = where; i < to->ngeoms; i++) 
        {
                geoms[i+1] = lwgeom_clone((LWGEOM *)to->geoms[i]);
        }

        if(TYPE_GETTYPE(what->type) == CURVETYPE) newtype = MULTICURVETYPE;
        else newtype = COLLECTIONTYPE;

        col = lwcollection_construct(newtype,
                to->SRID, NULL,
                to->ngeoms + 1, geoms);

        return (LWGEOM *)col;
}