Advertisement
Guest User

Untitled

a guest
Jan 24th, 2017
85
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.78 KB | None | 0 0
  1. func[var_] := Module[{x = var},
  2. gamma = 176;
  3. alphag = 0.01;
  4. alphajConstant = 0.00603;
  5. p = {Cos[Pi/6], 0, Sin[Pi/6]};
  6. current[t_] := Simplify`PWToUnitStep@Piecewise[{{3, t <= x}, {(5 - 3)/(150 - x) (t - x) + 3, x < t < 150}, {5, t >= 150}}];
  7. Beff[t_] := {0, 0, 1.5 - 0.8*(m[t].{0, 0, 1})};
  8. cons[t_] := -gamma*Cross[m[t], Beff[t]];
  9. tGilbdamp[t_] := alphag*Cross[m[t], cons[t]];
  10. tSlondamp[t_] :=
  11. current[t]*alphajConstant*gamma*Cross[m[t], Cross[m[t], p]];
  12. equ= {m'[t] == cons[t] + tGilbdamp[t] + tSlondamp[t], m[0] == {0, 0, 1}};
  13. sol1 = NDSolve[equ, {m}, {t, 0, 200}, MaxSteps -> [Infinity]];
  14. mm[t_] = m[t] /. sol1[[1]];
  15. Print[Plot[mm[t][[1]], {t, 0, 200}]];];
  16.  
  17. LaunchKernels[2]
  18. ParallelEvaluate[func[50], Kernels[][[1]]]
  19. ParallelEvaluate[func[100], Kernels[][[2]]]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement