double *FEM(double (*DENSITY),int sx,int sy,int sz,double ratio){
double ***STRAIN_ENERGY;
double MAX_STRAIN_ENERGY=0.;
double MIN_STRAIN_ENERGY=0.;
float angle = 0.0;
int rotate = 0;
//int mode = DEFAULT_LIGHT;
double ratio2;
////////////////////////////////double ratio; MOVED TO OPTIMIZATION_MAIN
////////////////////////////////int sx,sy,sz; MOVED TO OPTIMIZATION_MAIN
//double vx,vy,vz;
double K_matrix[ELDOF_SIZE][ELDOF_SIZE]={0.};
double (*K_matrix_temp)[ELDOF_SIZE];
//double deriv_matrix[DIMENSION][ELNODE_SIZE]={0.};
double (*deriv_matrix)[ELNODE_SIZE];
double (*deriv_matrix2)[ELNODE_SIZE];
double (*jacobian)[DIMENSION];
double (*Ijacobian)[DIMENSION];
double (*B_matrix)[ELDOF_SIZE];
double (*BT_matrix)[C_SIZE];
double (*C_matrix)[C_SIZE];
double (*CB_matrix)[ELDOF_SIZE];
double (*BTCB_matrix)[ELDOF_SIZE];
double **GLOBAL_K;
double *GLOBAL_U;
double *GLOBAL_R;
double *GLOBAL_U_ALLDOFS;
double *COORDS;
double *ELEM_DISP;
double *DISPLACEDCOORDS;
static double *FITNESS;
double *error;
double *YOUNGS_MODULUS;
double temp1;
double rmin=1.5;
int i,j,k,ii,jj,l,m,n;
int counter;
int i1;
int i2,j2,k2;
int start1,end1;
int start2,end2;
int start3,end3;
int ng=4;
int elem;
int temp;
int temp2;
int node[9]={0};
int LM[25]={0};//CONNECTIVITY MATRIX
int pinnodes[1000]={0};
int forcenodes[1000]={0};
int pinned=0;
int loaded=0;
int index;//size of the reduced K
int index2;//size of the reduced R
int index3;//used for U_ALLDOFS
int index4;//used for U_ALLDOFS
//for solver//
double old;
double sum;
double E=200.E6;
double v = 0.3;
double penal=1.;
double det=0;
double integration_factor=0;
double weight[5]={0,0.34785,0.65215,0.65215,0.34785};
printf("Starting FEM Analysis..\\n");
elem = (sx+1)*(sy+1)*(sz+1) * 3;
//DEFINES KU=R//
GLOBAL_K = (double **)malloc(sizeof(double *)*(elem+sizeof(double *)));
printf("1\\n");
GLOBAL_U = (double *)malloc(sizeof(double)*(elem+sizeof(double)));
printf("2\\n");
GLOBAL_R = (double *)malloc(sizeof(double)*(elem+sizeof(double)));
printf("3\\n");
GLOBAL_U_ALLDOFS = (double *)malloc(sizeof(double)*(elem+sizeof(double)));
COORDS = (double *)malloc(sizeof(double)*(elem+sizeof(double)));
DISPLACEDCOORDS = (double *)malloc(sizeof(double)*(elem+sizeof(double)));
ELEM_DISP = (double *)malloc(sizeof(double)*(ELDOF_SIZE));
printf("4\\n");
YOUNGS_MODULUS = (double *)malloc(sizeof(double)*((sz*sx*sy)+sizeof(double)));
FITNESS = (double *)malloc(sizeof(double)*((sz*sx*sy)+sizeof(double)));
//3차원배열 동적할당
STRAIN_ENERGY = (double***)malloc(sizeof(double**)* (sz+sizeof(double**)) );
STRAIN_ENERGY[0] = (double**)malloc( sizeof(double*)* (sz+sizeof(double*))*(sy+sizeof(double*)) );
STRAIN_ENERGY[0][0] = (double*)malloc(sizeof(double)*(sz+sizeof(double))*(sy+sizeof(double))*(sx+sizeof(double)));
//ASSIGNS POINTERS///
for(j=1; j <= sy; j++) {
STRAIN_ENERGY[0][j] = STRAIN_ENERGY[0][j-1] + sx;
}
for(i=1; i <= sz; i++) {
STRAIN_ENERGY[i] = STRAIN_ENERGY[i-1] + sy;
STRAIN_ENERGY[i][0] = STRAIN_ENERGY[i-1][0] + sy * sx;
for(j=1; j <= sy; j++) {
STRAIN_ENERGY[i][j] = STRAIN_ENERGY[i][j-1] + sx;
}
}
/////////////////////
error = (double *)malloc(sizeof(double)*(elem+sizeof(double)));
//MEMORY ALLOCATION ENDS//
확인해보니 제일 처음 나오는 malloc 부분에서 멈추는것 같습니다.
메모리 할당에 문제가 있는것같은데, 이것때문에 죽겠습니다. 도대체 제가 뭘 잘못짠건지 모르겠네요.
전공이 CS가 아니라 컴퓨터 내부에서 메모리를 가지고 어떻게 하는건지는 기본적인 지식밖에 없어서 코딩을 할때 뭐가 잘못한것같기도 한데 도대체 뭐가 잘못된건지 아무리 뚫어져라 쳐다봐도 모르겠네요. 해결법이 없을까요? 제가 뭘 잘못하고 있는건지 좀 알려주세요!
댓글 영역
획득법
① NFT 발행
작성한 게시물을 NFT로 발행하면 일주일 동안 사용할 수 있습니다. (최초 1회)
② NFT 구매
다른 이용자의 NFT를 구매하면 한 달 동안 사용할 수 있습니다. (구매 시마다 갱신)
사용법
디시콘에서지갑연결시 바로 사용 가능합니다.