From 0ee39cda7ef53b70ed466d0277d22913eddee689 Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@fe.uni-lj.si> Date: Wed, 23 Apr 2014 17:20:36 +0000 Subject: [PATCH] in debuggingg ptrocedure --- src/vertex.c | 27 +++++++++++++++++++-------- 1 files changed, 19 insertions(+), 8 deletions(-) diff --git a/src/vertex.c b/src/vertex.c index 694ee84..f733d49 100644 --- a/src/vertex.c +++ b/src/vertex.c @@ -6,6 +6,14 @@ #include "bond.h" #include<stdio.h> +ts_bool vertex_list_assign_id(ts_vertex_list *vlist, ts_uint id){ + ts_uint i; + for(i=0;i<vlist->n;i++){ + vlist->vtx[i]->id = id; + } + return TS_SUCCESS; +} + ts_vertex_list *init_vertex_list(ts_uint N){ ts_int i; ts_vertex_list *vlist=(ts_vertex_list *)malloc(sizeof(ts_vertex_list)); @@ -72,25 +80,25 @@ ts_bool vtx_remove_neighbour(ts_vertex *vtx, ts_vertex *nvtx){ /* find a neighbour */ /* remove it from the list while shifting remaining neighbours up */ - ts_int i,j=-1; + ts_uint i,j=0; for(i=0;i<vtx->neigh_no;i++){ +// fprintf(stderr,"neigh_addr=%ld\n", (long)vtx->neigh[i]); if(vtx->neigh[i]!=nvtx){ vtx->neigh[j]=vtx->neigh[i]; j++; } } - if(j==0) { - fatal("vtx_remove_neighbour: Error, vertices are not neighbours", 100); - } +// fprintf(stderr,"remove_neighbour: vtx1_addr=%ld, vtx2_addr=%ld\n",(long)vtx,(long)nvtx); /* resize memory. potentionally time consuming */ vtx->neigh_no--; vtx->neigh=(ts_vertex **)realloc(vtx->neigh,vtx->neigh_no*sizeof(ts_vertex *)); if(vtx->neigh == NULL && vtx->neigh_no!=0) - fatal("Reallocation of memory failed during removal of vertex neighbour in vtx_remove_neighbour",100); -fprintf(stderr,"first alloc"); + fatal("(1) Reallocation of memory failed during removal of vertex neighbour in vtx_remove_neighbour",100); +//fprintf(stderr,"first alloc"); /* repeat for the neighbour */ /* find a neighbour */ /* remove it from the list while shifting remaining neighbours up */ + j=0; for(i=0;i<nvtx->neigh_no;i++){ if(nvtx->neigh[i]!=vtx){ nvtx->neigh[j]=nvtx->neigh[i]; @@ -98,10 +106,12 @@ } } /* resize memory. potentionally time consuming. */ +// fprintf(stderr,"Neigbours=%d\n",nvtx->neigh_no); nvtx->neigh_no--; - nvtx->neigh=(ts_vertex **)realloc(nvtx->neigh,nvtx->neigh_no*sizeof(ts_vertex *)); + nvtx->neigh=(ts_vertex **)realloc(nvtx->neigh,nvtx->neigh_no*sizeof(ts_vertex *)); +// fprintf(stderr,"Neigbours=%d\n",nvtx->neigh_no); if(nvtx->neigh == NULL && nvtx->neigh_no!=0) - fatal("Reallocation of memory failed during removal of vertex neighbour in vtx_remove_neighbour",100); + fatal("(2) Reallocation of memory failed during removal of vertex neighbour in vtx_remove_neighbour",100); return TS_SUCCESS; } @@ -128,6 +138,7 @@ ts_bool vtx_add_cneighbour(ts_bond_list *blist, ts_vertex *vtx1, ts_vertex *vtx2){ ts_bool retval; retval=vtx_add_neighbour(vtx1,vtx2); + // retval=vtx_add_neighbour(vtx2,vtx1); if(retval==TS_SUCCESS) retval=vtx_add_bond(blist,vtx1,vtx2); return retval; -- Gitblit v1.9.3