| | |
| | | #include<math.h> |
| | | #include"energy.h" |
| | | |
| | | ts_bool poly_assign_filament_xi(ts_vesicle *vesicle, ts_tape *tape){ |
| | | ts_uint i; |
| | | |
| | | for(i=0;i<vesicle->filament_list->n;i++){ |
| | | vesicle->filament_list->poly[i]->k = tape->xi; |
| | | } |
| | | |
| | | return TS_SUCCESS; |
| | | } |
| | | |
| | | |
| | | ts_bool poly_assign_spring_const(ts_vesicle *vesicle){ |
| | | ts_uint i; |
| | | |
| | |
| | | else |
| | | { |
| | | /* Make filaments inside the vesicle. Helix with radius... Dist. between poly vertices put to 1*/ |
| | | dphi = 2.0*asin(1.0/2.0/vesicle->R_nucleus)*1.001; |
| | | dh = dphi/2.0/M_PI*1.001; |
| | | dphi = 2.0*asin(1.0/2.0/vesicle->R_nucleus)*1.01; |
| | | dh = dphi/2.0/M_PI*1.01; |
| | | for(i=0;i<poly_list->n;i++){ |
| | | for (j=0;j<poly_list->poly[i]->vlist->n;j++){ |
| | | ji = j + i*poly_list->poly[i]->vlist->n; |