<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:"Malgun Gothic";
panose-1:2 11 5 3 2 0 0 2 0 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@Malgun Gothic";
panose-1:2 11 5 3 2 0 0 2 0 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p
{mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
{mso-style-priority:34;
margin-top:0in;
margin-right:0in;
margin-bottom:0in;
margin-left:.5in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
/* List Definitions */
@list l0
{mso-list-id:216938851;
mso-list-type:hybrid;
mso-list-template-ids:-1572176310 67698705 67698713 67698715 67698703 67698713 67698715 67698703 67698713 67698715;}
@list l0:level1
{mso-level-text:"%1\)";
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-.25in;}
@list l0:level2
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-.25in;}
@list l0:level3
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
@list l0:level4
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-.25in;}
@list l0:level5
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-.25in;}
@list l0:level6
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
@list l0:level7
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-.25in;}
@list l0:level8
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-.25in;}
@list l0:level9
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
ol
{margin-bottom:0in;}
ul
{margin-bottom:0in;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>Hi Simon,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>That sounds great! Thanks a lot for your hard work.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>As you said, now I can reconstruct it successfully with --subsetsize option.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>I think the only concern left is the elongated reconstruction time.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>I’ve measured it for the same phantom image with different “subsetsize” values as well as w/o FFTW compiling.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>I just report it for your information.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><o:p> </o:p></span></p><p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo1'><![if !supportLists]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><span style='mso-list:Ignore'>1)<span style='font:7.0pt "Times New Roman"'> </span></span></span><![endif]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>w/o FFTW : 21.4 s<o:p></o:p></span></p><p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo1'><![if !supportLists]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><span style='mso-list:Ignore'>2)<span style='font:7.0pt "Times New Roman"'> </span></span></span><![endif]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>w/ FFTW & --subsetsize 1: 89 s<o:p></o:p></span></p><p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo1'><![if !supportLists]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><span style='mso-list:Ignore'>3)<span style='font:7.0pt "Times New Roman"'> </span></span></span><![endif]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>w/ FFTW & --subsetsize 2: 77 s<o:p></o:p></span></p><p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo1'><![if !supportLists]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><span style='mso-list:Ignore'>4)<span style='font:7.0pt "Times New Roman"'> </span></span></span><![endif]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>w/ FFTW & --subsetsize 4: 65 s<o:p></o:p></span></p><p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo1'><![if !supportLists]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><span style='mso-list:Ignore'>5)<span style='font:7.0pt "Times New Roman"'> </span></span></span><![endif]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>w/ FFTW & --subsetsize 6: 51 s<o:p></o:p></span></p><p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo1'><![if !supportLists]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><span style='mso-list:Ignore'>6)<span style='font:7.0pt "Times New Roman"'> </span></span></span><![endif]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>w/ FFTW & --subsetsize 7: (Error-graphic card seemed to be freezed)<o:p></o:p></span></p><p class=MsoListParagraph style='text-indent:-.25in;mso-list:l0 level1 lfo1'><![if !supportLists]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><span style='mso-list:Ignore'>7)<span style='font:7.0pt "Times New Roman"'> </span></span></span><![endif]><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>w/ FFTW & --subsetsize 8: (Error-“itk::ERROR: CUFFT ERROR #2 in rtkCudaFFTConvolutionImageFilter.cu”)<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>Of course, the difference was coming from the Ramp filter.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>Therefore, I’m going to use [--subsetsize 6] for my system and will update CUDA version soon!<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>Thanks.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>Yang<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>From:</span></b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> simon.rit@gmail.com [mailto:simon.rit@gmail.com] <b>On Behalf Of </b>Simon Rit<br><b>Sent:</b> Friday, November 20, 2015 11:28 AM<br><b>To:</b> Yang-Kyun Park <theday79@gmail.com><br><b>Cc:</b> rtk-users@openrtk.org<br><b>Subject:</b> Re: [Rtk-users] rtkfdk doesn't work correctly after a recent update<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><div><div><p class=MsoNormal>Hi Yang,<br>Thanks a lot for sharing the data. I have good news. I could reproduce the bug and I have solutions. First, about your dataset: it's a huge dataset and the hannY requires zero padding in both x and y. Because of this, cufft has to take the 2D FFT of projection subsets that are 4095x2048x16. Quite big!<o:p></o:p></p></div><div><p class=MsoNormal>When I ran it on my cuda 7.0, I immediately got:<br><br>itk::ExceptionObject (0x8b9cb30)<br>Location: "void CUDA_fft_convolution(const int3&, const int2&, float*, cufftComplex*)" <br>File: /home/srit/src/rtk/rtk/code/rtkCudaFFTConvolutionImageFilter.cu<br>Line: 97<br>Description: itk::ERROR: CUFFT ERROR #2<br><br>i.e., a GPU alloc issue, even if I have the largest possible memory on a nvidia desktop computer. Then, I ran it with cuda 6.5 and I got no error message but a result similar to yours.<o:p></o:p></p></div><div><p class=MsoNormal>Conclusions:<o:p></o:p></p></div><div><p class=MsoNormal>- if you add the option --subsetsize 1 to your rtkfdk command line, that should reconstruct correctly<o:p></o:p></p></div><div><p class=MsoNormal>- if you upgrade cuda, you should get an understandable error message,<o:p></o:p></p></div><div><p class=MsoNormal>- I don't think I'm going to look for a solution because the pb is probably in cufft but if someone is willing to, be my guest!<o:p></o:p></p></div><div><p class=MsoNormal>Simon<o:p></o:p></p></div></div><div><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>On Mon, Nov 16, 2015 at 4:47 PM, Yang-Kyun Park <<a href="mailto:theday79@gmail.com" target="_blank">theday79@gmail.com</a>> wrote:<o:p></o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in'><div><p class=MsoNormal>Hi Simon,<o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>For better understanding, I'm attaching the problematic case including projection files, geometry file(xml) and reconstructed images.(please find the GDrive links. The file size is ~800 MB)<o:p></o:p></p></div><div><p class=MsoNormal>This is phantom data and there is no patient information in it.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Regarding the fftw library, I'm using [fftw-3.3.4-dll64].<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>By the way, I found another clue. When I've tried it with a head phantom image (full-fan, ~ 190 projections), there was no problem. Is there any possibility that DDF and parker short filter play some roles on this issue?<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Feel free to let me know if further data is needed.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Thanks!<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Yang<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p></o:p></p><div style='border:solid #DDDDDD 1.0pt;padding:4.0pt 4.0pt 4.0pt 4.0pt;min-height:18px;max-height:18px'><p class=MsoNormal style='background:whitesmoke'><b><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222'><a href="https://drive.google.com/file/d/0BxuP3mON9V-qbVdQd0luc1lIb2M/view?usp=drive_web" target="_blank"><span style='border:none windowtext 1.0pt;padding:0in;text-decoration:none'><img border=0 id="_x0000_i1025" src="https://ssl.gstatic.com/docs/doclist/images/icon_9_archive_list.png"> </span><span style='color:#1155CC;border:none windowtext 1.0pt;padding:0in;text-decoration:none'>FFTW_issue.zip</span></a><o:p></o:p></span></b></p></div><p class=MsoNormal><o:p></o:p></p></div></div><div><div><div><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>On Mon, Nov 16, 2015 at 10:07 AM, Simon Rit <<a href="mailto:simon.rit@creatis.insa-lyon.fr" target="_blank">simon.rit@creatis.insa-lyon.fr</a>> wrote:<o:p></o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in'><div><div><div><div><p class=MsoNormal>Hi,<o:p></o:p></p></div><p class=MsoNormal>Thanks. It's really weird... When you run the tests with FFTW (compiling the "test" target), do they fail?<br>Can you provide your geometry file? It might be specific to this geometry.<o:p></o:p></p></div><p class=MsoNormal>Thanks again for looking into this issue, we'll find what's the problem.<o:p></o:p></p></div><p class=MsoNormal>Simon<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>On Mon, Nov 16, 2015 at 3:54 PM, Yang K Park <<a href="mailto:theday79@gmail.com" target="_blank">theday79@gmail.com</a>> wrote:<o:p></o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in'><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>Hi Simon,</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>Sorry for the late response.</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>As you suggested, I’ve tried rtkramp+CUDA with only one projection and it was successfully done for both w/ FFTW and w/o FFTW.</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>(Two results were exactly same.)</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>Regarding the CUDA version, I’m using v6.5.</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>Thanks a lot for the help!</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'>Yang</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D'> </span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>From:</span></b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> <a href="mailto:simon.rit@gmail.com" target="_blank">simon.rit@gmail.com</a> [mailto:<a href="mailto:simon.rit@gmail.com" target="_blank">simon.rit@gmail.com</a>] <b>On Behalf Of </b>Simon Rit<br><b>Sent:</b> Thursday, November 12, 2015 5:59 PM</span><o:p></o:p></p><div><div><p class=MsoNormal><br><b>To:</b> Yang K Park <<a href="mailto:theday79@gmail.com" target="_blank">theday79@gmail.com</a>><br><b>Cc:</b> <a href="mailto:rtk-users@openrtk.org" target="_blank">rtk-users@openrtk.org</a><br><b>Subject:</b> Re: [Rtk-users] rtkfdk doesn't work correctly after a recent update<o:p></o:p></p></div></div><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p><div><div><div><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Hi,<o:p></o:p></p></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Ok. Can you try rtkramp with only one projection? Then you shouldn't encounter this regions issue (that should be solved one day but that's a different story). <o:p></o:p></p></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Can you also indicate which version of CUDA you're using?<o:p></o:p></p></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>FYI, we use ITK for the FFT of the kernel of the ramp filter and cufft for the FFT of the projections. So if they are not in the same format, you're in trouble and that might be the issue here... The weird part is that I didn't reproduce this problem on linux (with and without FFTW) but it might be a CUDA version difference. I'll give it another try tomorrow when you provide the CUDA version.<o:p></o:p></p></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Thanks again for reporting the issue,<br>Simon<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>On Thu, Nov 12, 2015 at 11:39 PM, Yang K Park <<a href="mailto:theday79@gmail.com" target="_blank">theday79@gmail.com</a>> wrote:<o:p></o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0in;margin-bottom:5.0pt'><div><div><p>Hi Simon,<o:p></o:p></p><p> <o:p></o:p></p><p>Thanks for the advice. The current RTK was compiled with FFTW to use the phase extraction feature from Amsterdam shroud.<o:p></o:p></p><p> <o:p></o:p></p><p>Regarding your suggestions,<o:p></o:p></p><p> <o:p></o:p></p><p style='margin-left:.5in'>1)<span style='font-size:7.0pt'> </span>rtkramp w/ or w/o CUDA:<o:p></o:p></p><p>(command line: “ rtkramp --path projPath --regexp .*.his --hardware cuda --output Ramp_wCUDA.mha” )<o:p></o:p></p><p> <o:p></o:p></p><p>Then, I have encountered following error.<o:p></o:p></p><p>“CudaFFTRampImageFilter<..> CudaFFTConvolutionImageFilter assumes that input requested and buffered regions are equal.”<o:p></o:p></p><p> <o:p></o:p></p><p><img border=0 width=646 height=88 id="_x0000_i1026" src="cid:image001.png@01D1238E.C6AE3990"><o:p></o:p></p><p> <o:p></o:p></p><p>Without CUDA, it was successfully done.<o:p></o:p></p><p> <o:p></o:p></p><p style='margin-left:.5in'>2)<span style='font-size:7.0pt'> </span>I’ve re-compiled <u>ITK w/o FFTW</u> and built RTK again. r<u>tkfdk worked fine.</u><o:p></o:p></p><p style='margin-left:.75in'>-<span style='font-size:7.0pt'> </span>I’ve also repeated “rtkramp+CUDA” with this rtk version and the error was same as above.<o:p></o:p></p><p style='margin-left:.75in'> <o:p></o:p></p><p>In conclusion, my <u>FFTW seems to affect CUDA-based fdk reconstruction only</u>.<o:p></o:p></p><p>But I do need FFTW for the phase extraction. What should I do?<o:p></o:p></p><p> <o:p></o:p></p><p>Thanks.<o:p></o:p></p><p><span style='color:#888888'> </span><o:p></o:p></p><p><span style='color:#888888'>Yang</span><o:p></o:p></p><div><div><p> <o:p></o:p></p><p> <o:p></o:p></p><p>-----Original Message-----<br>From: <a href="mailto:simon.rit@gmail.com" target="_blank">simon.rit@gmail.com</a> [mailto:<a href="mailto:simon.rit@gmail.com" target="_blank">simon.rit@gmail.com</a>] On Behalf Of Simon Rit<br>Sent: Wednesday, November 11, 2015 4:16 PM<br>To: Yang K Park <<a href="mailto:theday79@gmail.com" target="_blank">theday79@gmail.com</a>><br>Cc: <a href="mailto:rtk-users@openrtk.org" target="_blank">rtk-users@openrtk.org</a><br>Subject: Re: [Rtk-users] rtkfdk doesn't work correctly after a recent update<o:p></o:p></p><p> <o:p></o:p></p><p>Hi,<o:p></o:p></p><p>I didn't manage to reproduce the two problems on linux with the same ITK version.<o:p></o:p></p><p>For point 1, did you copy itkMacro.h by hand? Because I don't see how replacing <itkMacro.h> by "itkMacro.h" can help if itkMacro.h is not in the same folder as rtkCudaUtilities.hcu For point 2, no clue. Do you use FFTW? It looks like a ramp filter issue. Can you try rtkramp --path projPath --regexp .*.his, with and without --hardware cuda? Maybe you can provide your geometry file in case it's a specific geometric issue.<o:p></o:p></p><p>Simon<o:p></o:p></p><p> <o:p></o:p></p><p>On Wed, Nov 11, 2015 at 7:50 PM, Yang K Park <<a href="mailto:theday79@gmail.com" target="_blank"><span style='color:windowtext;text-decoration:none'>theday79@gmail.com</span></a>> wrote:<o:p></o:p></p><p>> Hi all,<o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> I’ve updated rtk with the up-to-date revision and have encountered a <o:p></o:p></p><p>> couple of problems.<o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> 1. Compile error: RTK couldn’t find itkMacro.h while it compiled<o:p></o:p></p><p>> rtkCudaUtilities.hcu. It was solved by a minor fetching:<o:p></o:p></p><p>> <o:p></o:p></p><p>> #include <itkMacro.h> à #include "itkMacro.h" in rtkCudaUtilities.hcu<o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> 2. Output of rtkfdk with CUDA was not the same as before. It was very<o:p></o:p></p><p>> blurry as seen in the attached image.<o:p></o:p></p><p>> <o:p></o:p></p><p>> a. With “FDK with CPU” or “SART with CUDA”, there was no problem.<o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> I’m using VS2013 with win7, itk4.8.1.<o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> The command line I used was :<o:p></o:p></p><p>> <o:p></o:p></p><p>> rtkfdk --geometry 1111.xml --path projPath --verbos --regexp .*.his <o:p></o:p></p><p>> --output output.mha --spacing 1,1,1 --dimension 400,200,400 --hardware <o:p></o:p></p><p>> cuda --pad 1 --hann 5 --hannY 5<o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> Does anyone have any idea about this?<o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> Thanks.<o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> Yang<o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> <o:p></o:p></p><p>> _______________________________________________<o:p></o:p></p><p>> Rtk-users mailing list<o:p></o:p></p><p>> <a href="mailto:Rtk-users@public.kitware.com" target="_blank"><span style='color:windowtext;text-decoration:none'>Rtk-users@public.kitware.com</span></a><o:p></o:p></p><p>> <a href="http://public.kitware.com/mailman/listinfo/rtk-users" target="_blank"><span style='color:windowtext;text-decoration:none'>http://public.kitware.com/mailman/listinfo/rtk-users</span></a><o:p></o:p></p><p>> <o:p></o:p></p></div></div></div></div></blockquote></div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div></div></div></div></div></blockquote></div><p class=MsoNormal><o:p> </o:p></p></div></blockquote></div><p class=MsoNormal><o:p> </o:p></p></div></div></div></blockquote></div><p class=MsoNormal><o:p> </o:p></p></div></div></body></html>