| | |
|---|
| | | /* vim: set ts=4 sts=4 sw=4 noet : */ |
|---|
| | | #include<math.h> |
|---|
| | | #include<stdlib.h> |
|---|
| | | #include "general.h" |
|---|
| | | #include "sh.h" |
|---|
| | | |
|---|
| | | #include "io.h" |
|---|
| | | #include <string.h> |
|---|
| | | |
|---|
| | | |
|---|
| | | ts_spharm *sph_init(ts_vertex_list *vlist, ts_uint l){ |
|---|
| | |
|---|
| | | #ifdef TS_DOUBLE_DOUBLE |
|---|
| | | coord->e1=sqrt(x*x+y*y+z*z); |
|---|
| | | if(z==0) coord->e3=M_PI/2.0; |
|---|
| | | else coord->e3=atan(sqrt(x*x+y*y)/z); |
|---|
| | | else coord->e3=atan2(sqrt(x*x+y*y),z); |
|---|
| | | coord->e2=atan2(y,x); |
|---|
| | | #endif |
|---|
| | | #ifdef TS_DOUBLE_FLOAT |
|---|
| | |
|---|
| | | |
|---|
| | | return TS_SUCCESS; |
|---|
| | | } |
|---|
| | | |
|---|
| | | |
|---|
| | | ts_bool sph2cart(ts_coord *coord){ |
|---|
| | | coord->coord_type=TS_COORD_CARTESIAN; |
|---|
| | | ts_double x,y,z; |
|---|
| | | |
|---|
| | | x=coord->e1*cos(coord->e2)*sin(coord->e3); |
|---|
| | | y=coord->e1*sin(coord->e2)*sin(coord->e3); |
|---|
| | | z=coord->e1*cos(coord->e3); |
|---|
| | | |
|---|
| | | coord->e1=x; |
|---|
| | | coord->e2=y; |
|---|
| | | coord->e3=z; |
|---|
| | | |
|---|
| | | return TS_SUCCESS; |
|---|
| | | } |
|---|
| | | |
|---|
| | | |
|---|
| | | /* Function returns radius of the sphere with the same volume as vesicle (r0) */ |
|---|
| | | ts_double getR0(ts_vesicle *vesicle){ |
|---|
| | |
|---|
| | | ts_bool saveAvgUlm2(ts_vesicle *vesicle){ |
|---|
| | | |
|---|
| | | FILE *fh; |
|---|
| | | |
|---|
| | | fh=fopen("sph2out.dat", "w"); |
|---|
| | | char filename[10000]; |
|---|
| | | strcpy(filename, command_line_args.path); |
|---|
| | | strcat(filename, "sph2out.dat"); |
|---|
| | | fh=fopen(filename, "w"); |
|---|
| | | if(fh==NULL){ |
|---|
| | | err("Cannot open file %s for writing"); |
|---|
| | | return TS_FAIL; |
|---|
