Discussion :

[vasilpapa]

Bonjour,
J'ai un petit code qui inclu les fonctions pour nombres complexes. Pour Mandelbrot et Julia aucun problème.
(Voir les possibilités et les résultats sur http://haskellaskvasil.wordpress.com).
Par contre, dès que je modifie la dernière ligne d'itération pour dessiner Burning Ship, j'ai le message :

BurningShip.hs:77:26:
Couldn't match expected type `Complex' with actual type `Float'
In the return type of a call of `real'
In the first argument of `abs', namely `(real (z))'
In the first argument of `(+)', namely `abs (real (z))'

BurningShip.hs:77:41:
Couldn't match expected type `Complex' with actual type `Float'
In the return type of a call of `im'
In the first argument of `abs', namely `(im (z))'
In the second argument of `(+)', namely `abs (im (z))'


Voici le code en question :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
import Graphics.Rendering.OpenGL
import Graphics.UI.GLUT
import Data.IORef
 
data Complex = C (Float, Float) deriving (Show,Eq)
 
instance Num Complex where
    fromInteger n = C (fromIntegral n,0.0)
    C (x,y) * C (z,t) = C (z*x - y*t, y*z + x*t)
    C (x,y) + C (z,t) = C (x+z, y+t)
    abs (C (x,y))     = C (sqrt (x*x + y*y),0.0)
    signum (C (x,y))  = C (signum x , 0.0)
 
complex :: Float -> Float -> Complex
complex x y = C (x,y)
 
real :: Complex -> Float
real (C (x,y))    = x
 
im :: Complex -> Float
im   (C (x,y))    = y
 
magnitude :: Complex -> Float
magnitude = real.abs
 
main :: IO ()
main = do
  (progname,_) <- getArgsAndInitialize
  initialDisplayMode $= [DoubleBuffered]
  initialWindowSize $= Size 320 320
  createWindow "Mandelbrot Set"
  displayCallback $= display
  putStr "PATIENTEZ - La fenetre avec l'image sera affichee a la fin du traitement"
  display
  mainLoop
 
display = do
  clear [ColorBuffer] -- make the window black
  loadIdentity -- reset any transformation
  preservingMatrix drawMandelbrot
  swapBuffers -- refresh screen
 
 
drawMandelbrot =
  renderPrimitive Points $ do
    mapM_ drawColoredPoint allPoints
  where
      drawColoredPoint (x,y,c) = do
          color c -- set the current color to c
          vertex $ Vertex3 x y 0 
 
width = 320 :: GLfloat
height = 320 :: GLfloat
 
allPoints :: [(GLfloat,GLfloat,Color3 GLfloat)]  --  $ julia x y )
 
allPoints = [ ((x-100)/width,y/height,colorFromValue $ mandel x y) | 
                  x <- [-width+100..width+100], 
                  y <- [-height..height]]
 
colorFromValue n =
  let 
      t :: Int -> GLfloat
      t i = 0.5 + 0.5*sin( fromIntegral i / 1 )
  in
    Color3 (t n) (t (n+5)) (t (n+10))
 
mandel x y = 
  let r = realToFrac $ 0.5 * x / width 
      i = realToFrac $ 0.5 * y / height
  in
      f (complex r i) 0 64 
 
{-
julia x y =
      let r = realToFrac $ 2 * x / width 
          i = realToFrac $ 2 * y / height
	  s = -0.7
	  p = 0.2715
      in
          f (complex s p) (complex r i) 128-}
 
f :: Complex -> Complex -> Int -> Int
f c z 0 = 0
f c z n = if (magnitude z > 2 ) 
          then n
          else f c ((z*z) + c) (n-1)    -- Burning Ship ((abs(real(z)) + abs(im(z)))^2 + c) (n-1)

J'ai tout essayé, je n'y arrive pas. Pourrai-on svp me mettre sur la bonne voie ?

Si j'utilise module Complex (import Data.Complex) avec realPart et imagPart, j'ai un autre message :

BurningShipComplex.hs:57:26:
Expecting one more argument to `Complex'
In the type signature for `f':
f :: Double -> Double -> Complex -> Complex -> Int -> Int

Merci d'avance.

[stendhal666]

C'est curieux, chez moi le code compile tout à fait. En revanche j'ai lancé main à tout hasard, une fenêtre vide s'est affichée, 10, 20, 30 secondes ou plus et j'ai décidé d'interrompre les calculs...

[vasilpapa]

Vous avez lancé avec la ligne modifiée pour burining slip?
Seul ce cas ne marche pas.
Dans winCGHi ?
Il faut patienter un peu, comme tout les mandelbrot c'est un peu long.
Avec la ligne display avant mainLoop la fenêtre doit s'afficher a la fin du traitement.
Avons nous des versions différentes du CGHi ? J'ai installé la plateforme début mai 1914.

[stendhal666]

Avons nous des versions différentes du CGHi ? J'ai installé la plateforme début mai 1914.

mai 1914 alors oui, nous avons probablement des versions différentes... J'essaie de mettre ma plate-forme à jour régulièrement.

Je n'avais pas modifié la ligne pour Burning Ship, je crois voir maintenant de quoi il s'agit:
la signature de la fonction f est Complex -> Complex -> Int -> Int

Or : ((abs(real(z)) + abs(im(z)))^2 + c) qui correspond au deuxième argument, est un Float, pas un Complex. D'où la protestation du compilateur. J'imagine qu'en modifiant la signature vous devriez pouvoir compiler le code.

[vasilpapa]

J'ai essayé de modifier la signature, mais le compilateur protesté, c'est pire.
Si vous trouvez la bonne signature ca serait formidable.

[Jedai]

Il faut patienter un peu, comme tout les mandelbrot c'est un peu long.

Mon code final n'est pas particulièrement lent : il recalcule l'image plusieurs fois par seconde et j'obtiens un rendu fluide même sur une machine qui a bien 6,7 ans (testé pendant mes vacances).


Pour ton code, le problème ici c'est que "c" est un Complex. Comme tu essaies de l'ajouter à "(abs(real(z)) + abs(im(z)))^2", GHC s'attend ("expect") à ce que cette expression soit aussi un Complex. Mais il s'agit d'un Float car real et im renvoie des Floats et les opérations suivantes (abs, ^2) ne modifient pas le type. Je suppose que tu pourrais convertir de nouveau cette expression en Complex (la version de ce fichier) :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

let r = abs (real z) + abs (im z) in f c (complex (r^2) 0 + c) (n-1)

Il est également possible que tu ne souhaites pas que "c" soit un complexe, il faudra alors changer la signature et l'invocation de f.

--
Jedaï

[vasilpapa]

Ca ne marche pas; (complex (r^2)) a un seul argument alors qu'il en faut deux.
Mais vous m'avez mis sur la voie, voici la solution :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
f :: Complex -> Complex -> Int -> Int
f c z 0 = 0
f c z n = if (magnitude z > 2 ) 
          then n
          else let a = abs (real (z))
                   b = abs(im (z)) 
			in f c ( (complex a b)^2 + c) (n-1)

Merci.
Avant de fignoler la première image :

N'oubliez pas de changer y en -y!

[Jedai]

Ca ne marche pas; (complex (r^2)) a un seul argument alors qu'il en faut deux.

C'est pour ça qu'il y avait un 0 en second argument... Dans ton code d'origine tu mettais "|a|+|b|" au carré, pas "|a|+|b|i" mais il s'agissait je suppose d'une erreur de ta part si tu obtiens la bonne image avec ta nouvelle version.

[vasilpapa]

Oui, c'était une erreur, maintenant tout vas bien, merci. Même les zooms :