}
PrefixWrite( Commands[ RUN_SLOT ], "name ", 5 );
if( Commands[ MAP_SLOT ] != NULL ) {
PrefixWrite( Commands[ MAP_SLOT ], "option map=", 11 );
MapOption = FALSE; // make sure it isn't generated a second time.
}
PrefixWrite( Commands[ LIBRARY_SLOT ], "library ", 8 );
if( MapOption ) {
CommandOut( "option map" );
}
ListWrite( Commands[ OPTION_SLOT ] );
}
#define NL "\r\n"
static char TheHelp[] = {
banner1( "Microsoft to Watcom linker command translation utility ",
_MS2WLINK_VERSION_ ) NL
banner2( "1990" ) NL
banner3 NL
"usage: pass this program the same arguments that would be passed to the" NL
" Microsoft linker, and an equivalent set of wlink commands" NL
" will be written to standard output." NL
NL
"Usage:" NL
NL
"ms2wlink" NL
"ms2wlink @<response file>" NL
"ms2wlink <objs>,<exefile>,<mapfile>,<libs>,<deffile>" NL
NL
"Recognized options are:" NL
};
int main(){
int i, j, iter,isource,idir,step,irec;
clock_t start = clock();
setbuf(stdout,NULL);
/*** Initial Condition ***/
init();
printf("Start calculation\n");
/*** Calculating for iteratively ***/
for(iter=0;iter<MAXITR;iter++){
printf("%03d th iterative\n",iter);
media_coeff_3d();
init_iterate();
init_FILE3(iter);
//*********** Start backward calculation **********//
for(irec=0;irec<rec_num;irec++){
for(idir=0;idir<NDIRECT;idir++){
banner1(idir,irec);
set_zero_eh();
/*** Calculating for step ***/
for(step=0;step<it-1;step++){
backpropagation(idir, irec, step);
copytoEcal_b(EX,EY,EZ,step,idir, irec,ofe1,ofe2,ofe3);
if(step%(it/20)==it/20-1) printf("#");
}
printf("\n");
laplaceToFreq_back(irec,idir);
output_backwave(irec,idir);
} // END loop for direction
}// END loop for ireceiver
close_FILE3();
//*********** Start of FWD calculation **********//
/*** Calculating for each source ***/
for(isource=0;isource<shot_num;isource++){
printf("This is %2d loop for Transmitter \n",isource);
read_Eobs(isource); // Read observed value of Ex
//init_FILE2(isource, iter);
printf("[PROGRESS =>] ");
set_zero_eh();
media_coeff_3d();
/*** Calculating for step ***/
for(step=0;step<it-1;step++){
fwdpropagation(isource, step);
copytoEcal(EX,EY,EZ,step);
if(step%(it/20)==it/20-1) printf("#");
} // END loop for step
laplaceToFreq_2(isource);
output_fwdwave(isource);
//close_FILE2();
// Calculation of residual error
residualE(iter);
sensitivity();
// Convolution of Greenfunc and conjugated delE
conv_GdelE();
// Calculation of gamma
calc_gamma();
} // END loop for isource
// Calclulation of delta model
calc_phi();
//printf("sig2 %e\n",sig[0]);
media_coeff_sig_tmp();
//*********** again FWD calculation **********//
for(isource=0;isource<shot_num;isource++){
printf("This is %2d loop for AGAIN Transmitter \n",isource);
printf("[PROGRESS =>] ");
set_zero_eh();
/*** Calculating for step ***/
for(step=0;step<it-1;step++){
fwdpropagation(isource, step);
copytoEcal(EX,EY,EZ,step);
if(step%(it/20)==it/20-1) printf("#");
} // END loop for step
laplaceToFreq_3();
calc_alpha(isource);
} // END loop for isource
update_para2(iter);
show_error(iter);
} // END loop for iter
printf("\n");
printf("%f [s] \n",(double)(clock()-start)/CLOCKS_PER_SEC);
fclose(ofer);
}
