Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- CloseKernels[];
- LaunchKernels[8];
- xres2 = 8;
- yres2 = 8;
- l = 1;
- indTab = Flatten[Table[{ll, i, j}, {ll, Length[res1C[[l]]]}, {i, xres2 - 1}, {j, yres2 - 1}], 2];
- DistributeDefinitions["Global`*"];
- ParallelDo[
- ll = indTab[[ij, 1]];
- i = indTab[[ij, 2]];
- j = indTab[[ij, 3]];
- x0 = res1C[[l, ll, 1]];
- y0 = res1C[[l, ll, 2]];
- P1 = smplGrd[[x0 + 0, y0 + 0]];
- P2 = smplGrd[[x0 + 1, y0 + 1]];
- xint2 = Re[P1];
- xfin2 = Re[P2];
- xinc2 = (xfin2 - xint2)/(xres2 - 1);
- x2 = Table[xint2 + (i - 1) xinc2, {i, xres2}];
- yint2 = Im[P1];
- yfin2 = Im[P2];
- yinc2 = (yfin2 - yint2)/(yres2 - 1);
- y2 = Table[yint2 + (i - 1) yinc2, {i, yres2}];
- smplGrd2 = Table[x2[[ii]] + I y2[[jj]], {ii, xres2}, {jj, yres2}];
- OMint = Re[smplGrd2[[i + 0, j + 0]] nrm];
- OMfin = Re[smplGrd2[[i + 1, j + 0]] nrm];
- GMint = Im[smplGrd2[[i + 0, j + 0]] nrm];
- GMfin = Im[smplGrd2[[i + 0, j + 1]] nrm];
- Print[{OMint, OMfin, GMint, GMfin}]
- , {ij, (xres2 - 1) (yres2 - 1) Length[res1C[[l]]]}];
- Set::wrsym: Symbol smplGrd2 is Protected.
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement