ཻࢠ๏ʹ͍ͭͯ
ཻࢠ๏ͱ
ཻࢠ๏ͱ
ྲྀମ
ཻࢠ๏ͱ
ྲྀମ
ࢄԽʂ
ཻࢠ๏ͱ
ྲྀମ ཻࢠͷ͋ͭ·Γ
ࢄԽʂ
ཻࢠ๏ͱ
ʮཻࢠʯ
ɾҰൠͷཧ๏ଇʹै͏Θ͚Ͱͳ͍
ɾݶΒΕͨൣғͰ૬ޓ࡞༻͢Δ
ɾͱҐஔϕΫτϧΛͭ
ɾͦΕࣗମ͕ܭࢉج४ʹͳΔ
֨ࢠΛ͏߹ͱͷҧ͍
֨ࢠΛ͏߹ͱͷҧ͍
֨ࢠ๏
ɾ֨ࢠ͕มΛ࣋ͭ
ɾઢͰͭͳ͕ͬͨॴͰͷ૬ޓ࡞༻
ɾ֨ࢠੜ͕ܭࢉྔతʹॏ͍
ɾෳࡶͳࣗ༝ද໘ͷܭࢉ͕͍͠
ɾ࣍ݩ͕૿͑Δͱෳࡶʹ
֨ࢠΛ͏߹ͱͷҧ͍
֨ࢠ๏
ɾ֨ࢠ͕มΛ࣋ͭ
ɾઢͰͭͳ͕ͬͨॴͰͷ૬ޓ࡞༻
ɾ֨ࢠੜ͕ܭࢉྔతʹॏ͍
ɾෳࡶͳࣗ༝ද໘ͷܭࢉ͕͍͠
ɾ࣍ݩ͕૿͑Δͱෳࡶʹ
ཻࢠ๏
ɾཻࢠ͕มΛ࣋ͭ
ɾӨڹܘͰͷ૬ޓ࡞༻
ɾཻࢠΛฒΔ͚ͩ
ɾϥάϥϯδϡهड़Ͱྲྀମʹ࠷ద
ɾࢄԽ͕ࣜ୯७
༷ʑͳछྨ
.14๏
ʢ.PWJOH1BSUJDMF4JNVMBUJPOʣ
41)
41".๏
ʢ4NPPUIFE1BSUJDMF"QQMJFE.FDIBOJDTʣ
γϛϡϨʔγϣϯͷྲྀΕ
ཻࢠΛฒΔ
ԾΛܭࢉ͢Δ
ԾΛͬͯҠಈ͢Δ
ඇѹॖ݅Λຬͨ͢Α͏ʹ
ѹྗΛܭࢉ
ҐஔͱΛ࠶ܭࢉ
ྲྀମΛࢧ͢Δ̎ͭͷࣜ
φϏΤετʔΫεํఔࣜ
ɾྲྀମͷӡಈํఔࣜ
ɾV͕ɺH͕ॏྗՃ
ɾЛ͕ີɺӈลͷୈೋ߲೪ੑ߲
Du
Dt
=
1
⇢0
rP + ⌫r2u + g
AAACP3icbVBLSwMxGMz6rPVV9eglWAShWHaLoBehaBGPFewDum3Jptk2NMkuSVYoy/4j/4R/wZvoyYs38erNtLsHbR0IGWa+4UvGCxlV2rZfrKXlldW19dxGfnNre2e3sLffVEEkMWnggAWy7SFFGBWkoalmpB1KgrjHSMsbX0/91gORigbiXk9C0uVoKKhPMdJG6hduXF8iHNdcj8dRksQ1nVyeppqTxK4cBT07cQXyGIL1kiuilPcqME2Uptcw6ReKdtmeAS4SJyNFkKHeL7y7gwBHnAiNGVKq49ih7sZIaooZSfJupEiI8BgNScdQgThR3Xj23wQeG2UA/UCaIzScqb8TMeJKTbhnJjnSIzXvTcV/PY/Pbdb+RTemIow0EThd7EcM6gBOy4QDKgnWbGIIwpKat0M8QqY8bSrPm1Kc+QoWSbNSduyyc3dWrF5l9eTAITgCJ8AB56AKbkEdNAAGj+AZvII368n6sD6tr3R0ycoyB+APrO8fxL+wyg==
AAACP3icbVBLSwMxGMz6rPVV9eglWAShWHaLoBehaBGPFewDum3Jptk2NMkuSVYoy/4j/4R/wZvoyYs38erNtLsHbR0IGWa+4UvGCxlV2rZfrKXlldW19dxGfnNre2e3sLffVEEkMWnggAWy7SFFGBWkoalmpB1KgrjHSMsbX0/91gORigbiXk9C0uVoKKhPMdJG6hduXF8iHNdcj8dRksQ1nVyeppqTxK4cBT07cQXyGIL1kiuilPcqME2Uptcw6ReKdtmeAS4SJyNFkKHeL7y7gwBHnAiNGVKq49ih7sZIaooZSfJupEiI8BgNScdQgThR3Xj23wQeG2UA/UCaIzScqb8TMeJKTbhnJjnSIzXvTcV/PY/Pbdb+RTemIow0EThd7EcM6gBOy4QDKgnWbGIIwpKat0M8QqY8bSrPm1Kc+QoWSbNSduyyc3dWrF5l9eTAITgCJ8AB56AKbkEdNAAGj+AZvII368n6sD6tr3R0ycoyB+APrO8fxL+wyg==
AAACP3icbVBLSwMxGMz6rPVV9eglWAShWHaLoBehaBGPFewDum3Jptk2NMkuSVYoy/4j/4R/wZvoyYs38erNtLsHbR0IGWa+4UvGCxlV2rZfrKXlldW19dxGfnNre2e3sLffVEEkMWnggAWy7SFFGBWkoalmpB1KgrjHSMsbX0/91gORigbiXk9C0uVoKKhPMdJG6hduXF8iHNdcj8dRksQ1nVyeppqTxK4cBT07cQXyGIL1kiuilPcqME2Uptcw6ReKdtmeAS4SJyNFkKHeL7y7gwBHnAiNGVKq49ih7sZIaooZSfJupEiI8BgNScdQgThR3Xj23wQeG2UA/UCaIzScqb8TMeJKTbhnJjnSIzXvTcV/PY/Pbdb+RTemIow0EThd7EcM6gBOy4QDKgnWbGIIwpKat0M8QqY8bSrPm1Kc+QoWSbNSduyyc3dWrF5l9eTAITgCJ8AB56AKbkEdNAAGj+AZvII368n6sD6tr3R0ycoyB+APrO8fxL+wyg==
AAACP3icbVBLSwMxGMz6rPVV9eglWAShWHaLoBehaBGPFewDum3Jptk2NMkuSVYoy/4j/4R/wZvoyYs38erNtLsHbR0IGWa+4UvGCxlV2rZfrKXlldW19dxGfnNre2e3sLffVEEkMWnggAWy7SFFGBWkoalmpB1KgrjHSMsbX0/91gORigbiXk9C0uVoKKhPMdJG6hduXF8iHNdcj8dRksQ1nVyeppqTxK4cBT07cQXyGIL1kiuilPcqME2Uptcw6ReKdtmeAS4SJyNFkKHeL7y7gwBHnAiNGVKq49ih7sZIaooZSfJupEiI8BgNScdQgThR3Xj23wQeG2UA/UCaIzScqb8TMeJKTbhnJjnSIzXvTcV/PY/Pbdb+RTemIow0EThd7EcM6gBOy4QDKgnWbGIIwpKat0M8QqY8bSrPm1Kc+QoWSbNSduyyc3dWrF5l9eTAITgCJ8AB56AKbkEdNAAGj+AZvII368n6sD6tr3R0ycoyB+APrO8fxL+wyg==
࿈ଓͷࣜ
ɾྲྀମͷ࣭ྔอଘଇ
D⇢
Dt
+ ⇢r · u = 0
AAACIXicbZC7SgNBFIbPeo3xFrW0GRRBEMKujWkEUQtLBaOCG8LsZNYMzmWZOSuEZV/Bl/AVbLW3C3ZiIfgkThILjf4w8POfczhnviSTwmEYvgUTk1PTM7OVuer8wuLScm1l9cKZ3DLeZEYae5VQx6XQvIkCJb/KLKcqkfwyuT0a1C/vuHXC6HPsZbyl6I0WqWAUfdSubceppaw4jm3XlMUxljsDF2uaSBqzjsE4UUVe7oft2mZYD4cif030bTYPGh99AgCn7dpn3DEsV1wjk9S56yjMsFVQi4JJXlbj3PGMslt6w6+91VRx1yqGPyrJlk86JDXWP41kmP6cKKhyrqcS36kodt14bRD+W0vU2GZMG61C6CxHrtlocZpLgoYMcJGOsJyh7HlDmRX+dsK61CNDD7XqoUTjCP6ai916FNajM0/nEEaqwDpswDZEsAcHcAKn0AQG9/AIT/AcPAQvQT94HbVOBN8za/BLwfsXy0Gm4g==
AAACIXicbZC7SgNBFIZn4z3eVi1EbBZFEISwa2MaQTSFZQSTCG4Is5NZM2Quy8xZISz7Cr6Er2CrvV2wUwvBysdwcik0+sPAz3/O4Zz5ooQzA77/6hSmpmdm5+YXiotLyyur7tp63ahUE1ojiit9FWFDOZO0Bgw4vUo0xSLitBF1zwb1xi3Vhil5Cb2ENgW+kSxmBIONWu5+GGtMskqoOyrPKpAfDFwoccRxSNoKwkhkaX7st9xdv+QP5f01wdjsnpQ/+ptvX1vVlvsZthVJBZVAODbmOvATaGZYAyOc5sUwNTTBpItv6LW1EgtqmtnwR7m3Z5O2FyttnwRvmP6cyLAwpici2ykwdMxkbRD+W4vExGaIy82MySQFKslocZxyD5Q3wOW1maYEeM8aTDSzt3ukgy0ysFCLFkowieCvqR+WAr8UXFg6p2ikebSNdtA+CtAROkHnqIpqiKA79IAe0ZNz7zw7fedl1FpwxjMb6Jec929Xyai0
AAACIXicbZC7SgNBFIZn4z3eVi1EbBZFEISwa2MaQTSFZQSTCG4Is5NZM2Quy8xZISz7Cr6Er2CrvV2wUwvBysdwcik0+sPAz3/O4Zz5ooQzA77/6hSmpmdm5+YXiotLyyur7tp63ahUE1ojiit9FWFDOZO0Bgw4vUo0xSLitBF1zwb1xi3Vhil5Cb2ENgW+kSxmBIONWu5+GGtMskqoOyrPKpAfDFwoccRxSNoKwkhkaX7st9xdv+QP5f01wdjsnpQ/+ptvX1vVlvsZthVJBZVAODbmOvATaGZYAyOc5sUwNTTBpItv6LW1EgtqmtnwR7m3Z5O2FyttnwRvmP6cyLAwpici2ykwdMxkbRD+W4vExGaIy82MySQFKslocZxyD5Q3wOW1maYEeM8aTDSzt3ukgy0ysFCLFkowieCvqR+WAr8UXFg6p2ikebSNdtA+CtAROkHnqIpqiKA79IAe0ZNz7zw7fedl1FpwxjMb6Jec929Xyai0
AAACIXicbZC7SgNBFIZnvcZ4i1raLAYhIIRdG22EoCksI5gLZEOYncwmQ+ayzJwVwrKv4Ev4Crba24md2PkkziYpNPGHgZ//nMM584UxZwY879NZWV1b39gsbBW3d3b39ksHhy2jEk1okyiudCfEhnImaRMYcNqJNcUi5LQdjm/yevuBasOUvIdJTHsCDyWLGMFgo36pEkQak7Qe6JHK0jpkZ7kLJA45DshAQRCKNMmuvH6p7FW9qdxl489NGc3V6Je+g4EiiaASCMfGdH0vhl6KNTDCaVYMEkNjTMZ4SLvWSiyo6aXTH2XuqU0GbqS0fRLcafp7IsXCmIkIbafAMDKLtTz8txaKhc0QXfZSJuMEqCSzxVHCXVBujssdME0J8Ik1mGhmb3fJCFtkYKEWLRR/EcGyaZ1Xfa/q33nl2vUcTwEdoxNUQT66QDV0ixqoiQh6RM/oBb06T86b8+58zFpXnPnMEfoj5+sHkwakmw==
ϓϩάϥϜ্ͰͷࢄԽ
ཻࢠͷஔ
ཻࢠؒڑ l0
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZZBG8sEzQOSJc5OZpMh81hmZoWw5hNstbcTW3/GxsIvcfIoNPHAhcM593IuJ0o4M9b3P73cyura+kZ+s7C1vbO7V9w/aBiVakLrRHGlWxE2lDNJ65ZZTluJplhEnDaj4fXEbz5QbZiSd3aU0FDgvmQxI9g66ZZ3/W6x5Jf9KdAyCeakVDmqfd0DQLVb/O70FEkFlZZwbEw78BMbZlhbRjgdFzqpoQkmQ9ynbUclFtSE2fTVMTp1Sg/FSruRFk3V3xcZFsaMROQ2BbYDs+hNxH+9SCwk2/gyzJhMUkslmQXHKUdWoUkPqMc0JZaPHMFEM/c7IgOsMbGurYIrJVisYJk0zsuBXw5qrp0rmCEPx3ACZxDABVTgBqpQBwJ9eIJnePEevVfvzXufrea8+c0h/IH38QPkBJYa
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckS5idzCZDZmaXmVkhrPkEW+3tROz8BH/CxsIvcTZJoYkHLhzOuZdzOUHMmTau++nkVlbX1jfym4Wt7Z3dveL+QVNHiSK0QSIeqXaANeVM0oZhhtN2rCgWAaetYHSV+a07qjSL5K0Zx9QXeCBZyAg2VrrhPbdXLLlldwq0TLw5KVWO6l/srfpR6xW/u/2IJIJKQzjWuuO5sfFTrAwjnE4K3UTTGJMRHtCOpRILqv10+uoEnVqlj8JI2ZEGTdXfFykWWo9FYDcFNkO96GXiv14gFpJNeOmnTMaJoZLMgsOEIxOhrAfUZ4oSw8eWYKKY/R2RIVaYGNtWwZbiLVawTJrnZc8te3XbThVmyMMxnMAZeHABFbiGGjSAwAAe4BGenHvn2XlxXmerOWd+cwh/4Lz/AFEUl9Y=
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckS5idzCZDZmaXmVkhrPkEW+3tROz8BH/CxsIvcTZJoYkHLhzOuZdzOUHMmTau++nkVlbX1jfym4Wt7Z3dveL+QVNHiSK0QSIeqXaANeVM0oZhhtN2rCgWAaetYHSV+a07qjSL5K0Zx9QXeCBZyAg2VrrhPbdXLLlldwq0TLw5KVWO6l/srfpR6xW/u/2IJIJKQzjWuuO5sfFTrAwjnE4K3UTTGJMRHtCOpRILqv10+uoEnVqlj8JI2ZEGTdXfFykWWo9FYDcFNkO96GXiv14gFpJNeOmnTMaJoZLMgsOEIxOhrAfUZ4oSw8eWYKKY/R2RIVaYGNtWwZbiLVawTJrnZc8te3XbThVmyMMxnMAZeHABFbiGGjSAwAAe4BGenHvn2XlxXmerOWd+cwh/4Lz/AFEUl9Y=
AAAB+nicbVA9TwJBFHyHX4hfqKXNRmJiRfZstCTaWGL0gAQuZG/Zgw27e5fdPRNy8hNstbcztv4ZW3+JC1yh4CQvmcy8l3mZKBXcWIy/vNLa+sbmVnm7srO7t39QPTxqmSTTlAU0EYnuRMQwwRULLLeCdVLNiIwEa0fjm5nffmTa8EQ92EnKQkmGisecEuuke9HH/WoN1/EcaJX4BalBgWa/+t0bJDSTTFkqiDFdH6c2zIm2nAo2rfQyw1JCx2TIuo4qIpkJ8/mrU3TmlAGKE+1GWTRXf1/kRBozkZHblMSOzLI3E//1IrmUbOOrMOcqzSxTdBEcZwLZBM16QAOuGbVi4gihmrvfER0RTah1bVVcKf5yBaukdVH3cd2/w7XGdVFPGU7gFM7Bh0towC00IQAKQ3iGF3j1nrw37937WKyWvOLmGP7A+/wBNxyUNQ==
un = 0
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKsjTZC0MYyAfPAZImzk9lkyMzsMjMrhCWVv2CrvZ3Y+iU2Fn6Jk8RCEw9cOJxzL+dywkRwYzH+8HJLyyura/n1wsbm1vZOcXevYeJUU1ansYh1KySGCa5Y3XIrWCvRjMhQsGY4vJr4zXumDY/VjR0lLJCkr3jEKbFOuu2EMkvHXXWBu8USLuMp0CLxf0ipclD7vAOAarf41enFNJVMWSqIMW0fJzbIiLacCjYudFLDEkKHpM/ajioimQmy6cdjdOyUHopi7UZZNFV/X2REGjOSoduUxA7MvDcR//VCOZdso/Mg4ypJLVN0FhylAtkYTepAPa4ZtWLkCKGau98RHRBNqHWlFVwp/nwFi6RxWvZx2a+5di5hhjwcwhGcgA9nUIFrqEIdKCh4hCd49h68F+/Ve5ut5ryfm334A+/9G3QBmTc=
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKuTWyEEBvLBMwDkyXMTmaTITOzy8ysEJZU/oKt9haC2Nr6EzYWfomTR6GJBy4czrmXczlBzJk2rvvpZFZW19Y3spu5re2d3b38/kFDR4kitE4iHqlWgDXlTNK6YYbTVqwoFgGnzWB4NfGbd1RpFskbM4qpL3BfspARbKx02wlEmoy78tLt5gtu0Z0CLRNvTgrlo9oXe6l8VLv5704vIomg0hCOtW57bmz8FCvDCKfjXCfRNMZkiPu0banEgmo/nX48RqdW6aEwUnakQVP190WKhdYjEdhNgc1AL3oT8V8vEAvJJrzwUybjxFBJZsFhwpGJ0KQO1GOKEsNHlmCimP0dkQFWmBhbWs6W4i1WsEwa50XPLXo1204FZsjCMZzAGXhQgjJcQxXqQEDCAzzCk3PvPDuvzttsNePMbw7hD5z3H+ECmvM=
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKuTWyEEBvLBMwDkyXMTmaTITOzy8ysEJZU/oKt9haC2Nr6EzYWfomTR6GJBy4czrmXczlBzJk2rvvpZFZW19Y3spu5re2d3b38/kFDR4kitE4iHqlWgDXlTNK6YYbTVqwoFgGnzWB4NfGbd1RpFskbM4qpL3BfspARbKx02wlEmoy78tLt5gtu0Z0CLRNvTgrlo9oXe6l8VLv5704vIomg0hCOtW57bmz8FCvDCKfjXCfRNMZkiPu0banEgmo/nX48RqdW6aEwUnakQVP190WKhdYjEdhNgc1AL3oT8V8vEAvJJrzwUybjxFBJZsFhwpGJ0KQO1GOKEsNHlmCimP0dkQFWmBhbWs6W4i1WsEwa50XPLXo1204FZsjCMZzAGXhQgjJcQxXqQEDCAzzCk3PvPDuvzttsNePMbw7hD5z3H+ECmvM=
AAACAXicbVC7SgNBFL0bXzG+opY2g0GwCrM2phGCNpYRzAOTJcxOZpMhM7PLzKwQllT+gq32dmLrl9j6JU6SLTTxwIXDOfdyLidMBDcW4y+vsLa+sblV3C7t7O7tH5QPj1omTjVlTRqLWHdCYpjgijUtt4J1Es2IDAVrh+Obmd9+ZNrwWN3bScICSYaKR5wS66SHXiizdNpXV7hfruAqngOtEj8nFcjR6Je/e4OYppIpSwUxpuvjxAYZ0ZZTwaalXmpYQuiYDFnXUUUkM0E2/3iKzpwyQFGs3SiL5urvi4xIYyYydJuS2JFZ9mbiv14ol5JtVAsyrpLUMkUXwVEqkI3RrA404JpRKyaOEKq5+x3REdGEWldayZXiL1ewSloXVR9X/TtcqV/n9RThBE7hHHy4hDrcQgOaQEHBM7zAq/fkvXnv3sditeDlN8fwB97nD8cKl1I=
particleType[i] = WALL;
particleType[i] = FLUID;
ཻࢠͷஔ
ཻࢠؒڑ l0
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZZBG8sEzQOSJc5OZpMh81hmZoWw5hNstbcTW3/GxsIvcfIoNPHAhcM593IuJ0o4M9b3P73cyura+kZ+s7C1vbO7V9w/aBiVakLrRHGlWxE2lDNJ65ZZTluJplhEnDaj4fXEbz5QbZiSd3aU0FDgvmQxI9g66ZZ3/W6x5Jf9KdAyCeakVDmqfd0DQLVb/O70FEkFlZZwbEw78BMbZlhbRjgdFzqpoQkmQ9ynbUclFtSE2fTVMTp1Sg/FSruRFk3V3xcZFsaMROQ2BbYDs+hNxH+9SCwk2/gyzJhMUkslmQXHKUdWoUkPqMc0JZaPHMFEM/c7IgOsMbGurYIrJVisYJk0zsuBXw5qrp0rmCEPx3ACZxDABVTgBqpQBwJ9eIJnePEevVfvzXufrea8+c0h/IH38QPkBJYa
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckS5idzCZDZmaXmVkhrPkEW+3tROz8BH/CxsIvcTZJoYkHLhzOuZdzOUHMmTau++nkVlbX1jfym4Wt7Z3dveL+QVNHiSK0QSIeqXaANeVM0oZhhtN2rCgWAaetYHSV+a07qjSL5K0Zx9QXeCBZyAg2VrrhPbdXLLlldwq0TLw5KVWO6l/srfpR6xW/u/2IJIJKQzjWuuO5sfFTrAwjnE4K3UTTGJMRHtCOpRILqv10+uoEnVqlj8JI2ZEGTdXfFykWWo9FYDcFNkO96GXiv14gFpJNeOmnTMaJoZLMgsOEIxOhrAfUZ4oSw8eWYKKY/R2RIVaYGNtWwZbiLVawTJrnZc8te3XbThVmyMMxnMAZeHABFbiGGjSAwAAe4BGenHvn2XlxXmerOWd+cwh/4Lz/AFEUl9Y=
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckS5idzCZDZmaXmVkhrPkEW+3tROz8BH/CxsIvcTZJoYkHLhzOuZdzOUHMmTau++nkVlbX1jfym4Wt7Z3dveL+QVNHiSK0QSIeqXaANeVM0oZhhtN2rCgWAaetYHSV+a07qjSL5K0Zx9QXeCBZyAg2VrrhPbdXLLlldwq0TLw5KVWO6l/srfpR6xW/u/2IJIJKQzjWuuO5sfFTrAwjnE4K3UTTGJMRHtCOpRILqv10+uoEnVqlj8JI2ZEGTdXfFykWWo9FYDcFNkO96GXiv14gFpJNeOmnTMaJoZLMgsOEIxOhrAfUZ4oSw8eWYKKY/R2RIVaYGNtWwZbiLVawTJrnZc8te3XbThVmyMMxnMAZeHABFbiGGjSAwAAe4BGenHvn2XlxXmerOWd+cwh/4Lz/AFEUl9Y=
AAAB+nicbVA9TwJBFHyHX4hfqKXNRmJiRfZstCTaWGL0gAQuZG/Zgw27e5fdPRNy8hNstbcztv4ZW3+JC1yh4CQvmcy8l3mZKBXcWIy/vNLa+sbmVnm7srO7t39QPTxqmSTTlAU0EYnuRMQwwRULLLeCdVLNiIwEa0fjm5nffmTa8EQ92EnKQkmGisecEuuke9HH/WoN1/EcaJX4BalBgWa/+t0bJDSTTFkqiDFdH6c2zIm2nAo2rfQyw1JCx2TIuo4qIpkJ8/mrU3TmlAGKE+1GWTRXf1/kRBozkZHblMSOzLI3E//1IrmUbOOrMOcqzSxTdBEcZwLZBM16QAOuGbVi4gihmrvfER0RTah1bVVcKf5yBaukdVH3cd2/w7XGdVFPGU7gFM7Bh0towC00IQAKQ3iGF3j1nrw37937WKyWvOLmGP7A+/wBNxyUNQ==
un = 0
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKsjTZC0MYyAfPAZImzk9lkyMzsMjMrhCWVv2CrvZ3Y+iU2Fn6Jk8RCEw9cOJxzL+dywkRwYzH+8HJLyyura/n1wsbm1vZOcXevYeJUU1ansYh1KySGCa5Y3XIrWCvRjMhQsGY4vJr4zXumDY/VjR0lLJCkr3jEKbFOuu2EMkvHXXWBu8USLuMp0CLxf0ipclD7vAOAarf41enFNJVMWSqIMW0fJzbIiLacCjYudFLDEkKHpM/ajioimQmy6cdjdOyUHopi7UZZNFV/X2REGjOSoduUxA7MvDcR//VCOZdso/Mg4ypJLVN0FhylAtkYTepAPa4ZtWLkCKGau98RHRBNqHWlFVwp/nwFi6RxWvZx2a+5di5hhjwcwhGcgA9nUIFrqEIdKCh4hCd49h68F+/Ve5ut5ryfm334A+/9G3QBmTc=
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKuTWyEEBvLBMwDkyXMTmaTITOzy8ysEJZU/oKt9haC2Nr6EzYWfomTR6GJBy4czrmXczlBzJk2rvvpZFZW19Y3spu5re2d3b38/kFDR4kitE4iHqlWgDXlTNK6YYbTVqwoFgGnzWB4NfGbd1RpFskbM4qpL3BfspARbKx02wlEmoy78tLt5gtu0Z0CLRNvTgrlo9oXe6l8VLv5704vIomg0hCOtW57bmz8FCvDCKfjXCfRNMZkiPu0banEgmo/nX48RqdW6aEwUnakQVP190WKhdYjEdhNgc1AL3oT8V8vEAvJJrzwUybjxFBJZsFhwpGJ0KQO1GOKEsNHlmCimP0dkQFWmBhbWs6W4i1WsEwa50XPLXo1204FZsjCMZzAGXhQgjJcQxXqQEDCAzzCk3PvPDuvzttsNePMbw7hD5z3H+ECmvM=
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKuTWyEEBvLBMwDkyXMTmaTITOzy8ysEJZU/oKt9haC2Nr6EzYWfomTR6GJBy4czrmXczlBzJk2rvvpZFZW19Y3spu5re2d3b38/kFDR4kitE4iHqlWgDXlTNK6YYbTVqwoFgGnzWB4NfGbd1RpFskbM4qpL3BfspARbKx02wlEmoy78tLt5gtu0Z0CLRNvTgrlo9oXe6l8VLv5704vIomg0hCOtW57bmz8FCvDCKfjXCfRNMZkiPu0banEgmo/nX48RqdW6aEwUnakQVP190WKhdYjEdhNgc1AL3oT8V8vEAvJJrzwUybjxFBJZsFhwpGJ0KQO1GOKEsNHlmCimP0dkQFWmBhbWs6W4i1WsEwa50XPLXo1204FZsjCMZzAGXhQgjJcQxXqQEDCAzzCk3PvPDuvzttsNePMbw7hD5z3H+ECmvM=
AAACAXicbVC7SgNBFL0bXzG+opY2g0GwCrM2phGCNpYRzAOTJcxOZpMhM7PLzKwQllT+gq32dmLrl9j6JU6SLTTxwIXDOfdyLidMBDcW4y+vsLa+sblV3C7t7O7tH5QPj1omTjVlTRqLWHdCYpjgijUtt4J1Es2IDAVrh+Obmd9+ZNrwWN3bScICSYaKR5wS66SHXiizdNpXV7hfruAqngOtEj8nFcjR6Je/e4OYppIpSwUxpuvjxAYZ0ZZTwaalXmpYQuiYDFnXUUUkM0E2/3iKzpwyQFGs3SiL5urvi4xIYyYydJuS2JFZ9mbiv14ol5JtVAsyrpLUMkUXwVEqkI3RrA404JpRKyaOEKq5+x3REdGEWldayZXiL1ewSloXVR9X/TtcqV/n9RThBE7hHHy4hDrcQgOaQEHBM7zAq/fkvXnv3sditeDlN8fwB97nD8cKl1I=
࠷ॳͷͯ͢
ʹͦΖ͓͑ͯ͘
particleType[i] = WALL;
particleType[i] = FLUID;
ཻࢠͷஔ
ཻࢠؒڑ l0
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZZBG8sEzQOSJc5OZpMh81hmZoWw5hNstbcTW3/GxsIvcfIoNPHAhcM593IuJ0o4M9b3P73cyura+kZ+s7C1vbO7V9w/aBiVakLrRHGlWxE2lDNJ65ZZTluJplhEnDaj4fXEbz5QbZiSd3aU0FDgvmQxI9g66ZZ3/W6x5Jf9KdAyCeakVDmqfd0DQLVb/O70FEkFlZZwbEw78BMbZlhbRjgdFzqpoQkmQ9ynbUclFtSE2fTVMTp1Sg/FSruRFk3V3xcZFsaMROQ2BbYDs+hNxH+9SCwk2/gyzJhMUkslmQXHKUdWoUkPqMc0JZaPHMFEM/c7IgOsMbGurYIrJVisYJk0zsuBXw5qrp0rmCEPx3ACZxDABVTgBqpQBwJ9eIJnePEevVfvzXufrea8+c0h/IH38QPkBJYa
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckS5idzCZDZmaXmVkhrPkEW+3tROz8BH/CxsIvcTZJoYkHLhzOuZdzOUHMmTau++nkVlbX1jfym4Wt7Z3dveL+QVNHiSK0QSIeqXaANeVM0oZhhtN2rCgWAaetYHSV+a07qjSL5K0Zx9QXeCBZyAg2VrrhPbdXLLlldwq0TLw5KVWO6l/srfpR6xW/u/2IJIJKQzjWuuO5sfFTrAwjnE4K3UTTGJMRHtCOpRILqv10+uoEnVqlj8JI2ZEGTdXfFykWWo9FYDcFNkO96GXiv14gFpJNeOmnTMaJoZLMgsOEIxOhrAfUZ4oSw8eWYKKY/R2RIVaYGNtWwZbiLVawTJrnZc8te3XbThVmyMMxnMAZeHABFbiGGjSAwAAe4BGenHvn2XlxXmerOWd+cwh/4Lz/AFEUl9Y=
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckS5idzCZDZmaXmVkhrPkEW+3tROz8BH/CxsIvcTZJoYkHLhzOuZdzOUHMmTau++nkVlbX1jfym4Wt7Z3dveL+QVNHiSK0QSIeqXaANeVM0oZhhtN2rCgWAaetYHSV+a07qjSL5K0Zx9QXeCBZyAg2VrrhPbdXLLlldwq0TLw5KVWO6l/srfpR6xW/u/2IJIJKQzjWuuO5sfFTrAwjnE4K3UTTGJMRHtCOpRILqv10+uoEnVqlj8JI2ZEGTdXfFykWWo9FYDcFNkO96GXiv14gFpJNeOmnTMaJoZLMgsOEIxOhrAfUZ4oSw8eWYKKY/R2RIVaYGNtWwZbiLVawTJrnZc8te3XbThVmyMMxnMAZeHABFbiGGjSAwAAe4BGenHvn2XlxXmerOWd+cwh/4Lz/AFEUl9Y=
AAAB+nicbVA9TwJBFHyHX4hfqKXNRmJiRfZstCTaWGL0gAQuZG/Zgw27e5fdPRNy8hNstbcztv4ZW3+JC1yh4CQvmcy8l3mZKBXcWIy/vNLa+sbmVnm7srO7t39QPTxqmSTTlAU0EYnuRMQwwRULLLeCdVLNiIwEa0fjm5nffmTa8EQ92EnKQkmGisecEuuke9HH/WoN1/EcaJX4BalBgWa/+t0bJDSTTFkqiDFdH6c2zIm2nAo2rfQyw1JCx2TIuo4qIpkJ8/mrU3TmlAGKE+1GWTRXf1/kRBozkZHblMSOzLI3E//1IrmUbOOrMOcqzSxTdBEcZwLZBM16QAOuGbVi4gihmrvfER0RTah1bVVcKf5yBaukdVH3cd2/w7XGdVFPGU7gFM7Bh0towC00IQAKQ3iGF3j1nrw37937WKyWvOLmGP7A+/wBNxyUNQ==
un = 0
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKsjTZC0MYyAfPAZImzk9lkyMzsMjMrhCWVv2CrvZ3Y+iU2Fn6Jk8RCEw9cOJxzL+dywkRwYzH+8HJLyyura/n1wsbm1vZOcXevYeJUU1ansYh1KySGCa5Y3XIrWCvRjMhQsGY4vJr4zXumDY/VjR0lLJCkr3jEKbFOuu2EMkvHXXWBu8USLuMp0CLxf0ipclD7vAOAarf41enFNJVMWSqIMW0fJzbIiLacCjYudFLDEkKHpM/ajioimQmy6cdjdOyUHopi7UZZNFV/X2REGjOSoduUxA7MvDcR//VCOZdso/Mg4ypJLVN0FhylAtkYTepAPa4ZtWLkCKGau98RHRBNqHWlFVwp/nwFi6RxWvZx2a+5di5hhjwcwhGcgA9nUIFrqEIdKCh4hCd49h68F+/Ve5ut5ryfm334A+/9G3QBmTc=
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKuTWyEEBvLBMwDkyXMTmaTITOzy8ysEJZU/oKt9haC2Nr6EzYWfomTR6GJBy4czrmXczlBzJk2rvvpZFZW19Y3spu5re2d3b38/kFDR4kitE4iHqlWgDXlTNK6YYbTVqwoFgGnzWB4NfGbd1RpFskbM4qpL3BfspARbKx02wlEmoy78tLt5gtu0Z0CLRNvTgrlo9oXe6l8VLv5704vIomg0hCOtW57bmz8FCvDCKfjXCfRNMZkiPu0banEgmo/nX48RqdW6aEwUnakQVP190WKhdYjEdhNgc1AL3oT8V8vEAvJJrzwUybjxFBJZsFhwpGJ0KQO1GOKEsNHlmCimP0dkQFWmBhbWs6W4i1WsEwa50XPLXo1204FZsjCMZzAGXhQgjJcQxXqQEDCAzzCk3PvPDuvzttsNePMbw7hD5z3H+ECmvM=
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKuTWyEEBvLBMwDkyXMTmaTITOzy8ysEJZU/oKt9haC2Nr6EzYWfomTR6GJBy4czrmXczlBzJk2rvvpZFZW19Y3spu5re2d3b38/kFDR4kitE4iHqlWgDXlTNK6YYbTVqwoFgGnzWB4NfGbd1RpFskbM4qpL3BfspARbKx02wlEmoy78tLt5gtu0Z0CLRNvTgrlo9oXe6l8VLv5704vIomg0hCOtW57bmz8FCvDCKfjXCfRNMZkiPu0banEgmo/nX48RqdW6aEwUnakQVP190WKhdYjEdhNgc1AL3oT8V8vEAvJJrzwUybjxFBJZsFhwpGJ0KQO1GOKEsNHlmCimP0dkQFWmBhbWs6W4i1WsEwa50XPLXo1204FZsjCMZzAGXhQgjJcQxXqQEDCAzzCk3PvPDuvzttsNePMbw7hD5z3H+ECmvM=
AAACAXicbVC7SgNBFL0bXzG+opY2g0GwCrM2phGCNpYRzAOTJcxOZpMhM7PLzKwQllT+gq32dmLrl9j6JU6SLTTxwIXDOfdyLidMBDcW4y+vsLa+sblV3C7t7O7tH5QPj1omTjVlTRqLWHdCYpjgijUtt4J1Es2IDAVrh+Obmd9+ZNrwWN3bScICSYaKR5wS66SHXiizdNpXV7hfruAqngOtEj8nFcjR6Je/e4OYppIpSwUxpuvjxAYZ0ZZTwaalXmpYQuiYDFnXUUUkM0E2/3iKzpwyQFGs3SiL5urvi4xIYyYydJuS2JFZ9mbiv14ol5JtVAsyrpLUMkUXwVEqkI3RrA404JpRKyaOEKq5+x3REdGEWldayZXiL1ewSloXVR9X/TtcqV/n9RThBE7hHHy4hDrcQgOaQEHBM7zAq/fkvXnv3sditeDlN8fwB97nD8cKl1I=
࠷ॳͷͯ͢
ʹͦΖ͓͑ͯ͘
particleType[i] = WALL;
particleType[i] = FLUID;
นͳͷ͔ྲྀମͳͷ͔ͷ
λΠϓΛ͓ͬͯ͘
ཻࢠͷஔ
ཻࢠؒڑ l0
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZZBG8sEzQOSJc5OZpMh81hmZoWw5hNstbcTW3/GxsIvcfIoNPHAhcM593IuJ0o4M9b3P73cyura+kZ+s7C1vbO7V9w/aBiVakLrRHGlWxE2lDNJ65ZZTluJplhEnDaj4fXEbz5QbZiSd3aU0FDgvmQxI9g66ZZ3/W6x5Jf9KdAyCeakVDmqfd0DQLVb/O70FEkFlZZwbEw78BMbZlhbRjgdFzqpoQkmQ9ynbUclFtSE2fTVMTp1Sg/FSruRFk3V3xcZFsaMROQ2BbYDs+hNxH+9SCwk2/gyzJhMUkslmQXHKUdWoUkPqMc0JZaPHMFEM/c7IgOsMbGurYIrJVisYJk0zsuBXw5qrp0rmCEPx3ACZxDABVTgBqpQBwJ9eIJnePEevVfvzXufrea8+c0h/IH38QPkBJYa
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckS5idzCZDZmaXmVkhrPkEW+3tROz8BH/CxsIvcTZJoYkHLhzOuZdzOUHMmTau++nkVlbX1jfym4Wt7Z3dveL+QVNHiSK0QSIeqXaANeVM0oZhhtN2rCgWAaetYHSV+a07qjSL5K0Zx9QXeCBZyAg2VrrhPbdXLLlldwq0TLw5KVWO6l/srfpR6xW/u/2IJIJKQzjWuuO5sfFTrAwjnE4K3UTTGJMRHtCOpRILqv10+uoEnVqlj8JI2ZEGTdXfFykWWo9FYDcFNkO96GXiv14gFpJNeOmnTMaJoZLMgsOEIxOhrAfUZ4oSw8eWYKKY/R2RIVaYGNtWwZbiLVawTJrnZc8te3XbThVmyMMxnMAZeHABFbiGGjSAwAAe4BGenHvn2XlxXmerOWd+cwh/4Lz/AFEUl9Y=
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckS5idzCZDZmaXmVkhrPkEW+3tROz8BH/CxsIvcTZJoYkHLhzOuZdzOUHMmTau++nkVlbX1jfym4Wt7Z3dveL+QVNHiSK0QSIeqXaANeVM0oZhhtN2rCgWAaetYHSV+a07qjSL5K0Zx9QXeCBZyAg2VrrhPbdXLLlldwq0TLw5KVWO6l/srfpR6xW/u/2IJIJKQzjWuuO5sfFTrAwjnE4K3UTTGJMRHtCOpRILqv10+uoEnVqlj8JI2ZEGTdXfFykWWo9FYDcFNkO96GXiv14gFpJNeOmnTMaJoZLMgsOEIxOhrAfUZ4oSw8eWYKKY/R2RIVaYGNtWwZbiLVawTJrnZc8te3XbThVmyMMxnMAZeHABFbiGGjSAwAAe4BGenHvn2XlxXmerOWd+cwh/4Lz/AFEUl9Y=
AAAB+nicbVA9TwJBFHyHX4hfqKXNRmJiRfZstCTaWGL0gAQuZG/Zgw27e5fdPRNy8hNstbcztv4ZW3+JC1yh4CQvmcy8l3mZKBXcWIy/vNLa+sbmVnm7srO7t39QPTxqmSTTlAU0EYnuRMQwwRULLLeCdVLNiIwEa0fjm5nffmTa8EQ92EnKQkmGisecEuuke9HH/WoN1/EcaJX4BalBgWa/+t0bJDSTTFkqiDFdH6c2zIm2nAo2rfQyw1JCx2TIuo4qIpkJ8/mrU3TmlAGKE+1GWTRXf1/kRBozkZHblMSOzLI3E//1IrmUbOOrMOcqzSxTdBEcZwLZBM16QAOuGbVi4gihmrvfER0RTah1bVVcKf5yBaukdVH3cd2/w7XGdVFPGU7gFM7Bh0towC00IQAKQ3iGF3j1nrw37937WKyWvOLmGP7A+/wBNxyUNQ==
ཻࢠؒڑʹԠͯ͡
ྻͯ͠ฒ͍ͯ͘
un = 0
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKsjTZC0MYyAfPAZImzk9lkyMzsMjMrhCWVv2CrvZ3Y+iU2Fn6Jk8RCEw9cOJxzL+dywkRwYzH+8HJLyyura/n1wsbm1vZOcXevYeJUU1ansYh1KySGCa5Y3XIrWCvRjMhQsGY4vJr4zXumDY/VjR0lLJCkr3jEKbFOuu2EMkvHXXWBu8USLuMp0CLxf0ipclD7vAOAarf41enFNJVMWSqIMW0fJzbIiLacCjYudFLDEkKHpM/ajioimQmy6cdjdOyUHopi7UZZNFV/X2REGjOSoduUxA7MvDcR//VCOZdso/Mg4ypJLVN0FhylAtkYTepAPa4ZtWLkCKGau98RHRBNqHWlFVwp/nwFi6RxWvZx2a+5di5hhjwcwhGcgA9nUIFrqEIdKCh4hCd49h68F+/Ve5ut5ryfm334A+/9G3QBmTc=
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKuTWyEEBvLBMwDkyXMTmaTITOzy8ysEJZU/oKt9haC2Nr6EzYWfomTR6GJBy4czrmXczlBzJk2rvvpZFZW19Y3spu5re2d3b38/kFDR4kitE4iHqlWgDXlTNK6YYbTVqwoFgGnzWB4NfGbd1RpFskbM4qpL3BfspARbKx02wlEmoy78tLt5gtu0Z0CLRNvTgrlo9oXe6l8VLv5704vIomg0hCOtW57bmz8FCvDCKfjXCfRNMZkiPu0banEgmo/nX48RqdW6aEwUnakQVP190WKhdYjEdhNgc1AL3oT8V8vEAvJJrzwUybjxFBJZsFhwpGJ0KQO1GOKEsNHlmCimP0dkQFWmBhbWs6W4i1WsEwa50XPLXo1204FZsjCMZzAGXhQgjJcQxXqQEDCAzzCk3PvPDuvzttsNePMbw7hD5z3H+ECmvM=
AAACAXicbVC7SgNBFL0bXzG+ooKNzWAQrMKuTWyEEBvLBMwDkyXMTmaTITOzy8ysEJZU/oKt9haC2Nr6EzYWfomTR6GJBy4czrmXczlBzJk2rvvpZFZW19Y3spu5re2d3b38/kFDR4kitE4iHqlWgDXlTNK6YYbTVqwoFgGnzWB4NfGbd1RpFskbM4qpL3BfspARbKx02wlEmoy78tLt5gtu0Z0CLRNvTgrlo9oXe6l8VLv5704vIomg0hCOtW57bmz8FCvDCKfjXCfRNMZkiPu0banEgmo/nX48RqdW6aEwUnakQVP190WKhdYjEdhNgc1AL3oT8V8vEAvJJrzwUybjxFBJZsFhwpGJ0KQO1GOKEsNHlmCimP0dkQFWmBhbWs6W4i1WsEwa50XPLXo1204FZsjCMZzAGXhQgjJcQxXqQEDCAzzCk3PvPDuvzttsNePMbw7hD5z3H+ECmvM=
AAACAXicbVC7SgNBFL0bXzG+opY2g0GwCrM2phGCNpYRzAOTJcxOZpMhM7PLzKwQllT+gq32dmLrl9j6JU6SLTTxwIXDOfdyLidMBDcW4y+vsLa+sblV3C7t7O7tH5QPj1omTjVlTRqLWHdCYpjgijUtt4J1Es2IDAVrh+Obmd9+ZNrwWN3bScICSYaKR5wS66SHXiizdNpXV7hfruAqngOtEj8nFcjR6Je/e4OYppIpSwUxpuvjxAYZ0ZZTwaalXmpYQuiYDFnXUUUkM0E2/3iKzpwyQFGs3SiL5urvi4xIYyYydJuS2JFZ9mbiv14ol5JtVAsyrpLUMkUXwVEqkI3RrA404JpRKyaOEKq5+x3REdGEWldayZXiL1ewSloXVR9X/TtcqV/n9RThBE7hHHy4hDrcQgOaQEHBM7zAq/fkvXnv3sditeDlN8fwB97nD8cKl1I=
࠷ॳͷͯ͢
ʹͦΖ͓͑ͯ͘
particleType[i] = WALL;
particleType[i] = FLUID;
นͳͷ͔ྲྀମͳͷ͔ͷ
λΠϓΛ͓ͬͯ͘
ఆܭࢉ
n0 =
X
j6=i0
w |r0
j
r0
i0
|
AAACNHicbZBNSxxBEIZrzJfZfG3iMZdGCTGHLD1eIoTAopfkZiCrgrM79PTW7LZ294zdNcoy7o/xT+QveE2OAU+GXP0Hgr0fh2RNQcPD+1ZR1W9WauWJ81/R0r37Dx4+Wn7cePL02fMXzZevdn1ROYkdWejC7WfCo1YWO6RI437pUJhM4152tD3x907QeVXYbzQqsWvEwKpcSUFBSpsfbY9/Snxl0vowsXjM1NvxaaIxp/WzJDO1G6eHPf5+jnVwe/wscWowpHdpc423+LTYXYjnsNb+cnPSBoCdtHmV9AtZGbQktfD+IOYldWvhSEmN40ZSeSyFPBIDPAhohUHfraefHLM3QemzvHDhWWJT9e+JWhjvRyYLnUbQ0C96E/G/XmYWNlO+2a2VLStCK2eL80ozKtgkQdZXDiXpUQAhnQq3MzkUTkgKOTdCKPFiBHdhd6MV81b8NaSzBbNahtewCusQwwdow2fYgQ5IOIcL+AE/o+/RZfQ7+jNrXYrmMyvwT0XXt0tcrjY=
AAACNHicbZBNTxsxEIa9tLQ0/UrLsRerUdVwaOTtBaSqFSqXwgmkJiBlk5XXmU0cbO/WnqWKlvyY/okee+XaHpE4gbjyDyrhfBwgMJKlR+87oxm/Sa6kQ8ZOg6UHD5cfPV55Unn67PmLl9VXr1suK6yApshUZg8S7kBJA02UqOAgt8B1omA/Odya+PtHYJ3MzHcc5dDRvG9kKgVHL8XVT6bLPkeu0HE5jAz8oPL9+GekIMX6cZTo0o7jYZd9mGPp3S47jqzsD3AtrtZYg02L3oVwDrXN7f9HX3b+DHfj6nnUy0ShwaBQ3Ll2yHLslNyiFArGlahwkHNxyPvQ9mi4Btcpp58c03de6dE0s/4ZpFP15kTJtXMjnfhOzXHgFr2JeK+X6IXNmG50SmnyAsGI2eK0UBQzOkmQ9qQFgWrkgQsr/e1UDLjlAn3OFR9KuBjBXWh9bISsEe75dL6SWa2QN+QtqZOQrJNN8o3skiYR5Bc5IX/Jv+B3cBZcBJez1qVgPrNKblVwdQ1N9q+m
AAACNHicbZBNTxsxEIa9tLQ0/UrLsRerUdVwaOTtBaSqFSqXwgmkJiBlk5XXmU0cbO/WnqWKlvyY/okee+XaHpE4gbjyDyrhfBwgMJKlR+87oxm/Sa6kQ8ZOg6UHD5cfPV55Unn67PmLl9VXr1suK6yApshUZg8S7kBJA02UqOAgt8B1omA/Odya+PtHYJ3MzHcc5dDRvG9kKgVHL8XVT6bLPkeu0HE5jAz8oPL9+GekIMX6cZTo0o7jYZd9mGPp3S47jqzsD3AtrtZYg02L3oVwDrXN7f9HX3b+DHfj6nnUy0ShwaBQ3Ll2yHLslNyiFArGlahwkHNxyPvQ9mi4Btcpp58c03de6dE0s/4ZpFP15kTJtXMjnfhOzXHgFr2JeK+X6IXNmG50SmnyAsGI2eK0UBQzOkmQ9qQFgWrkgQsr/e1UDLjlAn3OFR9KuBjBXWh9bISsEe75dL6SWa2QN+QtqZOQrJNN8o3skiYR5Bc5IX/Jv+B3cBZcBJez1qVgPrNKblVwdQ1N9q+m
AAACNHicbZDPThsxEMa9UCiEP03psReLCAEHIm8vVEKVEFx6pFIDSGyy8jqzicH2LvYsKFryMLwEr8AVjkicqLjyDDhhD23oSJZ++r4ZzfhLciUdMvYQTE1/mJn9ODdfW1hcWv5U/7xy6LLCCmiJTGX2OOEOlDTQQokKjnMLXCcKjpKz/ZF/dAHWycz8xkEObc17RqZScPRSXN8xHfYjcoWOy9PIwDmV68PLSEGKG1dRoks7jE87bKvC0rsddhVZ2evjZlxvsCYbF30PYQUNUtVBXH+KupkoNBgUijt3ErIc2yW3KIWCYS0qHORcnPEenHg0XINrl+NPDumaV7o0zax/BulY/Xui5Nq5gU58p+bYd5PeSPyvl+iJzZh+b5fS5AWCEW+L00JRzOgoQdqVFgSqgQcurPS3U9Hnlgv0Odd8KOFkBO/h8FszZM3wF2vs7lXxzJGvZJVskJBsk13ykxyQFhHkmtySO3If3ASPwZ/g+a11KqhmvpB/Knh5BSKEq/o=
ཻࢠີͷج४
Өڹܘʹ͋Δཻۙࢠͱͷڑͷ̎ͷॏΈฏۉ
0 =
P
j6=i0
|r0
j
r0
i0
|2w |r0
j
r0
i0
|
P
j6=i0
w |r0
j
r0
i0
|
AAAClnicjVFNbxMxEPVu+SjhowEuSL1YVIj2QOQNB7iAqqKqPQaJtJXqZOV1ZhO3tnexZ0HRdv8EdyT+Vq/9JXU+Dm0KEiNZenrvzYz1Jiu18sjYZRSv3bv/4OH6o9bjJ0+fbbSfvzjyReUk9GWhC3eSCQ9aWeijQg0npQNhMg3H2fmXmX78A5xXhf2G0xIGRoytypUUGKi0/ZvrYB6JIfvEcydkzX1l0vqMW/hO1dvmgmemdk16NmTvlrAO9JBdDLs/uYYct/9t4U6NJ7jTrA7978a0vcU6bF70LkiWYGt3f+/XH0JIL21f8VEhKwMWpRbenyasxEEtHCqpoWnxykMp5LkYw2mAVhjwg3oeY0PfBGZE88KFZ5HO2ZsdtTDeT00WnEbgxK9qM/KvWmZWNmP+cVArW1YIVi4W55WmWNDZjehIOZCopwEI6VT4O5UTEa6D4ZKtEEqyGsFdcNTtJKyTfA3p7JFFrZNN8ppsk4R8ILvkkPRIn8hoLdqJutH7+FX8Od6PDxbWOFr2vCS3Ku5dAw9Xztg=
AAAClnicjVFNbxMxEPVuKZTwlcIFiYtFBbQHIm84wIWqKqrgmArSVqqTldeZTdza3q09SxVt909w58Sv4sovwfk4lBQkRrL09N6bGetNVmrlkbGfUbx2a/32nY27rXv3Hzx81N58fOSLyknoy0IX7iQTHrSy0EeFGk5KB8JkGo6z8w8z/fgrOK8K+wWnJQyMGFuVKykwUGn7O9fBPBJD9p7nTsia+8qk9Rm3cEHVq+aKZ6Z2TXo2ZK+XsA70kF0Nu5dcQ47b/7Zwp8YT3GlWh/53Y9reYh02L3oTJEuwtXew/+3H+ueXvbT9i48KWRmwKLXw/jRhJQ5q4VBJDU2LVx5KIc/FGE4DtMKAH9TzGBv6IjAjmhcuPIt0zl7vqIXxfmqy4DQCJ35Vm5F/1TKzshnzd4Na2bJCsHKxOK80xYLObkRHyoFEPQ1ASKfC36mciHAdDJdshVCS1QhugqNuJ2Gd5DCks08WtUGekedkmyTkLdkjn0iP9ImM1qKdqBu9iZ/Gu/FB/HFhjaNlzxPyR8W938Oiz1k=
AAAClnicjVFNbxMxEPVuKZTwlcIFiYtFBbQHIm84wIWqKqrgmArSVqqTldeZTdza3q09SxVt909w58Sv4sovwfk4lBQkRrL09N6bGetNVmrlkbGfUbx2a/32nY27rXv3Hzx81N58fOSLyknoy0IX7iQTHrSy0EeFGk5KB8JkGo6z8w8z/fgrOK8K+wWnJQyMGFuVKykwUGn7O9fBPBJD9p7nTsia+8qk9Rm3cEHVq+aKZ6Z2TXo2ZK+XsA70kF0Nu5dcQ47b/7Zwp8YT3GlWh/53Y9reYh02L3oTJEuwtXew/+3H+ueXvbT9i48KWRmwKLXw/jRhJQ5q4VBJDU2LVx5KIc/FGE4DtMKAH9TzGBv6IjAjmhcuPIt0zl7vqIXxfmqy4DQCJ35Vm5F/1TKzshnzd4Na2bJCsHKxOK80xYLObkRHyoFEPQ1ASKfC36mciHAdDJdshVCS1QhugqNuJ2Gd5DCks08WtUGekedkmyTkLdkjn0iP9ImM1qKdqBu9iZ/Gu/FB/HFhjaNlzxPyR8W938Oiz1k=
AAAClnicjVFNbxMxEPVuC5Tw0VAuSFwsIkR7IPKGA1yoKlBVjqlE2kp1svI6s4lb27vYs6Bouz+FH8a1v6ROugdIQWIkS0/vvZmx3mSlVh4Z+xXFG5v37j/Yeth59PjJ0+3us50TX1ROwkgWunBnmfCglYURKtRwVjoQJtNwml1+Xuqn38F5VdivuChhbMTMqlxJgYFKuz+5DuapmLCPPHdC1txXJq0vuIVvVL1prnhmatekFxP2toV1oCfsajL4wTXkuPtvC3dqNse9Zn3ofzem3R7rs1XRuyBpQY+0NUy713xayMqARamF9+cJK3FcC4dKamg6vPJQCnkpZnAeoBUG/LhexdjQ14GZ0rxw4VmkK/b3jloY7xcmC04jcO7XtSX5Vy0za5sx/zCulS0rBCtvF+eVpljQ5Y3oVDmQqBcBCOlU+DuVcxGug+GSnRBKsh7BXXAy6Cesnxyz3sGnNp4t8pK8IrskIe/JAflChmREZLQR7UWD6F38It6PD+OjW2sctT3PyR8VD28ARh3M3w==
ॏΈͱཻࢠີ
ॏΈͱཻࢠີ
re
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZZBG8sEzQOSJc5O7iZDZnaXmVkhrPkEW+3txNafsbHwS5w8Ck08cOFwzr2cywkSwbVx3U8nt7K6tr6R3yxsbe/s7hX3Dxo6ThXDOotFrFoB1Sh4hHXDjcBWopDKQGAzGF5P/OYDKs3j6M6MEvQl7Uc85IwaK92qLnaLJbfsTkGWiTcnpcpR7eseAKrd4nenF7NUYmSYoFq3PTcxfkaV4UzguNBJNSaUDWkf25ZGVKL2s+mrY3JqlR4JY2UnMmSq/r7IqNR6JAO7KakZ6EVvIv7rBXIh2YSXfsajJDUYsVlwmApiYjLpgfS4QmbEyBLKFLe/EzagijJj2yrYUrzFCpZJ47zsuWWvZtu5ghnycAwncAYeXEAFbqAKdWDQhyd4hhfn0Xl13pz32WrOmd8cwh84Hz9BO5ZV
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckIcxO7iZDZmeXmVkhrPkEW+3tROz8BH/CxsIvcfIoNPHAhcM593Iux48F18Z1P53Myura+kZ2M7e1vbO7l98/qOsoUQxrLBKRavpUo+ASa4Ybgc1YIQ19gQ1/eDXxG3eoNI/krRnF2AlpX/KAM2qsdKO62M0X3KI7BVkm3pwUSkfVL/5W/qh089/tXsSSEKVhgmrd8tzYdFKqDGcCx7l2ojGmbEj72LJU0hB1J52+OianVumRIFJ2pCFT9fdFSkOtR6FvN0NqBnrRm4j/en64kGyCy07KZZwYlGwWHCSCmIhMeiA9rpAZMbKEMsXt74QNqKLM2LZythRvsYJlUj8vem7Rq9p2yjBDFo7hBM7AgwsowTVUoAYM+vAAj/Dk3DvPzovzOlvNOPObQ/gD5/0HrjyYEQ==
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckIcxO7iZDZmeXmVkhrPkEW+3tROz8BH/CxsIvcfIoNPHAhcM593Iux48F18Z1P53Myura+kZ2M7e1vbO7l98/qOsoUQxrLBKRavpUo+ASa4Ybgc1YIQ19gQ1/eDXxG3eoNI/krRnF2AlpX/KAM2qsdKO62M0X3KI7BVkm3pwUSkfVL/5W/qh089/tXsSSEKVhgmrd8tzYdFKqDGcCx7l2ojGmbEj72LJU0hB1J52+OianVumRIFJ2pCFT9fdFSkOtR6FvN0NqBnrRm4j/en64kGyCy07KZZwYlGwWHCSCmIhMeiA9rpAZMbKEMsXt74QNqKLM2LZythRvsYJlUj8vem7Rq9p2yjBDFo7hBM7AgwsowTVUoAYM+vAAj/Dk3DvPzovzOlvNOPObQ/gD5/0HrjyYEQ==
AAAB+nicbVC7SgNBFL0bXzG+opY2g0GwCrs2WgZtLCOaByRLmJ3cTYbMzC4zs0KI+QRb7e3E1p+x9UucJFto4oELh3Pu5VxOlApurO9/eYW19Y3NreJ2aWd3b/+gfHjUNEmmGTZYIhLdjqhBwRU2LLcC26lGKiOBrWh0M/Nbj6gNT9SDHacYSjpQPOaMWifd6x72yhW/6s9BVkmQkwrkqPfK391+wjKJyjJBjekEfmrDCdWWM4HTUjczmFI2ogPsOKqoRBNO5q9OyZlT+iROtBtlyVz9fTGh0pixjNympHZolr2Z+K8XyaVkG1+FE67SzKJii+A4E8QmZNYD6XONzIqxI5Rp7n4nbEg1Zda1VXKlBMsVrJLmRTXwq8GdX6ld5/UU4QRO4RwCuIQa3EIdGsBgAM/wAq/ek/fmvXsfi9WCl98cwx94nz+URJRw
Өڹܘ
ॏΈͱཻࢠີ
re
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZZBG8sEzQOSJc5O7iZDZnaXmVkhrPkEW+3txNafsbHwS5w8Ck08cOFwzr2cywkSwbVx3U8nt7K6tr6R3yxsbe/s7hX3Dxo6ThXDOotFrFoB1Sh4hHXDjcBWopDKQGAzGF5P/OYDKs3j6M6MEvQl7Uc85IwaK92qLnaLJbfsTkGWiTcnpcpR7eseAKrd4nenF7NUYmSYoFq3PTcxfkaV4UzguNBJNSaUDWkf25ZGVKL2s+mrY3JqlR4JY2UnMmSq/r7IqNR6JAO7KakZ6EVvIv7rBXIh2YSXfsajJDUYsVlwmApiYjLpgfS4QmbEyBLKFLe/EzagijJj2yrYUrzFCpZJ47zsuWWvZtu5ghnycAwncAYeXEAFbqAKdWDQhyd4hhfn0Xl13pz32WrOmd8cwh84Hz9BO5ZV
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckIcxO7iZDZmeXmVkhrPkEW+3tROz8BH/CxsIvcfIoNPHAhcM593Iux48F18Z1P53Myura+kZ2M7e1vbO7l98/qOsoUQxrLBKRavpUo+ASa4Ybgc1YIQ19gQ1/eDXxG3eoNI/krRnF2AlpX/KAM2qsdKO62M0X3KI7BVkm3pwUSkfVL/5W/qh089/tXsSSEKVhgmrd8tzYdFKqDGcCx7l2ojGmbEj72LJU0hB1J52+OianVumRIFJ2pCFT9fdFSkOtR6FvN0NqBnrRm4j/en64kGyCy07KZZwYlGwWHCSCmIhMeiA9rpAZMbKEMsXt74QNqKLM2LZythRvsYJlUj8vem7Rq9p2yjBDFo7hBM7AgwsowTVUoAYM+vAAj/Dk3DvPzovzOlvNOPObQ/gD5/0HrjyYEQ==
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckIcxO7iZDZmeXmVkhrPkEW+3tROz8BH/CxsIvcfIoNPHAhcM593Iux48F18Z1P53Myura+kZ2M7e1vbO7l98/qOsoUQxrLBKRavpUo+ASa4Ybgc1YIQ19gQ1/eDXxG3eoNI/krRnF2AlpX/KAM2qsdKO62M0X3KI7BVkm3pwUSkfVL/5W/qh089/tXsSSEKVhgmrd8tzYdFKqDGcCx7l2ojGmbEj72LJU0hB1J52+OianVumRIFJ2pCFT9fdFSkOtR6FvN0NqBnrRm4j/en64kGyCy07KZZwYlGwWHCSCmIhMeiA9rpAZMbKEMsXt74QNqKLM2LZythRvsYJlUj8vem7Rq9p2yjBDFo7hBM7AgwsowTVUoAYM+vAAj/Dk3DvPzovzOlvNOPObQ/gD5/0HrjyYEQ==
AAAB+nicbVC7SgNBFL0bXzG+opY2g0GwCrs2WgZtLCOaByRLmJ3cTYbMzC4zs0KI+QRb7e3E1p+x9UucJFto4oELh3Pu5VxOlApurO9/eYW19Y3NreJ2aWd3b/+gfHjUNEmmGTZYIhLdjqhBwRU2LLcC26lGKiOBrWh0M/Nbj6gNT9SDHacYSjpQPOaMWifd6x72yhW/6s9BVkmQkwrkqPfK391+wjKJyjJBjekEfmrDCdWWM4HTUjczmFI2ogPsOKqoRBNO5q9OyZlT+iROtBtlyVz9fTGh0pixjNympHZolr2Z+K8XyaVkG1+FE67SzKJii+A4E8QmZNYD6XONzIqxI5Rp7n4nbEg1Zda1VXKlBMsVrJLmRTXwq8GdX6ld5/UU4QRO4RwCuIQa3EIdGsBgAM/wAq/ek/fmvXsfi9WCl98cwx94nz+URJRw
Өڹܘ
w(r) =
(
re
r
1 (r < re)
0 (r re)
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
ॏΈؔͷఆٛ
ॏΈͱཻࢠີ
re
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZZBG8sEzQOSJc5O7iZDZnaXmVkhrPkEW+3txNafsbHwS5w8Ck08cOFwzr2cywkSwbVx3U8nt7K6tr6R3yxsbe/s7hX3Dxo6ThXDOotFrFoB1Sh4hHXDjcBWopDKQGAzGF5P/OYDKs3j6M6MEvQl7Uc85IwaK92qLnaLJbfsTkGWiTcnpcpR7eseAKrd4nenF7NUYmSYoFq3PTcxfkaV4UzguNBJNSaUDWkf25ZGVKL2s+mrY3JqlR4JY2UnMmSq/r7IqNR6JAO7KakZ6EVvIv7rBXIh2YSXfsajJDUYsVlwmApiYjLpgfS4QmbEyBLKFLe/EzagijJj2yrYUrzFCpZJ47zsuWWvZtu5ghnycAwncAYeXEAFbqAKdWDQhyd4hhfn0Xl13pz32WrOmd8cwh84Hz9BO5ZV
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckIcxO7iZDZmeXmVkhrPkEW+3tROz8BH/CxsIvcfIoNPHAhcM593Iux48F18Z1P53Myura+kZ2M7e1vbO7l98/qOsoUQxrLBKRavpUo+ASa4Ybgc1YIQ19gQ1/eDXxG3eoNI/krRnF2AlpX/KAM2qsdKO62M0X3KI7BVkm3pwUSkfVL/5W/qh089/tXsSSEKVhgmrd8tzYdFKqDGcCx7l2ojGmbEj72LJU0hB1J52+OianVumRIFJ2pCFT9fdFSkOtR6FvN0NqBnrRm4j/en64kGyCy07KZZwYlGwWHCSCmIhMeiA9rpAZMbKEMsXt74QNqKLM2LZythRvsYJlUj8vem7Rq9p2yjBDFo7hBM7AgwsowTVUoAYM+vAAj/Dk3DvPzovzOlvNOPObQ/gD5/0HrjyYEQ==
AAAB+nicbVC7SgNBFL0bXzG+ooKNzWAQrMKujZYhNpYJmgckIcxO7iZDZmeXmVkhrPkEW+3tROz8BH/CxsIvcfIoNPHAhcM593Iux48F18Z1P53Myura+kZ2M7e1vbO7l98/qOsoUQxrLBKRavpUo+ASa4Ybgc1YIQ19gQ1/eDXxG3eoNI/krRnF2AlpX/KAM2qsdKO62M0X3KI7BVkm3pwUSkfVL/5W/qh089/tXsSSEKVhgmrd8tzYdFKqDGcCx7l2ojGmbEj72LJU0hB1J52+OianVumRIFJ2pCFT9fdFSkOtR6FvN0NqBnrRm4j/en64kGyCy07KZZwYlGwWHCSCmIhMeiA9rpAZMbKEMsXt74QNqKLM2LZythRvsYJlUj8vem7Rq9p2yjBDFo7hBM7AgwsowTVUoAYM+vAAj/Dk3DvPzovzOlvNOPObQ/gD5/0HrjyYEQ==
AAAB+nicbVC7SgNBFL0bXzG+opY2g0GwCrs2WgZtLCOaByRLmJ3cTYbMzC4zs0KI+QRb7e3E1p+x9UucJFto4oELh3Pu5VxOlApurO9/eYW19Y3NreJ2aWd3b/+gfHjUNEmmGTZYIhLdjqhBwRU2LLcC26lGKiOBrWh0M/Nbj6gNT9SDHacYSjpQPOaMWifd6x72yhW/6s9BVkmQkwrkqPfK391+wjKJyjJBjekEfmrDCdWWM4HTUjczmFI2ogPsOKqoRBNO5q9OyZlT+iROtBtlyVz9fTGh0pixjNympHZolr2Z+K8XyaVkG1+FE67SzKJii+A4E8QmZNYD6XONzIqxI5Rp7n4nbEg1Zda1VXKlBMsVrJLmRTXwq8GdX6ld5/UU4QRO4RwCuIQa3EIdGsBgAM/wAq/ek/fmvXsfi9WCl98cwx94nz+URJRw
Өڹܘ
w(r) =
(
re
r
1 (r < re)
0 (r re)
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
ॏΈؔͷఆٛ
ॏΈؔͷ૯ཻ͕ࢠີ
άϥσΟΤϯτϞσϧ
hr ii =
d
n0
X
j6=i
✓
j i
|rj ri
|
rj ri
|rj ri
|
w(|rj ri
|)
◆
AAACqHicbVFdaxQxFM2MX3X9Wis++RIswlZwmfHFghSWVkHfVnC7lWYdkuyd3bRJZppkKks6v8nf46v/QzAzsw+63Qshh3PO5V7OZaUU1iXJryi+dfvO3Xs793sPHj56/KT/dPfEFpXhMOGFLMwpoxak0DBxwkk4LQ1QxSRM2cVxo0+vwFhR6K9uVcJM0YUWueDUBSrr/ySS6oUEoimTlJRLQUxLZOKQ5IZyP6+9/p7UxFYq8+dEwyUWNZGQu0FnaJqy8zftJ2p/TZjypm6YDojrem3cFLabfwy2sPvEiMXS7Wf9vWSYtIVvgnQN9kaf/1yNEELjrP+bzAteKdCOS2rtWZqUbuapcYJLqHukslBSfkEXcBagpgrszLfB1vhVYOY4L0x42uGW/bfDU2XtSrHgVNQt7abWkFs1pjYmu/xg5oUuKwead4PzSmJX4OZqeC4McCdXAVBuRNgd8yUNobpw214IJd2M4CY4eTtMk2H6JaRzhLraQS/QSzRAKXqHRugTGqMJ4tHz6H30IfoYv47H8TT+1lnjaN3zDP1XMfsLcVvZ1Q==
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFda9swFJXdbuuyj2Yde9qLWBikgwV7LxuMQdkH7DGDpumoMiMp14laSfYkuSWo/k37PXvdL5ns+KFNc0HocM653Mu5rJTCuiT5G8U7u/fuP9h72Hv0+MnT/f6zgxNbVIbDhBeyMKeMWpBCw8QJJ+G0NEAVkzBlF18afXoJxopCH7tVCTNFF1rkglMXqKz/h0iqFxKIpkxSUi4FMS2RiU8kN5T7ee31r6QmtlKZPycafmNREwm5G64NTVN2/rb9RO2vCVPe1A2zBuK67oybwnbz1XALe0iMWCzdYdYfJKOkLXwXpB0YoK7GWf8fmRe8UqAdl9TaszQp3cxT4wSXUPdIZaGk/IIu4CxATRXYmW+DrfHrwMxxXpjwtMMte7PDU2XtSrHgVNQt7abWkFs1pjYmu/zDzAtdVg40Xw/OK4ldgZur4bkwwJ1cBUC5EWF3zJc0hOrCbXshlHQzgrvg5N0oTUbpj2Rw9LmLZw+9RK/QEKXoPTpC39EYTRCPXkQfo6/Rt/hNPI6n8c+1NY66nufoVsXsP0iD15k=
άϥσΟΤϯτϞσϧ
hr ii =
d
n0
X
j6=i
✓
j i
|rj ri
|
rj ri
|rj ri
|
w(|rj ri
|)
◆
AAACqHicbVFdaxQxFM2MX3X9Wis++RIswlZwmfHFghSWVkHfVnC7lWYdkuyd3bRJZppkKks6v8nf46v/QzAzsw+63Qshh3PO5V7OZaUU1iXJryi+dfvO3Xs793sPHj56/KT/dPfEFpXhMOGFLMwpoxak0DBxwkk4LQ1QxSRM2cVxo0+vwFhR6K9uVcJM0YUWueDUBSrr/ySS6oUEoimTlJRLQUxLZOKQ5IZyP6+9/p7UxFYq8+dEwyUWNZGQu0FnaJqy8zftJ2p/TZjypm6YDojrem3cFLabfwy2sPvEiMXS7Wf9vWSYtIVvgnQN9kaf/1yNEELjrP+bzAteKdCOS2rtWZqUbuapcYJLqHukslBSfkEXcBagpgrszLfB1vhVYOY4L0x42uGW/bfDU2XtSrHgVNQt7abWkFs1pjYmu/xg5oUuKwead4PzSmJX4OZqeC4McCdXAVBuRNgd8yUNobpw214IJd2M4CY4eTtMk2H6JaRzhLraQS/QSzRAKXqHRugTGqMJ4tHz6H30IfoYv47H8TT+1lnjaN3zDP1XMfsLcVvZ1Q==
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFda9swFJXdbuuyj2Yde9qLWBikgwV7LxuMQdkH7DGDpumoMiMp14laSfYkuSWo/k37PXvdL5ns+KFNc0HocM653Mu5rJTCuiT5G8U7u/fuP9h72Hv0+MnT/f6zgxNbVIbDhBeyMKeMWpBCw8QJJ+G0NEAVkzBlF18afXoJxopCH7tVCTNFF1rkglMXqKz/h0iqFxKIpkxSUi4FMS2RiU8kN5T7ee31r6QmtlKZPycafmNREwm5G64NTVN2/rb9RO2vCVPe1A2zBuK67oybwnbz1XALe0iMWCzdYdYfJKOkLXwXpB0YoK7GWf8fmRe8UqAdl9TaszQp3cxT4wSXUPdIZaGk/IIu4CxATRXYmW+DrfHrwMxxXpjwtMMte7PDU2XtSrHgVNQt7abWkFs1pjYmu/zDzAtdVg40Xw/OK4ldgZur4bkwwJ1cBUC5EWF3zJc0hOrCbXshlHQzgrvg5N0oTUbpj2Rw9LmLZw+9RK/QEKXoPTpC39EYTRCPXkQfo6/Rt/hNPI6n8c+1NY66nufoVsXsP0iD15k=
άϥσΟΤϯτϞσϧ
hr ii =
d
n0
X
j6=i
✓
j i
|rj ri
|
rj ri
|rj ri
|
w(|rj ri
|)
◆
AAACqHicbVFdaxQxFM2MX3X9Wis++RIswlZwmfHFghSWVkHfVnC7lWYdkuyd3bRJZppkKks6v8nf46v/QzAzsw+63Qshh3PO5V7OZaUU1iXJryi+dfvO3Xs793sPHj56/KT/dPfEFpXhMOGFLMwpoxak0DBxwkk4LQ1QxSRM2cVxo0+vwFhR6K9uVcJM0YUWueDUBSrr/ySS6oUEoimTlJRLQUxLZOKQ5IZyP6+9/p7UxFYq8+dEwyUWNZGQu0FnaJqy8zftJ2p/TZjypm6YDojrem3cFLabfwy2sPvEiMXS7Wf9vWSYtIVvgnQN9kaf/1yNEELjrP+bzAteKdCOS2rtWZqUbuapcYJLqHukslBSfkEXcBagpgrszLfB1vhVYOY4L0x42uGW/bfDU2XtSrHgVNQt7abWkFs1pjYmu/xg5oUuKwead4PzSmJX4OZqeC4McCdXAVBuRNgd8yUNobpw214IJd2M4CY4eTtMk2H6JaRzhLraQS/QSzRAKXqHRugTGqMJ4tHz6H30IfoYv47H8TT+1lnjaN3zDP1XMfsLcVvZ1Q==
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFda9swFJXdbuuyj2Yde9qLWBikgwV7LxuMQdkH7DGDpumoMiMp14laSfYkuSWo/k37PXvdL5ns+KFNc0HocM653Mu5rJTCuiT5G8U7u/fuP9h72Hv0+MnT/f6zgxNbVIbDhBeyMKeMWpBCw8QJJ+G0NEAVkzBlF18afXoJxopCH7tVCTNFF1rkglMXqKz/h0iqFxKIpkxSUi4FMS2RiU8kN5T7ee31r6QmtlKZPycafmNREwm5G64NTVN2/rb9RO2vCVPe1A2zBuK67oybwnbz1XALe0iMWCzdYdYfJKOkLXwXpB0YoK7GWf8fmRe8UqAdl9TaszQp3cxT4wSXUPdIZaGk/IIu4CxATRXYmW+DrfHrwMxxXpjwtMMte7PDU2XtSrHgVNQt7abWkFs1pjYmu/zDzAtdVg40Xw/OK4ldgZur4bkwwJ1cBUC5EWF3zJc0hOrCbXshlHQzgrvg5N0oTUbpj2Rw9LmLZw+9RK/QEKXoPTpC39EYTRCPXkQfo6/Rt/hNPI6n8c+1NY66nufoVsXsP0iD15k=
άϥσΟΤϯτϞσϧ
hr ii =
d
n0
X
j6=i
✓
j i
|rj ri
|
rj ri
|rj ri
|
w(|rj ri
|)
◆
AAACqHicbVFdaxQxFM2MX3X9Wis++RIswlZwmfHFghSWVkHfVnC7lWYdkuyd3bRJZppkKks6v8nf46v/QzAzsw+63Qshh3PO5V7OZaUU1iXJryi+dfvO3Xs793sPHj56/KT/dPfEFpXhMOGFLMwpoxak0DBxwkk4LQ1QxSRM2cVxo0+vwFhR6K9uVcJM0YUWueDUBSrr/ySS6oUEoimTlJRLQUxLZOKQ5IZyP6+9/p7UxFYq8+dEwyUWNZGQu0FnaJqy8zftJ2p/TZjypm6YDojrem3cFLabfwy2sPvEiMXS7Wf9vWSYtIVvgnQN9kaf/1yNEELjrP+bzAteKdCOS2rtWZqUbuapcYJLqHukslBSfkEXcBagpgrszLfB1vhVYOY4L0x42uGW/bfDU2XtSrHgVNQt7abWkFs1pjYmu/xg5oUuKwead4PzSmJX4OZqeC4McCdXAVBuRNgd8yUNobpw214IJd2M4CY4eTtMk2H6JaRzhLraQS/QSzRAKXqHRugTGqMJ4tHz6H30IfoYv47H8TT+1lnjaN3zDP1XMfsLcVvZ1Q==
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFda9swFJXdbuuyj2Yde9qLWBikgwV7LxuMQdkH7DGDpumoMiMp14laSfYkuSWo/k37PXvdL5ns+KFNc0HocM653Mu5rJTCuiT5G8U7u/fuP9h72Hv0+MnT/f6zgxNbVIbDhBeyMKeMWpBCw8QJJ+G0NEAVkzBlF18afXoJxopCH7tVCTNFF1rkglMXqKz/h0iqFxKIpkxSUi4FMS2RiU8kN5T7ee31r6QmtlKZPycafmNREwm5G64NTVN2/rb9RO2vCVPe1A2zBuK67oybwnbz1XALe0iMWCzdYdYfJKOkLXwXpB0YoK7GWf8fmRe8UqAdl9TaszQp3cxT4wSXUPdIZaGk/IIu4CxATRXYmW+DrfHrwMxxXpjwtMMte7PDU2XtSrHgVNQt7abWkFs1pjYmu/zDzAtdVg40Xw/OK4ldgZur4bkwwJ1cBUC5EWF3zJc0hOrCbXshlHQzgrvg5N0oTUbpj2Rw9LmLZw+9RK/QEKXoPTpC39EYTRCPXkQfo6/Rt/hNPI6n8c+1NY66nufoVsXsP0iD15k=
άϥσΟΤϯτϞσϧ
hr ii =
d
n0
X
j6=i
✓
j i
|rj ri
|
rj ri
|rj ri
|
w(|rj ri
|)
◆
AAACqHicbVFdaxQxFM2MX3X9Wis++RIswlZwmfHFghSWVkHfVnC7lWYdkuyd3bRJZppkKks6v8nf46v/QzAzsw+63Qshh3PO5V7OZaUU1iXJryi+dfvO3Xs793sPHj56/KT/dPfEFpXhMOGFLMwpoxak0DBxwkk4LQ1QxSRM2cVxo0+vwFhR6K9uVcJM0YUWueDUBSrr/ySS6oUEoimTlJRLQUxLZOKQ5IZyP6+9/p7UxFYq8+dEwyUWNZGQu0FnaJqy8zftJ2p/TZjypm6YDojrem3cFLabfwy2sPvEiMXS7Wf9vWSYtIVvgnQN9kaf/1yNEELjrP+bzAteKdCOS2rtWZqUbuapcYJLqHukslBSfkEXcBagpgrszLfB1vhVYOY4L0x42uGW/bfDU2XtSrHgVNQt7abWkFs1pjYmu/xg5oUuKwead4PzSmJX4OZqeC4McCdXAVBuRNgd8yUNobpw214IJd2M4CY4eTtMk2H6JaRzhLraQS/QSzRAKXqHRugTGqMJ4tHz6H30IfoYv47H8TT+1lnjaN3zDP1XMfsLcVvZ1Q==
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFda9swFJXdbuuyj2Yde9qLWBikgwV7LxuMQdkH7DGDpumoMiMp14laSfYkuSWo/k37PXvdL5ns+KFNc0HocM653Mu5rJTCuiT5G8U7u/fuP9h72Hv0+MnT/f6zgxNbVIbDhBeyMKeMWpBCw8QJJ+G0NEAVkzBlF18afXoJxopCH7tVCTNFF1rkglMXqKz/h0iqFxKIpkxSUi4FMS2RiU8kN5T7ee31r6QmtlKZPycafmNREwm5G64NTVN2/rb9RO2vCVPe1A2zBuK67oybwnbz1XALe0iMWCzdYdYfJKOkLXwXpB0YoK7GWf8fmRe8UqAdl9TaszQp3cxT4wSXUPdIZaGk/IIu4CxATRXYmW+DrfHrwMxxXpjwtMMte7PDU2XtSrHgVNQt7abWkFs1pjYmu/zDzAtdVg40Xw/OK4ldgZur4bkwwJ1cBUC5EWF3zJc0hOrCbXshlHQzgrvg5N0oTUbpj2Rw9LmLZw+9RK/QEKXoPTpC39EYTRCPXkQfo6/Rt/hNPI6n8c+1NY66nufoVsXsP0iD15k=
ཻࢠJͱཻۙࢠKͷؒͷޯϕΫτϧ
άϥσΟΤϯτϞσϧ
hr ii =
d
n0
X
j6=i
✓
j i
|rj ri
|
rj ri
|rj ri
|
w(|rj ri
|)
◆
AAACqHicbVFdaxQxFM2MX3X9Wis++RIswlZwmfHFghSWVkHfVnC7lWYdkuyd3bRJZppkKks6v8nf46v/QzAzsw+63Qshh3PO5V7OZaUU1iXJryi+dfvO3Xs793sPHj56/KT/dPfEFpXhMOGFLMwpoxak0DBxwkk4LQ1QxSRM2cVxo0+vwFhR6K9uVcJM0YUWueDUBSrr/ySS6oUEoimTlJRLQUxLZOKQ5IZyP6+9/p7UxFYq8+dEwyUWNZGQu0FnaJqy8zftJ2p/TZjypm6YDojrem3cFLabfwy2sPvEiMXS7Wf9vWSYtIVvgnQN9kaf/1yNEELjrP+bzAteKdCOS2rtWZqUbuapcYJLqHukslBSfkEXcBagpgrszLfB1vhVYOY4L0x42uGW/bfDU2XtSrHgVNQt7abWkFs1pjYmu/xg5oUuKwead4PzSmJX4OZqeC4McCdXAVBuRNgd8yUNobpw214IJd2M4CY4eTtMk2H6JaRzhLraQS/QSzRAKXqHRugTGqMJ4tHz6H30IfoYv47H8TT+1lnjaN3zDP1XMfsLcVvZ1Q==
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFda9swFJXdbuuyj2Yde9qLWBikgwV7LxuMQdkH7DGDpumoMiMp14laSfYkuSWo/k37PXvdL5ns+KFNc0HocM653Mu5rJTCuiT5G8U7u/fuP9h72Hv0+MnT/f6zgxNbVIbDhBeyMKeMWpBCw8QJJ+G0NEAVkzBlF18afXoJxopCH7tVCTNFF1rkglMXqKz/h0iqFxKIpkxSUi4FMS2RiU8kN5T7ee31r6QmtlKZPycafmNREwm5G64NTVN2/rb9RO2vCVPe1A2zBuK67oybwnbz1XALe0iMWCzdYdYfJKOkLXwXpB0YoK7GWf8fmRe8UqAdl9TaszQp3cxT4wSXUPdIZaGk/IIu4CxATRXYmW+DrfHrwMxxXpjwtMMte7PDU2XtSrHgVNQt7abWkFs1pjYmu/zDzAtdVg40Xw/OK4ldgZur4bkwwJ1cBUC5EWF3zJc0hOrCbXshlHQzgrvg5N0oTUbpj2Rw9LmLZw+9RK/QEKXoPTpC39EYTRCPXkQfo6/Rt/hNPI6n8c+1NY66nufoVsXsP0iD15k=
ཻࢠJͱཻۙࢠKͷؒͷޯϕΫτϧ
ޯϕΫτϧͷॏΈฏۉ
άϥσΟΤϯτϞσϧ
hr ii =
d
n0
X
j6=i
✓
j i
|rj ri
|
rj ri
|rj ri
|
w(|rj ri
|)
◆
AAACqHicbVFdaxQxFM2MX3X9Wis++RIswlZwmfHFghSWVkHfVnC7lWYdkuyd3bRJZppkKks6v8nf46v/QzAzsw+63Qshh3PO5V7OZaUU1iXJryi+dfvO3Xs793sPHj56/KT/dPfEFpXhMOGFLMwpoxak0DBxwkk4LQ1QxSRM2cVxo0+vwFhR6K9uVcJM0YUWueDUBSrr/ySS6oUEoimTlJRLQUxLZOKQ5IZyP6+9/p7UxFYq8+dEwyUWNZGQu0FnaJqy8zftJ2p/TZjypm6YDojrem3cFLabfwy2sPvEiMXS7Wf9vWSYtIVvgnQN9kaf/1yNEELjrP+bzAteKdCOS2rtWZqUbuapcYJLqHukslBSfkEXcBagpgrszLfB1vhVYOY4L0x42uGW/bfDU2XtSrHgVNQt7abWkFs1pjYmu/xg5oUuKwead4PzSmJX4OZqeC4McCdXAVBuRNgd8yUNobpw214IJd2M4CY4eTtMk2H6JaRzhLraQS/QSzRAKXqHRugTGqMJ4tHz6H30IfoYv47H8TT+1lnjaN3zDP1XMfsLcVvZ1Q==
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFda9swFJXdbuuyj2Yde9qLWBikgwV7LxuMQdkH7DGDpumoMiMp14laSfYkuSWo/k37PXvdL5ns+KFNc0HocM653Mu5rJTCuiT5G8U7u/fuP9h72Hv0+MnT/f6zgxNbVIbDhBeyMKeMWpBCw8QJJ+G0NEAVkzBlF18afXoJxopCH7tVCTNFF1rkglMXqKz/h0iqFxKIpkxSUi4FMS2RiU8kN5T7ee31r6QmtlKZPycafmNREwm5G64NTVN2/rb9RO2vCVPe1A2zBuK67oybwnbz1XALe0iMWCzdYdYfJKOkLXwXpB0YoK7GWf8fmRe8UqAdl9TaszQp3cxT4wSXUPdIZaGk/IIu4CxATRXYmW+DrfHrwMxxXpjwtMMte7PDU2XtSrHgVNQt7abWkFs1pjYmu/zDzAtdVg40Xw/OK4ldgZur4bkwwJ1cBUC5EWF3zJc0hOrCbXshlHQzgrvg5N0oTUbpj2Rw9LmLZw+9RK/QEKXoPTpC39EYTRCPXkQfo6/Rt/hNPI6n8c+1NY66nufoVsXsP0iD15k=
ॏΈฏۉͷن֨ԽΛ͍ͯ͠Δ
άϥσΟΤϯτϞσϧ
hr ii =
d
n0
X
j6=i
✓
j i
|rj ri
|
rj ri
|rj ri
|
w(|rj ri
|)
◆
AAACqHicbVFdaxQxFM2MX3X9Wis++RIswlZwmfHFghSWVkHfVnC7lWYdkuyd3bRJZppkKks6v8nf46v/QzAzsw+63Qshh3PO5V7OZaUU1iXJryi+dfvO3Xs793sPHj56/KT/dPfEFpXhMOGFLMwpoxak0DBxwkk4LQ1QxSRM2cVxo0+vwFhR6K9uVcJM0YUWueDUBSrr/ySS6oUEoimTlJRLQUxLZOKQ5IZyP6+9/p7UxFYq8+dEwyUWNZGQu0FnaJqy8zftJ2p/TZjypm6YDojrem3cFLabfwy2sPvEiMXS7Wf9vWSYtIVvgnQN9kaf/1yNEELjrP+bzAteKdCOS2rtWZqUbuapcYJLqHukslBSfkEXcBagpgrszLfB1vhVYOY4L0x42uGW/bfDU2XtSrHgVNQt7abWkFs1pjYmu/xg5oUuKwead4PzSmJX4OZqeC4McCdXAVBuRNgd8yUNobpw214IJd2M4CY4eTtMk2H6JaRzhLraQS/QSzRAKXqHRugTGqMJ4tHz6H30IfoYv47H8TT+1lnjaN3zDP1XMfsLcVvZ1Q==
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFdixMxFM2MX2v92K7iky/BRegKlhlfFBZl8QPUpwp2u7KpQ5LeadNNMmOSWSnZ+U3+A/+Hr/4PwcxMH7TbCyGHc87lXs5lpRTWJcmvKL5y9dr1Gzs3e7du37m729+7d2yLynAY80IW5oRRC1JoGDvhJJyUBqhiEibs7E2jT87BWFHoz25VwlTRuRa54NQFKuv/IJLquQSiKZOUlAtBTEtk4iXJDeV+Vnv9NamJrVTml0TDNyxqIiF3g87QNGXLp+0nan9BmPKmbpgOiIt6bdwUtpu/D7awB8SI+cIdZP39ZJi0hS+DdA32jz78OX/18edylPV/k1nBKwXacUmtPU2T0k09NU5wCXWPVBZKys/oHE4D1FSBnfo22Bo/DswM54UJTzvcsv92eKqsXSkWnIq6hd3UGnKrxtTGZJe/mHqhy8qB5t3gvJLYFbi5Gp4JA9zJVQCUGxF2x3xBQ6gu3LYXQkk3I7gMjp8N02SYfgrpvEZd7aCH6BEaoBQ9R0foPRqhMeLRg+gwehu9i5/Eo3gSf+mscbTuuY/+q5j9BXP120U=
AAACqHicbVFda9swFJXdbuuyj2Yde9qLWBikgwV7LxuMQdkH7DGDpumoMiMp14laSfYkuSWo/k37PXvdL5ns+KFNc0HocM653Mu5rJTCuiT5G8U7u/fuP9h72Hv0+MnT/f6zgxNbVIbDhBeyMKeMWpBCw8QJJ+G0NEAVkzBlF18afXoJxopCH7tVCTNFF1rkglMXqKz/h0iqFxKIpkxSUi4FMS2RiU8kN5T7ee31r6QmtlKZPycafmNREwm5G64NTVN2/rb9RO2vCVPe1A2zBuK67oybwnbz1XALe0iMWCzdYdYfJKOkLXwXpB0YoK7GWf8fmRe8UqAdl9TaszQp3cxT4wSXUPdIZaGk/IIu4CxATRXYmW+DrfHrwMxxXpjwtMMte7PDU2XtSrHgVNQt7abWkFs1pjYmu/zDzAtdVg40Xw/OK4ldgZur4bkwwJ1cBUC5EWF3zJc0hOrCbXshlHQzgrvg5N0oTUbpj2Rw9LmLZw+9RK/QEKXoPTpC39EYTRCPXkQfo6/Rt/hNPI6n8c+1NY66nufoVsXsP0iD15k=
Eʢ࣍ݩʣΛ͔͚ͯਖ਼͍͠ʹ
ϥϓϥγΞϯϞσϧ
hr2 ii =
2d
0n0
X
j6=i
( jw(|rj ri
|) iw(|rj ri
|))
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
ϥϓϥγΞϯϞσϧ
hr2 ii =
2d
0n0
X
j6=i
( jw(|rj ri
|) iw(|rj ri
|))
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
ϥϓϥγΞϯϞσϧ
hr2 ii =
2d
0n0
X
j6=i
( jw(|rj ri
|) iw(|rj ri
|))
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
ཻۙࢠʹПΛΔ
ϥϓϥγΞϯϞσϧ
hr2 ii =
2d
0n0
X
j6=i
( jw(|rj ri
|) iw(|rj ri
|))
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
ٯʹΒ͏
ϥϓϥγΞϯϞσϧ
hr2 ii =
2d
0n0
X
j6=i
( jw(|rj ri
|) iw(|rj ri
|))
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
AAACgXicbVFNaxsxFNRu+pG4H3Hb3EqLqCnYhxitKbRQCiG99JhCnQQie5G0b205knYraVvMZu/9i73mB/Q3VN7soXUyIBhm5vHEPF4q6Twhv6N45979Bw9393qPHj95ut9/9vzUFZUVMBWFKuw5Zw6UNDD10is4Ly0wzRWc8cvPG//sB1gnC/PNr0uYabYwMpeC+SCl/V9UMbNQQA3jis0ntFxKalsplZ9obpmoJ1lTh5jmGZsTMycNdZVO6xU18B3LhirI/XAzmK5+Dq8o17Vt0tVhR+TV6LA15d0mtXKx9KO0PyBj0gLfJklHBqjDSdq/plkhKg3GC8Wcu0hI6Wc1s14KBU2PVg5KJi7ZAi4CNUyDm9VtYw1+G5QM54UNz3jcqv9O1Ew7t9Y8JDXzS7ftbcQ7Pa63Nvv8w6yWpqw8GHGzOK8U9gXenANn0oLwah0IE1aGv2OxZKF1H47WC6Uk2xXcJqeTcULGydd3g6Pjrp5d9BK9QUOUoPfoCH1BJ2iKBPoTHUSvotfxTjyKSTy5icZRN/MC/Yf441+i8sVT
֦ࢄํఔࣜʹద༻͢ΔͳͲͰٻΊΒΕΔ
ఆܭࢉ
n0 =
X
j6=i0
w |r0
j
r0
i0
|
AAACNHicbZBNSxxBEIZrzJfZfG3iMZdGCTGHLD1eIoTAopfkZiCrgrM79PTW7LZ294zdNcoy7o/xT+QveE2OAU+GXP0Hgr0fh2RNQcPD+1ZR1W9WauWJ81/R0r37Dx4+Wn7cePL02fMXzZevdn1ROYkdWejC7WfCo1YWO6RI437pUJhM4152tD3x907QeVXYbzQqsWvEwKpcSUFBSpsfbY9/Snxl0vowsXjM1NvxaaIxp/WzJDO1G6eHPf5+jnVwe/wscWowpHdpc423+LTYXYjnsNb+cnPSBoCdtHmV9AtZGbQktfD+IOYldWvhSEmN40ZSeSyFPBIDPAhohUHfraefHLM3QemzvHDhWWJT9e+JWhjvRyYLnUbQ0C96E/G/XmYWNlO+2a2VLStCK2eL80ozKtgkQdZXDiXpUQAhnQq3MzkUTkgKOTdCKPFiBHdhd6MV81b8NaSzBbNahtewCusQwwdow2fYgQ5IOIcL+AE/o+/RZfQ7+jNrXYrmMyvwT0XXt0tcrjY=
AAACNHicbZBNTxsxEIa9tLQ0/UrLsRerUdVwaOTtBaSqFSqXwgmkJiBlk5XXmU0cbO/WnqWKlvyY/okee+XaHpE4gbjyDyrhfBwgMJKlR+87oxm/Sa6kQ8ZOg6UHD5cfPV55Unn67PmLl9VXr1suK6yApshUZg8S7kBJA02UqOAgt8B1omA/Odya+PtHYJ3MzHcc5dDRvG9kKgVHL8XVT6bLPkeu0HE5jAz8oPL9+GekIMX6cZTo0o7jYZd9mGPp3S47jqzsD3AtrtZYg02L3oVwDrXN7f9HX3b+DHfj6nnUy0ShwaBQ3Ll2yHLslNyiFArGlahwkHNxyPvQ9mi4Btcpp58c03de6dE0s/4ZpFP15kTJtXMjnfhOzXHgFr2JeK+X6IXNmG50SmnyAsGI2eK0UBQzOkmQ9qQFgWrkgQsr/e1UDLjlAn3OFR9KuBjBXWh9bISsEe75dL6SWa2QN+QtqZOQrJNN8o3skiYR5Bc5IX/Jv+B3cBZcBJez1qVgPrNKblVwdQ1N9q+m
AAACNHicbZBNTxsxEIa9tLQ0/UrLsRerUdVwaOTtBaSqFSqXwgmkJiBlk5XXmU0cbO/WnqWKlvyY/okee+XaHpE4gbjyDyrhfBwgMJKlR+87oxm/Sa6kQ8ZOg6UHD5cfPV55Unn67PmLl9VXr1suK6yApshUZg8S7kBJA02UqOAgt8B1omA/Odya+PtHYJ3MzHcc5dDRvG9kKgVHL8XVT6bLPkeu0HE5jAz8oPL9+GekIMX6cZTo0o7jYZd9mGPp3S47jqzsD3AtrtZYg02L3oVwDrXN7f9HX3b+DHfj6nnUy0ShwaBQ3Ll2yHLslNyiFArGlahwkHNxyPvQ9mi4Btcpp58c03de6dE0s/4ZpFP15kTJtXMjnfhOzXHgFr2JeK+X6IXNmG50SmnyAsGI2eK0UBQzOkmQ9qQFgWrkgQsr/e1UDLjlAn3OFR9KuBjBXWh9bISsEe75dL6SWa2QN+QtqZOQrJNN8o3skiYR5Bc5IX/Jv+B3cBZcBJez1qVgPrNKblVwdQ1N9q+m
AAACNHicbZDPThsxEMa9UCiEP03psReLCAEHIm8vVEKVEFx6pFIDSGyy8jqzicH2LvYsKFryMLwEr8AVjkicqLjyDDhhD23oSJZ++r4ZzfhLciUdMvYQTE1/mJn9ODdfW1hcWv5U/7xy6LLCCmiJTGX2OOEOlDTQQokKjnMLXCcKjpKz/ZF/dAHWycz8xkEObc17RqZScPRSXN8xHfYjcoWOy9PIwDmV68PLSEGKG1dRoks7jE87bKvC0rsddhVZ2evjZlxvsCYbF30PYQUNUtVBXH+KupkoNBgUijt3ErIc2yW3KIWCYS0qHORcnPEenHg0XINrl+NPDumaV7o0zax/BulY/Xui5Nq5gU58p+bYd5PeSPyvl+iJzZh+b5fS5AWCEW+L00JRzOgoQdqVFgSqgQcurPS3U9Hnlgv0Odd8KOFkBO/h8FszZM3wF2vs7lXxzJGvZJVskJBsk13ykxyQFhHkmtySO3If3ASPwZ/g+a11KqhmvpB/Knh5BSKEq/o=
ཻࢠີͷج४
Өڹܘʹ͋Δཻۙࢠͱͷڑͷ̎ͷॏΈฏۉ
0 =
P
j6=i0
|r0
j
r0
i0
|2w |r0
j
r0
i0
|
P
j6=i0
w |r0
j
r0
i0
|
AAAClnicjVFNbxMxEPVu+SjhowEuSL1YVIj2QOQNB7iAqqKqPQaJtJXqZOV1ZhO3tnexZ0HRdv8EdyT+Vq/9JXU+Dm0KEiNZenrvzYz1Jiu18sjYZRSv3bv/4OH6o9bjJ0+fbbSfvzjyReUk9GWhC3eSCQ9aWeijQg0npQNhMg3H2fmXmX78A5xXhf2G0xIGRoytypUUGKi0/ZvrYB6JIfvEcydkzX1l0vqMW/hO1dvmgmemdk16NmTvlrAO9JBdDLs/uYYct/9t4U6NJ7jTrA7978a0vcU6bF70LkiWYGt3f+/XH0JIL21f8VEhKwMWpRbenyasxEEtHCqpoWnxykMp5LkYw2mAVhjwg3oeY0PfBGZE88KFZ5HO2ZsdtTDeT00WnEbgxK9qM/KvWmZWNmP+cVArW1YIVi4W55WmWNDZjehIOZCopwEI6VT4O5UTEa6D4ZKtEEqyGsFdcNTtJKyTfA3p7JFFrZNN8ppsk4R8ILvkkPRIn8hoLdqJutH7+FX8Od6PDxbWOFr2vCS3Ku5dAw9Xztg=
AAAClnicjVFNbxMxEPVuKZTwlcIFiYtFBbQHIm84wIWqKqrgmArSVqqTldeZTdza3q09SxVt909w58Sv4sovwfk4lBQkRrL09N6bGetNVmrlkbGfUbx2a/32nY27rXv3Hzx81N58fOSLyknoy0IX7iQTHrSy0EeFGk5KB8JkGo6z8w8z/fgrOK8K+wWnJQyMGFuVKykwUGn7O9fBPBJD9p7nTsia+8qk9Rm3cEHVq+aKZ6Z2TXo2ZK+XsA70kF0Nu5dcQ47b/7Zwp8YT3GlWh/53Y9reYh02L3oTJEuwtXew/+3H+ueXvbT9i48KWRmwKLXw/jRhJQ5q4VBJDU2LVx5KIc/FGE4DtMKAH9TzGBv6IjAjmhcuPIt0zl7vqIXxfmqy4DQCJ35Vm5F/1TKzshnzd4Na2bJCsHKxOK80xYLObkRHyoFEPQ1ASKfC36mciHAdDJdshVCS1QhugqNuJ2Gd5DCks08WtUGekedkmyTkLdkjn0iP9ImM1qKdqBu9iZ/Gu/FB/HFhjaNlzxPyR8W938Oiz1k=
AAAClnicjVFNbxMxEPVuKZTwlcIFiYtFBbQHIm84wIWqKqrgmArSVqqTldeZTdza3q09SxVt909w58Sv4sovwfk4lBQkRrL09N6bGetNVmrlkbGfUbx2a/32nY27rXv3Hzx81N58fOSLyknoy0IX7iQTHrSy0EeFGk5KB8JkGo6z8w8z/fgrOK8K+wWnJQyMGFuVKykwUGn7O9fBPBJD9p7nTsia+8qk9Rm3cEHVq+aKZ6Z2TXo2ZK+XsA70kF0Nu5dcQ47b/7Zwp8YT3GlWh/53Y9reYh02L3oTJEuwtXew/+3H+ueXvbT9i48KWRmwKLXw/jRhJQ5q4VBJDU2LVx5KIc/FGE4DtMKAH9TzGBv6IjAjmhcuPIt0zl7vqIXxfmqy4DQCJ35Vm5F/1TKzshnzd4Na2bJCsHKxOK80xYLObkRHyoFEPQ1ASKfC36mciHAdDJdshVCS1QhugqNuJ2Gd5DCks08WtUGekedkmyTkLdkjn0iP9ImM1qKdqBu9iZ/Gu/FB/HFhjaNlzxPyR8W938Oiz1k=
AAAClnicjVFNbxMxEPVuC5Tw0VAuSFwsIkR7IPKGA1yoKlBVjqlE2kp1svI6s4lb27vYs6Bouz+FH8a1v6ROugdIQWIkS0/vvZmx3mSlVh4Z+xXFG5v37j/Yeth59PjJ0+3us50TX1ROwkgWunBnmfCglYURKtRwVjoQJtNwml1+Xuqn38F5VdivuChhbMTMqlxJgYFKuz+5DuapmLCPPHdC1txXJq0vuIVvVL1prnhmatekFxP2toV1oCfsajL4wTXkuPtvC3dqNse9Zn3ofzem3R7rs1XRuyBpQY+0NUy713xayMqARamF9+cJK3FcC4dKamg6vPJQCnkpZnAeoBUG/LhexdjQ14GZ0rxw4VmkK/b3jloY7xcmC04jcO7XtSX5Vy0za5sx/zCulS0rBCtvF+eVpljQ5Y3oVDmQqBcBCOlU+DuVcxGug+GSnRBKsh7BXXAy6Cesnxyz3sGnNp4t8pK8IrskIe/JAflChmREZLQR7UWD6F38It6PD+OjW2sctT3PyR8VD28ARh3M3w==
γϛϡϨʔγϣϯຊମ
ԾΛܭࢉ͢Δ
ԾΛͬͯҠಈ͢Δ
ඇѹॖ݅Λຬͨ͢Α͏ʹ
ѹྗΛܭࢉ
ҐஔͱΛ࠶ܭࢉ
void calGravity(void);
void calViscosity(void);
void moveParticle(void);
void collision(void);
void calPressure(void);
void calPressureGradient(void);
void moveParticleUsingPressureGradient(void);
γϛϡϨʔγϣϯຊମ
ԾΛܭࢉ͢Δ
ԾΛͬͯҠಈ͢Δ
ඇѹॖ݅Λຬͨ͢Α͏ʹ
ѹྗΛܭࢉ
ҐஔͱΛ࠶ܭࢉ
void calGravity(void);
void calViscosity(void);
void moveParticle(void);
void collision(void);
void calPressure(void);
void calPressureGradient(void);
void moveParticleUsingPressureGradient(void);
೪ੑ߲ void calViscosity(void);
⌫r2u ! ⌫
2d
0n0
X
j6=i
(uk
j
uk
i
)w(|rk
j
rk
i
|)
AAACbnicbVFdT9swFHXCBqx8lU3aC5qwQJPaB6qkL/BYjZfxxiQKSKSNbhyndWs7wXZaVSE/FPHGL+APIM1p0T5armzp+Jx7fK3jKONMG897dNy1Dx/XNzY/1ba2d3b36vufr3WaK0K7JOWpuo1AU84k7RpmOL3NFAURcXoTjc8r/WZClWapvDKzjPYEDCRLGAFjqbCeB9IuiDj020EkirwMFBsMDSiVTistUUCKdlwWAbe3xtD3ZN8rA52LsBgFkt5jVjYWzv44HJ38gaw5bTxUJ/VXmEP20Azrx17LmxdeBf4bOO5cvE46CKHLsP4cxCnJBZWGcND6zvcy0ytAGUY4LWtBrmkGZAwDemehBEF1r5jHU+Lvlolxkiq7pcFz9l9HAULrmYhspwAz1MtaRb6rRWJpsknOegWTWW6oJIvBSc6xSXGVPY6ZosTwmQVAFLNvx2QINmBjf6hmQ/GXI1gF1+2W77X8XzadH2hRm+gAHaEG8tEp6qCf6BJ1EUFPjuPUnC3nxf3qfnMPF62u8+b5gv4rt/EbUULBMg==
AAACbnicbVHdTtswGHWyDVj5K0PiBk2zQJPKBVXSG7hCCG4GV520AhJpoy+O07q1nWA7oCrNq/Ee0+72BHsBJJwW8VP2yZaOz/mOP+s4yjjTxvN+O+6Hj58WFpc+15ZXVtfW6xtfLnSaK0I7JOWpuopAU84k7RhmOL3KFAURcXoZjU4r/fKWKs1S+cuMM9oV0JcsYQSMpcJ6Hki7IOLQawWRKPIyUKw/MKBUeldpiQJStOKyCLi9NYaeJ3teGehchMUwkPQGs7Ixc/ZG4XD/GbK9u8akOqkXYQrZZC+s73pNb1r4PfCfwO7x2cPt0fn9sB3W/wZxSnJBpSEctL72vcx0C1CGEU7LWpBrmgEZQZ9eWyhBUN0tpvGU+LtlYpykym5p8JR97ShAaD0Wke0UYAZ6XqvI/2qRmJtsksNuwWSWGyrJbHCSc2xSXGWPY6YoMXxsARDF7NsxGYAN2NgfqtlQ/PkI3oOLVtP3mv5Pm84JmtUS2kY7qIF8dICO0Q/URh1E0B/HcWrOsvPP3XK/ut9mra7z5NlEb8ptPAJT3MKi
AAACbnicbVHdTtswGHWyDVj5K0PiBk2zQJPKBVXSG7hCCG4GV520AhJpoy+O07q1nWA7oCrNq/Ee0+72BHsBJJwW8VP2yZaOz/mOP+s4yjjTxvN+O+6Hj58WFpc+15ZXVtfW6xtfLnSaK0I7JOWpuopAU84k7RhmOL3KFAURcXoZjU4r/fKWKs1S+cuMM9oV0JcsYQSMpcJ6Hki7IOLQawWRKPIyUKw/MKBUeldpiQJStOKyCLi9NYaeJ3teGehchMUwkPQGs7Ixc/ZG4XD/GbK9u8akOqkXYQrZZC+s73pNb1r4PfCfwO7x2cPt0fn9sB3W/wZxSnJBpSEctL72vcx0C1CGEU7LWpBrmgEZQZ9eWyhBUN0tpvGU+LtlYpykym5p8JR97ShAaD0Wke0UYAZ6XqvI/2qRmJtsksNuwWSWGyrJbHCSc2xSXGWPY6YoMXxsARDF7NsxGYAN2NgfqtlQ/PkI3oOLVtP3mv5Pm84JmtUS2kY7qIF8dICO0Q/URh1E0B/HcWrOsvPP3XK/ut9mra7z5NlEb8ptPAJT3MKi
AAACbnicbVFdT9swFHXCYKV8FZB4maZZVEjlgcrpCzwieNkjk1ZAIm104zitqe1ktgOq0vxQxNt+AT9hTlvto+zKlo7PucfXOo5zwY0l5MXz1z6sb3xsbDa3tnd291r7B7cmKzRlfZqJTN/HYJjgivUtt4Ld55qBjAW7iyfXtX73xLThmfpupzkbSBgpnnIK1lFRqwiVWxALGPbCWJZFFWo+GlvQOnuutVQDLXtJVYbC3ZrAkKghqUJTyKh8DBX7gXnVWTiHk+jx7Dfkp8+dWX3Sf4Q55LPTqNUmXTIv/B4ES9BGy7qJWj/DJKOFZMpSAcY8BCS3gxK05VSwqhkWhuVAJzBiDw4qkMwMynk8FT5xTILTTLutLJ6zfztKkMZMZew6JdixWdVq8r9aLFcm2/RiUHKVF5YpuhicFgLbDNfZ44RrRq2YOgBUc/d2TMfgArbuh5oulGA1gvfgttcNSDf4RtqXV8t4GugTOkYdFKBzdIm+ohvURxS9ep7X9La8N//I/+x/WbT63tJziP4pv/MLKGq+9g==
ɾ೪ੑ߲ΛϥϓϥγΞϯϞσϧͰม
ɾ࣮ࡍʹܭࢉ͢ΔͷϕΫτϧͷ͝ͱ
ɾ࠷ޙʹͦΕΛՃʹͯ͋͛͠Δ
γϛϡϨʔγϣϯຊମ
ԾΛܭࢉ͢Δ
ԾΛͬͯҠಈ͢Δ
ඇѹॖ݅Λຬͨ͢Α͏ʹ
ѹྗΛܭࢉ
ҐஔͱΛ࠶ܭࢉ
void calGravity(void);
void calViscosity(void);
void moveParticle(void);
void collision(void);
void calPressure(void);
void calPressureGradient(void);
void moveParticleUsingPressureGradient(void);
ѹྗܭࢉ void calPressure(void);
void calPressure(void) {
calNumberDensity();
setBoundaryCondition();
setSourceTerm();
setMatrix();
solveSimultaniousEquationsByGaussEliminationMethod();
removeNegativePressure();
setMinimumPressure();
}
ѹྗܭࢉ void calPressure(void);
void calPressure(void) {
calNumberDensity();
setBoundaryCondition();
setSourceTerm();
setMatrix();
solveSimultaniousEquationsByGaussEliminationMethod();
removeNegativePressure();
setMinimumPressure();
}
w(r) =
(
re
r
1 (r < re)
0 (r re)
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
AAACUHicbVHLjtMwFL0pzIPMq8CSjUXFqF1MlcyGWTDSCDYsB4m2I9VR5bg3qTWOE+wbRlXUT+Mn2LGFLezZgftAgg5XsnzuOfdhHaeVVo6i6EvQevBwZ3dv/1F4cHh0fNJ+/GToytpKHMhSl/YmFQ61MjggRRpvKouiSDWO0ts3S330Ea1TpXlP8wqTQuRGZUoK8tSkPeLju67tXfIUc2Ua6Ue5BdeYUZdnVsjGTnDR2AW3Kp9R7yw+7dpXnutxHkYe8xw/sGUecjTTTX/Ik0m7E/WjVbD7IN6ADmzietL+yqelrAs0JLVwbhxHFSWNsKSkRj+ydlgJeStyHHtoRIEuaVYGLNgLz0xZVlp/DLEV+3dHIwrn5kXqKwtBM7etLcn/ammxtZmyi6RRpqoJjVwvzmrNqGRLd9lUWZSk5x4IaZV/O5Mz4X0k/wehNyXetuA+GJ7346gfvzvvXL3e2LMPz+A5dCGGl3AFb+EaBiDhE3yD7/Aj+Bz8DH61gnXpnxuewj/RCn8DX4OzmQ==
ѹྗܭࢉ void calPressure(void);
void calPressure(void) {
calNumberDensity();
setBoundaryCondition();
setSourceTerm();
setMatrix();
solveSimultaniousEquationsByGaussEliminationMethod();
removeNegativePressure();
setMinimumPressure();
}
ni < n0
AAACBHicbVC7SgNBFL0bXzG+opZaDAbBKuzaaGERtLFMwDwgWcPsZJIMmZ1dZ+4KYUnrL9hqbyeClf9hJ36Jk0ehiQcuHM65l3M5QSyFQdf9dDJLyyura9n13Mbm1vZOfnevZqJEM15lkYx0I6CGS6F4FQVK3og1p2EgeT0YXI39+j3XRkTqBocx90PaU6IrGEUr+aotLloBR0rUrdvOF9yiOwFZJN6MFEqH75UvACi389+tTsSSkCtkkhrT9NwY/ZRqFEzyUa6VGB5TNqA93rRU0ZAbP508PSLHVumQbqTtKCQT9fdFSkNjhmFgN0OKfTPvjcV/vSCcS8buuZ8KFSfIFZsGdxNJMCLjRkhHaM5QDi2hTAv7O2F9qilD21vOluLNV7BIaqdFzy16FdvOJUyRhQM4ghPw4AxKcA1lqAKDO3iEJ3h2HpwX59V5m65mnNnNPvyB8/EDkL+aaQ==
AAACBHicbVC7SgNBFL3rM8ZX1FKRwSBYhV0bLSyCNpYJmAcka5idzCZDZmbXmVkhLCn1F2y1txMhlf9ha+tPOHkUmnjgwuGcezmXE8ScaeO6n87C4tLyympmLbu+sbm1ndvZreooUYRWSMQjVQ+wppxJWjHMcFqPFcUi4LQW9K5Gfu2eKs0ieWP6MfUF7kgWMoKNlXzZYhfNgBqM5K3byuXdgjsGmifelOSLB8Py98PhsNTKfTXbEUkElYZwrHXDc2Pjp1gZRjgdZJuJpjEmPdyhDUslFlT76fjpATq2ShuFkbIjDRqrvy9SLLTui8BuCmy6etYbif96gZhJNuG5nzIZJ4ZKMgkOE45MhEaNoDZTlBjetwQTxezviHSxwsTY3rK2FG+2gnlSPS14bsEr23YuYYIM7MMRnIAHZ1CEayhBBQjcwRM8w4vz6Lw6b877ZHXBmd7swR84Hz+FZJvP
AAACBHicbVC7SgNBFL3rM8ZX1FKRwSBYhV0bLSyCNpYJmAcka5idzCZDZmbXmVkhLCn1F2y1txMhlf9ha+tPOHkUmnjgwuGcezmXE8ScaeO6n87C4tLyympmLbu+sbm1ndvZreooUYRWSMQjVQ+wppxJWjHMcFqPFcUi4LQW9K5Gfu2eKs0ieWP6MfUF7kgWMoKNlXzZYhfNgBqM5K3byuXdgjsGmifelOSLB8Py98PhsNTKfTXbEUkElYZwrHXDc2Pjp1gZRjgdZJuJpjEmPdyhDUslFlT76fjpATq2ShuFkbIjDRqrvy9SLLTui8BuCmy6etYbif96gZhJNuG5nzIZJ4ZKMgkOE45MhEaNoDZTlBjetwQTxezviHSxwsTY3rK2FG+2gnlSPS14bsEr23YuYYIM7MMRnIAHZ1CEayhBBQjcwRM8w4vz6Lw6b877ZHXBmd7swR84Hz+FZJvP
AAACBHicbVC7TgMxEPSFVwivACWNRYREFflooKCIoKEMEnlIyRH5nL3Eiu07bB9SdErLL9BCT4do+Q9avgQnuQISRlppNLOrWU2YCG4sIV9eYWV1bX2juFna2t7Z3SvvHzRNnGoGDRaLWLdDakBwBQ3LrYB2ooHKUEArHF1P/dYjaMNjdWfHCQSSDhSPOKPWSYHq8ctuCJZidU965QqpkhnwMvFzUkE56r3yd7cfs1SCskxQYzo+SWyQUW05EzApdVMDCWUjOoCOo4pKMEE2e3qCT5zSx1Gs3SiLZ+rvi4xKY8YydJuS2qFZ9Kbiv14oF5JtdBFkXCWpBcXmwVEqsI3xtBHc5xqYFWNHKNPc/Y7ZkGrKrOut5ErxFytYJs2zqk+q/i2p1K7yeoroCB2jU+Sjc1RDN6iOGoihB/SMXtCr9+S9ee/ex3y14OU3h+gPvM8fXHWYJA==
ѹྗܭࢉ void calPressure(void);
void calPressure(void) {
calNumberDensity();
setBoundaryCondition();
setSourceTerm();
setMatrix();
solveSimultaniousEquationsByGaussEliminationMethod();
removeNegativePressure();
setMinimumPressure();
}
1
⇢0
hr2Pik+1
i
=
1
t2
n⇤
i
n0
n0
AAACU3icbVFBaxNBFH7ZVq3RaqxHL4NFEKVhNxe9FEv14DGCaYvd7PJ28jYZMju7zLwthGWhP82/4MGDJ8GDV717cZIU1NQHM3zzfd+bN3yTVVo5DsMvnWBr+8bNWzu3u3fu7t6733uwd+LK2koayVKX9ixDR1oZGrFiTWeVJSwyTafZ/PVSP70g61Rp3vOionGBU6NyJZE9lfY+HMS5RdlEbRPbWZmEbazRTDXFBjONyWAY29U5VUkzfx61h3/8b0gzCk4G7ZozybNUHRh/R7Pc0t5+2A9XJa6D6ArsH736+ukSAIZp71s8KWVdkGGp0bnzKKx43KBlJTW13bh2VKGc45TOPTRYkBs3qwxa8cQzE5GX1i/DYsX+3dFg4dyiyLyzQJ65TW1J/lfLio3JnL8cN8pUNZOR68F5rQWXYhmwmChLkvXCA5RW+bcLOUOfD/tv6PpQos0IroOTQT8K+9E7n84xrGsHHsFjeAoRvIAjeAtDGIGEj/AdfsDPzufOryAIttfWoHPV8xD+qWD3N541uEo=
AAACU3icbVFBb9MwFHZDB6XboBtHLhbTpGmoVdILu6BNAwmOnUS3iaWJXtyX1qrjRPbLtCrKf+MvcODACYkDV7hzmdtOAjqeZOvz933Pz/qcFEpa8v2vDe9Bc+Pho9bj9ubW9pOnnZ3dc5uXRuBQ5Co3lwlYVFLjkCQpvCwMQpYovEhmbxb6xTUaK3P9geYFjjKYaJlKAeSouPOxG6YGRBXUVWimeeTXoQI9URhqSBRE/UFoludYRtXsZVC//uN/i4qAU9SvV5yODmPZ1e6OarHFnT2/5y+L3wfBHdg7Of72+aZ79G4Qd76H41yUGWoSCqy9CvyCRhUYkkJh3Q5LiwWIGUzwykENGdpRtcyg5vuOGfM0N25p4kv2744KMmvnWeKcGdDUrmsL8r9akq1NpvRoVEldlIRarAanpeKU80XAfCwNClJzB0AY6d7OxRRcPuS+oe1CCdYjuA/O+73A7wVnLp1TtqoWe85esAMWsFfshL1nAzZkgn1iP9hP9qvxpfHb87zmyuo17nqesX/K274Fg3K47g==
AAACU3icbVFBb9MwFHZDB6XboBtHLhbTpGmoVdILu6BNAwmOnUS3iaWJXtyX1qrjRPbLtCrKf+MvcODACYkDV7hzmdtOAjqeZOvz933Pz/qcFEpa8v2vDe9Bc+Pho9bj9ubW9pOnnZ3dc5uXRuBQ5Co3lwlYVFLjkCQpvCwMQpYovEhmbxb6xTUaK3P9geYFjjKYaJlKAeSouPOxG6YGRBXUVWimeeTXoQI9URhqSBRE/UFoludYRtXsZVC//uN/i4qAU9SvV5yODmPZ1e6OarHFnT2/5y+L3wfBHdg7Of72+aZ79G4Qd76H41yUGWoSCqy9CvyCRhUYkkJh3Q5LiwWIGUzwykENGdpRtcyg5vuOGfM0N25p4kv2744KMmvnWeKcGdDUrmsL8r9akq1NpvRoVEldlIRarAanpeKU80XAfCwNClJzB0AY6d7OxRRcPuS+oe1CCdYjuA/O+73A7wVnLp1TtqoWe85esAMWsFfshL1nAzZkgn1iP9hP9qvxpfHb87zmyuo17nqesX/K274Fg3K47g==
AAACU3icbVE9b9swEKXVfNVtUicduxAxChQNYkhemiVAkHTo6AK1HTSyhBN9sglTlECeAhiC/lv+RIauHbo2e5bSH0BbJweQeHzvHY94TAolLfn+j4b3Ymt7Z3fvZfPV6/2DN63Do4HNSyOwL3KVm+sELCqpsU+SFF4XBiFLFA6T2dVCH96isTLX32he4CiDiZapFECOilvfT8PUgKiCugrNNI/8OlSgJwpDDYmCqNsLzfIcy6ianQT1+V//Z1QEnKJuveJ09DGWp9rdUS22uNX2O/6y+FMQrEGbrasXt36G41yUGWoSCqy9CfyCRhUYkkJh3QxLiwWIGUzwxkENGdpRtcyg5u8dM+ZpbtzSxJfsvx0VZNbOs8Q5M6Cp3dQW5LNakm1MpvRsVEldlIRarAanpeKU80XAfCwNClJzB0AY6d7OxRRcPuS+oelCCTYjeAoG3U7gd4Kvfvvich3PHnvHjtkHFrBP7IJ9YT3WZ4LdsV/sN3to3DcePc/bWlm9xrrnLfuvvP0/05W1mw==
ѹྗܭࢉ void calPressure(void);
void calPressure(void) {
calNumberDensity();
setBoundaryCondition();
setSourceTerm();
setMatrix();
solveSimultaniousEquationsByGaussEliminationMethod();
removeNegativePressure();
setMinimumPressure();
}
ai,j =
(
1
⇢0
2d
0n0
w(|r⇤
j
r⇤
i
|) (j 6= i)
1
⇢0
2d
0n0
P
j06=i
w(|r⇤
j0
r⇤
i
|) (j = i)
AAACyXicjVFNb9NAEF2br9Z8BTj2YhEBSUUjOwfgUqmCCxKXIpG2Utax1utxsunu2t1dA2HrE3+w1175JWwSI5WEAyOt9PbNzJvRm6ziTJsouvb8W7fv3L23sxvcf/Dw0ePOk6cnuqwVhREteanOMqKBMwkjwwyHs0oBERmH0+z8wzJ/+hWUZqX8YhYVJIJMJSsYJcZRaecKjwOSWvZ63hziDKZMWurkdBMc4EIRauPGYjUrJ1Gz/g9zR3A3ICeTSDr6W+8SZ8KqJp1P9g9ayCb7l/2XvTmWcBGyPsbBf6phXYvUzl+tG2+IO67Z1j9k/QCDzP8sjZO0040G0SrCbRC3oIvaOE47v3Be0lqANJQTrcdxVJnEEmUY5eAkaw0VoedkCmMHJRGgE7tyvglfOCYPi1K5J024Ym92WCK0XojMVQpiZnoztyT/mcvExmRTvEssk1VtQNL14KLmoSnD5VnDnCmghi8cIFQxt3tIZ8RZbNzxA2dKvGnBNjgZDuJoEH8edo/et/bsoD30HPVQjN6iI/QRHaMRot4bD3vgFf4n/8L/7v9Yl/pe2/MM/RX+z98yu975
AAACyXicjVFNb9NAEF2br9Z8BTj2YhEBSUUjOwfgUqmCCxKXIpG2Utax1utxsunu2t1dA2HrE3+w1175JWwSI5WEAyOt9PbNzJvRm6ziTJsouvb8W7fv3L23sxvcf/Dw0ePOk6cnuqwVhREteanOMqKBMwkjwwyHs0oBERmH0+z8wzJ/+hWUZqX8YhYVJIJMJSsYJcZRaecKjwOSWvZ63hziDKZMWurkdBMc4EIRauPGYjUrJ1Gz/g9zR3A3ICeTSDr6W+8SZ8KqJp1P9g9ayCb7l/2XvTmWcBGyPsbBf6phXYvUzl+tG2+IO67Z1j9k/QCDzP8sjZO0040G0SrCbRC3oIvaOE47v3Be0lqANJQTrcdxVJnEEmUY5eAkaw0VoedkCmMHJRGgE7tyvglfOCYPi1K5J024Ym92WCK0XojMVQpiZnoztyT/mcvExmRTvEssk1VtQNL14KLmoSnD5VnDnCmghi8cIFQxt3tIZ8RZbNzxA2dKvGnBNjgZDuJoEH8edo/et/bsoD30HPVQjN6iI/QRHaMRot4bD3vgFf4n/8L/7v9Yl/pe2/MM/RX+z98yu975
AAACyXicjVFNb9NAEF2br9Z8BTj2YhEBSUUjOwfgUqmCCxKXIpG2Utax1utxsunu2t1dA2HrE3+w1175JWwSI5WEAyOt9PbNzJvRm6ziTJsouvb8W7fv3L23sxvcf/Dw0ePOk6cnuqwVhREteanOMqKBMwkjwwyHs0oBERmH0+z8wzJ/+hWUZqX8YhYVJIJMJSsYJcZRaecKjwOSWvZ63hziDKZMWurkdBMc4EIRauPGYjUrJ1Gz/g9zR3A3ICeTSDr6W+8SZ8KqJp1P9g9ayCb7l/2XvTmWcBGyPsbBf6phXYvUzl+tG2+IO67Z1j9k/QCDzP8sjZO0040G0SrCbRC3oIvaOE47v3Be0lqANJQTrcdxVJnEEmUY5eAkaw0VoedkCmMHJRGgE7tyvglfOCYPi1K5J024Ym92WCK0XojMVQpiZnoztyT/mcvExmRTvEssk1VtQNL14KLmoSnD5VnDnCmghi8cIFQxt3tIZ8RZbNzxA2dKvGnBNjgZDuJoEH8edo/et/bsoD30HPVQjN6iI/QRHaMRot4bD3vgFf4n/8L/7v9Yl/pe2/MM/RX+z98yu975
AAACyXicjVFNb9NAEF2br9Z8BTj2YhEBSUUjOwfgUqmCCxKXIpG2Utax1utxsunu2t1dA2HrE3+w1175JWwSI5WEAyOt9PbNzJvRm6ziTJsouvb8W7fv3L23sxvcf/Dw0ePOk6cnuqwVhREteanOMqKBMwkjwwyHs0oBERmH0+z8wzJ/+hWUZqX8YhYVJIJMJSsYJcZRaecKjwOSWvZ63hziDKZMWurkdBMc4EIRauPGYjUrJ1Gz/g9zR3A3ICeTSDr6W+8SZ8KqJp1P9g9ayCb7l/2XvTmWcBGyPsbBf6phXYvUzl+tG2+IO67Z1j9k/QCDzP8sjZO0040G0SrCbRC3oIvaOE47v3Be0lqANJQTrcdxVJnEEmUY5eAkaw0VoedkCmMHJRGgE7tyvglfOCYPi1K5J024Ym92WCK0XojMVQpiZnoztyT/mcvExmRTvEssk1VtQNL14KLmoSnD5VnDnCmghi8cIFQxt3tIZ8RZbNzxA2dKvGnBNjgZDuJoEH8edo/et/bsoD30HPVQjN6iI/QRHaMRot4bD3vgFf4n/8L/7v9Yl/pe2/MM/RX+z98yu975
ѹྗܭࢉ void calPressure(void);
void calPressure(void) {
calNumberDensity();
setBoundaryCondition();
setSourceTerm();
setMatrix();
solveSimultaniousEquationsByGaussEliminationMethod();
removeNegativePressure();
setMinimumPressure();
}
Ax = b
AAACBXicbVC7TsNAEFyHVzCvACWNRYREFdk00EQEaCiDRB5SYqLz5Zyccne27s6IyLJExy/QQk+HaCn5Blq+hHNCAQkjrXY0s6tdTRAzqrTrflqFhcWl5ZXiqr22vrG5VdreaaookZg0cMQi2Q6QIowK0tBUM9KOJUE8YKQVjC5yv3VLpKKRuNbjmPgcDQQNKUbaSDfdgKdnd1k170HWK5XdijuBM0+8H1I+/bCr9wBQ75W+uv0IJ5wIjRlSquO5sfZTJDXFjGR2N1EkRniEBqRjqECcKD+dfJ05B0bpO2EkTQntTNTfGyniSo15YCY50kM16+Xiv17AZy7r8MRPqYgTTQSeHg4T5ujIySNx+lQSrNnYEIQlNb87eIgkwtoEZ5tQvNkI5knzqOK5Fe/KLdfOYYoi7ME+HIIHx1CDS6hDAzBIeIQneLYerBfr1Xqbjhasn51d+APr/RvOXJsT
AAACBXicbVC7TgMxEPSFVzheAUoaiwiJKrqjgSYiQEMZJPKQkhDZji+xYvtOtg8Rna7mF2iBmg7RUvIJiJYvwZekgISRVjua2dWuBkecaeN5X05uYXFpeSW/6q6tb2xuFbZ36jqMFaE1EvJQNTHSlDNJa4YZTpuRokhgTht4eJH5jVuqNAvltRlFtCNQX7KAEWSsdNPGIjm7S8tZx2m3UPRK3hhwnvhTUjz9cMvR86db7Ra+272QxIJKQzjSuuV7kekkSBlGOE3ddqxphMgQ9WnLUokE1Z1k/HUKD6zSg0GobEkDx+rvjQQJrUcC20mBzEDPepn4r4fFzGUTnHQSJqPYUEkmh4OYQxPCLBLYY4oSw0eWIKKY/R2SAVKIGBuca0PxZyOYJ/Wjku+V/CuvWDkHE+TBHtgHh8AHx6ACLkEV1AABCjyAR/Dk3DsvzqvzNhnNOdOdXfAHzvsP1yuchw==
AAACBXicbVC7TgMxEPSFVzheAUoaiwiJKrqjgSYiQEMZJPKQkhDZji+xYvtOtg8Rna7mF2iBmg7RUvIJiJYvwZekgISRVjua2dWuBkecaeN5X05uYXFpeSW/6q6tb2xuFbZ36jqMFaE1EvJQNTHSlDNJa4YZTpuRokhgTht4eJH5jVuqNAvltRlFtCNQX7KAEWSsdNPGIjm7S8tZx2m3UPRK3hhwnvhTUjz9cMvR86db7Ra+272QxIJKQzjSuuV7kekkSBlGOE3ddqxphMgQ9WnLUokE1Z1k/HUKD6zSg0GobEkDx+rvjQQJrUcC20mBzEDPepn4r4fFzGUTnHQSJqPYUEkmh4OYQxPCLBLYY4oSw0eWIKKY/R2SAVKIGBuca0PxZyOYJ/Wjku+V/CuvWDkHE+TBHtgHh8AHx6ACLkEV1AABCjyAR/Dk3DsvzqvzNhnNOdOdXfAHzvsP1yuchw==
AAACBXicbVA9T8MwEL3wWcpXgZHFokJiqhIWWJAKLIxFoh9SGyrbdVqrdhLZDqKKMvMXWGFnQ6z8DlZ+CU6bAVqedLqn9+50p0diwbVx3S9naXlldW29tFHe3Nre2a3s7bd0lCjKmjQSkeoQrJngIWsabgTrxIphSQRrk/F17rcfmNI8Cu/MJGa+xMOQB5xiY6X7HpHp5WN2kXeS9StVt+ZOgRaJV5AqFGj0K9+9QUQTyUJDBda667mx8VOsDKeCZeVeolmM6RgPWdfSEEum/XT6dYaOrTJAQaRshQZN1d8bKZZaTySxkxKbkZ73cvFfj8i5yyY491MexolhIZ0dDhKBTITySNCAK0aNmFiCqeL2d0RHWGFqbHBlG4o3H8EiaZ3WPLfm3brV+lURTwkO4QhOwIMzqMMNNKAJFBQ8wwu8Ok/Om/PufMxGl5xi5wD+wPn8AUJUmUY=
γϛϡϨʔγϣϯຊମ
ԾΛܭࢉ͢Δ
ԾΛͬͯҠಈ͢Δ
ඇѹॖ݅Λຬͨ͢Α͏ʹ
ѹྗΛܭࢉ
ҐஔͱΛ࠶ܭࢉ
void calGravity(void);
void calViscosity(void);
void moveParticle(void);
void collision(void);
void calPressure(void);
void calPressureGradient(void);
void moveParticleUsingPressureGradient(void);
ѹྗޯ߲
ɾ࠷ѹྗΛ͏͜ͱͰޯ͕ৗʹਖ਼ʹ
ɾੱྗ͕ಇ͖ඇѹॖ݅Λຬͨͨ͠ஔʹ
ɾਖ਼͍͠ʹͳΔͨΊͷՃ͕ग़Δɹɹɹɹɹɹɹ
void calPressureGradient(void);
uk+1
i
= u⇤
i
1
⇢0
hrPik+1
i
t
1
⇢0
hrPik+1
i
=
1
⇢0
d
n0
X
j6=i
Pk+1
j
ˆ
Pk+1
i
|r⇤
j
r⇤
i
|2
(r⇤
j
r⇤
i
)w(|r⇤
j
r⇤
i
|)
!
AAADSniclVJLbxMxEPZuC5TllcKRi0VElRQa7eYCl0oVcECcgkTaSnWy8jqzGzde79b2giJ3fxT/gzviBjfO3BAXnAcFknJgJEvfzHyf7XkkpeDahOEnz9/YvHL12tb14MbNW7fvNLbvHuqiUgz6rBCFOk6oBsEl9A03Ao5LBTRPBBwlk+ez/NFbUJoX8o2ZljDIaSZ5yhk1LhRve69IAhmXNhVU8Ezu1sEOSXJb1TEf2smjqN6/cHf3SKoos1FtiRoXw7AmgspMAJE0ERT3iJq7v5TkBQhDsSEk2PlP6f46f+6Paitnjq7y2J4SCWeYu6sgNa0FoRefLm7YI2NqbO+ijNqezwpRtSO4QhbQ1XQ+7NatSzPtd61/SNpE8Wxs2gEBOfrdubjRDDvh3PA6iJagiZbWixtfyahgVQ7SMEG1PonC0jxOBZzJgaXKcCagDkiloaRsQjM4cVDSHPTAzkdf44cuMsJpodyRBs+jfyoszbWe5olj5tSM9WpuFrw0l+QrL5v06cByWVYGJFs8nFYCmwLP9gqPuAJmxNQByhR3f8dsTN1IjNu+wLUmWm3EOjjsdqKwE73uNg+eLZu0he6jB6iFIvQEHaCXqIf6iHnvvY/eZ++L/8H/5n/3fyyovrfU3EN/2cbmTw1HE+I=
AAADSniclVJLbxMxEPZuC5TllcKRi0VElRQa7eYCl0oVcECcgkTaSnWy8jqzGzde79b2giJ3fxT/gzviBjfO3BAXnAcFknJgJEvfzHyf7XkkpeDahOEnz9/YvHL12tb14MbNW7fvNLbvHuqiUgz6rBCFOk6oBsEl9A03Ao5LBTRPBBwlk+ez/NFbUJoX8o2ZljDIaSZ5yhk1LhRve69IAhmXNhVU8Ezu1sEOSXJb1TEf2smjqN6/cHf3SKoos1FtiRoXw7AmgspMAJE0ERT3iJq7v5TkBQhDsSEk2PlP6f46f+6Paitnjq7y2J4SCWeYu6sgNa0FoRefLm7YI2NqbO+ijNqezwpRtSO4QhbQ1XQ+7NatSzPtd61/SNpE8Wxs2gEBOfrdubjRDDvh3PA6iJagiZbWixtfyahgVQ7SMEG1PonC0jxOBZzJgaXKcCagDkiloaRsQjM4cVDSHPTAzkdf44cuMsJpodyRBs+jfyoszbWe5olj5tSM9WpuFrw0l+QrL5v06cByWVYGJFs8nFYCmwLP9gqPuAJmxNQByhR3f8dsTN1IjNu+wLUmWm3EOjjsdqKwE73uNg+eLZu0he6jB6iFIvQEHaCXqIf6iHnvvY/eZ++L/8H/5n/3fyyovrfU3EN/2cbmTw1HE+I=
AAADSniclVJLbxMxEPZuC5TllcKRi0VElRQa7eYCl0oVcECcgkTaSnWy8jqzGzde79b2giJ3fxT/gzviBjfO3BAXnAcFknJgJEvfzHyf7XkkpeDahOEnz9/YvHL12tb14MbNW7fvNLbvHuqiUgz6rBCFOk6oBsEl9A03Ao5LBTRPBBwlk+ez/NFbUJoX8o2ZljDIaSZ5yhk1LhRve69IAhmXNhVU8Ezu1sEOSXJb1TEf2smjqN6/cHf3SKoos1FtiRoXw7AmgspMAJE0ERT3iJq7v5TkBQhDsSEk2PlP6f46f+6Paitnjq7y2J4SCWeYu6sgNa0FoRefLm7YI2NqbO+ijNqezwpRtSO4QhbQ1XQ+7NatSzPtd61/SNpE8Wxs2gEBOfrdubjRDDvh3PA6iJagiZbWixtfyahgVQ7SMEG1PonC0jxOBZzJgaXKcCagDkiloaRsQjM4cVDSHPTAzkdf44cuMsJpodyRBs+jfyoszbWe5olj5tSM9WpuFrw0l+QrL5v06cByWVYGJFs8nFYCmwLP9gqPuAJmxNQByhR3f8dsTN1IjNu+wLUmWm3EOjjsdqKwE73uNg+eLZu0he6jB6iFIvQEHaCXqIf6iHnvvY/eZ++L/8H/5n/3fyyovrfU3EN/2cbmTw1HE+I=
AAADSniclVJLbxMxEPZuC5TllcKRi0VElRQa7eYCl0oVcECcgkTaSnWy8jqzGzde79b2giJ3fxT/gzviBjfO3BAXnAcFknJgJEvfzHyf7XkkpeDahOEnz9/YvHL12tb14MbNW7fvNLbvHuqiUgz6rBCFOk6oBsEl9A03Ao5LBTRPBBwlk+ez/NFbUJoX8o2ZljDIaSZ5yhk1LhRve69IAhmXNhVU8Ezu1sEOSXJb1TEf2smjqN6/cHf3SKoos1FtiRoXw7AmgspMAJE0ERT3iJq7v5TkBQhDsSEk2PlP6f46f+6Paitnjq7y2J4SCWeYu6sgNa0FoRefLm7YI2NqbO+ijNqezwpRtSO4QhbQ1XQ+7NatSzPtd61/SNpE8Wxs2gEBOfrdubjRDDvh3PA6iJagiZbWixtfyahgVQ7SMEG1PonC0jxOBZzJgaXKcCagDkiloaRsQjM4cVDSHPTAzkdf44cuMsJpodyRBs+jfyoszbWe5olj5tSM9WpuFrw0l+QrL5v06cByWVYGJFs8nFYCmwLP9gqPuAJmxNQByhR3f8dsTN1IjNu+wLUmWm3EOjjsdqKwE73uNg+eLZu0he6jB6iFIvQEHaCXqIf6iHnvvY/eZ++L/8H/5n/3fyyovrfU3EN/2cbmTw1HE+I=
γϛϡϨʔγϣϯຊମ
ԾΛܭࢉ͢Δ
ԾΛͬͯҠಈ͢Δ
ඇѹॖ݅Λຬͨ͢Α͏ʹ
ѹྗΛܭࢉ
ҐஔͱΛ࠶ܭࢉ
void calGravity(void);
void calViscosity(void);
void moveParticle(void);
void collision(void);
void calPressure(void);
void calPressureGradient(void);
void moveParticleUsingPressureGradient(void);
·ͱΊ
γϛϡϨʔγϣϯͷྲྀΕ
ཻࢠΛฒΔ
ԾΛܭࢉ͢Δ
ԾΛͬͯҠಈ͢Δ
ඇѹॖ݅Λຬͨ͢Α͏ʹ
ѹྗΛܭࢉ
ҐஔͱΛ࠶ܭࢉ
αϯϓϧͷಈ࡞ྫ
࣍ʹ͚ͯ
ɾ$ʹҠ২ˍվΛऴΘΒͤΔ
ɾ࢛ষͳͲΛࢀߟʹฒྻԽʢҰ୴0QFO.1*ʣ