void main(int argc, char *argv[]){
int t;
char fname1[100], fname2[100];
int m = sprintf(fname1,"phi");
int n = sprintf(fname2,"velocities");
MPI_Init(&argc,&argv);
MPI_Comm_size(MPI_COMM_WORLD,&numtasks);
MPI_Comm_rank(MPI_COMM_WORLD,&taskid);
numworkers = numtasks - 1;
allocate_memory_phasefields(MESHX);
#ifdef FLUID
allocate_memory_FluidVariables(MESHX,pmesh);
#endif
if(taskid == MASTER) {
initialize_phasefields();
#ifdef FLUID
initialize_velocities(MESHX);
#endif
mpi_distribute(MESHX);
for (t = 0; t <= phi_timesteps; t++){
#ifdef FLUID
gauss_siedel();
#endif
if (t%save_phi == 0) {
receivefrmworker(phi_old);
receivefrmworker(u_old);
receivefrmworker(v_old);
write2file( t, MESHX, phi_old, fname1);
write2file1( t, MESHX, u_old, v_old, fname2);
}
}
}else {
mpi_distribute(MESHX);
for (t = 0; t <= phi_timesteps; t++){
mpiexchange(taskid, phi_old, MESHX);
mpiexchange(taskid, mu_old, MESHX);
boundary_mpi(taskid, phi_old, MESHX);
boundary_mpi(taskid, mu_old, MESHX);
laplacian(phi_old, lap_phi, MESHX);
laplacian(mu_old, lap_mu, MESHX);
#ifdef FLUID
computeH(u_old,v_old,Hx,Hy);
RHS_fn(Hx,Hy,rhs_fn,MESHX, pmesh);
LHS_fn(MESHX, pmesh);
boundary_pressure(taskid);
gauss_siedel();
ns_solver(start, end, phi_old);
update_velocities(MESHX);
#endif
solverloop();
phi_update();
if (t%save_phi == 0) {
sendtomaster(taskid, phi_old);
}
}
}
}
void
initialize()
{
set_limits();
initialize_velocities();
randomize_particles();
initialize_pbest();
intialize_gbest();
}
