Codice: Seleziona tutto
typedef double Ele;
typedef Ele* PEle;
typedef PEle* PPEle;
int row = 5;
int col =3;
typedef lapack_complex_double* complessa;
void DataAllocation(PPEle A, PPEle B, PPEle U, PPEle V, complessa AC, complessa BC, complessa UC, complessa VC){
int i,j,k;
A=(PPEle)malloc(sizeof(PEle)*row+1);
for (i=0; i <row; i++){
A[i] = (PEle)malloc(sizeof(Ele)*col+1);
for(j=0;j<col;j++){
A[i][j]=0;
}
}
B=(PPEle)malloc(sizeof(PEle)*row+1);
for (i=0; i <row; i++){
B[i] = (PEle)malloc(sizeof(Ele)*col+1);
for(j=0;j<col;j++){
B[i][j]=0;
}
}
U=(PPEle)malloc(sizeof(PEle)*row+1);
for (i=0; i <row; i++){
U[i] = (PEle)malloc(sizeof(Ele)*row+1);
for(j=0;j<col;j++){
B[i][j]=0;
}
}
V=(PPEle)malloc(sizeof(PEle)*col+1);
for (i=0; i <col; i++){
V[i] = (PEle)malloc(sizeof(Ele)*col+1);
for(j=0;j<col;j++){
V[i][j]=0;
}
}
AC=(complessa)malloc(sizeof(lapack_complex_double)*col*row+1);
BC=(complessa)malloc(sizeof(lapack_complex_double)*col*row+1);
UC=(complessa)malloc(sizeof(lapack_complex_double)*row*row+1);
VC=(complessa)malloc(sizeof(lapack_complex_double)*col*col+1);
return;
}