FRACTAL DE MANDELBROT

Commentaire de jesusonline le 11/02/2007 01:21:03 administrateur CS

Voici une version non managé :

    ''' ------------------------------------------------------------------------
    ''' <summary>
    ''' Draw the fractal.
    ''' </summary>
    ''' <remarks> For more performance, use unmanaged code [unsafe]. </remarks>
    ''' ------------------------------------------------------------------------
    Public Sub DraMandelbrotFractal()
        Dim gfx As Graphics = Me.CreateGraphics()
        Me._bmpBuffer = New Bitmap(Me.Width, Me.Height, gfx)

        Dim halfX As Integer = 2 * Me._bmpBuffer.Width / 3
        Dim halfY As Integer = Me._bmpBuffer.Height / 2
        Dim i As Integer, j As Integer = 0

        Dim bmpData As Imaging.BitmapData = Me._bmpBuffer.LockBits(New Rectangle(0, 0, width, Height), System.Drawing.Imaging.ImageLockMode.ReadWrite, System.Drawing.Imaging.PixelFormat.Format32bppArgb)
        Dim newPixel(width * height - 1) As Integer
        Marshal.Copy(bmpData.Scan0, newPixel, 0, newPixel.Length)
        For j = 0 To Me.Height - 1
            For i = 0 To Me.Width - 1
                Dim complex As Complex = New Complex((i - halfX) / 150D, (j - halfY) / 150D)

                Dim temp As Integer = newPixel(Me.Width * j + i)
                newPixel(Me.Width * j + i) = Me.GetColor(Me.GetStep(complex)).ToArgb()

            Next
        Next

        Marshal.Copy(newPixel, 0, bmpData.Scan0, newPixel.Length)
        Me._bmpBuffer.UnlockBits(bmpData)
        Graphics.FromImage(Me._bmpBuffer).DrawString("Mandelbrot's Fractal", Me._font, Me._fontBrush, Point.Empty)
        gfx.DrawImage(Me._bmpBuffer, Point.Empty)
    End Sub


Pour en savoir plus il faut lire http://www.techheadbrothers.com/Articles.aspx?Id=8d3eb481-8a98-42a6-8033-e851c797aa60&p=3

Commentaire de Bidou le 11/02/2007 09:46:39 administrateur CS

Cyril> Oui, en fait je pensais à un bloc unsafe (que j'ai écrit et qui marche) mais les couleurs ne joue plus tout à fait. J'ai pas encore regardé en détail pourquoi... mais je vais certainement trouvé la réponse dans ce très bon article dont tu as donné le lien.

Badr> Aucun risque, après utilisation et comparaison des deux, je me sens nettement plus à l'aise avec C# (et il n'a jamais été question d'arrêter C#) :-p

Commentaire de Bidou le 11/02/2007 17:49:08 administrateur CS

C'est pas de "l'unmanaged", mais de "l'unmanaged .NET" là est certainement la différence...
http://vivekthangaswamy.blogspot.com/2006/11/unsafe-code-in-c.html

Le code de JesusOnline n'est pas directement unmanaged .NET à mon sens car il ne fait pas intervenir de bloc unsafe.
Voici le code que j'ai fait en C# que j'appelerais unmanaged (.NET) :

unsafe
{
  byte* newPixel = (byte*)(void*)bmpData.Scan0;
  for (int y = 0; y < height; ++y)
  {
     for (int x = 0; x < width; x += 3)
     {
        Color color = XXX;
        newPixel[0] = color.B;
        newPixel[1] = color.G;
        newPixel[2] = color.R;
        newPixel += 3;
     }
     newPixel += offset;
   }
}

Commentaire de jesusonline le 11/02/2007 17:58:47 administrateur CS

Marshal copy et lock bits sont des méthodes unmanaged (au sens .net) par conséquent mon code est unmanaged :-)

A ma connaissane on ne peut pas utiliser les pointeurs en VB ... :-( donc on peut pas faire la meme chose mais je suis curieux de voir la différence de perf entre C# et VB.net et la méthode MSIL que TKfé a posté sur son blog ...

Commentaire de BruNews le 11/02/2007 19:26:46 administrateur CS

oh la misère, faut avoir tué père et mère pour se taper cette syntaxe.
C'est tellement plus simple à écrire en vrai C avec en plus de bien meilleures perfs et aucun setup, j'espère que le blaireau qui a incité MS à prendre le terme "C++" en dotnetterie aura depuis été affecté à tondre la pelouse.
Merci pour les infos.

Commentaire de BruNews le 11/02/2007 19:41:46 administrateur CS

"assembleur" a déjà une signification pour le commun dev dont je fais partie:
- mise au format humain des opcodes directement exécutables par le cpu.

Toutes ces usurpations de langage ne devraient pas berner grand monde.

"l'assembleur des langages interprété", voila bien une suite de mots dont chaque suivant contredit le précédent.

Commentaire de BruNews le 11/02/2007 19:50:21 administrateur CS

ok Bidou, j'avais bien noté avec ta réponse plus haut au vu de data1->Scan0.

Commentaire de EBArtSoft le 20/05/2008 22:40:33 administrateur CS

Voila ce que je propose :

    Private Structure ARGB
        Dim r As Byte
        Dim g As Byte
        Dim b As Byte
        Dim a As Byte
    End Structure

    Private Declare Ansi Sub Copy Lib "kernel32" Alias "RtlMoveMemory" (ByVal memTo As Integer, ByRef memFrom As ARGB, ByVal Length As Integer)

    Public Sub DraMandelbrotFractal()

        'Me.Width = ((Me.Width + 3) And Not 3) 'Seulement pour 24Bits RGB

        Static _Width As Integer = Me.Width
        Static _Height As Integer = Me.Height
        Static halfX As Integer = _Width * 2 / 3
        Static halfY As Integer = _Height * 0.5
        Static siz As Integer = _Width * _Height * 4
        Static z1 As Single = 1 / Me._depth * 64
        Static z2 As Single = 1 / Me._depth * 255
        Static StepX As Integer = 0

        Dim Pixels(_Height - 1, _Width - 1) As ARGB
        For j As Integer = 0 To (_Height - 1)
            For i As Integer = 0 To (_Width - 1)
                stepx = Me.GetStep(New Complex((i - halfX) / 150D, (j - halfY) / 150D))
                If (stepx > 0) Then
                    With Pixels(j, i)
                        .a = 255
                        .r = z1 * stepx
                        .b = z2 * stepx
                    End With
                End If
            Next
        Next

        Me._bmpBuffer = New Bitmap(_Width, _Height, System.Drawing.Imaging.PixelFormat.Format32bppArgb)
        Dim bmpData As Imaging.BitmapData = Me._bmpBuffer.LockBits(New Rectangle(0, 0, _Width, _Height), System.Drawing.Imaging.ImageLockMode.ReadWrite, System.Drawing.Imaging.PixelFormat.Format32bppArgb)
        'Marshal.Copy(NewPixel, 0, bmpData.Scan0, NewPixel.Length)
        Copy(bmpData.Scan0, Pixels(0, 0), siz)
        Me._bmpBuffer.UnlockBits(bmpData)

        Graphics.FromImage(Me._bmpBuffer).DrawString("Mandelbrot's Fractal", Me._font, Me._fontBrush, Point.Empty)
        Me.CreateGraphics().DrawImage(Me._bmpBuffer, Point.Empty)

    End Sub

Apres difficile de faire mieu sans changer le comportement des autres methodes dont celles qui prennent le plus de temps : Complexe et GetStep

Et puis on pourrais creer le bitmap au moment de l'initialisation et pas a chaque fois, bref l'optimisation peut ce faire un peu partout.

@+