JSL pour la vidéo.
picpath = ""; // "D:\Chales\pic\m";
//http://www.fractalforums.com/mandelbrot-and-julia-set/good-coordinates-for-mandelbrot-zooms/
//
//Points on the boundary of the Mandelbrot set generally have the most interesting orbits. The easiest boundary points to compute are:
//
//* the spike along the negative real axis
//* the boundary of the main cardioid: r = (1 - cos(theta))/2, x = r*cos(theta)+0.25, y = r*sin(theta)
//* the boundary of the period 2 disk: r = 0.25, x = r*cos(theta)-1, y = r*sin(theta)
eps = 1.6; // slightly outside
deep = 10;
steps = 8000;
angle = 3;
scale = 2; // <<<<<<< video made with 20, very slow
xsize = 192 * scale;
ysize = 108 * scale;
diag = Sqrt( xsize ^ 2 + ysize ^ 2 );
mat = J( ysize, xsize, . );
New Window( "mandelbrot", bb = Border Box( Text Box( "wait" ) ) );
Wait( .2 );
num = 1e6;
For( t = 0, t < steps, t += 1,
If( eps > 1e-9, eps *= .997 );
deep += .3;
angle += eps / 50;
radius = (1 - Cos( angle )) / 2 + eps;
x = radius * Cos( angle ) + 0.25;
y = radius * Sin( angle );
top = y - eps;
bottom = y + eps;
left = x - eps * xsize / ysize;
right = x + eps * xsize / ysize;
d = Floor( Min( 1e5, deep ) );
Parallel Assign( {d = d}, mat[r( top, bottom ), c( left, right )] = Mandelbrot( d, 4, c, r ) );
mat = mat / d;
red = mat ^ .5;
green = mat ^ 9;
both = Loc( red :*</img> green == 1 );
red[both] = .1;
green[both] = .4;
blue = 1 - mat ^ .2;
blue[both] = .1;
img = New Image( "rgb", {red, green, blue} );
(bb << child) << delete;
bb << append( img );
num += 1;
If( picpath != "", img << saveimage( picpath || Char( num ) || ".png" ) );
bb << updatewindow;
Wait( 0 );
);
Il semble que je n'ai jamais documenté cela, je vais voir ce que je peux faire.
Oui, il effectue une interpolation linéaire de 1..nrows en haut..en bas. (similaire pour c, ncols, gauche, droite).
Si vous omettez le bit (haut, bas), la valeur r dans l'expression sera 1..nrows.
La partie {d = d} doit également être expliquée. Chaque thread obtient sa propre copie de tout ce qu'il doit référencer ; la variable du côté gauche est le nom du thread de travail et le côté droit est le nom du thread principal. Sur une machine à 4 processeurs, cette affectation est effectuée quatre fois. S'il s'agit d'une matrice énorme dans laquelle vous copiez, vous ne constaterez peut-être pas d'amélioration des performances.
Mandelbrot n’a besoin d’aucune autre donnée et évolue bien. La plupart des problèmes réels doivent être lus à partir d'un autre tableau et peut-être également faire une comptabilité.
Si vous voulez jeter un œil sous les couvertures, essayez ceci :
m=J(10,10,0);
x=0;
parallelassign({x=x++},m[r,c]=x);
show(x,m);
Ma machine à 8 processeurs produit ceci :
x = 8 ;
m =
[ 0 1 2 3 4 5 6 7 0 1,
2 3 4 5 6 7 0 1 2 3,
4 5 6 7 0 1 2 3 4 5,
6 7 0 1 2 3 4 5 6 7,
0 1 2 3 4 5 6 7 0 1,
2 3 4 5 6 7 0 1 2 3,
4 5 6 7 0 1 2 3 4 5,
6 7 0 1 2 3 4 5 6 7,
0 1 2 3 4 5 6 7 0 1,
2 3 4 5 6 7 0 1 2 3];
Bien entendu, vous ne devez pas écrire de code qui dépend de ce modèle.
Notez également que le côté droit de l’affectation de la matrice peut être plus compliqué.
m=J(400,600,0);
parallelassign({},m[r(-1,1),c(-1,1)]=(
z=sqrt(r^2+c^2);
z=if(z<.707,.707,z)-.5;
1-z; // last item in ; sequence is result
));
m=m/max(m); // for heat color
newwindow("example",newimage(heatcolor(m,"spectral")));
"z" n'existe qu'à l'intérieur de chacun des threads, une copie distincte pour chacun.
That flat area, dark red, was created by the if statementThat flat area, dark red, was created by the if statementThat flat area, dark red, was created by the if statementThat flat area, dark red, was created by the if statementCette zone plate, rouge foncé, a été créée par l'instruction if
Et oui, vous pouvez utiliser des déclarations. Vous ne devez essayer d'effectuer aucune forme d'E/S, sauf écrire dans le journal. Ne créez pas de boîtes d'affichage, ne mettez pas à jour les boîtes d'affichage, ne lisez pas et n'écrivez pas de fichiers disque. S'il ne plante pas, ce sera toujours un problème.
Vous pouvez utiliser des globaux, mais ils utilisent un mutex pour sérialiser l'accès et le code ne peut pas être plus rapide qu'un code non threadé.
m=J(3,5,0);
global:words="";
parallelassign({},m[r(-1,1),c(-1,1)]=(
z=sqrt(r^2+c^2);
show(r+0,c+0,r+c,z); // r and c are special, use an expression
global:words||=char(z);
z=if(z<.707,.707,z)-.5;
1-z; // last item in ; sequence is result
));
show(words);
global:words est construit dans un ordre imprévisible, déterminé par la vitesse d'exécution des différents threads.
JSL est compliqué (en interne), tout comme le threading. La fonction parallelassign essaie de garder tout sain (ces copies, par exemple). Si vous trouvez un moyen de contourner quelque chose qu'il essaie de vous empêcher de faire, vous risquez d'endommager les structures de données internes lorsque deux threads y accèdent en même temps, ce qui entraînera le blocage de JMP à un moment donné. Économisez souvent.
Oh... une dernière remarque sur les matrices 1D : utilisez deux indices ! Il existe des situations déroutantes si vous avez un vecteur de ligne 1xN et essayez d'utiliser un seul indice. Utilisez R,C et ignorez la valeur constante R==1.
Je me suis amusé avec ce script, en apportant quelques modifications et en créant la série d'images à convertir en vidéo pour réaliser mes propres visites de Mandelbrot.
Question : La fonction Parallel Assign() -- Je recherche une référence de syntaxe dessus ?Je ne le vois pas dans les livres de script JMP, mais il y a une entrée dans le Guide de script JMP sous Aide.Je ne comprends pas ce que vous faites avec la fonction Parallel Assign dans votre script.
Parallel Assign( {d = d}, mat[r( top, bottom ), c( left, right )] = Mandelbrot( d, 4, c, r ) );que font les r(haut, bas) et c(gauche, droite) pour vous ?le haut, le bas, la gauche et la droite sont déterminés juste avant dans le script pour déterminer les limites des coordonnées du cadre que vous créez.Je suppose qu'il effectue une sorte d'interpolation entre le haut, le bas et la gauche, la droite pour obtenir les coordonnées y et x à utiliser ?Existe-t-il une description générale de cette syntaxe quelque part ?
