| | |
|---|
| | | #include "sh.h" |
|---|
| | | #include "vesicle.h" |
|---|
| | | #include "initial_distribution.h" |
|---|
| | | #include <math.h> |
|---|
| | | |
|---|
| | | int main(){ |
|---|
| | | ts_vesicle *vesicle=initial_distribution_dipyramid(4,60,60,60,0.15); |
|---|
| | |
|---|
| | | fprintf(stderr,"co(%d,%d)=%e\n",i,j,vesicle->sphHarmonics->co[i][j]); |
|---|
| | | } |
|---|
| | | } |
|---|
| | | ts_double r0; |
|---|
| | | vesicle_volume(vesicle); |
|---|
| | | fprintf(stderr,"Volume=%e\n",vesicle->volume); |
|---|
| | | r0=getR0(vesicle); |
|---|
| | | fprintf(stderr,"r0=%e\n",r0); |
|---|
| | | preparationSh(vesicle,r0); |
|---|
| | | calculateYlmi(vesicle); |
|---|
| | | calculateUlm(vesicle); |
|---|
| | | |
|---|
| | | for(i=0;i<vesicle->sphHarmonics->l;i++){ |
|---|
| | | for(j=0;j<2*i+1;j++){ |
|---|
| | | fprintf(stderr,"ulm(%d,%d)=%e\n",i,j+1,vesicle->sphHarmonics->ulm[i][j]); |
|---|
| | | } |
|---|
| | | } |
|---|
| | | |
|---|
| | | |
|---|
| | | |
|---|
| | | sph_free(vesicle->sphHarmonics); |
|---|
| | | //vesicle_free(vesicle); |
|---|
| | | vesicle_free(vesicle); |
|---|
| | | return 0; |
|---|
| | | } |
|---|
