[c] mi servono numeri con piu cifre...

[mamo139]

purtroppo costanti e variabili di questo programma (ad esempio const double Y1 = 0.6501800130615235 stanno iniziando ad essere troppo poco precise e il programma funziona peggio...
devo usare numeri con piu cifre, ma da come potete vedere devo usare e poter lavorare senza problemi con i numeri complessi...
come posso fare??

grazie mille per ogni aiuto!!!

#include "ppm_image.h" 
#include <complex> 
#include <math.h> 
#include <windows.h> 

using namespace std; 


#define WIDTH  2000
#define HEIGHT 2000

#define DEPTH 4000
               
                 const double Y1 = 0.6501800130615235;
const double X0 = -0.0984813251342773;
                                 const double X1 = -0.0984810866699219;
                 const double Y0 = 0.6501797746582032;



unsigned int * hsv_rgb(unsigned int *);

int main() 
   { 

   complex<double> value, point; 

   double r, i, modulo; 
   long cnt; 
   
   //dati utili per colorare
   long cnt_tot=0, media_cnt=0, media_cnt_n=0, c_cnt=0;

    
   // Creo un'immagine vuota di dimensione WIDTH*HEIGHT 
   PpmImage img(WIDTH, HEIGHT); 

   for(int x=0; x<WIDTH; x++){ 
      r = X0+((double(x)*(X1-X0))/ double(WIDTH));       
      for(int y=0; y<HEIGHT; y++){ 
         i = Y0+((double(y)*(Y1-Y0))/double(HEIGHT)); 
         
         point = complex<double>(r,i); 
         value = point; 
          
         // iniziamo il ciclo 
         
         for(cnt=0; cnt<DEPTH; cnt++){ 
            value = value*value+point; 
            modulo = abs(value); 
            if(modulo>2)break; 
         } 
         

         if(modulo<=2){ 
            img.SetPixel(x,(HEIGHT-y),0,0,0); 
         } 
         else{ 
            // invece di disegnare un pixel bianco ne disegno uno colorato che dipende da cnt 
            unsigned int red, green, blue; //colore, contrasto, luminosità
               
            unsigned int *col;
            col = new unsigned int[3]; 
            unsigned int *col2;

            
            //cnt_tot=cnt_tot+cnt;
            //media_cnt_n++;
            c_cnt = cnt/4;
            
            red   = (c_cnt)%360; 
            green = 250; 
            blue  = 250;   

            if(red>359)   red=359; 
            if(green>255) green=255; 
            if(blue>255)  blue=255; 
            
            col[0]=red;
            col[1]=green;
            col[2]=blue;

            col2 = hsv_rgb(col);
            delete []col;            

            img.SetPixel(x, (HEIGHT-y), col2[0], col2[1], col2[2]); 
            delete []col2;
         } 
      } 
      printf("%d - %d%%\n",x,x*100/WIDTH);
   } 

   // Salvo l'immagine 
   img.SaveToFile("mandelbrot.ppm"); 
   
   //media_cnt=cnt_tot/media_cnt_n;
   //printf("cnt medio: %d\n",media_cnt);
   
   printf("zoom consigliati:\n\n");
   printf("(Y)\n");
   printf("1: %.16f\n",(double)Y1);
   printf("2: %.16f\n",(double)(((Y1+Y0)/2)+Y1)/2 );   
   printf("3: %.16f\n",(double)(Y1+Y0)/2);
   printf("4: %.16f\n",(double)(((Y1+Y0)/2)+Y0)/2 );      
   printf("5: %.16f\n",(double)Y0);   
   
   printf("\n(X)\n");
   printf("1: %.16f\n",(double)X0);
   printf("2: %.16f\n",(double)(((X0+X1)/2)+X0)/2 );   
   printf("3: %.16f\n",(double)(X0+X1)/2);
   printf("4: %.16f\n",(double)(((X0+X1)/2)+X1)/2 );      
   printf("5: %.16f\n",(double)X1);   
   
   system("pause");
   return 0;
}









unsigned int * hsv_rgb(unsigned int hsv[3]){
         
    //printf("REVERSE:\nH:%d\nS:%d\nV:%d\n",hsv[0],hsv[1],hsv[2]);          
    double *dRGB;
    dRGB = new double[3];    
    
    double H = ((double)hsv[0]);
    double S = ((double)hsv[1])/255;
    double V = ((double)hsv[2])/255;
    
    double hi = ((int)(H/60))%6;
    double f = (H/60) - (double)hi;
    double p = V*(1-S);
    double q = V*(1-(f*S));
    double t = V*(1-(1-f)*S);
    
    if (hi == 0){ dRGB[0]=V;dRGB[1]=t;dRGB[2]=p;}
    else if (hi == 1){ dRGB[0]=q;dRGB[1]=V;dRGB[2]=p;}
    else if (hi == 2){ dRGB[0]=p;dRGB[1]=V;dRGB[2]=t;}
    else if (hi == 3){ dRGB[0]=p;dRGB[1]=q;dRGB[2]=V;}
    else if (hi == 4){ dRGB[0]=t;dRGB[1]=p;dRGB[2]=V;}
    else if (hi == 5){ dRGB[0]=V;dRGB[1]=p;dRGB[2]=q;}
    else printf("ERRORE!");

    //printf("R:%f\nG:%f\nB:%f\n",dRGB[0],dRGB[1],dRGB[2]);
    
    unsigned int *RGB;
    RGB = new unsigned int[3];    
    
    RGB[0] = (int)(dRGB[0]*255);
    RGB[1] = (int)(dRGB[1]*255);
    RGB[2] = (int)(dRGB[2]*255);
    
    delete[]dRGB; 
    
    //printf("R:%d\nG:%d\nB:%d\n",RGB[0],RGB[1],RGB[2]);  
    return RGB;   
}

[cionci]

Se ti servono più cifre dovresti usare una qualche libreria che ti permetta di utilizzare numeri con precisione arbitraria. Per il C librerie di questo tipo sono un po' scomode da usare.