123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410 |
- /*
- ****************************************************************************
- *
- * MODULE: Vector library
- *
- * AUTHOR(S): Dave Gerdes, CERL.
- * Update to GRASS 5.7 Radim Blazek.
- *
- * PURPOSE: Lower level functions for reading/writing/manipulating vectors.
- *
- * COPYRIGHT: (C) 2001 by the GRASS Development Team
- *
- * This program is free software under the GNU General Public
- * License (>=v2). Read the file COPYING that comes with GRASS
- * for details.
- *
- *****************************************************************************/
- #include <stdlib.h>
- #include <grass/vector.h>
- /* These routines all eventually call calloc() to allocate and zero
- ** the new space. BUT It is not neccessarily safe to assume that
- ** the memory will be zero. The next memory location asked for could
- ** have been previously used and not zeroed. (e.g. compress())
- */
- /* alloc_node (map, add)
- ** alloc_line (map, add)
- ** alloc_area (map, add)
- ** alloc_points (map, num)
- ** node_alloc_line (node, add)
- ** area_alloc_line (node, add)
- **
- ** Allocate array space to add 'add' elements
- */
- /* allocate new node structure */
- struct P_node *dig_alloc_node()
- {
- struct P_node *Node;
- Node = (struct P_node *) G_malloc(sizeof(struct P_node));
- if (Node == NULL)
- return NULL;
- Node->n_lines = 0;
- Node->alloc_lines = 0;
- Node->lines = NULL;
- Node->angles = NULL;
- return (Node);
- }
- /* free node structure */
- void dig_free_node(struct P_node *Node)
- {
- if (Node->alloc_lines > 0) {
- G_free(Node->lines);
- G_free(Node->angles);
- }
- G_free(Node);
- }
- /* dig_node_alloc_line (node, add)
- ** allocate space in P_node, lines and angles arrays to add 'add' more
- ** lines
- **
- ** Returns 0 ok or -1 on error
- */
- int dig_node_alloc_line(struct P_node * node, int add)
- {
- int num;
- char *p;
- G_debug(3, "dig_node_alloc_line(): add = %d", add);
- num = node->n_lines + add;
- p = G_realloc(node->lines, num * sizeof(plus_t));
- if (p == NULL)
- return -1;
- node->lines = (plus_t *) p;
- p = G_realloc(node->angles, num * sizeof(float));
- if (p == NULL)
- return -1;
- node->angles = (float *)p;
- node->alloc_lines = num;
- return (0);
- }
- /* Reallocate array of pointers to nodes.
- * Space for 'add' number of nodes is added.
- *
- * Returns: 0 success
- * -1 error
- */
- int dig_alloc_nodes(struct Plus_head *Plus, int add)
- {
- int size;
- char *p;
- size = Plus->alloc_nodes + 1 + add;
- p = G_realloc(Plus->Node, size * sizeof(struct P_node *));
- if (p == NULL)
- return -1;
- Plus->Node = (struct P_node **) p;
- Plus->alloc_nodes = size - 1;
- return (0);
- }
- /* allocate new line structure */
- struct P_line *dig_alloc_line()
- {
- struct P_line *Line;
- Line = (struct P_line *) G_malloc(sizeof(struct P_line));
- if (Line == NULL)
- return NULL;
- return (Line);
- }
- /* free line structure */
- void dig_free_line(struct P_line *Line)
- {
- G_free(Line);
- }
- /* Reallocate array of pointers to lines.
- * Space for 'add' number of lines is added.
- *
- * Returns: 0 success
- * -1 error
- */
- int dig_alloc_lines(struct Plus_head *Plus, int add)
- {
- int size;
- char *p;
- size = Plus->alloc_lines + 1 + add;
- p = G_realloc(Plus->Line, size * sizeof(struct P_line *));
- if (p == NULL)
- return -1;
- Plus->Line = (struct P_line **) p;
- Plus->alloc_lines = size - 1;
- return (0);
- }
- /* Reallocate array of pointers to areas.
- * Space for 'add' number of areas is added.
- *
- * Returns: 0 success
- * -1 error
- */
- int dig_alloc_areas(struct Plus_head *Plus, int add)
- {
- int size;
- char *p;
- size = Plus->alloc_areas + 1 + add;
- p = G_realloc(Plus->Area, size * sizeof(struct P_area *));
- if (p == NULL)
- return -1;
- Plus->Area = (struct P_area **) p;
- Plus->alloc_areas = size - 1;
- return (0);
- }
- /* Reallocate array of pointers to isles.
- * Space for 'add' number of isles is added.
- *
- * Returns: 0 success
- * -1 error
- */
- int dig_alloc_isles(struct Plus_head *Plus, int add)
- {
- int size;
- char *p;
- G_debug(3, "dig_alloc_isle():");
- size = Plus->alloc_isles + 1 + add;
- p = G_realloc(Plus->Isle, size * sizeof(struct P_isle *));
- if (p == NULL)
- return -1;
- Plus->Isle = (struct P_isle **) p;
- Plus->alloc_isles = size - 1;
- return (0);
- }
- /* allocate new area structure */
- struct P_area *dig_alloc_area()
- {
- struct P_area *Area;
- Area = (struct P_area *) G_malloc(sizeof(struct P_area));
- if (Area == NULL)
- return NULL;
- Area->n_lines = 0;
- Area->alloc_lines = 0;
- Area->lines = NULL;
- Area->alloc_isles = 0;
- Area->n_isles = 0;
- Area->isles = NULL;
- Area->centroid = 0;
- return (Area);
- }
- /* free area structure */
- void dig_free_area(struct P_area *Area)
- {
- if (Area->alloc_lines > 0)
- free(Area->lines);
- if (Area->alloc_isles > 0)
- free(Area->isles);
- G_free(Area);
- }
- /* alloc new isle structure */
- struct P_isle *dig_alloc_isle()
- {
- struct P_isle *Isle;
- Isle = (struct P_isle *) G_malloc(sizeof(struct P_isle));
- if (Isle == NULL)
- return NULL;
- Isle->n_lines = 0;
- Isle->alloc_lines = 0;
- Isle->lines = NULL;
- Isle->area = 0;
- return (Isle);
- }
- /* free isle structure */
- void dig_free_isle(struct P_isle *Isle)
- {
- if (Isle->alloc_lines > 0)
- G_free(Isle->lines);
- G_free(Isle);
- }
- /* allocate room for 'num' X and Y arrays in struct line_pnts
- ** returns -1 on out of memory
- */
- int dig_alloc_points(struct line_pnts *points, int num)
- {
- int alloced;
- char *p;
- alloced = points->alloc_points;
- /* alloc_space will just return if no space is needed */
- if (!(p =
- dig__alloc_space(num, &alloced, 50, (char *)points->x,
- sizeof(double)))) {
- return (dig_out_of_memory());
- }
- points->x = (double *)p;
- alloced = points->alloc_points;
- /* alloc_space will just return if no space is needed */
- if (!(p =
- dig__alloc_space(num, &alloced, 50, (char *)points->y,
- sizeof(double)))) {
- return (dig_out_of_memory());
- }
- points->y = (double *)p;
- alloced = points->alloc_points;
- /* alloc_space will just return if no space is needed */
- if (!(p =
- dig__alloc_space(num, &alloced, 50, (char *)points->z,
- sizeof(double)))) {
- return (dig_out_of_memory());
- }
- points->z = (double *)p;
- points->alloc_points = alloced;
- return (0);
- }
- /* allocate room for 'num' fields and category arrays
- ** in struct line_cats
- ** returns -1 on out of memory
- */
- int dig_alloc_cats(struct line_cats *cats, int num)
- {
- int alloced;
- char *p;
- /* alloc_space will just return if no space is needed */
- alloced = cats->alloc_cats;
- if (!(p =
- dig__alloc_space(num, &alloced, 1, (int *)cats->field,
- sizeof(int)))) {
- return (dig_out_of_memory());
- }
- cats->field = (int *)p;
- alloced = cats->alloc_cats;
- if (!(p =
- dig__alloc_space(num, &alloced, 1, (int *)cats->cat,
- sizeof(int)))) {
- return (dig_out_of_memory());
- }
- cats->cat = (int *)p;
- cats->alloc_cats = alloced;
- return (0);
- }
- /* area_alloc_line (area, add)
- ** allocate space in P_area for add new lines
- **
- ** Returns 0 ok or -1 on error
- */
- int dig_area_alloc_line(struct P_area * area, int add)
- {
- int num;
- char *p;
- num = area->alloc_lines + add;
- p = G_realloc(area->lines, num * sizeof(plus_t));
- if (p == NULL)
- return -1;
- area->lines = (plus_t *) p;
- area->alloc_lines = num;
- return (0);
- }
- /* area_alloc_isle (area, add)
- ** allocate space in P_area for add new isles
- **
- ** Returns 0 ok or -1 on error
- */
- int dig_area_alloc_isle(struct P_area * area, int add)
- {
- int num;
- char *p;
- G_debug(5, "dig_area_alloc_isle(): add = %d", add);
- num = area->alloc_isles + add;
- p = G_realloc(area->isles, num * sizeof(plus_t));
- if (p == NULL)
- return -1;
- area->isles = (plus_t *) p;
- area->alloc_isles = num;
- return (0);
- }
- /* dig_isle_alloc_line (isle, add)
- ** allocate space in P_isle for add new lines
- **
- ** Returns 0 ok or -1 on error
- */
- int dig_isle_alloc_line(struct P_isle * isle, int add)
- {
- int num;
- char *p;
- G_debug(3, "dig_isle_alloc_line():");
- num = isle->alloc_lines + add;
- p = G_realloc(isle->lines, num * sizeof(plus_t));
- if (p == NULL)
- return -1;
- isle->lines = (plus_t *) p;
- isle->alloc_lines = num;
- return (0);
- }
- /* for now just print message and return error code */
- int dig_out_of_memory()
- {
- fprintf(stderr, "OUT OF MEMORY!\n");
- return (-1);
- }
|