Numerical Solution of Diffusion Models  in Biomedical Imaging on Multicore Processors

<table>Algorithm SI. Behavior of the relative residual <svg height="17.4125" id="M319" style="vertical-align:-3.3907pt" version="1.1" viewbox="0 0 73.412498 17.4125" width="73.412498" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.625,12.575)"><path d="M162 -163h-60v866h60v-866zM309 -163h-60v866h60v-866z" id="x2016"></path></g><g transform="matrix(.017,-0,0,-.017,7.611,12.575)"><path d="M393 379q-9 -16 -28 -29q-15 -10 -23 -2q-19 19 -36 19q-21 0 -52 -38q-57 -72 -82 -126l-40 -197q-23 -3 -75 -18l-7 7q49 196 74 335q7 43 -2 43q-7 0 -30 -14.5t-47 -37.5l-16 23q37 42 82 73t67 31q41 0 15 -113l-11 -50h4q41 71 85 117t77 46q29 0 45 -26
q13 -21 0 -43z" id="x1D45F"></path></g>
<g transform="matrix(.012,-0,0,-.012,14.787,16.65)"><path d="M244 607q0 -25 -15.5 -43t-37.5 -18q-19 0 -32 13t-13 35q0 21 15 41t39 20q20 0 32 -14t12 -34zM222 91q-29 -33 -79 -68t-75 -35q-13 0 -19 7.5t-6 31t10 65.5l62 253q5 26 -1 26q-21 0 -72 -43l-13 24q43 40 91 68t71 28q30 0 10 -78l-71 -274q-8 -30 3 -30
q16 0 76 48z" id="x1D456"></path></g>
<g transform="matrix(.017,-0,0,-.017,18.65,12.575)"><path d="M162 -163h-60v866h60v-866zM309 -163h-60v866h60v-866z" id="x2016"></path></g>
<g transform="matrix(.012,-0,0,-.012,25.637,16.65)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g>
<g transform="matrix(.017,-0,0,-.017,31.975,12.575)"><path d="M368 703l-264 -866h-60l265 866h59z" id="x2F"></path></g><g transform="matrix(.017,-0,0,-.017,38.978,12.575)"><path d="M162 -163h-60v866h60v-866zM309 -163h-60v866h60v-866z" id="x2016"></path></g><g transform="matrix(.017,-0,0,-.017,45.965,12.575)"><path d="M393 379q-9 -16 -28 -29q-15 -10 -23 -2q-19 19 -36 19q-21 0 -52 -38q-57 -72 -82 -126l-40 -197q-23 -3 -75 -18l-7 7q49 196 74 335q7 43 -2 43q-7 0 -30 -14.5t-47 -37.5l-16 23q37 42 82 73t67 31q41 0 15 -113l-11 -50h4q41 71 85 117t77 46q29 0 45 -26
q13 -21 0 -43z" id="x1D45F"></path></g>
<g transform="matrix(.012,-0,0,-.012,53.138,16.65)"><path d="M241 635q53 0 94 -28.5t63.5 -76t33.5 -102.5t11 -116q0 -58 -11 -112.5t-34 -103.5t-63.5 -78.5t-94.5 -29.5t-95 28t-64.5 75t-34.5 102.5t-11 118.5q0 58 11.5 112.5t34.5 103t64.5 78t95.5 29.5zM238 602q-32 0 -55.5 -25t-35.5 -68t-17.5 -91t-5.5 -105
q0 -76 10 -138.5t37 -107.5t69 -45q32 0 55.5 25t35.5 68.5t17.5 91.5t5.5 105t-5.5 105.5t-18 92t-36 68t-56.5 24.5z" id="x30"></path></g>
<g transform="matrix(.017,-0,0,-.017,59.475,12.575)"><path d="M162 -163h-60v866h60v-866zM309 -163h-60v866h60v-866z" id="x2016"></path></g>
<g transform="matrix(.012,-0,0,-.012,66.462,16.65)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g>
</svg> versus 3 scale steps. <svg height="16.8125" id="M320" style="vertical-align:-3.41936pt" version="1.1" viewbox="0 0 78.300003 16.8125" width="78.300003" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.625,11.925)"><path d="M600 0h-557v24l268 633l28 8l261 -641v-24zM497 50l-194 489l-196 -489h390z" id="x394"></path></g><g transform="matrix(.017,-0,0,-.017,11.521,11.925)"><path d="M471 456q-22 -60 -37 -85q-20 -3 -84 -3q-9 0 -29 2t-29 2q-38 -145 -54 -234q-15 -80 14 -80q23 0 68 38l14 -25q-78 -83 -137 -83q-50 0 -50 64q0 25 7 57l67 266q-57 4 -98 -9t-80 -49l-20 21q14 18 27 31t36.5 30t56.5 26t73 9q26 0 82.5 -2t89.5 -2q21 0 32 5t24 24
z" id="x1D70F"></path></g>
<g transform="matrix(.012,-0,0,-.012,19.725,16.012)"><path d="M409 504l-29 -5q-16 60 -44.5 96t-86.5 36q-54 0 -82.5 -32t-28.5 -77q0 -51 30.5 -82.5t96.5 -65.5l35.5 -18t33 -19.5t33.5 -23.5l27 -25.5t24.5 -32t13.5 -36.5t6 -45q0 -80 -62 -134.5t-160 -54.5q-47 0 -98 15q-22 7 -50 21q-8 23 -27 155l30 7q7 -27 18 -52
t30 -51.5t49 -42.5t67 -16q56 0 88 32.5t32 87.5q0 51 -31.5 82t-98.5 67q-80 44 -110 73q-55 53 -55 124q0 75 56 126.5t150 51.5q53 0 126 -23z" id="x53"></path></g><g transform="matrix(.012,-0,0,-.012,25.234,16.012)"><path d="M303 0h-265v28q62 5 76 19.5t14 77.5v400q0 63 -14 77.5t-76 19.5v28h265v-28q-62 -5 -76 -19.5t-14 -77.5v-400q0 -64 14 -78t76 -19v-28z" id="x49"></path></g>
<g transform="matrix(.017,-0,0,-.017,34.625,11.925)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,49.312,11.925)"><path d="M241 635q53 0 94 -28.5t63.5 -76t33.5 -102.5t11 -116q0 -58 -11 -112.5t-34 -103.5t-63.5 -78.5t-94.5 -29.5t-95 28t-64.5 75t-34.5 102.5t-11 118.5q0 58 11.5 112.5t34.5 103t64.5 78t95.5 29.5zM238 602q-32 0 -55.5 -25t-35.5 -68t-17.5 -91t-5.5 -105
q0 -76 10 -138.5t37 -107.5t69 -45q32 0 55.5 25t35.5 68.5t17.5 91.5t5.5 105t-5.5 105.5t-18 92t-36 68t-56.5 24.5z" id="x30"></path></g><g transform="matrix(.017,-0,0,-.017,57.471,11.925)"><path d="M113 -12q-24 0 -39.5 16t-15.5 42q0 24 16 40.5t40 16.5t40 -16.5t16 -40.5q0 -26 -16 -42t-41 -16z" id="x2E"></path></g><g transform="matrix(.017,-0,0,-.017,61.347,11.925)"><path d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" id="x31"></path></g><g transform="matrix(.017,-0,0,-.017,69.523,11.925)"><path d="M137 343l67 33q37 17 63 17q79 0 129.5 -53t50.5 -131q0 -92 -58 -156.5t-147 -64.5t-147 68t-58 182q0 63 17 119t43 95.5t61.5 72t69 52t67.5 31.5q62 22 128 33l6 -32q-56 -11 -108 -35q-149 -71 -184 -231zM227 337q-47 0 -95 -27q-6 -23 -6 -70q0 -93 36 -155.5
t96 -62.5q53 0 78 45.5t25 105.5q0 68 -35 116t-99 48z" id="x36"></path></g>
</svg>.</table>

International Journal of Biomedical Imaging

fig4

Figure 4

Figure 4: Numerical Solution of Diffusion Models  in Biomedical Imaging on Multicore Processors 