En we are back in business!!

Via een tip van mijn assistent, heb ik gevonden hoe ik mijn linear system of equations kan opstellen. Deze was namelijk te vinden in het PHD thesis van Kwatra! Na een week van studie op die thesis, kon ik dus deze week beginnen aan de implementatie. JOY!!

Eerst moest ik een matrix library vinden. Na heel wat zoeken, ben ik op een fantastische vergelijkingspagina terecht gekomen. Het is een echte aanrader! Eerst wat zitten knoeien met MTL en ITL, en hoewel deze erg goed leken, wou deze zelfs na verscheidene malen proberen, niet compileren in VS.NET!

Uiteindelijk ben ik dus terecht gekomen bij de newmat library! Alles ging vlot tot ik een matrix vermenigvuldiging probeerde! Deze duurde voor 256*4096 (en dat was dan een kleine afbeelding) anderhalve minuut!!! En als je dan weet dat deze 256 keer of meer in de eerste iteratie van het algoritme zou voorkomen, dan kan je al 6 en half uur wachten voor een iteratie!! Net iets te lang dus! Na bespreking met de auteur, mijn begeleider en enkele vrienden, heb ik de methode wat aangepast, zonder matrix vermingvuldigingen en nu duurt een iteratie een fractie van een seconde!!! Als output krijg ik wel nog rubbish:

Bron                           Output
  

Maar ik zit in alle geval terug in de goede richting!! HALLELUJA!

Nog enkele dankwoordjes:

• Dank aan Björn-Ove Heimsund voor de vergelijkingspagina tip!
• Dank aan de MTL en ITL mailing list!
• Dank aan Robert Davies, de maker van newmat, voor de tips en het snel antwoorden op mijn emails!
• Dank aan Nele en Chris voor het testen van matrices in Java en Fortran! :-p
• En zoals altijd thx Ares voor de tips!

This entry was posted on Thursday, March 23rd, 2006 at 17:14 and is filed under Algemeen, C++ / OpenGL / Cg, Dutch. You can follow any responses to this entry through the RSS 2.0 feed. You can leave a response, or trackback from your own site.