<div dir="ltr"><span style="text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Hamid, Scott,</span><br style="text-decoration-style:initial;text-decoration-color:initial"><br style="text-decoration-style:initial;text-decoration-color:initial"><span style="text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Same as previously: looks like a NIRFAST-Matlab question. Would you like to have this filed in our NIRFAST issue tracker?</span><div style="text-decoration-style:initial;text-decoration-color:initial"><a href="https://github.com/nirfast-admin/NIRFAST/issues">https://github.com/nirfast-admin/NIRFAST/issues</a></div><div><br></div>Best,<div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><span style="color:rgb(136,136,136);font-size:12.8000001907349px">Alexis Girault</span><br style="color:rgb(136,136,136);font-size:12.8000001907349px"><span style="color:rgb(136,136,136);font-size:12.8000001907349px">R&D Engineer in Medical Computing</span><br style="color:rgb(136,136,136);font-size:12.8000001907349px"><span style="color:rgb(136,136,136);font-size:12.8000001907349px">Kitware, Inc.</span><br style="color:rgb(136,136,136);font-size:12.8000001907349px"><br style="color:rgb(136,136,136);font-size:12.8000001907349px"><a href="http://www.kitware.com/" rel="noreferrer" style="color:rgb(17,85,204);font-size:12.8000001907349px" target="_blank">http://www.kitware.com</a><br style="color:rgb(136,136,136);font-size:12.8000001907349px"><font color="#999999"><a href="javascript:void(0);" target="_blank"><span style="font-size:12.8000001907349px">(919) 969-6990 x3</span>25</a></font><br></div></div></div></div></div></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, May 15, 2018 at 5:16 AM Stanislaw Wojtkiewicz <<a href="mailto:s.wojtkiewicz@cs.bham.ac.uk">s.wojtkiewicz@cs.bham.ac.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-GB" link="blue" vlink="purple"><div class="m_-1140050059574826654WordSection1"><p class="MsoNormal">Hi All<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">One more problem to report. <u></u><u></u></p><p class="MsoNormal" style="text-autospace:none">In time-resolved modality (file ‘<span style="font-size:13.0pt;font-family:"Courier New";color:black">femdata_stnd_tr.m</span>’) the incomplete Cholesky factorisation is used to generate a preconditioner:<u></u><u></u></p><p class="MsoNormal" style="text-autospace:none"><span style="font-size:13.0pt;font-family:"Courier New";color:forestgreen">% R = cholinc(MASS1,1e-3);</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p><p class="MsoNormal" style="text-autospace:none">However, Mathworks removed this function completely since it was not a true incomplete Cholesky (it utilized incomplete LU with pivoting).<u></u><u></u></p><p class="MsoNormal" style="text-autospace:none">I suggest to replace the missing function by:<u></u><u></u></p><p class="MsoNormal" style="text-autospace:none"><span style="font-size:13.0pt;font-family:"Courier New";color:black">R = ichol(MASS1,struct(</span><span style="font-size:13.0pt;font-family:"Courier New";color:#a020f0">'type'</span><span style="font-size:13.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:13.0pt;font-family:"Courier New";color:#a020f0">'ict'</span><span style="font-size:13.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:13.0pt;font-family:"Courier New";color:#a020f0">'droptol'</span><span style="font-size:13.0pt;font-family:"Courier New";color:black">,1e-3,</span><span style="font-size:13.0pt;font-family:"Courier New";color:#a020f0">'shape'</span><span style="font-size:13.0pt;font-family:"Courier New";color:black">,</span><span style="font-size:13.0pt;font-family:"Courier New";color:#a020f0">'upper'</span><span style="font-size:13.0pt;font-family:"Courier New";color:black">));</span><span style="font-size:12.0pt;font-family:"Courier New""><u></u><u></u></span></p><p class="MsoNormal" style="text-autospace:none">Works in my case.<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><span lang="EN-US" style="font-family:Consolas">Regards<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US" style="font-family:Consolas">Stanislaw Wojtkiewicz, PhD<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US" style="font-family:Consolas"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US" style="font-family:Consolas">---<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US" style="font-family:Consolas">School of Computer Science<u></u><u></u></span></p><p class="MsoNormal"><span lang="PL" style="font-family:Consolas">University of Birmingham<u></u><u></u></span></p><p class="MsoNormal"><span lang="PL" style="font-family:Consolas">B15 2TT<u></u><u></u></span></p><p class="MsoNormal"><u></u> <u></u></p></div><div id="m_-1140050059574826654DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br>
<table style="border-top:1px solid #d3d4de">
<tr>
<td style="width:55px;padding-top:13px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient" target="_blank"><img src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" alt="" width="46" height="29" style="width:46px;height:29px"></a></td>
<td style="width:470px;padding-top:12px;color:#41424e;font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Virus-free. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient" style="color:#4453ea" target="_blank">www.avast.com</a>
</td>
</tr>
</table><a href="#m_-1140050059574826654_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1"> </a></div></div>_______________________________________________<br>
Nirfast mailing list<br>
<a href="mailto:Nirfast@public.kitware.com" target="_blank">Nirfast@public.kitware.com</a><br>
<a href="https://public.kitware.com/mailman/listinfo/nirfast" rel="noreferrer" target="_blank">https://public.kitware.com/mailman/listinfo/nirfast</a><br>
</blockquote></div>