An Efficient Algorithm for EM Scattering from Anatomically Realistic Human Head Model Using Parallel CG-FFT Method

Algorithm 1

/*𝑆𝑢𝑚𝑛𝑜𝑟𝑚𝑜𝑓𝑒𝑎𝑐ℎ𝑃𝑟𝑜𝑐𝑒𝑠𝑠*/

𝑀𝑃𝐼_𝐴𝑙𝑙𝑟𝑒𝑑𝑢𝑐𝑒(&𝑟𝑛𝑜𝑟𝑚,&𝑟𝑛𝑜𝑟𝑚1,1,𝑀𝑃𝐼_𝐹𝐿𝑂𝐴𝑇,𝑀𝑃𝐼_𝑆𝑈𝑀,𝑤𝑜𝑟𝑙𝑑);

/*𝑣𝑒𝑐𝑡𝑜𝑟-𝑣𝑒𝑐𝑡𝑜𝑟𝑃𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑛𝑒𝑎𝑐ℎ𝑝𝑟𝑜𝑐𝑒𝑠𝑠*/

   𝑓𝑜𝑟(𝑖𝑛𝑡𝑖=0;𝑖<𝑎𝑙𝑙;𝑖++)

   {

      𝑐𝑗𝑥[0][0][𝑖]=𝑐𝑗𝑥[0][0][𝑖]+𝑎𝑚∗𝑐𝑝𝑥[0][0][𝑖];

      𝑐𝑗𝑦[0][0][𝑖]=𝑐𝑗𝑦[0][0][𝑖]+𝑎𝑚∗𝑐𝑝𝑦[0][0][𝑖];

      𝑐𝑗𝑧[0][0][𝑖]=𝑐𝑗𝑧[0][0][𝑖]+𝑎𝑚∗𝑐𝑝𝑧[0][0][𝑖];

      𝑐𝑟𝑥[0][0][𝑖]=𝑐𝑟𝑥[0][0][𝑖]−𝑎𝑚∗𝑐𝑡𝑥[0][0][𝑖];

      𝑐𝑟𝑦[0][0][𝑖]=𝑐𝑟𝑦[0][0][𝑖]−𝑎𝑚∗𝑐𝑡𝑦[0][0][𝑖];

      𝑐𝑟𝑧[0][0][𝑖]=𝑐𝑟𝑧[0][0][𝑖]−𝑎𝑚∗𝑐𝑡𝑧[0][0][𝑖];

      𝑟𝑒𝑎𝑥=𝑎𝑏𝑠(𝑐𝑟𝑥[0][0][𝑖]);

      𝑟𝑒𝑎𝑦=𝑎𝑏𝑠(𝑐𝑟𝑦[0][0][𝑖]);

      𝑟𝑒𝑎𝑧=𝑎𝑏𝑠(𝑐𝑟𝑧[0][0][𝑖]);

      𝑟𝑛𝑜𝑟𝑚=𝑟𝑛𝑜𝑟𝑚+𝑟𝑒𝑎𝑥∗𝑟𝑒𝑎𝑥+𝑟𝑒𝑎𝑦∗𝑟𝑒𝑎𝑦+𝑟𝑒𝑎𝑧∗𝑟𝑒𝑎𝑧;

      }