<HTML><BODY style="word-wrap: break-word; -khtml-nbsp-mode: space; -khtml-line-break: after-white-space; ">Hi again,<DIV><BR class="khtml-block-placeholder"></DIV><DIV>I believe, that I might have found the issue myself possibly. When I use 10% of all voxels for sampling such as in line 335 in DeformableRegistration8.cxx of release 3.2.0, which is</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" face="Monaco"> </FONT><FONT class="Apple-style-span" color="#760F50" face="Monaco">const</FONT><FONT class="Apple-style-span" face="Monaco"> </FONT><FONT class="Apple-style-span" color="#760F50" face="Monaco">unsigned</FONT><FONT class="Apple-style-span" face="Monaco"> </FONT><FONT class="Apple-style-span" color="#760F50" face="Monaco">int</FONT><FONT class="Apple-style-span" face="Monaco"> numberOfSamples = fixedRegion.GetNumberOfPixels() / </FONT><FONT class="Apple-style-span" color="#0000FF" face="Monaco">10</FONT><FONT class="Apple-style-span" face="Monaco">;</FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" face="Monaco"><BR class="khtml-block-placeholder"></FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" face="Monaco">then I get the bad_alloc as described. Changing the default example to only use 5% of the voxels i.e. </FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" face="Monaco"><BR class="khtml-block-placeholder"></FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" color="#760F50" face="Monaco">const</FONT><FONT class="Apple-style-span" face="Monaco"> </FONT><FONT class="Apple-style-span" color="#760F50" face="Monaco">unsigned</FONT><FONT class="Apple-style-span" face="Monaco"> </FONT><FONT class="Apple-style-span" color="#760F50" face="Monaco">int</FONT><FONT class="Apple-style-span" face="Monaco"> numberOfSamples = fixedRegion.GetNumberOfPixels() / <FONT class="Apple-style-span" color="#0000FF">20</FONT></FONT><FONT class="Apple-style-span" face="Monaco">;</FONT><FONT class="Apple-style-span" face="Monaco"></FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" face="Monaco"><BR class="khtml-block-placeholder"></FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" face="Monaco">solved the memory bug, my app is up and running now using about 3 GBytes of memory. Usually I would blame it to 32bit applications now, but mine (and I also checked part of the ITK libraries for it) are 64 bit (as is the system, where I run it). Thus I would assume now, that maybe either the const unsigned int above or some other int pointer would be guilty of throwing a bad_alloc, because the pointer is not able managing addesse space requiring more than 32 bits. Do you agree to this? To what degree is ITK capable (I mean of course capable from scratch, without modifications to the code!) of handling more than 32bit memory? Would I need to change some pointers of ITK to handle some large volumes? My images have up to 90,000,000 voxels.</FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" face="Monaco"><BR class="khtml-block-placeholder"></FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" face="Monaco">Best</FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" face="Monaco">Ruben</FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" face="Monaco"><BR class="khtml-block-placeholder"></FONT></DIV><DIV><BR><DIV><BR><DIV>Anfang der weitergeleiteten E-Mail:</DIV><BR class="Apple-interchange-newline"><BLOCKQUOTE type="cite"><P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">Hi Luis,</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px; font: 12.0px Helvetica; min-height: 14.0px"><BR></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">I just came back to the registration issue. I ran DeformableRegistration8 from gdb by typing</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px; font: 12.0px Helvetica; min-height: 14.0px"><BR></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">gdb DeformableRegistration8</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px; font: 12.0px Helvetica; min-height: 14.0px"><BR></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">afterwards I ran the program from the gdb prompt with my inputfiles like this</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px; font: 12.0px Helvetica; min-height: 14.0px"><BR></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">>> run fixedImage.tif movingImage.tif outputImage.tif</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px; font: 12.0px Helvetica; min-height: 14.0px"><BR></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">I don't know, if there is a way to get more details from gdb about the running program, but so far I got this from gdb:</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px; font: 12.0px Helvetica; min-height: 14.0px"><BR></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">[Thread debugging using libthread_db enabled]</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">[New Thread 47202034567056 (LWP 24727)]</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px; font: 12.0px Helvetica; min-height: 14.0px"><BR></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">Starting Registration</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">terminate called after throwing an instance of 'std::bad_alloc'</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica"><SPAN class="Apple-converted-space"> </SPAN>what():<SPAN class="Apple-converted-space"> </SPAN>St9bad_alloc</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px; font: 12.0px Helvetica; min-height: 14.0px"><BR></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">Program received signal SIGABRT, Aborted.</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">[Switching to Thread 47202034567056 (LWP 24727)]</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">0x00002aee143e55ed in *__GI_raise (sig=24727) at ../nptl/sysdeps/unix/sysv/linux/raise.c:67</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica">67<SPAN class="Apple-converted-space"> </SPAN>../nptl/sysdeps/unix/sysv/linux/raise.c: No such file or directory.</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px"><FONT face="Helvetica" size="3" style="font: 12.0px Helvetica"><SPAN class="Apple-converted-space"> </SPAN>in ../nptl/sysdeps/unix/sysv/linux/raise.c</FONT></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px; font: 12.0px Helvetica; min-height: 14.0px"><BR></P> <P style="margin: 0.0px 0.0px 0.0px 0.0px; font: 12.0px Helvetica; min-height: 14.0px"><BR></P> </BLOCKQUOTE></DIV><BR></DIV></BODY></HTML>