| | |
| | | } |
| | | xmlFree(values); |
| | | } |
| | | if(!xmlStrcmp(property_name,(const xmlChar *)"direct_interaction_force")){ |
| | | values=xmlNodeListGetString(doc,child->xmlChildrenNode,1); |
| | | vals=(char *)values; |
| | | token=strtok(vals," "); |
| | | idx=0; |
| | | while(token!=NULL){ |
| | | if(idx<vesicle->vlist->n){ |
| | | vesicle->vlist->vtx[idx]->direct_interaction_force=atof(token); |
| | | } else if(vesicle->tape->nmono && vesicle->tape->npoly && idx<vesicle->vlist->n+vesicle->tape->nmono*vesicle->tape->npoly) { |
| | | polyidx=(idx-vesicle->vlist->n)/vesicle->tape->nmono; |
| | | monoidx=(idx-vesicle->vlist->n)%vesicle->tape->nmono; |
| | | vesicle->poly_list->poly[polyidx]->vlist->vtx[monoidx]->direct_interaction_force=atof(token); |
| | | } else { |
| | | filidx=(idx-vesicle->vlist->n-vesicle->tape->nmono*vesicle->tape->npoly)/vesicle->tape->nfono; |
| | | fonoidx=(idx-vesicle->vlist->n-vesicle->tape->nmono*vesicle->tape->npoly)%vesicle->tape->nfono; |
| | | //fprintf(stderr,"filidx=%d, fonoidx=%d, coord=%s,%s,%s\n",filidx,fonoidx,token[0],token[1],token[2]); |
| | | vesicle->filament_list->poly[filidx]->vlist->vtx[fonoidx]->direct_interaction_force=atof(token); |
| | | } |
| | | idx++; |
| | | token=strtok(NULL," "); |
| | | } |
| | | xmlFree(values); |
| | | } |
| | | xmlFree(property_name); |
| | | } |
| | | |