| | |
| | | /* vim: set ts=4 sts=4 sw=4 noet : */ |
| | | #include "general.h" |
| | | #include "sh.h" |
| | | #include <stdlib.h> |
| | |
| | | |
| | | preparationSh(vesicle,r0); |
| | | calculateYlmi(vesicle); |
| | | /* |
| | | ts_coord *coord=(ts_coord *)malloc(sizeof(ts_coord)); |
| | | ts_double fi, theta; |
| | | |
| | | for(i=0;i<vesicle->vlist->n;i++){ |
| | | |
| | | cart2sph(coord,vesicle->vlist->vtx[i]->x, vesicle->vlist->vtx[i]->y, vesicle->vlist->vtx[i]->z); |
| | |
| | | } |
| | | |
| | | free(coord); |
| | | */ |
| | | 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]); |
| | | } |
| | | } |
| | | |
| | | |
| | | storeUlm2(vesicle); |
| | | |
| | | sph_free(vesicle->sphHarmonics); |
| | | vesicle_free(vesicle); |