<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:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Verdana;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:Aptos;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:12.0pt;
font-family:"Aptos",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Aptos",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></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" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Dear Simon and other RTK developers,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Thanks a lot for your suggestions. So far, I have not been able to find what causes the problem, but I found out a few things which may help.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">The size of the reconstructed volume is 400x400x1313 which with 16-bit values is about 400 megabytes.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">The size of the projections is 1285x1944x1536 (where 1285 is the number of projections), which with 16-bit values is 7.14 gigabytes.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">The maximum number of projections I can use with GPU-based conjugate gradient is 719, which is 3.999 GB. When I use 720 projections (4.005 GB) or more, I get the illegal memory
access error from Cuda as in the previous mail. <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">I have tried to use the nocudacg option, but the threshold remains exactly same. Also, resizing the reconstruction volume has no effect on it.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">I have also tried to run the same script on another machine with a GPU with 8 GB of memory instead of 20 GB, but the limit of the projections that I can load is still 719 (3.999
GB) on both machines. This makes me think the limit is somehow imposed by the RTK software (perhaps a 32-bit limit, since 2^32 bits is equal to 4 gigabytes).
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">The error only happens in the GPU-based version of the conjugate gradient code; the GPU-based SART code works fine (with nprojpersubset=1285), which I think should rule out
problems with my ITK setup. Could there be anything wrong with my RTK setup (I can include the CMakeCache file if it is helpful), or is it more likely that the issue is caused by a bug in the conjugategradient script (or its include files)?<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Kind regards,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Mathé<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"><img width="157" height="76" style="width:1.6354in;height:.7916in" id="Picture_x0020_1400708847" src="cid:image001.png@01DD362D.26DF2790"></span><span style="font-family:"Calibri",sans-serif"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"> </span><span style="font-family:"Calibri",sans-serif"><o:p></o:p></span></p>
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Mathé Zeegers</span></b><span style="font-family:"Calibri",sans-serif"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Preclinical X-ray CT Specialist
</span><span style="font-family:"Calibri",sans-serif"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:#006E77"><a href="mailto:m.zeegers@milabs.com"><span style="color:#006E77">m.zeegers@milabs.com</span></a></span><span style="font-family:"Calibri",sans-serif"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Verdana",sans-serif"> </span><span style="font-family:"Calibri",sans-serif"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:8.0pt;font-family:"Verdana",sans-serif;color:#666666">This message may contain confidential information and is intended exclusively for the addressee. If you receive this message unintentionally, please do not use
the contents but notify the sender immediately by return e-mail. MILabs B.V. is a legal person by public law and is registered at the Chamber of Commerce for Midden-Nederland under no. 302.075.15.</span><span style="font-family:"Calibri",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><span style="font-family:"Calibri",sans-serif"> </span><o:p></o:p></p>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<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 <simon.rit@creatis.insa-lyon.fr>
<br>
<b>Sent:</b> Wednesday, August 26, 2026 2:04 PM<br>
<b>To:</b> Mathé Zeegers <M.Zeegers@milabs.com><br>
<b>Cc:</b> rtk-users@openrtk.org<br>
<b>Subject:</b> Re: Artefacts with GPU-based iterative reconstruction with RTK<o:p></o:p></span></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Dear Mathé,<o:p></o:p></p>
<p class="MsoNormal">Good news! I'm glad it helped.<o:p></o:p></p>
<p class="MsoNormal">I'm not sure about the memory issue. Can you indicate the size of the reconstructed volume and projections? To reduce GPU memory consumption, you can try --nocudacg which avoid computnig the CG step on the GPU but still does the forward
and backprojections on the GPU. You can also monitor the memory consumption if you activate RTK_PROBE_EACH_FILTER to track which filter eats up all the memory (display the report with --verbose).<br>
Cheers,<br>
Simon<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">On Mon, Aug 24, 2026 at 10:54<span style="font-family:"Arial",sans-serif"> </span>AM Mathé Zeegers <<a href="mailto:M.Zeegers@milabs.com">M.Zeegers@milabs.com</a>> wrote:<o:p></o:p></p>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Dear Simon and other RTK developers,</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Thank you for replying so fast! And also for the suggestion to CC the RTK user mailing list.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Your suggestion of changing the step size was exactly the solution that I needed: the reconstructions now look as
they should. </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Setting the step value equal to the voxel size made the conjugate gradient reconstructions initially take very long,
until I set the backprojector from CudaRayCast to CudaVoxelBased. Now I am very pleased with the conjugate gradient performance.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">The only issue that is coming up now, though, is that I get issues with GPU memory when I load too many projections
(720 to be specific, which in our case amounts to 4.00 GB). I get the following error:</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">itk::ExceptionObject (0000000408EFE618)</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Location: "unknown"</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">File: ..\include\itkCudaDataManager.h</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Line: 54</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Description: Cuda Error : an illegal memory access was encountered</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">I am working with a GPU that has 20 GB of memory, so I am a bit surprised that I get this problem. Also, I do not
have this issue with the SART/OSSART/SIRT script, and neither do I have it when I use the CPU-only option. Are you familiar with this issue? Do I perhaps need to set a certain (large-memory) option that I might have missed? The usage of the binning option
solves the problem, but I would like to use the full resolution of the projections. I also tried reducing the reconstruction volume size, but this does not change anything. Any pointers on this would be very helpful.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Regarding the discontinuation of the config option: it is not optimal for us, but we can work around it by putting
back the code for the config option in in our local versions or by writing a small wrapper for it.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Kind regards,</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Mathé</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"><img border="0" width="157" height="76" style="width:1.6354in;height:.7916in" id="m_5930278996077743136Picture_x0020_112524925" src="cid:image001.png@01DD362D.26DF2790"></span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"> </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:10.0pt;font-family:"Verdana",sans-serif">Mathé Zeegers</span></b><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Preclinical X-ray CT Specialist
</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><a href="mailto:m.zeegers@milabs.com" target="_blank"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:#006E77">m.zeegers@milabs.com</span></a><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-family:"Verdana",sans-serif"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:8.0pt;font-family:"Verdana",sans-serif;color:#666666">This message may contain confidential information and is intended exclusively for the addressee. If
you receive this message unintentionally, please do not use the contents but notify the sender immediately by return e-mail. MILabs B.V. is a legal person by public law and is registered at the Chamber of Commerce for Midden-Nederland under no. 302.075.15.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-family:"Calibri",sans-serif"> </span><o:p></o:p></p>
<div style="border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0cm 0cm 0cm;border-color:currentcolor currentcolor">
<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"> Simon Rit <<a href="mailto:simon.rit@creatis.insa-lyon.fr" target="_blank">simon.rit@creatis.insa-lyon.fr</a>>
<br>
<b>Sent:</b> Monday, August 17, 2026 2:57 PM<br>
<b>To:</b> Mathé Zeegers <<a href="mailto:M.Zeegers@milabs.com" target="_blank">M.Zeegers@milabs.com</a>><br>
<b>Cc:</b> rtk-users <<a href="mailto:rtk-users@openrtk.org" target="_blank">rtk-users@openrtk.org</a>><br>
<b>Subject:</b> Re: Artefacts with GPU-based iterative reconstruction with RTK</span><o:p></o:p></p>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Hi Mathé,<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Nice to hear that RTK is useful. I have CCed the RTK user mailing list, you'll probably have faster answers with it in the future. Luckily this time, I'm back from my annual leave
today!<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I bet the problem is the hard coded step length in the GPU ray casting. There's been an update since v2.7.0 to fix this issue<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><a href="https://github.com/RTKConsortium/RTK/commit/990ae96992529ca8f54abe4c1ab48e40950a8092" target="_blank">https://github.com/RTKConsortium/RTK/commit/990ae96992529ca8f54abe4c1ab48e40950a8092</a><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">See <a href="https://docs.openrtk.org/en/latest/documentation/docs/rtk_3_migration_guide.html#cuda-forward-ray-cast-step-size-default" target="_blank">https://docs.openrtk.org/en/latest/documentation/docs/rtk_3_migration_guide.html#cuda-forward-ray-cast-step-size-default</a>.<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Can you try setting the --step option to 0.08 to see if it helps? This would be the new default of the new version.<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I hope it helps,<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Simon<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">PS: sorry about the --config, it's already been merged <a href="https://github.com/RTKConsortium/RTK/commit/d0f15386cbc1c89376753469eea23debcdf78a3f" target="_blank">https://github.com/RTKConsortium/RTK/commit/d0f15386cbc1c89376753469eea23debcdf78a3f</a>.
The goal was to uniformize the Python and C++ applications. Will that be a big issue for you? If yes, please open an issue to trigger a conversation.<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">On Mon, Aug 17, 2026 at 1:51<span style="font-family:"Arial",sans-serif"> </span>PM Mathé Zeegers <<a href="mailto:M.Zeegers@milabs.com" target="_blank">M.Zeegers@milabs.com</a>>
wrote:<o:p></o:p></p>
<blockquote style="border:none;border-left:solid windowtext 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt;border-color:currentcolor currentcolor currentcolor rgb(204,204,204)">
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">Dear Simon Rit and/or other RTK developers,</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">I am writing you to ask for your suggestions as I am stuck with an RTK-related problem for a long time
now. I am not sure if I have used the right e-mail address, so my apologies if I have not.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">Since more than a year, I am working for MILabs as the X-ray CT developer.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">For a long time, we have used RTK for GPU-based FDK reconstructions, and we have only used an SART in
its CPU-based version. We are now looking into integrating (more) GPU-based iterative methods.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">I have implemented a GPU-based version of CGLS using the templates from the RTK github repository (version
2.7.0, commit eb3fd5b). I use this in combination with ITK 5.7.0 and ITKCudaCommon 2.1.0. The problem that keeps coming up is that my reconstructions do not look at all like what they are supposed to look like with the GPU-based reconstruction. I have attached
an image to this mail to clarify this.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">In my reconstruction, I use 720 images with the default settings from the rtkprojectors_section.ggo file.
The number of CGLS iterations is set to 30 and the voxel size is set to 0.08.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">When I run CGLS with forward projector CudaRayCast and backprojector CudaRayCast, the result is very
grainy and unsharp.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">When I change the backprojector to CudaVoxelBased, the result initially looks better with a very low
number of iterations, but with more iterations the reconstruction volume rapidly becomes very stripy and the details disappear as well. As a reference, I also ran the CPU reconstruction and the results look totally fine in this case.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">I have tried to switch to the RTK implementation of SART/OSSART/SIRT, but I encounter the same problem
there as well.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">To double check, I have also used the ASTRA and TIGRE GPU-reconstructions and these give correct results.
Still, I would like to get this GPU-reconstruction to run with RTK since in my opinion it is a more comprehensive toolbox, and the GPU reconstruction is much faster than the CPU reconstruction.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">Do you know what the problem is here and how this can be resolved?</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">Please let me know if something is unclear or if additional information is required.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">Thank you in advance.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">Kind regards,</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">Mathé Zeegers</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:black">PS: In the RTK v3 migration guide it is stated that the -config command-line option will be removed because
it was not known to be used in practice. I would like to let you know that we actually use that feature very frequently.</span><o:p></o:p></p>
<div id="m_5930278996077743136m_-6505260438397447093Signature">
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;color:black"> </span><o:p></o:p></p>
<p style="margin:0cm"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"><img border="0" width="157" height="76" style="width:1.6354in;height:.7916in" id="m_5930278996077743136Picture_x0020_1" src="cid:image001.png@01DD362D.26DF2790"></span><o:p></o:p></p>
<p style="margin:0cm"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif"> </span><o:p></o:p></p>
<p style="margin:0cm"><b><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Mathé Zeegers</span></b><o:p></o:p></p>
<p style="margin:0cm"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif">Preclinical X-ray CT Specialist
</span><o:p></o:p></p>
<p style="margin:0cm"><a href="mailto:m.zeegers@milabs.com" target="_blank"><span style="font-size:10.0pt;font-family:"Verdana",sans-serif;color:#006E77">m.zeegers@milabs.com</span></a><o:p></o:p></p>
<p style="margin:0cm"><span style="font-family:"Verdana",sans-serif"> </span><o:p></o:p></p>
<p style="margin:0cm"><span style="font-size:8.0pt;font-family:"Verdana",sans-serif;color:#666666">This message may contain confidential information and is intended exclusively for the addressee. If you receive this message unintentionally, please do not use
the contents but notify the sender immediately by return e-mail. MILabs B.V. is a legal person by public law and is registered at the Chamber of Commerce for Midden-Nederland under no. 302.075.15.</span><o:p></o:p></p>
<p style="margin:0cm"> <o:p></o:p></p>
</div>
</blockquote>
</div>
</blockquote>
</div>
</body>
</html>