ix: código parámetros fractrace

Enviado por mfg el Mar, 07/11/2017 - 10:50

El proceso iterativo en el programa 'fractrace' se muestra a continuación.

Fue necesario añadir un número pequeño pero finito en los denominadores para evitar divergencias en los cálculos numéricos. Ésta no es una solución satisfactoria pues introduce una pequeña deformación que no debe estar allí. Recordemos que el producto de escatores es cerrado en el conjunto $\mathbb{E}_{-}^{1+n}$, de manera que se debe poder resolver este problema numérico de una forma mas satisfactoria.

/* ix - fractrace - tomotrace */
    aar = ar * ar; aai = ai * ai; aaj = aj * aj;
    double r2 = aar + aai + aaj + aai * aaj /(aar+.001);
  for (n = 1; n < zkl; n++) {   
    r2 = aar + aai + aaj + aai * aaj /(aar+.001);
    tar =  aar - aai - aaj + aai * aaj /(aar+.001) +br;
    tai = 2 * ar * ai - 2 * ai * aaj / (ar+.001) +bi;
    taj = 2 * ar * aj - 2 * aj * aai / (ar+.001) +bj;
    ar=tar; ai=tai; aj=taj;
    aar = ar * ar; aai = ai * ai; aaj = aj * aj;
    r2 = aar + aai + aaj + aai * aaj /(aar+.001);
    if (r2 > 400) { tw = n; break; }

    }
 
Ejemplo de archivo .xml para bajar y procesar en 'fractrace': fractrace-formula-example.zip