home - stfc computational mathematics groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö...

26

Upload: others

Post on 30-Jul-2020

4 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

�������������� � � �������� �������������! "$#!���&% '(�&)*�� � � ,+

- .0/2143�576 89/ :7;2.0/21489/<1 =?>4@93�6 89/ AB;<3 5C.DAB89>4/<[email protected]/<:IHJ>46 K 3IAB;2>ML NO>45 AB;<3 P 3�57K 8QAB89.E/

3�89143�/SRT.E@9U<3 =<5C>4V<@93�K

WMXZYQ[7\^]�]OY`_!]OYbac]ed

f�g�hji�kJf�l7imonpnpqsrutwvyxzru{b|~}�}��~t���v�}�xz�����,{9��v�q!v�|�v�qs�9}zru}�{b�B}zxz�z��t�xz������}�v�|��Sv�t�{b�~}zxzr�t�tyv���rsxzr���},��������r�����v|`�p�B�����J{b����r�����|`��vyq��~��}Jvy|~����r�����|`����t�x�{b��}�{��&q�v����Q��}�twv�q��������z��rsxzr vy|S��ru�Q��|`��v�qs�~��np��{b�pqu����}�¡¢ ���!��{�}zx£n�{bnp�pq�v��¤����xz��{9�Jr�}¥np��{b�¦vy�pqs��xz�~�B}�n���t�xz�§v�q�xz�§vy|~}z�¨{b�z��v�xzr�{�|�©jvy|~t�ª�{Q}c����xz�~{9��¡«mo|rs��n�{b�zx§v�|`xj�9r�}§vb�9��v�|`x§vb���¬{b�9xz�pru}��,��xz�~{9�&r�}Oxz�¦vyx�v¥t­�¦vy|~���Z{b�`n�{bqu�®�����Q�pr�����}¬v¥t�{���npq���x��¬����}�x§v��zx�¡¯ |°xz��ru}on¦vyn�����±�²¥�³r�|`����}zxzr��`v�x���xz�~���~}­��{��®xz�~���§vyxzr�{�|¦vyq�´&�z�`q�{w�J����xz�~{9���¨{��¤xz�pr�}�vynpnpqsr�tyvyxzr�{�|e¡¢ ��ru}¬����xz�~{9���p{9��}µ|�{�x¬��������r����&v!t�{���npq���x��¶����}zx§vy�zxcv���x����cv!t��¦vy|~�Q�¶{b�¦n�{bqu��¡ ¢ ���¶t�{�|`xz�zrs�p�pxzr�{�|{b�cxz�pru}³n¦v�n����!ru}�xz�p�������¨{�q���¡�·jrs��}zx�±er�x�ru}�}��~{w²¶|Dxz�¦v�x!xz����t��¦vy|~�Q��{b�Zn�{bq���twvy|�����t�{b|~}zr��p�������vb}�v¸t­�¦vy|~���,{��!©jvy|~t�ª�{Q}���v�}zr�}�¡E¹9��t�{b|~��±µ�,{���rs|~�¸xz����n�{�q���|~�wvy�,v�q�{9t�ºQ���4��r��Q��|`��vyq�������v��np��������|`x�{bxz�~���¥��ru�Q��|`��v�qs�~��}cx�{C���&t�{b��n��px����Jv�t�t��p�§vyx���q���¡ ¢ �pr�����±�²¥�³}��~{w²Sxz��v�x�v�n�{�q��&t��~{Q}���|t�qu{Q}���x�{�vy|»��r�����|`��vyq��~���Iv��Tq��yvb��x�{�vJq�{Q}�}7{b�¥|`�p�����zrutwvyq¶}zx§vy�prsq�rsx¼��¡�½o}z�¦vyq�qs��±exz�pr�}7²¤r�qsqc|~{bxq��yvb�Cx�{��9r�}§vb}zx����µ�p��xc�������~t��¥xz�~�&v�t�t��p�§vbt��7{b�¦xz�~�ot�{���np��x�������r�����|`��vyq��~��}cv�²¶v��¾����{��¿xz�~�£n�{bqu��¡À �p�����zrutwvyqe�§Á�vy��n�qu��}¤r�qsq���}zxz�§v�x��³xz�~�!xz�~��{b�z�`¡

Â�Ã�Ä�ÅCÆ~Ç�È®ÉQʾ�o���z��r�xzr�v�|���r�����|`��vyq��~��np��{b�pqu���˱O©ev�|~t�ª�{Q}¶����xz�~{9��±~�§v�xzr�{�|�v�q�´&�z�`q�{w��}���������|�t���±}zn���t�xz�§vyq�xz�§vy|~}z�¨{b�z��v�xzr�{b|

f�ÌÍhMÎÏÌÍÐ�hµÑ³ÉwÒjÓOÔ�Ã`Õ�Ö�Õ`×�ØpÉ�ÉwÙsÚµÕQØ9Ö�Ù�ƦÛjÉ�Ê�Ü�Ýb·¥ÞwÝ

ßc{b��np�px§vyxzru{b|¦vyqj¹9t�ru��|~t��Bvy|~�Jà¬|~�brs|~�����zrs|~��á&��n¦vy�zxz����|`xm¥xzq�vb}³ßc��|`xz���â �pxz�~���z�¨{�����monpnpq���x�{�|¸©jvy��{b�§v�x�{b�z�ã¶Á�{b|Dã£ä�Þ�Þ!å�æ¶äÀ {w�����¾�����³çQÝ9±OÞwèQè�è

Page 2: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

é #�ê,ëµ�&ê�ë¬�ì í Û�ÖwÇ�Æ�ȵÒeÕ�Ö�Ù�Æ�Û ìî iBïjÃ�É�ðZÃ�ÕbÖwÇ�Ø~×ZÖwÇ�Ø�ÛjÉ�ñ¼Æ�Ç�ò�Ø�ÖbÙ�Æ�ÛSócØ�ÛjÕQô�Æ�É�ò�ÃQÖ�ïjÆ�È îõ iBïjÃ�Ç�Ø9Ö�Ù�Æ�ÛjØ~×£ócØ~ÛjÕQô�Æ�É�ò�ÃQÖ�ïjÆ�È õö ÷ ÒjÇ�ø�Ù�ÛjøOù¬Ù�ò�ðµ×�Ù�Õ`ÙsÖËÇyÃ�É�ÖwØpÇwÖwÉ�Ø~ÛeÈS×�Æ�Õ�ú�Ù�Ûjø ûü iBïjÃ�ÛOÒ�òDÃ`Ç�Ù�Õ�Ø�×cÉ�ÖyØ�ÓµÙ�×�ÙsÖ§ÄýÆ~ñ¤Ö�ïjÃ�Ç�Ø9Ö�Ù�ƦÛeØ~×¥ócØ~ÛeÕ�ôQÆ~É�òDÃQÖbïeÆ�È ì�îþ ÿ Ò®ò�Ã`Ç�Ù�ÕQØ~×cÃ��OØ~ò�ð¬×�à ì�þ� lBÆ�ÛjÕ`×�ÒjÉwÙ�ƦÛeÉ î��f k�Ã`ȵÒjÕbÖ�Ù�ƦÛ4Æ�ñ�Ø�Ç�Ã`Õ�ÖwØ~Ûeø¦Ò®×�ØpÇJò�Ø�ÖyÇbÙ��4ÖyÆDÖwÇ�Ù�ȵÙ�Øpø�ƦÛeØ~×¥ñ Æ~Ç�ò î9ö

r

Page 3: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

� �Oê,ëµ)Z#Ë ��£ë«��#�ê¢ ���7}z�����z��t�x&{��®xz�pr�}&n¦v�n�����r�}&xz�~�¾t�{b��np�px§vyxzru{b|�{b�£v�|`�p�¾�����!{b�«��r�����|`��v�qs�~��}!v�|���xz�~�Bt�{b�z��� �}zn�{�|��9r�|�����r�����|`����t�x�{b��}¶{b�jxz�~�!��vyxz�zrsÁJn���|~t�rsq

����������� � ����ýå�� �zÞ�¡�Þ��²¶��������� vy|~� � vy���¶q�vy�����¶vy|~��}zn�v���}��¤�����z��rsxzr vy|���v�xz�zr�t���}cvy|~�!��r�}¬n�{Q}zrsxzr��Q���¨}�����r"�¬�p�$#p|�r�x���¡% �&twvyq�q&��v�|���r��Q��|`��vyq�����vy|~���°v¾t�{��z����}zn�{�|��9r�|��7��r�����|`����t�x�{��(')�*�+�,�-�Zr�}¶twvyq�q����Jvy|J��ru�Q��|pn¦vyr���¡monpnpqsrutwvyxzru{b|~}Z{b��xz��ru}Z�¨{b�z� vy�zr�}��¶rs|�}zxz�z�~t�xz�p�§vyq¦��|~�brs|~�����zrs|~�.�*/o�zr�����}�±�©O��²¶r�}¥v�|��,¹�r���{b|�Þyè�0QÜ1�vy|~�Tvbt�{b�~}zxzr�t¾�,{9��v�qZv�|¦vyqs��}�ru}2�43¬r�����t���Þwè10�Þ��§¡ ¢ �~����{Q}zx!t�{b����{�|pqs���~}����¸����xz�~{9�p}Bvy���Cxz�~�}zn���t�xz�§vyqOxz�§vy|~}z�¨{b�z�Ivyxzr�{�|�©jv�|�t�ª�{Q}c����xz�~{9���¨à¬�zr�t�}�}�{b|Dvy|~� â �p����Þyè�0Qå�± À {��p�5�zão��ru��±�3µvy�zqu��xzx�±à¬�zrut�}�}­{�|�vy|~�76���|~}­��|4Þwè�098:�¶vy|~��rsx�}���qu{9t­ºË�Q����}�ru{b|;�*/o�zrs�,��}�±e©O��²¶r�}!v�|��D¹�r���{b|4ÞyèQè=<>�§¡ ¢ �~�©jvy|~t�ª�{�}O����xz�~{9���p�prsqu�p}Zv&´&�z�`qu{��B}zn¦vbt��µ�¨{b�jxz�~�£}zn���t�xz�§vyq�xz�§vy|~}z�¨{b�z��v�xzr�{b|��4�!?A@)�B�DCFE,�ý}zx§v��zx5�rs|~������{b� v³|~{b|~ª�����{³����t�x�{��HG E ¡ ¢ ���¶����xz�~{9��x ��n�rutwvyq�qs��t�{b|`�������Q��}Z�Q�pr�t�º`qs���¨{��Zxz�~�o��ru�Q��|���vyqs�~��}{b�I�zÞ�¡�Þ��£|~�wv��oxz�~��n�{bqu�J@¬¡LKc���z��{���x���|e±¦v�����q vyxzr��Q��q���q vy���Q��|`�p�B�����!{��Z��ru�Q��|���vyqs�~��}�r�}o²¶v�|`x�����±}§v��DÞyå9±�v�ÞyåQåB{b�¶��{b����±�v�|~�Ixz�~��|�r�xo��v������Bv�����v�|`x§vb����{b�~}¥x�{�t���v�|~�Q�!xz�~�³n�{bqu�³x�{�}�n������M� ��nxz�~��t�{b|`����������|~t���¡D½¶|��¨{��zxz�p|�v�x���qs�`±jxz�pru}C�����Q�prs����}¾����}zx§vy�zxzr�|~��xz��������xz�~{9�T����{b��xz�~�I}§v�����{b�v�|~��²NG E ¡ ¢ �~�I�9ru} v��9��vy|`x§v��Q��r�}Cxz�¦v�xCxz�~��´&�z�`qu{w���¦vb}zru}¾r�}7t�{b��n�qu��x���q���xz�p��{w²¤|Sv�²¶v���¡ % �np��{bn�{Q}��C�~}zrs|~��xz�~���§vyxzr�{�|¦vyq¥´&�z�`qu{�������xz��{9��±®�p������q�{�n����T�`� â �p�~���zÞwè10:<>�§±e²¶�pr�t��(�p{9��}!|~{bx�����Q�pr����CvCt�{b��npq���x��&����}zx§vy�zx¶²¤�~��|O@Dt��¦vy|~�Q��}�¡

¯ |�xz�~�C©jv�|�t�ª�{Q}¥����xz�~{9��±�xz���7vynpnp��{�Á�rs�Ivyx�����r�����|pn¦vyrs��}�±�twv�qsq���� â r�x�ª�n¦v�rs��}�±�vy���Bt�{b��n��px��������{b� vI�����z��rsxzr vy|¸xz�zru��r vb��{b|¦vyqj��v�xz�zr Á�v�|��Ërs|�xz�~�¾�§vyxzru{b|¦vyq«´&�z�`q�{w�Ë����xz�~{9������{b� vy|¸�pnpn�����o��}�}���|`�������¸��vyxz�zrsÁ»n¦vyr���¡ ¯ |�xz�pr�}�n¦vyn�����±«²¥���~}���v�|Svyqsx����z|�v�xzrs���I�¨{���xz�~���§vyxzru{b|¦vyq³´&�z�`q�{w�����xz�~{9��±P#p��}zx��¨{��z�¾�pq vyx����¸rs|Q� â �p�~��Þwè�è�01�§±�}�{Ixz�¦vyx�xz�~� â r�x�ªC��vyq��~��}Cv�����t�{���np��x����D����{���v�o���z��rsxzr�v�|���v�xz�zr Á�¡o·~{b�oxz�pru}o���wvb}�{�|O±�²¥�Btwv�qsqjxz�~������xz��{9�Ë�§vyxzr�{�|¦vyq�©jvy|~t�ª�{�}�¡ ¯ x³��rSR�����}�����{b��§vyxzru{b|¦vyq�´&�z�`q�{w��v�}®©jv�|�t�ª�{Q}T� �����z��rsxzr vy|Bnp��{b�pq����,}U���9rSR�����}O����{b� m��z|�{�q��9rV�¨m��z|�{�q��9r�ÞyèQÝpÞ�±w¹�v�vb�Þwè�èQç1����|�{�|��Ï�o���z��rsxzr vy|Jn���{��pq�����}W�§¡ ¢ �~�&�§v�xzr�{b|¦v�qj´&�z�`q�{w������xz�~{9�Jr�}¥|�{�x¥�Q���z�J{���x���|I�~}������¨{b�}�{bq��`rs|~��v�npn�q�r�twv�xzr�{�|�}�¡�ão|��7���wvb}�{b|�ru}³xz�¦v�x!xz�����,��xz�~{9�Dru}³|~{�x!�Q���z��²¥��q�q¬�p|~�p����}zx�{9{9��¡ ¯ x�r�}�p���zrs������rs|Ëv¾�9rSR�������|�x£²¶v���v�}Zxz�~��©ev�|~t�ª�{Q}µ����xz�~{9�Jvy|~��xz�������&r�}¶v³q vy�����!vy��{��p|`x¥{b�����������p{b�rs|°}���xzxzrs|~��xz�~�³��v��zr�{���}¶n¦vy�§vy�,��x�����}�¡ ¯ |°xz��ru}¤n¦vyn�����±�²¥��q�{9{bºJv�x¤xz�~���§vyxzru{b|¦vyqj´&�z�`qu{w�I����xz��{9�����{b�GvC�9rSR�������|�x!vy|~��q��!v�}¤m¶Áp��q â �p�~��¡ ¢ �~�!{b|pqs������������{��?rs|�xz�~�!����xz�~{9�J²¥�&n���{�n�{�}­�&ru}¤xz�~�n�{�q���¡ ¢ ���!t�{b|`xz�zr����pxzru{b|¸{��exz�pru}¤n¦vyn����¥r�}¤xz�p�������¨{bq���¡Z·er���}zx�±�rsx&}z�~{w²o}cxz�¦vyx¶xz�~��t­�¦vy|~���!{b��n�{bqu�twv�|����Bt�{b|~}zr��p�������¸vb}!v�t��¦vy|~�Q��{b�£©jvy|~t�ª�{�}¤�¦vb}zr�}�¡&¹9��t�{�|���±���{w�`rs|~��xz�~�¾n�{bq���|~�wvy��v�q�{9t�ºQ���â r�x�ª���vyqs�~����v��Tn������Q��|�x�{bxz�~���7��ru�Q��|`��v�qs�~��}Bx�{�����t�{b��n��px����»vbt�t����§v�x���qs�`¡ ¢ �prs����±�²¥��}z�~{w²xz�¦vyx¥v!n�{�q��¶t­�~{�}­��|�t�q�{�}��£x�{7v â rsx�ª¶��vyqs�~�¶��v��Cq��yvb��x�{7v�qu{Q}�}Z{���|`�p�����zrutwvyq�}zx§vy�prsq�rsx �`¡µ½�}��¦v�qsqs�`±xz�pr�}¶²¤r�qsqO|~{bx¶q��wv��Ix�{��9r�}§v�}�x����£�p�pxo�������~t��!xz�~�Bvbt�t����§v�t��J{b�jxz�~� â rsx�ª!��v�qs�~��}&vC�pr�x&v�²¶v�������{b�xz�~�¤n�{bq��Q¡ ¯ |�xz�pr�}c}zrsxz�¦vyxzru{b|e±Q²¥�o}z�~�Q�Q��}�x«v�t�{b��npqu��x��£����}�x§v��zxc{b�¦xz�~�!©jvy|~t�ª�{�}®����xz�~{9�C²¶rsxz��xz�~�|~��²�n�{bq��Q¡

¢ �~�Inpq�v�|S{b�&xz�~��n¦vyn�����ru}Iv�}C�¨{bq�q�{w²�}�¡ ·jrs��}zx�±¤²¥�I�����`r���² xz�~��}zn���t�xz�§vyqoxz�§v�|~}��¨{��z��vyxzru{b|©jvy|~t�ª�{�}c����xz�~{9��rs|YX�ç9¡ % �&r�|`xz��{9�9��t���xz�~�!�§vyxzru{b|¦vyq�©jvy|~t�ª�{Q}c����xz�~{9�Ir�|ZX$[�v�|��Ë}z��{w² }�{����{b��rsx�}�n���{�n����zxzru��}�¡ ¯ |ZXU<�±p²¥�B�§Á�npq vyrs|°��{w² t�{�|`�������Q��� â rsx�ª�n¦vyr���}otwv�|����!q�{9t�º����¸vy|~���~{w² ²¥�twv�|Inp�p���Q�!�p|`²¶vy|�x���� â r�x�ª&�Q��t�x�{b��}�¡ ¢ �pr�}¥r�}¶|�{�x¤|~��² ²¥{b�zº¦±��p�pxo²¥��}z�~{w²ýxz�¦vyx�t­�¦v�|���rs|~�Cxz�~�n�{�q��Bvy��x����¤qu{9t­º`r�|~� â rsx�ª�n¦vy��}¶��v��I�����9�~t���xz�~�³qu������q�{b�«vbt�t��p�§v�t��J{b��xz��������xz�~{9��¡ ¯ |\X�Ý9±�r�xor�}}z�~{�²¶|¸xz�¦vyx!xz�~�C����xz��{9���Iv��°q�{Q}��B��v�t­º�²¶vy����}zx§vy�prsq�rsx �°²¤�~��|¸xz�~�Cn�{bq��Br�}!t��~{Q}���|�t�qu{Q}��¾x�{Ëvâ r�x�ª£��v�qs�~��¡ ¯ |]X�Ü�±�²¥�¥np����}���|`x¥v&|`�p�����zr�twv�q�����}��pq�xZ�¨{b�£vy|�v�n�npq�r�twv�xzr�{b|�t�{b��rs|~�!����{b� vbt�{���}zxzrut�}�¡% �!t�{b|~t�q����p�&xz�~�!n¦vyn����¤²¶rsxz�Ëxz���!��v�rs|°t�{b|~t�q���}zru{b|~}crs|�X^8�¡

Þ

Page 4: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

¢ �~�o�9rs�,��|~}�ru{b|J{b��xz����n���{��pq����_�zÞ�¡�Þ`�¬r�}c�p��|~{bx�������� a�¡cbZ�!dDe£²¥�&�p��|~{bx��¶xz�~�HfQxz�It�{bqs�p��|{b�¦xz�~�or��p��|`xzr�x ���Ivyxz�zr Á�gp¡ ¢ �p��{b�~����{��pxZxz�~�¤n¦vyn�����±1�-hDiC�p��|~{bx���}µxz����}zx§vy|~�~vy����rs|p|~���cnp��{9�9�~t�x�±vy|~�kj �lj m�xz���Ërs|~�9��t����(x ²¥{=� |~{b�z��±n� h ��iExz�~�7��rs|p|~����np��{9�9��t�x�vy|~�kj$�oj$pq�sr � h ��� xz�~�rs|~�9�~t����O�2|~{b�z�Ë¡ ¢ �~�B�2|~{b�z� {��Oxz���!��v�xz�zr Á7tvuxw�yMz`y�±��p��|~{�x����,���Yj`t{j pµ±�r�}¶�p�$#p|����I�`���vyÁ+|"}~|��M� E j`t��oj$pµ¡ ¢ �pr�}or���npqsr���}oxz�¦vyx2j`t��oj$p\��j�t j$p�j �lj pE�¨{b��vy|`�]�7uYw�y�vy|~�¸v�|`�7t�uw�yMz�y�¡ % �¥v�q�}�{o�p�$#p|~�cxz�~����·p��{b����|prs�~}j|�{��z�<{b��tvuxw�y�z:�cv�}�j`t{j pV�O���,� �e � E j�t�dDe:j mp+� E*� m ¡bZ�ËÒZ±�²¥���p��|~{bx��&xz�~�!��vbt��prs|~�!np����t�r�}zru{b|��p�$#p|~���°vb}¥xz�����9r�R�������|~t��!����x ²¥����|0ÞBvy|~��xz���!|~�§Á�x�~{`v�xzrs|~�In�{brs|�x¾|`�p�B������� �¤ru�b�¦vy� ÞyèQè�Ü��§¡�bZ�Z���4�A�U�J��²¥���p��|~{bx���xz�~��}���xB{b�¶��r�����|`��v�qs�~��}B{b��zÞ�¡�Þ��§¡ % �Bv�q�}�{¾�~}��B�9�+���>�*t��¤v�|��]�M�)�����*t��Zx�{��p��|~{bx��!xz�~�³��rs|prs�B�p��v�|��J��vwÁ9rs�B���G}zrs|~�b�pq vy���vyq��~��{���tM����}zn���t�xzrs����qs�`¡ ¢ �~��t�{b|~�9rsxzr�{�|I|����B�����or�}¶�p�$#p|~���I�`�O��m`�*t��c���M�)���=�*t��D�:�M�+� �>�*t��§¡� ���¸� �P »�(£ëµ)c¡£¢(ë¬)c¡Bê¸�o¤Q#�)Z% ¡�ë«��#�ê �¥¡¾êL¦o#�� % �oë§�¸#� ¯ |Txz�pr�}B}���t�xzr�{�|e±�²¥���p��}�t��zrs���Cxz�~��}�n���t�xz�§v�q¬xz�§vy|~}z�¨{b�z�Ivyxzr�{�|0©ev�|~t�ª�{Q}&�,��xz�~{9�°�¨{��7�o���z��r�xzr�vy|��ru�Q��|`��v�qs�~��np��{b�pq����,}�¡�·~{b�Bxz����{b����xzrutwvyq¶v�|���rs��n�qu������|�x§vyxzr�{�|4�p��x§vyr�q�}�±®²¥��}��¦v�qsqZ�����¨���¾x�{�xz�~�qsr�x����§vyxz�p����¡

¢ �~��vyr�� ru}£x�{�t�{b��np�px���±��¨{��o��rs����|�G E v�|~��@�±�vB��v�}zr�}LG E �^¨^¨~¨^�,G �D© E �¨{b�¶xz���B´&�z�`qu{��J}zn¦vbt��ª �D© E � }zn¦vy|�«`G E �~�4�Q?Y@)�B� C�E �¬G E �^¨^¨^¨^���5�4�Q?Z@P�J� CFE �B� � G EW­ ¨ � ç9¡�Þ��¢ ���7��vyxz�zrsÁ®�4�k?;@)�J�DCFEU�?r�}�twv�qsq�����xz����}zn���t�xz�§vyqµxz�§vy|~}z�¨{b�z��v�xzr�{�|Ev�|~�x@¯u®°$ru}³xz�~�Cn�{bqu��¡¢ ���B|¦vy���B}zn���t�xz�§v�qjxz�§vy|~}z�¨{b�z��v�xzr�{b|Dt�{b����}o����{b�Gxz���B� vbt�x�xz�¦vyx&r����*�+�D�V�&}§vyxzru}5#���}J�zÞ�¡�Þ��¥xz����|�4±o�x�*�{?Z@)�DCFEW�,�-�£r�}³vy|°��r�����|pn�v�rs�¶{b�L�4�Q?Y@)�B�DC�EW��¡ ¢ �~��}zn���t�xz�§v�q�xz�§vy|~}z�¨{b�z�Ivyxzr�{�|¸©jvy|~t�ª�{Q}����xz�~{9�St�{b��n��px���}���ru�Q��|pn¦vyr���}7�4±M�,�-��{b�!xz�~��}�n���t�xz�§v�qoxz�§vy|~}z�¨{b�z�Ivyxzr�{�|e¡ ¢ �~������q vyxzru{b|��²�@!³Y±9C�E¶vyq�q�{w²o}¥�~}£x�{�t�{b��np�px��!xz�~��t�{b�z����}zn�{�|~��r�|~�!�+´ }¬�¨à¬�zrut�}�}�{b|¸v�|~� â ���~��Þwè10Qå1�§¡

¢ �~�x� �����z��rsxzr vy|-��©jvy|~t�ª�{Q}�����xz�~{9��� ©jvy|~t�ª�{�}IÞyèQÝQå�±P3µvy�zq���xzx�Þwè10Qå9±jà¬�zr�t�}�}�{b| vy|~� â �p�~�Þwè10Qå�± À {��p�5�zão��ru�,��x¶vyq ¡OÞwè�098:�®�p�prsqu�p}¤v�´��z��q�{w��}zn¦vbt��o�¨{��¶�o���z��r�xzr�v�|I��v�xz�zr�t���}�±9���px¥����twvy�~}��xz�~�Z}zn���t�xz�§vyqQxz�§v�|~}��¨{��z��vyxzru{b|!ru}�|~{�|����o���z��rsxzr�v�|e±�xz�~�¬����xz�~{9��twvy|p|~{bx����¬�~}����!�¨{b�c�4�£?(@)�B� C�E �rs|�r�x�}�}�x§v�|~��v����(�¨{b�z�Ë¡ý¹`r�|�t��Y�4�µ?�@)�B� C�E �^ru}�}­��qs�¶��vb���z{br�|`xC²¶rsxz������}zn���t�x�x�{�xz�~�O�^r�|�|~���np��{9�9��t�x�±jxz�~������xz�~{9�Ttwvy|T�����~}����T²¤�~��|\� {b�zxz�~{��Q{b|¦v�qsr�ªyvyxzr�{�|�r�}B����npq�{w������¡�� ¹9����� À {��p�5�ão��r��Ë��x³vyq ¡OÞwè�098:�¬�¨{��£xz�~���p��x§vyr�q�}�¡·�7mo|¸v�q��Q{��zrsxz�p�r�}¶�br��Q��|�����q�{w²�¡f�×�ø�Æ~Ç�ÙsÖ�ï�ò îV¸ ì Î�hOðZÃ`ÕbÖyÇyØ~׬ÖwÇ�Ø~ÛjÉ�ñ Æ~Ç�òDØ9Ö�Ù�Æ¦Û ócØ~ÛjÕQôQÆ~É�ÑÞ�¡^/or��Q��|OG E uxw y }��~t���xz�¦vyx£j$G E j$p� Þ�¡©O��xn¹+ºL�ýå�v�|��]G�ºL�ýå9¡ç9¡ ·~{b��f � Þ��^¨�¨S¨��U»ç�¡�ÞQ¡ ¢ �§vy|~}z�¨{b�z��v�xzr�{�|7¼c½§e��¾�4��?Z@)�B�DCFE,��G e�¡ç�¡ ç�¡eßc{b��np�px§vyxzru{b|Ë{b��t�{9�$¿,t�r���|�x�¼�ÀFe��¯G9he ��½§eb¡ç�¡·[�¡��¿{��zxz�~{Q�Q{�|¦vyqsruªwvyxzru{b|x¼HÁ,eL�¯½§eH?;À�eWG$eH?¹�e CFE G e CFE ¡ç�¡Â<�¡eßc{b��np�px§vyxzru{b|Ë{b��t�{9�$¿,t�r���|�x�¼o¹�eL�²j`Á,e�j pµ¡ç�¡ Ý�¡ À {��z��vyq�r�ªwv�xzr�{�|7¼cG e © E �ÃÁDe^�~¹�e�¡àµ|��9�¨{��

·p��{b�¹`x���n~}¥ç�¡�ÞWÄpç9¡ Ýorsxc�¨{bq�q�{w²�}¬xz�¦vyxHÅ �D© E �¾ÆÂG E �^¨~¨^¨��,G �D© EÈÇ �¨{b�z��}¥v�� {b�zxz�~{Q��{b|¦vyqp�¦vb}zr�}«{b��xz�~�´&�z�`q�{w��}�n¦v�t��J� ç�¡�Þ��§¡ ¯ |�� vbt�x�±9¹`x���nËç9¡·[&ru}¥vJ/o�§v��2�Ϲ9t��p��r��9x¥{b�zxz�~{��Q{b|¦v�qsr�ªyvyxzr�{�|�}zx���ne±�²¶�~������ÀFevy|~�x¹�e CFE v�����/o�§vy�2��¹9t­�p��ru��xBt�{9� ¿�t�r���|`x�}�¡ ¢ �pr�}B{b�zxz�~{��Q{b|¦v�qsr�ªyvyxzr�{�|�}zx���nEtyvy|�q�{�}­��}�x§v��prsqsr�x �

ç

Page 5: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

vy|~�Cxz�~�������¨{b���o{�|~�¤{���x���|�����npqu{���}Zn�v��zxzr�v�q�v�|~���§Á�x����z|¦vyq¦}���q���t�xzrs���£n¦v��zxzr�vyq~����{��zxz��{��Q{�|�v�qsruªwvyxzru{b|�*/o�zrs����}c��x¥vyq ¡~ÞyèQè=<>�j{b�c������|����pqsq~����{��zxz��{��Q{�|�v�qsruªwvyxzru{b|É� á!v�|pr���q ±�/o�§vb���p±9´�vy�p����vy|,vy|~��¹`x���²L�vy�zx³Þwè>8�Ü9±Q¹9{b����|�}���|¸ÞyèQè�ç��§¡eà¬q�rs��r�|¦vyxzr�{�|�{b�)½§eov�|~�!ÁDe«����{b� moq���{b�zr�xz��� ç9¡�Þ£qu�wvb�p}¬x�{�xz�~�¤xz�p�����x����z�?����t��p�z����|~t�������q vyxzr�{�|

�4�¯?Z@)�B� CFE ��G e��¯G e © E ¹�e�³YG$e^ÀFec³¥G e C�E ¹�e CFE ¨ � ç9¡ ç��ßc{bqsqu��t�xzr�|��7xz���!xz�p�����!x����z� ����t��p�z����|~t��!����q vyxzr�{�|~}£�¨{��cf � Þ=�^¨^¨~¨^�U»�q��wv���}¥x�{

�4�Q?Y@)�B� C�E �BÅ � �ÊÅ �D© EÌË � � � ç9¡Í[��²¶�������

Ë � �ÎÏÏÏÏÏÏÏÐ

À E ¹ E¹ E ¡ ¡ ¡ ¡ ¡ ¡¡ ¡ ¡ ¡ ¡ ¡ ¹ � C�E¹ � CFE À �¹ �

ÑÓÒÒÒÒÒÒÒÔ

r�},vx»!³2Þ�Õ¥»�xz�zr��9r vb�Q{�|¦vyq¥��vyxz�zrsÁ�±µ�¨{��z�����(�`�0xz�~�It�{9� ¿�t�r���|`x�} ÀFe�±T¹�ew¡ À {bx���xz�¦vyx Ë � �ÅBh�D© E �{�4��?�@)�B�DCFE,�BÅ � ¡S©���x Ë � ����xz�~�7»�Õ\»Tq��wv���r�|~��}z�p�p��vyxz�zrsÁ»{b� Ë � ¡ ¢ ����|Svy|~{bxz�~���|~{bx§v�xzr�{b|J�¨{b�¥xz�~�!����t��p�z����|~t�������q�v�xzr�{�|Iru}�4��?�@P�J� CFE �BÅ � ��Å � Ë � ³ZG �D© E ¹ � d h� ¨

¢ �~�A�{��zxz��{��Q{�|�v�qOnp��{:�z��t�xzr�{b|D{b���4�Ê?¥@)�B�DC�ED������±=�¸{b|`x�{�xz�~��´&�z�`qu{w��}z�p��}zn¦vbt��³np��{=��9�~t���}�v�npn���{�Á�r���vyx�����ru�Q��|pn¦vyr���}�±�tyvyqsqu��� â rsx�ª�n¦vyr���}�±µv�|~�Etwv�|E����t�{���np�px����(v�}��¨{�qsq�{w²�}�¡°©O��x�4±��,Ö>�Z}§vyxzru}z��� Ë � Ö{�¯±:Ö�±�xz�~��|�4±��,iM�®²¶rsxz�!i.��Å � Ö�ru}Ztwv�qsq�����v â rsx�ª¶n�v�rs�«{b�T�4�×?É@)�J�DCFE5��²¤r�xz�����}zr��9�¦vyq

Ø � �4��?Z@)�B� CFE ��iB?¥±:i � ç9¡Ó<>�� ¹ � G �D© E d h� Övy|~������}zr��9�¦vyqO|~{b�z�ÚÙ{�²j Ø j p�;¹ ��Û d�h� Ö Û ¡¯ �¬¹�e¥� å»vyx˹`x���n ç�¡Â<E{��Bmoq���{b�zrsxz�p� ç�¡�Þ�±¤xz�~��|M}�n¦v�|-«`G E �~¨^¨^¨��,G e ­ �¨{��z��}�v�| r�|`��vy�zr vy|`x}z�p��}zn¦vbt���{����zÞQ¡�Þ��§¡ ¢ ��ru}o���wv�|~}oxz�¦vyxoxz�~�7t�{b�z����}zn�{b|~�9rs|~� â r�x�ª¾n¦vyr���}!vy���B��r�����|pn¦vyrs��}�¡ ¢ ��ru}or�}�����z���p|pqsrsº���q���x�{J{9t�t��p�or�|¸n��§v�t�xzr�t��Cv�|���xz�~�������¨{b����±�²¥�Cv�}�}��p����xz�¦vyxI¹ E �^¨^¨^¨~�ȹ �£Ü å9¡ ¯ |¸xz�pr�}twv�}­�Q± Ë � ru}¤tyvyqsqu���I�p|p�����9�~t������*/�{�qs�p�¸v�|~�.K¥v�|�©O{�vy|�Þwè�èQÜ9±`n¦vb�Q��[=<`Ü1�cvy|~�Iru}¤{��e���pq�qO�§vy|pº¦¡ÝUÞß=à5á�â^à£ãÌß�ä=å�ß:æ4áÉçWß:Þ`èUéMêWæ4ß:ÞëDâ�ã"ì�âUâ^Þ]í�æ4ãïî�ð�ä:æ¶àWê�è$ß:à �4�Q?Z@P�J�DCFE,� ä:Þ-á �����v�F�¬�Pñ ìcâá�â^Þ-ß�ãÌâ ë$ò ± ä í�æ4ãïîAå�ä�óÍé�âHè$ß:à �4�Ê?Z@P�J�DCFE,� ä:Þ-áÉë$òBô¥ä.í�æ4ãõîAå�ä1ó·é�âHè$ß=à ���7�v����� á1âÌö§Þ-âWáä:êLô �Ã@!³¥±9C�E�÷!ø�ù â2å�ä�óÍé�â^êHôYä=à5â2ä1ó·ê ß�çWä�ó¶óúâWá ù ä:àWû2ß=Þ�æ4ç�í�æ4ãïîAå�ä�óÍé�â^ê(ìTæ4ã ù ãÌä:à�ü�â�ã @vý4þ ä:æSü�â ñþ ä=à$óúâ�ã�ã�ä=ÞVá å�ä=Þá�â^àOÿ�ß=àUê ã�������� ñ�� óúâ^æ �ð�â^Þ7ä:Þ-á å�ä:Þ�á�â^àOÿ�ß=àUê ã������� ñ í�é ù â ��������� ÷�� ò � ìcâá�â^Þ-ß�ãÌâAâ^æSü�â^ÞFå�ä�óÍé�â^ê�ß,è �����k�F�¬�§÷� ���¸� )c¡�ëZ��#�ê7¡£¢��¥¡¾êL¦o#�� % �oë§�¸#� ¢ �����§v�xzr�{�|�v�q!©jvy|~t�ª�{Q}�����xz�~{9��r�}Jvy| �§Á�x���|~}zr�{�|ý{��!xz�~�D}zn���t�xz�§vyq�xz�§vy|~}z�¨{b�z��v�xzr�{�|�©jvy|~t�ª�{Q}����xz�~{9��±�xz�¦v�xZ��vyº���}¥v�t­�¦v�|����¤{���n�{bq��I@Ix�{���uO°?n�{�}­}zr���qu�£²¶rsxz�~{���x«����}zx§vy�zxzrs|~��xz�~�o����xz�~{9��¡

[

Page 6: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

¯ xµr�}«v&}�n���t�r vyq9twv�}��c{b�~xz���¥�§vyxzru{b|¦vyq�´��z��q�{w�¾�,��xz�~{9��±b²¶�pr�t���r�}µ�§Á�x���|~}zrs����qs�7�9r�}�t��~}�}����¾�`��� â �p�~�Þwè10:<�± â �p�~��Þwè�è:<p± â �p�~��Þwè�è�0��§¡bZ���B��q�xzrsnpqs��rs|~��� ç9¡Í[��c�`���Q?Y@)�vy|~������{b����vy|pr�}zrs|~�Bxz���!x����z�,}�±�²¥�!�¦v��Q�

�(Å �D© EÈË � ���BÅ �D© E �Ìg � ³¥@ Ë � ��¨ �Ì[9¡�Þ��¢ ��ru}�r�}�xz�����§v�xzr�{�|�v�q!´&�z�`qu{��4����t��p�z����|~t�������q�v�xzr�{�| {�� â �p�~�¥�zÞwè10:<>�§¡ � ¢ ���¸t�{bqs�p��|~}�{b� Ë ��¨{bq�q�{w²S����{b� xz�~�!{b�zxz�~{��Q{b|¦v�qsr�ªyvyxzr�{�|�np��{9t���}�}�±�²¶�prsq���xz���!t�{bqs�p��|~}¥{b�og � vy���oxz�~��� t�{b|`xzr�|`�¦vyxzru{b|P´����t�x�{b��}�¡Í��¹`��npn�{Q}��&²¥�!²¶vy|`x¶x�{�t­�¦v�|����&xz�~�³n�{bqu�&x�{����¯@�¡ â ��²¶�zrsx�� �Ì[�¡�Þ��¥v�}

�4�Q?����J�5Å �D© E5Ë � ���BÅ �D© E �Ìg � ³¯�4@�?���� Ë � ��¨ �Ì[9¡ ç��©O��x�� � �Êg � ³¯�4@�?���� Ë � v�|��Jt�{�|�}zru�����£xz�~�7æ â �p��t�{���n�{�}�r�xzr�{�|� � ��� � ���"! � å$#²¶�������%�Íru}Cv]» ³MÞ{Õ×».³MÞC�p|prsx§v��z����vyxz�zrsÁEvy|~�&� ru}³�pnpn�����xz�zr vy|~�b�pq vy��¡ ¯ x¾ru}�r���n�{��zx§vy|`xx�{�|~{bx���xz�¦v�x'� � r�}!�p|��������~t�����xz�zr��9r vb�Q{�|¦vyqcv�|~�Dxz�`�~}B{b�¥���pqsqZ�§v�|pº¦¡ ¢ �pru}�r���npqsr���}�xz�¦vyx�� r�}rs|��Q���zxzr��pq���¡4©O��x)( �D© E � Å �D© E �Gvy|~��( � ����|�{�x��,xz�~�!#p��}�x.»Et�{�qs�p��|~}C{��*( �D© E ¡ ¢ �~��|(�`��B��q�xzrsnpqs��rs|~�7�Ì[9¡ ç1�c{b|�xz�~�!q�����x£���×�4�¯?+���B�DC�Eov�|~�I{�|�xz�~�!�zr��b��x¤�`�,�BCFE�±�²¥�&��v����

�4�Q?����B� CFE ��( � � ( �D© E.- �- � � � h Ë � � CFE �Ì[9¡Í[��( h�D© E ��( �D© E � g( h�D© E � �4�Ê?����B� CFE ��( � � - � ¨¢ ����}­��������v�xzr�{�|�}³����|�{�x��CvC����q�vyxzru{b|¸}zr���rsq vy�ox�{�xz�~��©ev�|~t�ª�{Q}¥����t��p�z����|~t��B����q�v�xzr�{�|O¡ ¢ �~���Ivyrs|�9r�R�������|~t��!r�}¶xz��v�x - � ru}£|~{bx³vCxz�zr��9r vb�Q{�|¦vyq���v�xz�zr Á�¡ À {bx��&xz�¦vyx - � r�}³�o���z��r�xzr�vy|°}zrs|~t�� - � �( h� �{�4�Q?����J� CFE �/( � ¡ ¢ �~��´&�z�`q�{w�J}z�p��}zn¦vbt��!twvy|J|~{�²M���!�§Á�n¦v�|��p���J�`�Ëvb�p�9rs|~��4�Q?����B� CFE ��½ �D© E �~�5�4�Ê?����B� CFE �B� m ½ �D© E �^¨S¨�¨x�{ â vy|~�Q�=�0( � �§±�v�}c�¨{b�¶xz���B©jvy|~t�ª�{Q}c����xz�~{9��¡¢ �~�Zxz�zr��9r�v��Q{�|�v�qsruªwvyxzru{b|B{b� - � ru}®}zxz�§vyru�b�`xz�¨{��z²¶vy���B�`��xz���¥v�npn�q�r�twv�xzr�{�|C{���{��zxz�~{Q�Q{�|¦vyq�xz�§vy|~}5��¨{b�z�Ivyxzr�{�|~}£{�|���{bxz�J}zr��p��}¥{�� - � vy|~��{b|�xz�~�&�zr����`x¶{b�1( �D© E ¡�� ¹����oxz�~�!��ru}�t��~}­}zru{b|J{b|�xz�~��æ â����xz�~{9��r�|Y�*/&{bq��p��vy|~��K¤vy|¸©O{�vy|TÞwèQè�Ü��¥vy|~�Ëmonpn���|~�9r Á�m(�§¡ ¯ |°xz��v�x&twv�}���±��Ì[9¡Í[��¶twvy|Ë���wvyq�qs�����t�{�|~}�ru�p�������Ëvb}�v�©jvy|~t�ª�{�}«xz�p����� � x����z�?����t����z����|�t�������q vyxzru{b|e¡% �¬|~{w²E�¦v��Q�¶v¤np��{9t������p���Z�¨{��¬t��¦vy|~�br�|��¶xz�~�«n�{bqu�c{b�9xz�~�¥}�n���t�xz�§v�q`xz�§vy|~}z�¨{b�z��v�xzr�{�|�©jvy|~t�ª�{Q}����xz�~{9�B²¤r�xz��{��pxZv&t�{���npq���x��Z����}�x§v��zx�¡ À {bx��Zxz�¦vyxµrs��npq�r�t�r�x®����}zx§v��zx�}�� ¹9{b����|~}���|ËÞwèQè�ç��§±y�p�$�¦v�xzr�{b|vy|~��np�p���br�|~�x� ©O���~{b�~t��Ëvy|~��¹9{b����|�}���|»ÞwèQè�Ü��§±��p������q�{�n������¨{��&xz�~��©jvy|~t�ª�{Q}¤�,��xz�~{9��±Ovy����|~{w²n�{�}�}�r��pq��&rs|°t�{b�B�prs|¦vyxzru{b|J²¤r�xz�¸vCt��¦vy|~�Q�!{��jn�{bqu��¡ ¢ �pr�}�vynpnp��{`v�t­�Jr�}¶��vyxz�~���Ivyxzr�tyvyqsq��J���Q�prs��v��q���|`x�±¦�p�px��Q�prsx����9rSR�������|�x&rs|�np�§vbt�xzrut���±�����{b�Gxz���Brs��npq�r�t�r�xzqs�7#pqsx����������§vyxzr�{�|¦vyq«´��z��q�{w������xz��{9�� á&��¹pv��¾�pq�v�|�Á�±32������z����������|�vy|~�7bZ�pqsxz�~����qoÞwè�è>8=�«xz�¦vyx³²¥{b�zº9}&r��������9r�v�x���qs�D{b|�xz�~�C�§v�xzr�{b|¦v�q´&�z�`q�{w������t��p�z����|~t��!����q vyxzru{b|Y�Ì[�¡�Þ`�§¡�m�|¸vyqu�Q{��zrsxz�p��¨{b�¶xz�~��t��¦vy|~�Q�!{��jn�{bqu�&r�}¶�br��Q��|Ë����qu{w²�¡f�×�ø�Æ~Ç�ÙsÖ�ï�ò õ+¸ ì Î�lCïjØ~Ûeø~Ã�Æ~ñ ÷ Ʀ×�Ã~ÑÞ�¡ ·~{b�z�4� � ��g � ³¯�4@�?���� Ë � ¡

<

Page 7: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

ç9¡ ·¦vbt�x�{b�zr�ª��5� � �6� ! � å # ¡[9¡�ßc{b��n��px�� - � ���¬h Ë � �BCFEc� E7 C98�: g¬?;��h<! �BCFEå #�= ¡<p¡ m�n�npq�� xz�~��{b�zxz�~{��Q{b|¦v�qoxz�§vy|~}z�¨{b�z�Ivyxzr�{�|~}?>*vy|~� @> � ! > å

å ÞA# }�{�xz�¦v�x - � ¼ú�@>�h - � > r�}¶xz�zr��9r�v��Q{b|¦v�q ¡<p¡�½¶n��~vyx��5( �D© E �ÊÅ �D© E � @>7¡% �!�p�$#p|~��xz�~���§vyxzr�{�|¦vyqj©jvy|~t�ª�{Q}¥����xz�~{9�I�`�°v�t�{��¾�prs|¦v�xzr�{b|°{b�µmoq���{b�zr�xz���,}&ç9¡�ÞBvy|~��[�¡�ÞQ±r ¡ ��¡¬vy��x�����¹�x���nCç�¡ Ýcrs|Bmoqu�Q{b�zr�xz�p� ç9¡�Þ�±yv¶t­�¦v�|����¬{��`n�{�q��¬ru}�n�{Q}�}zrs�pqu��¡emo|7vyqu�Q{��zrsxz�p�ý�¨{���v£}zn���t�r�#~tvynpnpqsrutwvyxzru{b|�r�}���rs����|�r�|Em�q��Q{��zrsxz�p��Ü�¡�Þ�¡ ¢ �~���§v�xzr�{b|¦v�qc©jvy|~t�ª�{�}&����xz�~{9�¸r�}7vy|T�§Á�x���|�}zru{b|0{b�xz�~��©jvy|~t�ª�{�}¬����xz�~{9��±�}�{�r�x¥rs|p�~���zrsx�}c��v�|`��np��{bn����zxzr���}¥{b�Oxz�~�oq vyxzx�����¡ ¢ ���!{�|pqs���9r�R�������|~t��oq�r���}rs|Ëxz�~��t��¦vy|~�Q�!{��jn�{bqu��¡ ¢ �pru}¤�¦vb}¥rs��npq�r�twv�xzr�{b|~}¶{b|Jxz�~�¾´&�z�`qu{w�I}z�p�~}zn�v�t��£¼iBïeÃ`Æ~Ç�Ã9ò õ+¸ ìCB â�ã Å �D© E ä:Þ-á Ë � ëDâ¥çWß:ûIðVé�ãÌâUá�ë$ò¯ã ù âZê�ð�âWç�ã"à,ä�ó�ã"à,ä:ÞFê¶è$ß=àUû{ä�ã4æ4ß:Þ B ä=ÞVçWî�ß:êû2â�ã ù ß�á ì§æ4ã ù ð�ß�óúâ @ ä:Þ-á!ó â�ã ( �D© E ä:Þ-á - � ëDâAçWß:ûIðVé�ãÌâUá�ë$ò�D�óÂü�ß:àWæ4ã ù ûFE ÷ � ñ ã ù â^Þ

â vy|~���:�0( �D© E � � â v�|����:�ÌÅ �D© E ��� �Ì[9¡Ó<>�â vy|~���:�0( � � � â v�|����:�5�4�Q?¥@)�B� C�E �4�Q?����J�5Å � �(¨ �Ì[9¡ Ý��

÷ Ç�Æ�Æ~ñ*¹�r�|~t��*( �D© E �ÊÅ �D© E � @>7±-�Ì[�¡Â<>���¨{bqsqu{w²o}�¡®àZ�Q�¦v�xzr�{b|�Ì[�¡ Ý���r�}c}z�~{�²¶|Jvb}Z�¨{�qsqu{�²�}�¡¬¹`rs|~t��( �D© E �¾Å �D© E � @>¿vy|~�G� � �H�5� ���I� @>��$� @>¬hJ� �§±�²¥���¦v��Q�A( �D© E @>�h.� �K( � >�h.���vÅ �D© E � � ¡·~{bqsqu{w²¤r�|��7xz�����p�$#p|prsxzru{b|�{��L� � vy|~�×� ç�¡·[1�§±�²¥�!��v����Å �D© E � � � Å � ³��4@�?����5Å �D© E*Ë �� Å � ³��4@�?����$�4�Ê?Y@)�J� CFE �JÅ �� �4��?�@P�J� CFE �4��?+���B�5Å � ¨¢ ���!np��{9{b�j{����Ì[�¡ Ý1�µ�¨{�qsqu{�²�}c����{b�?xz�~��{b�~}����z��vyxzru{b|Jxz��v�x â vy|~�Q�=�0( � ��� â vy|~�Q�=�ÌÅ �D© E � � �§¡ M¢ ����t­�¦v�|����¤{���n�{bq��¥r�}Z�����z��}�r���rsq vy�Zx�{Bv�|�rs��npqsrut�rsx«����}zx§vy�zxZr�|�xz���³©jvy|~t�ª�{�}®����xz�~{9��� ¹9{b����|�}���|Þwè�èQç1�§¡ ¯ |(xz�pr�}7twvb}���±¬v�©jv�|�t�ª�{Q}�np��{9t���}�}!r�}¾r���npqsrut�rsxzq��(v�n�npq�r����Tx�{�xz����|~��²}zx§v��zxzrs|~���Q��t�x�{b�½ E �N( � d E ±�r ¡ �Q¡�²¤r�xz��{��px!�§Á�npqsrut�rsxzq��¸t�{���np��xzr�|~��½Hm`�^¨^¨^¨��,½ �D© E ±jvy|~��xz�~��´&�z�`qu{��°}z���~}zn¦vbt���r�}#pqsx��������°����v��Ivyxz�zr ÁDn�{bq��`|~{b��r vyqH�02������z�������Q��|Tv�|~��¹`n���|~t��IÞyèQè>8�±�©�����{��~t���ÞyèQèQÜ�±�©O���~{b�~t��Þwè�èQè1�§¡T·~{bq�q�{w²¶rs|~�;�Ì[9¡ Ý1�§±µxz�~����r��Q��|`��vyq�����}C|~�wvy�{@�vy���J��|p�¦vy|~t����0rs| â v�|~�Q�:�0( � ��vy|~�Exz�~{Q}��|~�wvy���Tvy������v���n����D�`�¸xz�~�����p|~t�xzr�{�| �*��?O���D�M�*�!?®@)�§¡O� À {bx��Cxz�¦vyxBmoqu�Q{��zrsxz�p� [9¡�Þ�twv�|����t�{�|�}zru���������Jvb}¥v�xz�§v�|~}�r�xzr�{�|�{����¦vy�z��{�|�rut â r�x�ª&��vyq�����}«²¶rsxz�Jx§vy������xH@�x�{B��v��z��{b|prut â rsx�ª���vyqs�~��}²¶rsxz��x§v����Q��x��x� ¹`qu��r �Ïn���|¸v�|�����vy|°�p����K£{b��}zx�ÞwèQè�Ü��§¡Í�% �~��|�xz�~�¤n�{bq��¥r�}ct��~{Q}���|����Q�¦vyq~x�{Cv!�¦vy�z��{�|�rut â rsx�ª¤��v�qs�~��±`xz�����$R���t�xc{b�ev�t��¦vy|~�Q�¶{b��n�{bqu�r�}¶�br��Q��|��`�Jxz�~�³�¨{�qsqu{�²¶rs|~�Bxz����{b�����Ë¡iBïeÃ`Æ~Ç�Ã9ò õ+¸�îKB â�ã Å �D© E ä:Þ-á Ë � ëDâ¥çWß:ûIðVé�ãÌâUá�ë$ò¯ã ù âZê�ð�âWç�ã"à,ä�ó�ã"à,ä:ÞFê¶è$ß=àUû{ä�ã4æ4ß:Þ B ä=ÞVçWî�ß:êû2â�ã ù ß�áZìTæ4ã ù ð�ß�óúâ @�÷ B â�ãLô E ëDâ7ä ù ä=àUû{ß:ÞFæ4çÉí�æ4ãïî�å�ä�óÍé�â ÷QP&ù ß�ß=ê$âã ù â�Þ-â^ì�ð�ß1ó â �k� ô E ÷B â�ã ( �D© E ä=ÞVá - � ëDâ2çWß:ûIðVé�ãÌâUá�ë$òAD�ó ü�ß=àUæ4ã ù û4E ÷ � ñ ã ù â^Þ�ä,è$ã�â^àBã ù â2ç ù ä:Þ�ü�â2ß5è�ð�ß1ó â ñ ã ù â ù ä:àJRû2ß=Þ�æ4ç.í�æ4ãõîÉå�ä1ó·é�â£ô E á=æ¶ê$äWð�ð�âWä:àWê�ä=ÞVá @ ë,âWçUß=û2â^ê ä ù ä=àUû{ß:ÞFæ4ç�í�æ4ãõî�å�ä1ó·é�â�ìTæ4ã ù í�æ4ãïî�å�âWç�ã�ß=àG �D© E ÷2ø�ù â2ß1ã ù â^à ù ä=àUû{ß:ÞFæ4ç�í�æ4ãõîJå�ä1ó·é�â^ê¬ä:Þ-áAí�æ4ãïîJå�âWç�ã�ß=àUê(à5â^û{ä:æ¶Þ7éMÞ-ç ù ä=ÞMü�âWá ÷

Ý

Page 8: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

÷ Ç�Æ�Æ~ñ ·er���}�x¶|~{bx��&xz�¦vyx�S_¼ú� â v�|����:�0( � ��� â vy|~���:�5�4� ?×@P�J�DCFE �4�¯?&���J�UT � �Z²¶�~�����VT �r�}¥xz�~�!��vyxz�zrsÁË{�� â rsx�ª!����t�x�{b��}�¡©O��x£�4± E �D� E �¤²¶rsxz�± E � � ô E ?;@)� CFE vy|~�7� E �²Å � Ö E ����v â rsx�ª¾n¦v�rs��v�|�� Ø E �¾G �D© E ¹ � d h� Ö Exz�~��t�{��z����}zn�{b|~�9rs|~�¸����}zr��9�¦vyq ¡ % �Ë}��¦v�qsq¶}z��{w²?xz�¦vyx7�5�4@Z?K��� C�E � Ø E �Cr�}�v â r�x�ª�n¦vyr����¨{b��xz�~�n�{�q�� �¦¡�·er���}zx�±j²¥��np��{w����xz�¦v�x Ø E uWS°¡Ëá&� #�|~� i�� �4�k?¯@P�J� CFE �4�¾?W���B�È� E uXS°¡�¹`rs|~t���4�Q?Y@)�B� C�E ��� E �¯± E � E ³ Ø E ±�²¥�!�¦v��Q�i E � � E ³¯�4@�?����ȱ E � E ³¯�4@É?+��� Ø E� ô E ?��ô E ?Z@ � E ³¯�4@É?Y��� Ø E �¾�4@É?+��� Ø E ¨% �&|~{w²ýnp��{w�Q��xz�¦vyxB�5�4@]?Y��� CFE � Ø E �cr�}�v â r�x�ª!n�v�rs��{b�H�4�Q?¥@)�B� C�E � vy��x����¤xz�~��t­�¦v�|�����{b�n�{�q���¡ ¢ �~�������¨{b���!²¥��}z�~{�² xz�¦vyx

�4�Q?����J� CFE � Ø E ? Þ@�?�� Ø E[Z p,S ¨¢ ��ru}£ru}o������r���vyqu��|`xox�{

� �4��?Z@)�B� CFE �4�¯?+���B�5Å � � h � : �4�¯?+���B� C�E � Ø E ? Þ@�?�� Ø E = �ýå�¨à¬q vy��{��§vyxzr�|��7xz���!qu����x5� �¦v�|��J}zr��p���br��Q��}

Å h� �{�4�Ê?Y@)�B� CFE � Ø E ? Þ@�?�� Å h� � Ø E ? @É?+�@É?+� Å h� � �4�Ê?�@P�J� CFE � Ø E²¶��rut­�Ër�}¶ª�����{�}zrs|~t���ÅBh� � Ø E �ýå�¡©O����� ô �� ô E �,�V�¶���7vy|~{bxz�~������v��z��{b|prut â rsx�ª¾n¦v�rs�o�¨{��&n�{�q���@�¡ % ��|~{�² np��{w���¾xz�¦vyx2� ô �,�V�r�}�vC�¦vy�z��{�|pr�t â r�x�ª!n�v�rs�¥�¨{b�¶n�{�q��5�¦¡c©���x�i{�¾�4�Q?¥@)�B� C�E �4� ?����J�È�xu,S°¡ ¢ �~��|

i{� ô ?��ô ?�@ �{³¯�4@�?���� Ø ¨¹`rs|~t�� Ø �¯G �D© E ¹ � d�h� Ö�vy|~��i�q�r��!rs|\S�±�v�q�}�{A�Ëqsr���}£r�|]S°¡ % �!}zxzrsq�qO��v����!x�{�}��~{w²ýxz�¦v�x�4� ?����B� CFE ���!? Þô ?�� � Z p,S ¨¢ ��ru}£ru}o������r���vyqu��|`xox�{

� �4�Q?¥@)�B� C�E �4� ?����J�5Å � � h � : �4�¯?+���B� C�E �¬��? Þô ?�� � = �ýå�¨â ��twvyq�qsrs|~��xz�¦vyx(�4�Z?O@P�J�DCFE5�¬�]� Þ:�9� ô ?O@P�È�(³ Ø ±`��q vy��{��§vyxzru{b|�{���xz�~�ox����z��}crs|�xz�~�oqu����x5� �¦v�|��}zr��p�!�br��Q��}

Å h� �{�4��?Z@)�B� CFE ��� ? Þô ?�� Å h� ��� ? @�?��ô ?�� Å h� � �4�Ê?Y@)�B� C�E �¬�� ô ?Y@ô ?+� Å h� �{�4��?Z@)�B� CFE ��� ? Þô ?¥@ Å h� �¬�

� ô ?Z@ô ?�� : Þô ?¥@ Å h� ��� ³¥Å h� � Ø = ? Þô ?�� Å h� �¬���²¶��rut­�Ër�}¶ª�����{�}zrs|~t���ÅBh� � Ø �ýå�¡ M

Ü

Page 9: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

Þwå ºÞwå C mÞwå9C�^

_ `a bc _d _e fg h ijå ç <k h _ f bl` dnmpoqdsrUt1u h j g�awv

·jr��b�p���][9¡�Þ=¼ â ��}zr��9�¦vyq¥|~{b�z��}B�Q����}��~}C�¦v��z��{b|pr�t â r�x�ª,��v�qs�~��}�¡ ¢ ���J}��Q�¦vy����}7}z��{w² xz�~��n¦vyrs��}� ô �,Ù����¨}­���¬� ç�¡Â<>�5���¨{b�c<�å�©ev�|~t�ª�{Q}µ}zx���n~}¬²¤r�xz��n�{bqu��@O�k?CÞ�¡ ¢ ���¥�p�pqsq���x�}c}��~{w²�xz�~�¤n¦vyr���}�� ô �,ÙM�vy��x����ot��¦vy|~�br�|��Bxz�~�³n�{bqu�&rs|`x�{ �!�ýç�¡¢ �pr�}¥r�}¥|~{�²Mrsqsq��~}�xz�§v�x����I����vC|����,���zr�twv�qj�§Á�v���npq���¡x��OØ~ò�ð¬×�à õ�¸ ì ¢ �~���Ivyxz�zr Á¯�Gr�}����wv�q&}z�`������xz�zr�t�vy|~�4xz�zr��9r vb�Q{�|¦vyq¶²¤r�xz�ýç�´ }�{�| xz�~���Ivyrs|

�9r�v��Q{b|¦v�qov�|~�Ê?BÞ�´ }�{�|(xz�~��x ²¥{T{:R-� �9r�v��Q{�|�v�q�}�vy|~�;��r�}�xz�~��r��p��|�xzrsx �»��vyxz�zrsÁ�¡Ãb«{�xz�4�¦v�����9rs����|�}zru{b|MçQåQå×Õ çQå�å9¡ ¢ �~�°����}��pq�x�}�{b��xz�pr�}Ë�§Á�v���npq���vy��������|~���§vyx������~}zrs|~��2¸vyxzq vy� {b|Mvá&à¥ß�moqsnp�¦v9¡ % �on����z�¨{��z�����]»!�¯<`å�}zx���n~}¤{��Oxz�~�!}�n���t�xz�§v�q�xz�§vy|~}z�¨{b�z��v�xzr�{�|�©jv�|�t�ª�{Q}Z����xz��{9�²¶rsxz��@x�k?CÞ�vy|~�Ëvy|Jrs|prsxzr vyqO����t�x�{b�¶²¶rsxz�����Q�¦vyqj��qu������|`x�}�¡

¢ �~�¬#p���C}z��v�qsq���}�x��¦vy�z��{�|�rut â rsx�ªB��vyqs�~��}�vy��� ô E<y çM¨Ó<�<`è=<oÞwå9C�^ ± ô m y çM¨ ç:<`Ü�å¤Þwå9C{z­± ô z yÜM¨Ó<�ÞyÝQè¤Þyå>C9z�± ô ^|y Þ=¨ çQè:<�8£Þwå9C m ±�vy|~� ô�} y ç�¨�Þwè10Qè¤Þyå>C m ¡ ¢ ��� â rsx�ª!��vyq�����t�{��z����}zn�{�|��9r�|��Bx�{ ô Er�}L± E[y åM¨ èQèQè10¾v�|~�,�¦vb}¥����}zr��9�¦vyqO|~{��z� ç�¨ ܤÞyå>C�^�¡m���x�����vyq�x����zrs|~�Jxz���7n�{bqu�¾r�|`x�{]�7� ô E ±¦xz�~�A#p����}z��vyq�q���}zx&�¦vy�z�,{b|pr�t â rsx�ªB��vyqs�~��}Bv���� @ô EVy?CÞ=¨ å�åQå�å9± @ô m y ç�¨ ç=<`ÜQå¶Þwå>C9z�± @ô z y Ü�¨Â<�ÞyÝQè¤Þyå>C9zz± @ô ^%y Þ�¨ çQè=<M8£Þyå>C m ±Zvy|~� @ô�} y çM¨�Þwè�0�è¤Þwå9C m ¡

¢ ����}z��v�qsq���}�x¶��vyq��~� ô E ��çM¨Â<1<`ç�è¤Þwå9C�^c�¦vb}��9r�}§vynpn��yvy������¡¤¹`rs|~t����Iru}&v â rsx�ª���vyqs�~��±�xz�~��t�{b�z��� �}zn�{�|��9r�|�� â r�x�ª¬����t�x�{b�Oru})#�q�x��������B{b�pxj{b��xz���£©ev�|~t�ª�{Q}��¦v�}�ru}�²¤r�xz��xz���§#�q�x����H�4�£?�@)�B�DCFE �4�£?5���B�§¡·~{bqsqu{w²¤r�|�� ¢ �~��{������Ú[9¡ ç9±�xz�~�!�¦vy�z��{�|�rut â rsx�ª&��v�qs�~� ô E r�}¶��vynpn����I{�|`x�{ @��¾?BÞQ¡m���x����o{�|~�¾v��p��r�xzr�{�|�v�q�©jvy|~t�ª�{�}£}zx���n˲¤r�xz��n�{�q����¦±�xz���!qu{Q}zx â rsx�ª���vyq����!ru}¤����t�{w����������±�}zrs|~t��{b|~�B©jvy|~t�ª�{�}¤}zx���n~}!v��p��}I½ �D© E x�{�xz�~�B}�n¦v�t��!rs|�²¶�pr�t���xz�~� â rsx�ª��Q��t�x�{b��}!v����¾t�{���np�px�����¡ ¢ �~�}z��vyq�q���}zx¥�¦vy�z��{�|pr�t â r�x�ª&��vyq����!ru}ov��`vyr�|ËçM¨Â<1<`ç�è¤Þwå C�^ ±�vy|~�,xz�~�!{�xz�����o}z��v�qsq���}�x¥�¦vy�z�,{b|pr�t â rsx�ª��vyq��~��}&v����&|�{w² ç�¨�Þyè�0Qè¶Þwå>C9z­±�ÜM¨�Þ�Þwè1[¤Þwå9C{z ±�Þ=¨ çQåQÜpÞZÞwå9C m ±�vy|~�ËçM¨ åpÞ^<�8£Þwå9C m ¡% �!�9r��Jxz����}§vy�,�!�§Á�n����zrs����|�xo�p�px¤|~{w² t���v�|~�Q����xz�~�³n�{bqu�³r�|`x�{ �!�ýç�¡µ·jr��������B[�¡�Þ!}z�~{w²o}xz�~�,�¦v��z��{b|pr�t â r�x�ª,��v�qs�~��}É� ô �C²¤r�xz�(xz�~��rs������}�ru�9��v�q¥|�{��z��}��4Ù��¾�����¨{����7�¨}­����v�����}W�Cv�|�� vy��x�������p�pqsq���x�}U��xz�~�It��¦vy|~�Q��{��¤n�{bq��Q¡D¹`rs|~t��,��r�}B|~{bx�v���v��z��{b|prut â rsx�ª���vyq��~��± ¢ �~��{b����� [�¡ ç��p{9��}|~{bx�vynpnpqs�`¡ ¢ �~��{b�����Ú[9¡�Þ³}zxzrsqsq®vynpnpqsr���}�±`r ¡ ��¡�xz�~�!}z�p�~}�n¦v�t��or�}H#pqsx��������,�`��xz�~�&�§vyxzru{b|¦vyq¦���p|�t�xzr�{b|�*��?~���D�M�*��?O@)�§¡ ¢ �pru}£}z�~�Q����}zx�}Zxz�¦vyxcxz�~� â r�x�ª&��vyq�����}«|~�wv�����vy���onp�~}z�����Iv�²¶v���v�|���xz�¦vyx¥xz�~�t�{��z����}zn�{b|~�9rs|~�B����}zr��9�¦vyq�|~{��z��}¥�b��{w²ý��{����oxz�¦vy|Jxz�~�!{bxz�~���¥����}zr��9�¦vyqu}�¡ À {�x��oxz�¦vyx¥rsx¥ru}«x�{�x§vyq�qs�vbt�t���npx§vy�pq��&xz�¦v�x£����}�ru�9��v�qO|~{b�z��}¶�b��{w²ý²¶rsxz�Dv¾� vbt�x�{b��� ô ?Y@)�D�9� ô ?Y���§±�}zrs|~t��!xz�~� â r�x�ª!��vyqs�~��}vyqu}�{C�b��{w²�²¶rsxz�Ëxz�pr�}¤� v�t�x�{b��¡¯ | xz�~�¸rs��n�q�r�t�r�xzqs������}zx§vy�zx�����©jvy|~t�ª�{�}�����xz�~{9�¾� ¹9{�����|~}­��|2ÞyèQèQç�±£ß¥vyq��Q��xzxzr ± â ��r�t��~��q�v�|��

¹9{b����|~}���|�Þwè�è:<9�exz�~�&|~��² xz�zru�9r�vb��{b|¦vyq¦��vyxz�zrsÁ,twv�|,���&t�{���np��x����,�`��xz����v�t�xzr�{�|\� h Ë � � ²¶�~������»ru}�v¤�p|�r�x§vy�z����v�xz�zr Á�¡ ¢ �pr�}jr�}µvo}z�`������xz�zrutZ�¨{b�z�B�pq�v9¡ ¢ �~�Z|~��²Exz�zru��r vb��{b|¦vyq`�Ivyxz�zr Á�ru}e�~}z�¦vyq�qs�8

Page 10: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

t�{���np��x����»²¶rsxz���p�pq����$� t��¦vb}zrs|~�;�*/&{bq���� v�|��;K¥v�|�©O{`v�| ÞwèQè�Ü��§¡ m�q��Q{��zrsxz�p� [�¡�ÞJn�����}­��|`x�}�v|~{b|~}z�`������xz�zr�t!�¨{b�z�Ë¡Zm }z�`������xz�zrut!�¨{b�z� twvy|J���!{��px§vyrs|~���°vb}¥�¨{bq�q�{w²o}�¡ â ��²¤�zr�x�� �Ì[9¡ ç1�¤v�}

�4�Q?+���B�5Å � Ë � ?¥�JÅ � � � �k?A�4�Q?Y@)�J�ÈG �D© E ¹ � d�� � �Ì[9¡ Ü��²¶rsxz��� � ��g�³¯�4@�?���� Ë � ¡Z·�v�t�x�{b�zruª��5� � ���5�C±`xz�~��|��BCFET�X� CFE� �C¡Zmoqu}­{7�p�$#p|���( � ��Å � �vy|~� - � ����h Ë � �JCFE�±�xz�~��|- � � Þ@�?�� �Ìg�?�� h � CFE ��� Þ@�?�� �Ìg¬?;� h � C�E� ���I�²¶��rut­�¸ru}!v�}���������xz�zrut¾�¨{b�z�B�pq�v�¡�½¤|p�¨{b�zxz�p|¦vyx���qs��±�rsx!tyvy|p|~{bxo���Ct�{���np�px������`�¸vJæ â }zx���n;�¨{b��p�pq��Q� � t��¦vb}zrs|~�>�!{�| Ë � ¡\2�{�����{w�Q����±j²¥��}zxzrsq�qc�¦v�����|�{�x��Q��x�t�{b��n��px���� - � ¡¸¹�r���rsq vy�zq���x�{�xz�~�rs��n�q�r�t�r�xzqs������}zx§vy�zx�����©jv�|�t�ª�{Q}&����xz��{9�;� ¹9{b����|�}���|�ÞyèQèQç1�§±¦xz�pr�}�twv�|�����vbt��pr��������¸�`�¸��v�|prsnp���q�v�xzrs|~���Ì[9¡ Ü��§¡ % ���p{�|~{bx�t�{b|~}zr��p���¶xz�pr�}&v�npn���{�vbt��¸v�|`�������zxz�~����±�}zrs|~t��!xz�~�����!r�}¶|�{�x!v�����v�|`x§vb���t�{���n¦vy������x�{�moq���{b�zr�xz���Ú[�¡�ÞQ¡� � �¸)�����ê��\����%  �¢w��o��ë )Z�&�jëT¡B)¬ëµ� ¡Bê� ¢�#]��°��ê��% ����|�»�r�}c�����z��q�v����Q��±`xz���!}zx�{��§vb�Q��{b�Oxz�~�&�¦vb}zr�}Z����t�x�{���}£����t�{�����}cnp��{b�prs�prsxzr��Q��v�|~��r�x¥��v������|~��t���}�}§vy�z�³x�{o����}zx§vy�zx�xz�~�¥©ev�|~t�ª�{Q}O����xz��{9��¡®mo|7��qu���`v�|`xj����}zx§vy�zxzrs|~�!v�q���{b�zrsxz�p��ru}exz�~�crs��n�q�r�t�r�xzqs�����}zx§vy�zx����C©jvy|~t�ª�{Q}c�¨{b�jm��z|�{�q��9r¶������xz��{9�!� ¹9{b����|~}���|�Þwè�èQç�±yߥvyq��Q��xzxzr���xµvyq ¡`Þwè�è:<>�§¡ ¢ �pr�}O����xz��{9�t�{���np����}�}­��}¤xz�~��©jvy|~t�ª�{Q}¥��v�}zr�}¤r�|`x�{J{b|~��{b�Z}z��v�qsqu���&�9rs�,��|~}�ru{b|�����xz�p��{w²¤r�|��Jv�²¶v���v�n¦vy�zx�{b�xz�~�!}z���~}zn¦vbt���xz��v�x¤ru}£�p|pqsr�ºQ��qs��x�{C�¦v��Q��vC}zru�b|pr�#~twv�|`x¥t�{b|�xz�zrs�p�pxzr�{b|Ëx�{Cxz�~�!t�{b|��Q��������|~t��!{b�jxz�~�²¶vy|�x�������r��Q��|`��vyq�����}�¡�mo��x����!xz�~��t�{b��n�����}­}zru{b|e±�xz�~��©jvy|~t�ª�{�}&����xz�~{9��twvy|»vb�p�¸|~��²©jvy|~t�ª�{Q}����t�x�{b��}�x�{¸xz�~�I�¦v�}�ru}�¡TmG}zrs��r�q�vy�7����}zx§vy�zx�²¶v�}��p������q�{�n����E�¨{��Cxz�~�I�§v�xzr�{�|�v�q�´��z��q�{w�E����xz��{9�� â �p�~��Þwè�è�09±pá&��¹pvy�B��q vy|�Á���x³vyq ¡OÞwèQè98:�§¡% �¬tyvy|��~}­�£vy|!rs��n�q�r�t�r�xj����}zx§vy�zxj²¤r�xz�B�§Á�vbt�xj}��pr���x�}�� ¹9{�����|~}­��|,ÞyèQè�ç���r�|�xz�~�¥©jvy|~t�ª�{Q}�����xz��{9�{b�j���Q�prs��vyqu��|`xzq���±b²¥�ctwvy|�np�p���Q�¬�p|`²¶v�|`x���� â rsx�ª¬����t�x�{���}n� ©O���~{���t��!vy|~�7¹�{�����|~}���|�ÞwèQè�Ü9±�2�{b����vy|Þwè�èQÜ1�§¡�©���x Ë ���T� � �T����� ²¶rsxz� �I� ��9r�v����4± E �^¨^¨^¨^�D± � �³����xz�~�,��r��Q��|��p��t�{b��n�{Q}zrsxzru{b|T{b� Ë � ¡2��pqsxzr�n�q��`rs|~��� ç�¡·[��c{�|�xz�~�!�zr��b��x¤�`� �T� ±�²¥�³�¦v�����4�Q?Y@)�J� CFE �JÅ ���T� �ÊÅ �D© EÌË � �T� ¨% rsxz� T �D© E �¾Æ T � � �D© EÈÇ ��Å �D© E ! �T� Þ #vy|~� �

� ��! �T� Þ�# h Ë � �T� ��! ���¹ � d�h� �T� # �²¥���Q��x¥xz�~�!����q�v�xzr�{�| �4�Q?Y@)�B� C�E ��T � �;T �D© E � � ¨À {�x��¤xz�¦vyx

�� r�}¶v!��vyxz�zrsÁ,²¤r�xz�É»I³(Þ¤��{w²�}¥vy|~��»�t�{bqs�p��|~}Zxz�¦vyx¥t�{b|~}zr�}zx�}c{b�jvB»{Õ.»��9r�v��Q{b|¦v�q��v�xz�zr ÁCrs|�xz�~�I»J#p��}�xµ��{w²o}�²¤r�xz��xz��� â r�x�ªc��vyqs�~��}�± E �^¨~¨^¨^�,± � {b|7xz���¥��v�rs|��9r�v��Q{�|�v�qpvy|~�J¹ � d�h� �T�rs|�xz���I»c³¸Þ¶}�xµ��{w²�¡ ¢ �~�£q vb}zx¬��{w²�t�{�|`x§vyr�|�}�xz�~�£����}zru���¦v�q�r�|p�¨{b�z��v�xzr�{b|��¨{b��xz��� â r�x�ª¤n¦vyr���}�¡ ¢ �~�

0

Page 11: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

t�{�qs�p��|~}¤{��pT � vy���&xz�~���¿{��zxz�~{Q�Q{�|¦vyq�}���x¶{�� â rsx�ª!����t�x�{���}ov�|���xz�~�&q�v�}�x¶t�{�qs�p��|]� �D© E �¯G �D© Er�}¥xz�~����r�����t�xzru{b|�{��jxz���!����}zru���¦v�q�� ç9¡Â<9�§¡¹`�pnpn�{�}­�Jxz�¦vyx�xz�~��t�{�qs�p��|~}�{b�*T � vy|~� � � vy����{����p��������}�{Dxz�¦vyx�xz�~��� â r�x�ª���vyq��~��}�{b�rs|�x�������}zx�vy���&rs|�xz�~�&q��yvb�9rs|~�Ct�{bqs�p��|~}¥{b�

�� vy|~��xz�~��»�??�,qu��}�}crs|`x�������}�xzr�|~� â r�x�ªo��vyq�����}ov����!vyxxz�~����|~��¡ % ����|�²¥�!x§vyº��&xz�~��#p��}zx���t�{bq�����|�}¶{b� � � r�|`x�{ � � ±�²¥�!�¦v����

�4�Q?Y@)�B� C�E ��T � �vÆ T � � �D© E�Ç � � ²¶rsxz�� � � ÎÏÏÏÏÐ

± E ¡ ¡ ¡ ±�

¹ � d�h� � �Ñ ÒÒÒÒÔ ¨ �4<p¡�Þ��

¢ �§vy|~}z�¨{b�z��r�|���xz�~�|�2³�ÞBÕ�����vyxz�zrsÁ � � x�{�xz�zru�9r�vb��{b|¦vyqO�¨{��z�x�{ @Ë � �`��{��zxz�~{Q�Q{�|¦vyqOxz�§vy|~}z�¨{b�5���v�xzr�{b|~}�vy|~�Ëvynpnpqs�`r�|���xz�~��}��!xz�§vy|~}z�¨{b�z�Ivyxzr�{�|~}«x�{)T � ±���r���qu��}³v¾|~��²�����t����z����|�t��³����q�vyxzru{b|�4�Q?Y@)�B� C�E � @Å � ? @Å � © E @Ë � �ýå�¨ �4<p¡ ç��

% �¬twv�qsq�xz��ru}Oxz�z��|~twv�xzr�{�|�npq���}Oxz�§v�|�}z�¨{��z��vyxzru{b|!x�{¤xz�zru�9r�vb��{b|¦vyq��¨{b�z�Mn��p����rs|~�p¡ ¯ xjr�}jvyqu}�{ov�t­�pr����Q����`��»!?��Drs��n�q�r�t�r�x�����}zx§vy�zx�}!²¤r�xz�E��Á�vbt�xB}z�prs��x�} � ¹9{b����|~}���|SÞwè�èQç9±�2�{����`vy|4ÞyèQè�Ü��§¡�·���{�����������±»£?]�¸vb�p�9rsxzr�{�|¦vyqj©jvy|~t�ª�{Q}crsx����§v�xzr�{�|�}¥np��{9�9��t���v�©ev�|~t�ª�{Q}c�¦vb}zru}£{�����r�����|~}zr�{�|7»B³ Þ�¡

â rsx�ª�n¦vyr���}¾²¶rsxz�4v�����}�ru�9��v�qc|~{b�z� ¹ �9Û d�h� ÖWe Û }z��v�qsq����¾xz�¦vy| v��br��Q��|0x�{bqu���§vy|~t��Jvy����qu{9t­º�����¡¢ �����Jvy���&t�{�|~}�ru�p�������Ëvb}c�§Á�v�t�x¶��r�����|pn¦vyrs��}�¡cbZ����{w�`rs|~�Bxz���!t�{b�z����}zn�{�|~��r�|~�¾����}zr��9�¦vyq¦x����z��}¥x�{xz�~�!�zr��b��x5� �¦vy|~�Ë}�ru�p��±��4<p¡�Þ��c����t�{b����}

�4�Q?Z@P�J� CFE �VT � ?¯Æ T � � �D© E�Ç Î ÏÐ �*� åå � � C �å ¹ � d�h� ��� C �

ÑÓÒÔ �¾Æ � �D© E ¹ � d h� ��� å Ç

²¶rsxz� �*� � �9r�v����4± E �^¨^¨^¨~�,± � �Zvy|~� � � �¾Æ ����� � C � Ç ¡ ¢ �§v�|�}z�¨{��z��vyxzru{b|Cx�{�xz�zr��9r vb�Q{�|¦vyqp�¨{b�z�Í�br��Q��}�~} �4�Q?Y@)�B� C�E �"� T � @Å � C �5� ?"� T � @Å � C �È© E � ! �*� å

å @Ë � C � # �¾Æs� � å Ç �4<p¡Í[��²¶rsxz�+@Ë � C � v��§?'�M³,ÞPÕ��§?'�®xz�zru��r vb��{b|¦vyq���vyxz�zrsÁ7v�|~�/� � �¯� �D© E ¹ � d�h� ��� ¡ ¢ ���¥©jv�|�t�ª�{Q}�����xz��{9�twv�|�����t�{�|`xzrs|������I����{b�?xz�pr�}�n�{br�|`x�¡£©O{9t�º`rs|~� â rsx�ª�n¦vyrs��}¥���wv�|~}¤xz�¦vyx¶xz�~����v�����t�{b|~}zr��p�������°vb}�§Á�v�t�x�¡ % �&xz�~��|¸v�}�}��p���&xz�~���zr����`x5� �¦v�|��Ë}zr��p�!{b�L�4<p¡Í[��cª�����{�¡£©O{9t­º�rs|~�Crs|�xz��{9���~t���}&vy|Ë���z��{��ors|xz�~�!����t��p�z����|~t��!����q�v�xzr�{�|�{b�jxz�~��{b���p���£j$� �D© E ¹ � d�h� ��� j pV�O� Û ¹ �>Û j^d�h� ��� j mw¡¢ �~�£|~��²4xz�zru�9r�vb��{b|¦vyqp��v�xz�zr Á @Ë � C � tyvy|����¤t�{�|~}�ru�p��������v�}¬������|�t�{b��n��px��������Cxz�~�&©jvy|~t�ª�{Q}����xz�~{9�°vynpnpqsr����,x�{ �Ìg�?�T � T h� �B�$�4�Ê?�@P�J� CFE � �ÌgB?�T � T h� �J�¬¨ �4<p¡Ó<>�À {�x��¤xz�¦vyx¥å!r�}£vy|���ru�Q��|pn¦vyr��Z²¤r�xz���¾�pq�xzrsnpqsrut�rsx¼���!v�|~��xz����t�{bq��p��|~}c{b��T � vy���o��r��Q��|`�Q��t�x�{b��}�¡ ¯ ��0�����å��,iM�cr�}³vy|���r��Q��|�n¦v�rs��{b�L�4<p¡Â<9�§±�xz�~��|]T h� ��i���å9±�}�{É�Ìg�?&T � T h� �J�$�4�Ê?¥@)�B� C�E ��i!���:ivy|~�

�4�Q?¥@)�B� C�E �¬iA?+�=i � T � T h� �{�4�Ê?Y@)�B� C�E ��i� T � � �|� T h� ��iB³�� h� �¬iM�¬¨m&}�v¾����}z�pqsx�± j:�4�Q?Y@)�B� C�E ��i£?Y�=i+j$p��µjlT � � h� �¬i�j$pZ�²j�� � j pV�Z¨

è

Page 12: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

¹`rs|~t���xz�~�,��r��Q��|�n¦v�rs��}C{��J�4�¾?�@)�B�DC�ED� xz�¦vyx�}zxzr�qsqc�¦v��Q��x�{�����t�{���np��x����4vy����t�{���np��x����4vb}��ru�Q��|pn¦vyr���}£{��H�4<p¡Â<9�§±�xz�~�!����}zr��9�¦vyq¦|�{��z� {b�Oxz�~�&q�{9t�ºQ��� â rsx�ª&n¦vyr���}£q�rs��r�x�}£xz�~�Bvbt��pr�����v���qu��vbt�t�����§vbt���¡% �~��|É�4±��,�-�®r�}£voq�{9t�ºQ��� â rsx�ª¥n¦vyr��®�¨{����4�O?.@)�B� C�E ��±�²¥�¶t�{b|~}zr��p���L� ô �,�-�®²¶rsxz� ô �¯@I³.± C�Ev,qu{9t­º������¦vy�z��{�|�rut â r�x�ª�n¦vyrs���¨{������¥�N�F�¬�e¡ % �~��|�xz�~��n�{�q���t��¦vy|~�Q��}³x�{~��±e²¥��t�{b|~}zr��p����5� ô ?���� CFE �,�-�cv³qu{9t­º���� â r�x�ª¤n¦vyr��c�¨{b���4��?����J� CFE ��¡ % �~��|�xz�~�on�{bqu�ot���v�|~�Q��}�±Qxz�������or�}¥v!�zr�}zºxz�¦vyx!xz�~�C����}zru���¦v�q¬|~{b�z������{w²o}�¡�©O{9t­º`r�|~� â rsx�ªCn¦vyr���}!²¤r�xz�Ev�����q�v�xzrs����qs�¸q vy�����¾����}zr��9�¦vyqµ|~{b�z�np��������|`x�}¤xz�~�Bvbt�t��p�§vyx��!t�{���np��x§v�xzr�{�|�{��®{�xz�~���o��r��Q��|�n¦v�rs��}�¡ ¯ |�{����p���¥x�{����wvb}z�p���&xz�~�����z��{b�¶{b|xz�~������t��p�z����|~t��¾����q vyxzru{b|¸v���x����oqu{9t­º�rs|~��{b��xz�~� â r�x�ª�����t�x�{���}¤²¥�Btwvy|�vb}¤²¥��qsqj�~}­�!xz�~� â v��`q���r��b��Q�~{bxzru��|`x��!�v�VhD� �4�Ã?$���J�DCFE,����rs|~}zx��wvb��{��oÞ:�9� ô ?;���§±�²¤�pr�t����br��Q��}&xz���7}z��vyq�q���}zxo����}�ru�9��v�q|~{b�z�Í�¨{b�¥v!�br��Q��| â rsx�ª¶�Q��t�x�{b��¡ % �¶�p{!|�{�xc�p{�xz�pru}Zrs|�np�§vbt�xzrut���±`}zrs|~t��¤xz�pru}Z�����Q�prs����}¶v����9r�xzr�{b|¦v�qt�{���np��x§v�xzr�{�|�v�q�²¥{b�zº�¡

¢ �~��{b����� <p¡ ç,��rs����}Bvy|��pn�n�������{���|~��{b|�xz�~�C����}zru���¦v�q¬|~{b�z��v���x����Bt­�¦v�|����C{��cn�{bqu��¡ ¢ �~�xz�~��{b����� �~}���}¥xz�~�³�¨{�qsqu{�²¶rs|~�Cqu������v�¡ócÃ9òTòDØ öP¸¨ì j:�4�Q?����J� CFE �4��?�@P�J�^j p�� ��vyÁ�������¡ £¢ p¥¤£¦¦¦¦¦

�2?Y@�2?�� ¦¦¦¦¦ ¨÷ Ç�Æ�Æ~ñ ©O��x�§ ���!}z�~t­�Ëxz��v�x��µ�X§kh¨§<¡Z·er���}zx¤|~{bx��!xz�¦vyxj:�4�Q?����J� CFE �!j$p � ��vyÁ|"}~|��M� E j:�4��?+���B� C�E �¬�oj$p� �IvwÁ|ª©�}~|¬« � E j=�0§ �4�Q?����J� CFE § h �$�0§Ã�-�^j$m� ��vyÁ�­�����¡ £¢ p¥¤ ¦¦¦¦ Þ�2?�� ¦¦¦¦ ¨¹9��t�{�|���±�}zr�|�t��{�4��?Y���B�DC�E^�4�Q?Y@)�B�H��g(³¯�0�{?Y@)�$�4�Q?����B�DCFE,��±�²¥�!��v����j:�4��?+���B� C�E �4�Q?Y@)�B�^j$p � ��vwÁ�­�����s £¢ p£¤ ¦¦¦¦ Þ�³ �£?�@�A?�� ¦¦¦¦ �²¶��rut­�Ënp��{�����}¶xz���!qu������v�¡ M

iBïeÃ`Æ~Ç�Ã9ò ö)¸�îKB â�ã �5� ô ?�@P�DC�EW�,�-� ëDâAä£í�æ4ãõî�ð�ä=æ¶à�è$ß:à �4� ?¥@)�B�DC�ED� ìTæ4ã ù à,â^êUæ4á=é�ä�ó+Þ-ß:àWûÙ2�²j=�4�Q?Y@)�B� C�E �¬��?¯� ô ?�@P� C�E �oj$pQ¨DTè$ãÌâ^à�ã ù âAç ù ä:Þ�ü�âAß,è�ð�ß�óúâ ñ ã ù âBÞVâ^ì à5â^êWæ4á:é�ä1ó+ÞVß=àUûNê ä�ã4æ¶ê¶ö�â^ê® �²j:�4��?Y���B� C�E �¬�.?¯� ô ?���� CFE �lj p�� ¦¦¦¦¦

ô ?¥@ô ?Y� ¦¦¦¦¦ : ��vyÁ�������¡ £¢ p¥¤ ¦¦¦¦¦�2?Y@�2?�� ¦¦¦¦¦ = Ùx¨÷ Ç�Æ�Æ~ñ bZ�I�B�pqsxzr�n�q��`rs|~�7�4�Q?¥@)�B�DC�ED�¬��?Z±=�]� Ø ²¶rsxz�7±£�²� ô ?¥@)�DCFE£�`��4�Ê?Y@)�J�§±�²¥��¦v��Q� ���&� ?�±��&³;���P�zÞ�³Y@+±��n�¾�4�Q?Y@)�B� Ø ¨bZ������{b����vy|pr�}zrs|~�Bxz����������v�xzr�{�|�rs|�x�{

�4�Q?����J�È�&� ?�±��§³¥�¬�&�zÞH³¯�4@]?����ȱ1�H�¾�4� ?¥@)�B� ØÞyå

Page 13: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

vy|~���`�J�¾�pqsxzr�npqs�`r�|���²¶rsxz�Z�4�Q?����B�DCFE�±�²¥�!�¦v��Q��4�¯?+���B� C�E ��� ? ±

Þn³¯�4@�?����ȱ � � ÞÞ�³��4@�?����ȱ �4� ?����J� CFE �4�Q?Y@)�J� Ø �

� ô ?¥@ô ?Y� �4�Q?+���B� C�E �4�Q?Y@)�B� Ø ¨¢ ���!np��{9{b�O�¨{�qsqu{�²�}¥����{b�

j:�4�Q?����J� CFE �4�Q?Y@)�J� Ø j p¥�²j=�4�Q?����J� CFE �4�¯?Z@)�B�^j$pLj Ø j p;¨ M©O{9t�º`rs|~� â rsx�ª�n¦v�rs��}B�Q��|~���§v�x���}Cv�|E���z��{b��rs|0xz��������t��p�z����|~t�������q�v�xzr�{b|0{b�¥xz�~��{b���p���B{b�¥xz�~�

����}zr��9�¦vyq�|~{b�z��}�¡ % �~��|�xz���Zn�{bqu�¬t���v�|~�Q��}�±�xz�pr�}O|~{b�z� t­�¦vy|~����}jvb}�²¥��q�q9v�|���}�{o�p{9��}�xz�~�Z���z��{b�j{b|xz�~�¤|~��²S����t��p�z����|~t��¶����q�v�xzr�{b|e¡¬m&}Zqu{b|~�Bvb}Zxz�~�¤����}zru���¦v�q�|~{b�z�Ír�}Z�B�~t­�J}z��vyq�q����«xz��v�|�xz�~�o|~{b�z�{b�jxz�~�!xz�zr��9r�v��Q{�|�v�q���v�xz�zr Á�±�xz�~�����!r�}¶|~{C���wvyq�n���{��pq����Ë¡ % �~��|Jxz�~�!n�{bqu�!t­�¦v�|�����}�±�xz�~�!|�{��z� {b�xz�~�¤xz�zru��r vb��{b|¦vyqp��v�xz�zr Á7ru}Z��{b�~�b�pq��C}zn��yvyº`r�|��!�B�pqsxzrsnpq�r������`�Oj:�4�?]���B� C�E �4�?É@)�B�^j$p«±9vy|~��}�{rsx�r�}«vbt�t���npx§vy�pq��Zxz�¦vyx�xz�~�c����}zr��9�¦vyq�}��b��{w²0²¤r�xz�Cxz�pr�}�� vbt�x�{���¡¬�o{w²¥���Q����±�xz�~�c� vbt�x�{b� Û ô ?2@ Û � Û ô ?A� Û}z�~{b�pq���|~{bx¶�b��{w²Sx�{9{�q vy������¡ ¢ �pr�}¥r�}c|~{bxcxz�~�!twv�}��o�¨{b�cxz�~�&q�{9t�º���� â rsx�ª&����t�x�{���}c� vy��v�²¶v�������{b�xz�~��|~��²2n�{bqu��¡�m&}Bv�t�{b|~t�qs�~}zr�{�|e±¦²¥������xzx�����n��px�xz�~��|~��²2n�{bqu��v�²¶v��¸����{���xz�~��q�{9t�ºQ��� â rsx�ªn¦vyr���}�¡ ¯ |�n��§v�t�xzr�t��£npr�t�º`r�|���xz�~�£n�{bqu�£|~�wv��¬q�{9t�º���� â r�x�ª£����t�x�{���}Zr�}¬|~{bxµ�Q���z��rs|`x�������}zxzrs|~�Cv�|`�`²¶v��}zrs|~t��&xz�~��}��!�¦v����Bvyq����wvb�9��������|�t�{���np�px�����¡x��OØ~ò�ð¬×�à öP¸ ì â ��tyvyqsq�xz�~�¾np��{b�pqu��� ����{���à�Á�v���npq��{[9¡�ÞQ¡ ¢ �~�¾����}z�pqsx�}&{��Zxz�pr�}!�§Á~vy��npq��B²¥������Q��|����§vyx����4�~}zrs|~��2¸v�xzq�vy�e¡ % �In����z�¨{b�z�,����»v�*ÝQå�}zx���n~}�{b�!xz�~�¸}zn���t�xz�§vyq&xz�§v�|~}��¨{��z��vyxzru{b|©jvy|~t�ª�{�}c����xz�~{9�I²¶rsxz�7@x�v?BÞ¾vy|~��v�|�rs|pr�xzr�vyqO����t�x�{��¤²¶rsxz�D���Q�¦vyqj��qu������|�x�}�¡ ¢ �~� â r�x�ª���v�qs�~�²¶rsxz�J}z��vyq�q���}zxZ����}zru���¦v�q�|~{b�z�Ír�}�±£�ýå�¨ èQè98bÜ&{�� ô �¯@J³]±9CFE��ýçM¨Ó<�<oÞyå>C�^§¡ ¢ �~�ot�{b�z����}zn�{�|~��r�|~�����}zr��9�¦vyq¦|~{b�z�r�}¶Ý�¨·[¤Þyå>C°¯�¡§3µ������rs|~�¾xz�~�!t�{��z����}zn�{�|��9r�|�� â r�x�ª&n¦vyrs�¥r�|,xz�~�&����t����z����|�t��&����q�v�xzr�{b|t����wv�x���}�vy|Ë���z��{b��{b�jxz�~�¾{����p���¶{b�µÝ¤Þwå C±¯ ¡�bc���¨{b���&np�p���br�|~�p±�²¥�&�¦v��Q�

j­².j~���²j=�4�Q?Z@P�J� CFE �BÅ � ?ZÅ �D© EÈË � j$mL� å!¨¢ �����¨{��p�¬#p��}zxB��{�²�}B{b� ¢ vy�pq�� <�¡�Þ�}z�~{w²¿xz�~�����z��{��¾²¶�~��|�xz�~��n�{�q���t��¦vy|~�Q��}Bv���x�����q�{9t�º`rs|~��±�¡¢ ���(#p��}zx¥n�{bqu�5�!�¯<Cq�r���}¥x�{�xz�~�³�zru�b�`x�{b�§�)�4�B�U�B�§¡ ¢ ���!}���t�{b|~�,n�{bqu�&r�}¥xz���!���yvy|Ë{b�jxz�~�³x¼²¥{}z��vyq�q���}zx â r�x�ª!��vyqs�~��}�¡

¯ |I��{bxz�Ëtwvb}���}�±`xz���!����q vyxzr��Q�!���z��{��¤������vyr�|�}¶{b�Oxz�~�!{��������o{��cÞwå>C9³�¡ ¯ |�xz�~�!}���t�{�|~�,tyvb}���±�xz�~�vy�~}�{bq��px������z��{b�Br�}��B��t��Eq vy��������²¤�pr�t��T���wvy|~}�xz�¦vyxBxz����}z��v�qsqc��r��Q��|`��vyq�����}¾{��(�4�k?X���B�DC�ED�twv�|�|~{�x¾���,�¨{���|~�E�����z�»vbt�t����§v�x���qs�`¡ À {bx����~{w²¥���Q����xz�¦vyxCxz�~���pn�n����C��{b�p|~�({b�¶xz�~�,����}�ru�9��v�q|~{b�z� rs| ¢ �~��{b����� <�¡ çCr�}¶�§vyxz�~����n���}�}zrs��ru}zxzr�t�¡¤·~{b��xz��ru}onp��{b�pq����xz�~����npn����o��{b�p|~���¨{bq�q�{w²¤r�|~�¢ ����{b�����k<�¡ çZr�}cÞ�¨ Ü�0�±§²¶�pr�t���r�}O�¾�~t�����ru�Q������xz�¦vy|�ÞZÞyå C9z ¡ ¢ �~�¬���yvb}�{b|!ru}�xz��v�xHj:�4�B?V���J� CFE �4�A?@)�B� Ø j$p�r�}µ�¾�~t��I}z��v�qsqu���¬xz�¦vy|xj:�4�x?����B� CFE �4�x?�@)�J�^j pHj Ø j$p«¡ ¢ �~�¶����}zr��9�¦vyq Ø r�}Zn�{9{b��rs|�xz�~��p{b��r�|¦vy|`x¶��ru�Q��|`����t�x�{���}¥{b���4�®?����J� CFE �4�\?×@)�B�«²¶�pr�t���v����o²¥��q�q�����np����}���|`x������`��xz�~��´&�z�`q�{w�}zn¦vbt���¡®m�q�}�{!|�{�x��¤xz�¦vyx�±Qrs|�xz�pru}ctwvb}���±bxz�~� â v��`qu��r������Q�~{bxzru��|`x���t�{b�z����}�n�{b|~�p}¬²¥��q�q�x�{�Þ��M� ô ?]���§¡¢ ���!���wv�}�{b|�r�}¶xz��v�x¤xz�~�!|~��²ý����}�ru�9��v�qOr�}�vyq����wvb�9��� vyr��zqs��{b�zxz�~{��Q{b|¦v�q�x�{Cxz�~� â rsx�ª!����t�x�{���¡

·jrs|¦v�qsqs�`±�²¥��rsq�qs�~}zxz�§vyx���xz�~�I�9r�R�������|~t�������x¼²¥����|0q�{9t�º`rs|~�������¨{b���Jvy|~�4v���x����¾xz�~�Jt­�¦v�|�����{b�n�{�q���¡ ¢ �~�o���z��{���}c{b|�xz�~�¤����t��p�z����|~t��o����q vyxzr�{�|Jvy���¶�br��Q��|���� ¢ v��pq��L<p¡�ÞQ¡ â ��tyvyqsqpxz�¦vyxcxz�~�¶q�{9t�ºQ���

Þ�Þ

Page 14: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

¢ v���qu�(<p¡�Þ=¼¬à¬�z��{b�¶rs|Ëxz�~�³����t��p�z����|~t�������q�vyxzru{b|°vy��x����¶q�{9t�º`r�|���{b�µv â r�x�ª!n�v�rs��¡� <�¨ å Þ=¨ ç�ç¤Þwå9C{z©O{9t­º�rs|~�C�����¨{b���!t���v�|~�Q�!{b�jn�{bq���¼j�²!j^� ݤÞyå>C } ÞZÞyå>C9zj�²!j^�)��j - � j m ÞZÞyå>C9³ ÞZÞyå>C9³©O{9t­º`r�|~��vy��x����ot��¦vy|~�Q�!{��en�{bqu�=¼j�²!j^� ç¤Þyå>C°¯ ݤÞyå>CFEj�²!j^�)��j - � j m ݤÞyå>C9´ ݤÞyå>C�^

�¦vy�z��{�|pr�t â r�x�ª���vyq�����ru} ô �çM¨Ó<�<oÞyå>C�^�¡ % ����|�xz�~� â r�x�ª���vyq�����}Cv���������t�{b��n��px����»vy��x�����xz�~�t���v�|~�Q�¶{b��n�{bqu�ov�|��7xz����|�xz����}z��vyq�q���}zx¬�¦v��z��{b|pr�t â r�x�ª¤��v�qs�~�¤ru}¬q�{9t�ºQ����±`xz�~�¤���z��{���}¥vy���¶}z��vyq�q����²¶����|4xz���Jn�{�q���r�}C� v���v�²¶v����0�Ê�_<9�§±£v�|��0q�vy��������²¶�~��|(xz�~�In�{bq���r�}�|~�wv��Cxz�~�Iqu{9t­º���� â rsx�ª��vyq��~�2�0�!� Þ=¨ çQç¤Þyå C9z �§¡ ¢ �~�o���wv�}�{b|�r�}cxz�¦vyxcxz�~��´&�z�`q�{w��}z�p�~}�n¦v�t��¤ru}�#pq�x���������²¤r�xz�Ixz�~�o��v�xz�zr Á�4�x?�@)�J� CFE �4�7?����B�¬²¤�pr�t���rs��n���{w����} â r�x�ª£��v�qs�~��}Z|~�wv��c@�v�|~�C� vy�£v�²¶v��7����{b�µ�C���pxZ������{w�Q��}xz�~{Q}��!t�qu{Q}��&x�{��¦¡� ���¸� ê��¸% �&)c��L¡A¢ �jëT¡�'T��¢w��ë�¶ #2¤ ë§�¸� )c¡�ë«��#�ê7¡A¢ �¥¡Bê�¦o#��

% �oë§�¸#� m?������vyr�|prs|~���Q�~��}zxzr�{b|4ru}Cxz�~�,|����,���zr�twv�qo}zx§v���r�qsr�x �({��oxz�~�I�§v�xzr�{�|�v�q¶©jvy|~t�ª�{�}¾����xz�~{9��¡ ¢ �~�|`�p�����zrutwvyqµ�����¦v��`ru{b�p�B{b�Zxz�~��©jvy|~t�ª�{Q}&����xz�~{9�¸r�}!²¥��q�q¬�p|~�p����}zx�{9{9��¡ ¢ �`nprutwvyq�qs��±exz�~�C����xz��{9�������v�����}j|`�p�����zr�tyvyqsq���²¥��qsq`²¶�~��|Cr�xjr�}�twvy�������pq�qs��rs��npqu������|`x����!¼�xz�pr�}�rs|~t�q����p��}jn¦vy�zxzr vyq`����{��zxz�~{Q�:�{b|¦v�qsr�ªyvyxzr�{�|×�*/o�zr�����}&��x�vyq ¡jÞyèQè=<>�c{b�¶���pqsqj����{b�zxz�~{��Q{b|¦v�qsr�ªyvyxzr�{�|×� á!v�|pr���qj��x�vyq ¡eÞwè>8�Ü��§¡cmo|�rs�{�npqsrut�rsx&����}�x§v��zxo��v��Jq�{�}­�B}zx§vy�prsq�rsx¼�I²¶�~��|¸xz�~��r���npqsr�t�rsx!}z�prs��x�r�}&t�qu{Q}��!x�{Jv â rsx�ª���vyq����Q¡I3¬�p����rs|~�rs|~}zx��wv��,{��Ors��npq�r�t�r�x£����}�x§v��zxzrs|~�¾ru}£np�����¨���z�����,r�|�xz��v�x¶}�r�xz�¦vyxzr�{�|� ©�����{��~t���vy|~�J¹9{b����|~}���|TÞwè�èQÜ��§¡moqu}­{�}­���2�02������z�������Q��|°vy|~��¹`n���|~t���Þwè�è>8:�¬�¨{��&v�|~{bxz�~���¥|`�p�����zrutwvyqj}zx§vy�pr�qsrsx¼�I����}��pq�x�¡ ¢ �~��{�|pqs��9r�R�������|~t��&����x ²¥����|Jxz�~�B©ev�|~t�ª�{Q}Z�,��xz�~{9��v�|���xz�~�&�§vyxzru{b|¦vyq®´��z��q�{w������xz�~{9��r�}cxz�~�!t­�¦vy|~���!{b�n�{�q���±p}­{7rsx¶r�}¶}���¿�t�r���|`x�x�{�}zxz�~����xz�~��}zx§vy�prsq�rsx �Ë{b�jxz�pr�}¶}zx���ne¡% �¥t�{b|~}zr��p���µxz�p�����o}�{b�p��t���}¬{b�����z��{b��}µrs|�xz�~�¤�§v�xzr�{�|�v�q¦©ev�|~t�ª�{Q}j����xz�~{9�]¼exz�~�¶����vyq���v�xzr�{�|�{b�xz�~�C©jv�|�t�ª�{Q}¥����t��p�z����|~t�������q�v�xzr�{�|�'�xz�~�¾t��¦vy|~�Q��{b��n�{�q��{'�vy|~�Ixz�~�B}�{bqs�pxzru{b|�{��®q�rs|~�wv��o}z�9}zx�����}²¶rsxz����?�@P� {b�L�Q?�����¡iBïeÃ7ócØ�ÛjÕQô�Æ�ɳÇyÃ�Õ`ÒjÇyÇ�Ã9ÛjÕQÃ�Ç�Ã9×�Ø9Ö�Ù�Æ¦Û ¢ �~��#p��}zxµ}­{��p��t��c{b�~���z��{���}eru}®rs|Bxz�~�¥����vyq���v�xzr�{�|�{b�pxz�~�©jvy|~t�ª�{�}¤����t����z����|�t��¾����q vyxzr�{�|e¡�¹`�pn�n�{Q}���xz��v�x&xz�~�C}zn���t�xz�§vyqjxz�§vy|~}z�¨{b�z�Ivyxzr�{�|�©jvy|~t�ª�{�}¤����xz��{9�np��{9�9��t���}ov�|�v�t�t��p�§vyx��!����t����z����|�t��³����q�vyxzru{b|e±�r ¡ ��¡

�4��?Z@)�B� CFE �BÅ � ?¥Å �D© E�Ë � ��²²¶������� ² �¦vb}�v�}z��v�qsqµ|�{��z��t�{b��n�v������¸x�{Ixz�~��{bxz�~����x����z�,}�r�|�xz�~�����Q�¦v�xzr�{b|e¡ % ��}§v��¸xz�¦vyxxz�~�!����xz�~{9�Iru}¤�¦vbt�º`²¶vy���Ë}zx§vy�pq���¡ % ��v�}�}��p���&xz�¦v�x¤²¶�~��|/o�§vy�{�Ϲ9t�����r��9x&{��zxz��{��Q{�|�v�qsruªwvyxzru{b|��²¶rsxz������{b�zxz�~{��Q{b|¦v�qsr�ªyvyxzr�{�|-��ru}Z��}������¨{��Z��v�º`rs|~� �4�?É@)�B� C�E �¬G eo{��zxz��{��Q{�|�v�q9�Q����}��~}�G e C�E v�|��

Þyç

Page 15: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

G eors|�m�q��Q{��zrsxz�p�Gç9¡�Þ�±�²¥�&�¦v��Q�j­².j pV� · j Ë � j$m§Ò � Ý9¡�Þ��� j:�4�Q?¥@)�B� C�E �!j$p�Ò� ��vyÁ�������¡ £¢ p¥¤ Û �{?Z@ Û C�E Ò®¨

¹9�����Ìb���¸{���t­ºJÞwè�èQÜ�±MX�ç9¡Â<�¡ Ý1�ev�|~��� �¤ru�b�¦vy� Þwè�èQÜ9± ¢ �~��{������Þ`09¡�Þwç����¨{��¬��{9�9r�#~���./o�§vy�2��¹9t­�p��ru��x²¶����|É����gp¡ ¢ �~�&���z��{b�|² twvy|�����t�{b�,�oq�v����Q�¤²¶�~��|]@�r�}¥t�qu{Q}��ox�{�v�|J��r�����|`��vyq��~��¡cã¶���`r�{���}zq���±xz�~�&����q�v�xzrs���&���z��{���±�xz�pr�}¥r�}(j�²!j$p-�o��j Ë � j$m¥r�}¶vyq�²¶v���}c}��Ivyqsq ±�r ¡ ��¡¬xz�~�oq�v����Q�&��ru�Q��|���vyqs�~��}¥{��c�4�\?@)�B�DC�ED� twvy|,vyq�²¶v���}¬���¶t�{b��n��px����7����q�v�xzrs����qs�,vbt�t��p�§v�x���qs��¡ % �~��|�n¦vy�zxzr vyqp����{��zxz��{��Q{�|�v�qsruªwvyxzru{b|r�}B�~}­���Ex�{¸��vyr�|`x§vyr�|({b�zxz�~{��Q{b|¦v�qsrsx¼�E{��¤xz�~�Ë©jvy|~t�ª�{�}�����t�x�{���}�±cvy| vb�p�9rsxzr�{�|¦vyq¥���z��{��C�B��}zx7���vb�p�p����x�{�xz�~�¾�zru�b�`x5� �¦vy|~��}zr��p�Bt�{b�z����}�n�{b|~�9rs|~��x�{�x����z��}&²¶�pr�t��¸��v�����������|¸t�{b��np�px����¸�9�p�zrs|~�xz�~�!����{b�zxz�~{Q��{b|¦vyq�r�ªwv�xzr�{�|O¡ ¢ �~��}���vy���!x ��n�rutwvyq�qs�J{b�jxz�~��{b���p����{b� r Ò�j Ë � j$m�¡lCïjØ~Ûeø~Ã�Æ�ñ�ðZƦ×�à ¹`��npn�{Q}���xz�¦v�x!moq���{b�zr�xz��� [9¡�Þ�ru}!�§Á���t���x�����rs|T�§Á�v�t�x�vy�zr�xz�p����xzr�tQ¡ % �B�p{vb}�}z�p����±��~{w²¥���Q����±�v�|¸���z��{b��{b|�� v�|��~�C¡o·~{bq�q�{w²¤r�|~�O� �¶r��b�¦v�� ÞyèQèQÜ�±�ßZ��v�npx�����Þ�01�§±�²¥�¾��v��vb}�}z�p���&xz�¦vyx¶xz�~�³�¦v�t­º`²¶v����I���z��{b�£j�¹ � j m¤rs|Ëxz�~�Cæ â � v�t�x�{b�zruªwvyxzru{b|�5� �X� � ³&¹ �r�}¥��{��p|��p���J����{b� v���{w�Q�!�`��j­� � j$m ® ÒD²¶������� ® ����n���|~�p}¤{�|7»¦¡¢ �~��|;�¬h¨� � �BCFE�� g � ?���hº¹ � �BC�E�±¬}�{�²¶rsxz�;( �D© E �NÅ �D© E �vy|~� - � �F��h Ë � �BC�E�±®²¥��¦v��Q�

�4�¯?+���B� C�E ��( �D© E � h � � � CFE ?�( �D© E¨- � � �4�Q?����B� CFE �4� ?¥@)�B�º²/� CFE�4�Q?����J� CFE ��( � ?�( �D© E¨- � � �4�Q?����B� CFE �4� ?¥@)�B�º²/� CFE � Ý9¡ ç��³ �4�Q?����B� CFE ��( �D© E � h ¹ � � CFE ¨¢ ���!|~{b�z�?{��®xz�~�(#p��}zx¤x����z�?r�|�xz�~�!�zr��b��x5� �¦vy|~�Ë}�ru�p�&r�}¥��{b�~�b�pqs�J}zn��yvyº`r�|��7��{b�p|~�p�����`�

j=�4�Q?����B� CFE �4�Q?Y@)�B�^j$pLj�� CFE j$m=j Ë � j$m§Ò®� � Ý9¡Í[��²¶�������.j­�BC�E$j mL� Þ��=�9�+�����0� � �§¡ ¢ ���!|~{b�z� {��exz�~��}���t�{�|���x����z�?r�}¶��{b�~�b�pq���}�n��wvyº�rs|~�C��{��p|��p����`� j:�4�Q?����J� CFE �.j pT��m��0� � �zÒ � Ý9¡Ó<>�²¶����������m`�0� � ��� �M�)���1�0� � �D�=�9�+���>�0� � �Cru}Cxz�~��t�{�|~��r�xzr�{�|4|`�p�¾������{b� � � ¡ ¢ �~�I� vbt�x�{b��}É�4�²?���J� CFE �4�Ê?�@P�J�&v�|~�×�4�Q?����J� CFE �¿�~}z�¦vyqsq�����v������,{9����}�x��¿|~{b�z��}¶�p|�qu��}�}*�Iru}ot�qu{Q}��&x�{�vy|��ru�Q��|`��v�qs�~��¡ ¢ �~��}zrs|~�b�pq vy�¥��vyq�����}¤{���� � npq v��Ëv¾��rsx§vyqO��{�q���¡ócÃ9òTòDØ ü-¸¨ìC» â^ÞVß1ã�â ã ù â â^æSü�â^Þ=ð�ä:æ¶àWê2ß,è Ë � ë$ò �4±We��,ÖWe^� ìTæ4ã ùxj$ÖWe�j mB� Þ è$ß=à f�� Þ��^¨�¨S¨��U» ä:Þ-á±We��v� ô eH?Y@)�DCFE�÷ » âÌö§Þ-âBã ù âJà5â^êWæ4á:é�ä1ó+ÞVß=àUû Ù�eI�²j:�4�Q?Y@)�B�DCFE,�BÅ � ÖWeH?Y±We^Å � ÖWe�j$pH÷2ø�ù â^Þ

�M�+� �9�0� � �s� ��rs|e ¼½ : ô e�?��ô en?�@ = m ³¯�4@]?���� m Ù meJ¾¿ E�� m ��M�)���1�0� � �ÁÀ ��vwÁe ¼½ : ô e�?��ô eH?Y@ = m ³¯�4@�?+��� m Ù me�¾¿ E�� m ¨

Þ`[

Page 16: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

÷ Ç�Æ�Æ~ñ ·~{bqsqu{w²¤r�|��Ëxz����rs���p�$#p|prsxzru{b|~}�±�²¥���¦v��Q�]�M�+� �>�0� � � � j­� � ÖVj m�� �M�)���=�0� � ���¨{���v�|`�Ö.uOwÉ�&²¶rsxz�\j$Ö-j$mI� Þ�¡¬·p��{b�?xz�~�Bvb}�}z����n�xzru{b|~}¶{b�jxz�~�³qu������v�±�²¥���p���zrs���

� m�+��� �0� � �H�²j­� � ÖWe:j mm �  ! �Ìg(³¯�4@�?���� Ë � �ÈÖWe�4@]?����̹ � d�h� ÖWe # ÂÂÂÂÂmm

� : ô en?Y�ô eH?Y@ = m ³¯�4@É?+��� m Ù me �vy|~��²¥�!��v����BvC}zrs��r�q�vy�¶rs|~������v�qsr�x �J�¨{b�¥xz�~�³q vy������}zx¤}zr�|����pq�vy�¥��v�qs�~��¡ M% ����| �¾r�}µt�q�{�}��«x�{ ô vy|~�Cxz�~�¥����}zr��9�¦vyq9|~{b�z��Ù�ru}¬}z��v�qsq ±bxz�~��|��M�+� �9�0� � �eru}¬}z��v�qsq ¡ % �~��|Cxz�~�����r�}�v¾�¦vy�z�,{b|pr�t â rsx�ª!��vyq����&|~�wv���@¬±��9�P�����0� � �µru}£q vy������¡ â {b�~����q���}zn��wv�º`rs|~�p±9²¤�~��|Jxz�~�³����}�ru�9��v�q|~{b�z��}�vy����}z��vyq�q ±�xz�~��|~{b�z�G{b��xz�~�¾���z��{b�&r�|�xz�~������t����z����|�t�������q vyxzru{b|�v���x����&t��¦vy|~�Q��{b��n�{�q��!r�}��{��p|��p���J�`� : �IvwÁ�­�����s £¢ p£¤ Û �2?¥@ Û CFE =Yà : ��vwÁ�­�����s £¢ p£¤ Û �2?�� Û CFE = m ¨iBïeÃ�É�Ʀ×�ÒOÖ�Ù�Æ�Û(Æ~ñ¥×�Ù�ÛjÃ`ØpÇCÉ�ÄOÉ�ÖwÃ9ò�É ¢ �~�o}�{bq��pxzr�{b|,{b��xz�~�¶qsrs|~�wv��c}z�9}zx����Ú�4�?É@)�B�ȽTeI����G$e�¦vb}�v¸����}zr��9�¦vyq�ÄWe,}�{¸xz�¦vyx,ÄWe7� ��G e�?²�4�µ?�@)�J�Ƚ§e�¡ ¢ �~�I�¦vbt�º`²¶vy���4���z��{b� Å ru}���� #�|~����`�v�4��?�@P��³ÆÅ��Ƚ§e;� ��G$e�±o}�{�Å(½§e¥� ÄUeI²¶rsxz�qj�Å j�Ç · ÒHj$��?�@P�.j­Ç �*/&{bq��p��v�|��K¥v�|T©�{�vy|(ÞyèQè�Ü��§¡ ¢ �pr�}&rs��n�q�r���}!xz�¦vyx!xz�~��©jv�|�t�ª�{Q}o�,��xz�~{9��r�}!t�{���np��xzr�|~�Ixz�~����ru�Q��|pn¦vyr���}!{b��4�O?.@)�7³%Å��DC�E5��rs|~}zx��wvb�7{b�o�4�O?.@)�B�DC�ED��¡ % �~��| �O?.@)�ýru}®rsq�q��¼t�{b|~�9rsxzru{b|~����±Qxz���¶�9r�R�������|~t������x ²¥����|Z�4�Q?Y@)�B�DC�ED�vy|~��4�Q?Y@)�ó�Å��DCFE,�¿tyvy|J�����Q�prsx��!�pr����O¡ % �&��v����!xz��v�x

�4�Q?¥@)�ʳ+Å�� CFE �Ã?��4�Q?Y@)�B� C�E �µ�k?B�4��?�@P�J� CFE Å2�4�Q?¥@)�ʳ�Ŭ� CFE �²¶��rut­��ru}Z|�{��z�¾²¶r�}��&{���xz�~�&{b���p���¥{��&��m~�4�Z?7@)�B�^j:�4�Z?O@P�¥³GŬ�DCFE$j$m�Ò«¡ ¢ �~�!©jvy|~t�ª�{Q}¬����xz��{9�np��{9�9��t���} Ë � vy|~�OÅ �D© E }�{Cxz�¦vyx�4��?�@P�ó�Å�� C�E �BÅ � ?;Å �D© E�Ë � �X²²¶rsxz�;j­².j pV�T}��Ivyqsqjt�{b��n¦vy�����Ix�{7j:�4�Ê?¥@)��³+Å�� CFE �.j pc¡£m³}­}z�p��r�|���xz�¦v�xo|~{����z��{b��r�}¶��vb�p�²¶rsxz��xz�~��t���v�|~�Q�!{b�jn�{bq��Q±�²¥�!��v����

�4�Q?����B� C�E �/( � ?�( �D© E¨- � ��4�Q?����B� CFE �4�Q?¥@)�ʳ�Ŭ�º²�� C�E ³¯�4��?Y���B� C�E Å�Å �D© EÌË � � CFE ¨¢ ���!q vb}zx¥x����z�?r�}&v�q�}�{��4��?Y���B� C�E Å5( �D© E¨- � vy|~�,ru}o{��exz�~��{b���p����{b�j:�4�Q?����B� C�E j$m�j$�Q?Y@)�!j$m1j - � j$m�Ò;� � Ý9¡ Ý��{b�¥��{������pqs�J}zn��wv�º`rs|~�p± : ��vyÁ�������¡ £¢ p¥¤ Û �{?Y� Û C�E = m ¨

Þ~<

Page 17: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

¢ vy�pq���Ý9¡�Þ�¼ ¯ qsqs�~}zxz�§vyxzru{b|0{b�c����{w²¤xz�T{b�«��{b�p|~��r�|~�����z��{���}³�¨{���à�Á�v���npq���Ý�¡�Þ�¡ % �C}z�~{�² ��{b�p|~��}�¨{b�¥xz�~���9r�R�������|`x&���z��{b��}¤r�|Ixz�~�!np�����`ru{b�~}&v�|�v�qs�9}zru}

bc{b�p|~�p} mo�~}­{�qs�px�� â ��q vyxzr��Q�� � Ý9¡·[1� � Ý�¡Â<>� � Ý9¡ Ý�� ���z��{�� ���z��{��å�¨ åQåpÞwç èM¨Í[¤Þwå9C�E º 89¨ ݤÞyå>CFE º Þ=¨Ó<oÞwå9C�E E Þ=¨Í0¤Þwå9C�E º Þ=¨ú8£Þwå9C�E0zçM¨ åQå Þ=¨ ܤÞwå9C�E m Þ�¨ å¤Þyå>CFE m 0M¨ ܤÞwå9C�E0z <�¨ ܤÞwå9C�E0z [M¨Í0¤Þwå9C�EÈ^

ÿ Ò�ò�Ã`Ç�Ù�ÕQØ~×£Ã1�OØ~òDðµ×�à ¢ �~�!�¨{bq�q�{w²¶rs|~���§Á�v���npq��&r�qsqs�~}zxz�§vyx���}¤xz�~�!xz����{b�z��rs|Ëxz��ru}o}���t�xzru{b|e¡x��OØ~ò�ð¬×�à ü-¸ ì â ��tyvyqsq!xz�~����ru�Q��|`��v�qs�~�¸np��{b�pq���� ����{�� à¬Á�vy��npqu�×[�¡�ÞQ¡ ¢ �~�°����}��pq�x�}�{b��xz�pr�}�§Á�v���npq��³vy���¶�Q��|~���§v�x������~}�r�|~�/2¸v�xzq�v���{b|Jv�á&à¥ß(moq�n��¦v�¡ % �¶n����z�¨{��z����� »!�ýÝ�å!}�x���n�}«{b��xz�~�}zn���t�xz�§vyq¬xz�§v�|~}��¨{��z��vyxzru{b|0©jvy|~t�ª�{�}&����xz�~{9�¸²¤r�xz�Y@¯�N?BÞ,vy|~�Tvy|�rs|prsxzr vyqZ����t�x�{���²¤r�xz�E������v�q��qu������|`x�}�¡ ¢ �~���§Á�npq�r�t�r�xzqs�Ët�{b��np�px����Ë���z��{b��{b|Jxz�~�³����t��p�z����|~t��!����q�v�xzr�{b|Ër�}

j�²!j^�O�²j:�4�Q?Y@)�B� CFE �BÅ � ?;Å �D© E�Ë � j^� y Ý�¨ 8£Þyå CFEª³²¶��rut­�J}§vyxzrS#~��}�� Ý9¡�Þ��§¡ ¢ ���¶����q�v�xzrs���o���z��{���j�²!j^�)��j Ë � j$mZru}c{b��xz�~��} v����¶{b���p���¥}zrs|~t��Bj Ë � j$m y ÞQ¡¢ vy�pq��¥Ý�¡�ÞZ}z�~{w²o}Oxz�~�¥��}zxzrs��v�x���}H� Ý9¡·[1�§±�� Ý9¡Ó<>�jvy|~�!� Ý9¡ Ý1���¨{b�jxz�~�£���z��{b�µ{b|Bxz���Z����t��p�z����|~t��c����q�v�xzr�{b|vy��x����!xz����t­�¦vy|~���C{b�«n�{bqu�¾�¨{��!x ²¥{�t��~{brut���}!{b�É��¡�·jrs��}zx�±w�� åM¨ åQåpÞwçCr�}!xz�~�C���yvy|�{b�cxz�~�Cx¼²¥{}z��vyq�q���}zx��¦vy�z��{�|pr�t â r�x�ª���vyq�����}�± ô E �?çM¨Ó<�<�çQè¤Þyå>C�^!vy|~� ô m.� ç�¨�Þyè�0Qç¶Þwå>C9z­¡Ë¹9��t�{b|~��±��¥�çqsru��}&rs|�xz�~�C��ru�p��qu��{b�����4�B�U�B�§¡ ¢ �~�C����q�vyxzr��Q�����z��{b���p{9��}o|~{�x��b��{w² rs|�xz�~�C}§vy���B²¶v���vb}&xz�~�vy�~}�{bq��px��Z���z��{b��±b}zr�|�t���j - � j meru}O��}z�¦vyq�qs�!{��9xz�~�c{b���p���j{b�9�IvwÁ ���¡ £¢ p¥¤W�zÞ�� Û �&?/� Û �§¡ % �Ztwv�|¾}����µ����{b�¢ vy�pq��«Ý9¡�Þ�xz�¦vyx)#p|~�9rs|~�oxz�~�¬q�v����Q�µ��r�����|`��vyq��~��}j{b�+�4�£?5���B� CFE ��r�}O|~{£np��{b�pqu��� ²¤�¦v�x�}­{9���Q����¡ ¢ �~�vbt�t��p�§vyx��!t�{���np��x§v�xzr�{�|�{��Oxz����}z��v�qsq����¤��r��Q��|`��vyq�����}£��ru�b�`x¶���&�,{b���!�9r�¿�t���q�x�¡ ¯ xoru}£r���n�{��zx§vy|`xx�{C|~{bx��Q±�xz��{��~�b�e±�xz�¦vyx¥��{�xz��xz�~�Bvy�~}�{bqs�px���v�|�������q vyxzrs���B���z��{b��}¥�¦v������b��{w²¶|O¡

¢ �~�!�¨{bqsqu{w²¤r�|����§Á�n����zrs����|�x&r�}³vy|�r�qsq���}zxz�§v�xzr�{b|D{b�jxz�~�¾�~vy|~�����o{��®nprut­º`r�|~��xz����n�{�q��!|~�wvy�³vy|��ru�Q��|`��v�qs�~��¡h�ÖwÃ`ðSf % �&�§vy|7».�ýÝQå¾r�x����§vyxzru{b|~}¶{b�jxz�~�¾©jv�|�t�ª�{Q}c����xz�~{9�J²¤r�xz��n�{�q��(@x�k?BÞQ¡h�ÖwÃ`ð�g % �cv�qsx���������xz�~�cn�{�q��¬r�|`x�{V�!� ô E ²¶�~����� ô E r�}Oxz�~�¥}��Ivyqsqu��}zxO�¦vy�z��{�|�rut â r�x�ªZ��vyqs�~��¡ % �n����z�¨{b�z�,����Ý�vb�p�9rsxzr�{�|¦vyq¦©jvy|~t�ª�{�}¬rsx����§v�xzr�{�|�}�¡ % �o��|~�p�����pn�²¶rsxz�Ëv�©ev�|~t�ª�{Q}�����t��p�z����|~t��

����q�vyxzru{b|Ë{b��{b���p����» �ýÝ�Ý9¡h�ÖwÃ`ð l % �¥vyq�x��������Bxz���cn�{bq��¥v��`v�rs|Brs|`x�{(@O�k?CÞ¥vy|~�¾n����z�¨{b�z�,����Ý&v��p��r�xzr�{�|�v�qp©jvy|~t�ª�{�}j}�x���n�}

}�{Bxz��v�x�»!� ÜQå9¡·jr����p����Ý9¡�ÞJ}��~{w²�}Cxz�~�I|~{b�z��}�{��oxz�~�����z��{b��}�{b��xz�~�I����t��p�z����|~t��Ë����q�v�xzr�{b|~}É� ç�¡ ç1�¾�¨{b�Jf®�

Þ=�~¨^¨^¨^�U»��¨{b�¥xz�~�!xz��������}zx���n~}�¡Zm&}£²¥��twv�|�}����Q±�²¥�!��v����Bvy��x�����» �ýÝ�å�©jvy|~t�ª�{�}«r�x����§vyxzru{b|~}�4��?Z@)�B� CFE �BÅ } º�?¥Å } EÌË } º �X² } º

²¶rsxz�¥j­² } º�j~� y ÝM¨ ݤÞwå9C�E0³§±�²¶��rut­�D} v�xzr�}5#~��}(� Ý9¡�Þ��§¡¯ |�¹`x���n¾b�±«xz�~��n�{bqu��r�}Ëvyq�x���������rs|`x�{;�²� ô EGy çM¨Â<1<oÞwå9C�^�¡ % �¸��Á�n���t�x�vT}zr��b|prS#�tyvy|`x�b��{w²¶xz�S{b�&xz�~�¸v���}�{�qs�px������z��{���rs|�xz�~������t��p�z����|~t��Ë����q�v�xzr�{b|S}zrs|~t��O�9��� �>�0� } º � y ÝM¨wÞwå>C°¯�v�|��

ÞyÝ

Page 18: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

��m`�0� } º � y Þ=¨ ܤÞwå�³§¡,mo��x������¦v��`rs|~�Ët­�¦vy|~�����¸xz����n�{bqu��±exz�~������t��p�z����|~t�������q vyxzr�{�|���v�}Bvy|0���z��{b�{b�¥xz�~�J{b���p���7{b�³ç¤Þwå m �¨}�����xz�~��n��wv�ºEv�x(fZ�ÝpÞ��§¡ ¢ �~����{b�p|~�¯� Ý�¡·[1�³ru}¾{��¤xz�pr�}7{b���p���7}zrs|~t����vyÁ �������¡ £¢ p¥¤ Û �2?¥@ Û � Û � ?�� Û y Þ=¨ ݶÞwå E m v�|��O�9��� �>�0� � � y ݤÞyå C°¯ ¡oão|�xz�~�7Ý�Þ�}zx¶rsx����§v�xzr�{�|O±�xz�~����z��{��Z�9��{bn~}®x�{�Þyå C�^ ²¶�pr�t���r�}Z{���xz�~�¶{b���p���«{b�~��vwÁ �­�����s £¢ p£¤ Û �H?%� Û C�E Ò«±`vy|~���9��{bn~}¬���Q��|������zxz�~����¨{b��f Ü Ý�Þ�¡m���x����jxz���c}���t�{b|~��t­�¦vy|~���c{b�9n�{�q���±yxz�~�c���z��{���}®�9��{�n¾�p�px®|~{�xex�{�xz���Zqu������qp{b�9��vyÁ �­�����¡ £¢ p¥¤ Û �o?@ Û C�EzÒ«¡ ¢ �~�¶Ý!vb�p�9rsxzr�{�|¦vyq~©jvy|~t�ª�{�}er�x����§vyxzru{b|~}�±b��{w²¥��������±Qnp��{9���~t��¥����t��p�z����|~t��¥����q vyxzr�{�|~}¬{b�~xz�¦vyx{b���p����¡

¢ �~�o�¨{bq�q�{w²¶rs|~�¾{��~}­���z��vyxzr�{�|�r�}¶vyqu}­{Brs|`x�������}zxzrs|~�p¡Z¹`�pn�n�{Q}��oxz�¦vyx¥²¥�&�p{B|�{�x¶��{�ÝCv����9r�xzr�{b|¦v�q©jvy|~t�ª�{�}®r�x����§vyxzru{b|~}Zrs|J¹`x���n~}�b4v�|��Jß�¡`mo��x����¥¹`x���n�b�±Q²¥�¤�¦v��Q�B� Ý9¡ ç1�§¡ ¯ |J¹`x���nËß�±�²¥��t�{b��np�px��� � mʤ� � gA³Ã�0�!?®@P� - � � � � ²¶rsxz� � �p|prsx§v��z�Tvy|~� � �pn�n�����xz�zr�v�|����pq�vy��¡�¹`rs|~t��C�¨{��B¹`x���n¥b- � � E7 C{8 : g�?���h ! �BCFEå #�= ±�²¥�¬#p|~��xz�¦vyx � �Ë��h&vy|~� � � ! �BC�Eå # ¡&¹9{p±�v���x����!¹`x���n0ß�±²¥���Q��x

�4��?Z@)�B� CFE ��( �Ì� ?�( �D© E �V� h - � � CFE �X²Ê³¯�4�Q?¥@)�J� CFE �/( �D© E � h ¹ � ¨À {�x��&xz��v�x¤xz�~��x ²¥{�x����z�,}¤rs|Ëxz�~�³�zru�b�`x5� �¦vy|~��}zru���!�¦v��Q��}z��v�qsq�|~{b�z�,}o{b�jxz�~��{b���p����{b�µÒ«¡ ¯ � �vy|~�

�vy��� çWß:û�ð-é�ã�âWá ����{��Í� � mʤ� ±�vy|Ë���z��{��¤r�|�xz�~�!����t��p�z����|~t��!����q�v�xzr�{�|�r�}¥rs|�xz��{9���~t����Ixz�¦vyx�twvy|���B{b�µxz�~�C{b���p����{b�cÒo��m~�0� � mʤ� � y Þ=¨ ܤÞwå�³�ÒZ¡ À {bx���xz�¦vyx���m��0� � mʤ} º �I�¾��m��0� } º �§¡ ¯ |¸� vbt�x�±�rs|��§Á~vbt�x

vy�zr�xz�p����xzr�t�� � mʤ} º ²¥{b�pqu�I���!xz�~�B�Q��|~���§v�qsr�ª����Ir�|`�Q����}­�B{b�L� } º ¡ ¯ |�np�§vbt�xzrut���±��~{w²¥��������±~²¥���p{�|~{bxvyx�v�qsqOq�{�}���ÜC�9r��br�x�}£�p�px�{b|pqs�ËçC}zrs|~t��

j=�4�Q?Y@)�B� C�E �/( �Ì� ?�( �D© E �V� h - � � C�E j^� y çM¨ ݤÞwå C�E¬^ ¨·p��{b�^xz��ru}��§Á�vy��n�qu��±�²¥��t�{b|~t�qs�~�p�¾xz�¦v�x�nprut­º`r�|~��v�n�{bqu�C|~�wvy�Bv â r�x�ªC��vyq��~�C��v��¸���7��v�|��

�Q����{b�~}�¡IbZ�px&npr�t�º`r�|��Jv�n�{�q���²¤�~��|°xz���7{bqu��n�{�q���r�}o|~�wv��!v â r�x�ª���v�qs�~�¾��v��J���7��v�|~�Q����{���}³vb}²¥��q�q ¡% ��vyqu}�{Ct�{�|�t�qs�~�p�³xz�¦v�x¤r�|¸v¾np�§vbt�xzr�twv�qjvyqu�Q{b�zr�xz�p� r�x¤ru}¤}z��¿�t�ru��|`x�x�{C��{b|pr�x�{b���9�+���>�0� � �§¡ ¯ �vyx!xz�~�A#p��}zx�t���v�|~�Q�C{��Zn�{�q���±&�9��� �>�0� � E ¤� �or�}!}z��vyq�q ±�²¥�C��v���rs|�xz��{9���~t��,v,q�v����Q��vy�~}�{bqs�px��7���z��{b�rs|Jxz�~�&����t��p�z����|~t��!����q vyxzr�{�|e¡ ¯ |�xz�~�&|~�§Á�x¶t���v�|~�Q�!{b�On�{bqu�oxz�pr�}¥��v���q��wv��,x�{�v¾q�v����Q�(�M�)���=�0� � mʤ� �vy|~�Ëvb��vyrs|Jq�v����Q��v���}�{�qs�px��!���z��{b��}¥��v�������}z�pqsx�¡Î " �¸% �&)c��L¡A¢4�<Ï]¡B%  �¢��¢ ��ru}��§Á�v���npq��Iru}�����q�v�x����(x�{�xz���°vbt�{��~}�xzrutI}zr��¾�pq�v�xzr�{�|�{b�!vTåM¨Ó<�Ð Õ4åM¨Ó<�Ð Õ(åM¨ å�Ü�Ð }§v���npq����v����!{��µv¾n�{���{:�¼��q�v�}�xzrut¤��v�x����zr�v�q ¡ ¢ �~�&��v�x����zr�v�qOr�}¥��{9�p��qsqu������}zr�|��,v¾x ²¥{=� np�¦vb}��(bZr�{�x¥��{9�p��qvbt�t�{b�p|`xzr�|��7�¨{b�¥º`rs|~���Ivyxzr�tBv�|�������t��¦vy|pr�tyvyqOrs|�x����§vbt�xzru{b|~}¥����x¼²¥����|�xz�~� �¨��q vb}zxzr�t��¥}zºQ��q���x�{�|�v�|��xz�~��n�{����×� v�t�{b�~}zxzr�t����p��ru�k� ¹�v�|~�����ýv�|~�;3¬r�}zx����¸Þyè>8bå�±¬¹`rs��{�|O± % �e±ÉÑ�ru��|pº`r���²¶r�t�ª¸vy|~�;3µvy�pqÞwè10QÜ1�§¡ ¢ �~�¶�¨{bqsqu{w²¤r�|��!��v�x����zr�v�q�np��{�n����zxzru��}¬�¦v��Q�¶������|�}���qu��t�x����]¼e�¨{��Zxz�~�o}zº���q���x�{�|e±�xz�~�ÉÒ£{b�p|~���{9�9�pqs�~}oru}�Þ^<`å�åQå�å�Ó �ÌÐ m ±¦xz�~�A3�{bru}�}­{�|��§vyxzru{JåM¨·[�Ý9±�vy|~��xz�~�C�p��|~}zrsx¼�(Þ�[�åQå1»�Ô���Ð%z§¡ ¢ �~�¾n�{b����p�pr��I�¦v�}¤�p��|~}zrsx¼�DÞ�¨ ç�çQÝ�»�Ô���Ð%z§±�xz�~��}�{���|~�J}zn�������r�}I[:<�å�ÐO��Á�±Qxz���!n�{b��{Q}zr�x �ËåM¨ èQÝ9±`xz�~���~{w²ý��� �}zr�}zxzrs��rsx �7r�}cå�±bxz���IbZr�{�x¬� vbt�x�{b��r�}oÞQ±�xz�~�H�p��ru���p��q�º¾�,{9���pq���}�Þ^<�ÞyÜQå�å�Ó �ÌÐ�ç�±�v�|��7xz���¥x�{b�zxz�~{�}�r�x �r�}BÞ=¨ ç9¡ ¢ �~���9r�}�t�����x��(#p|prsx��$�¼��q�������|`x���{9�p��qO����qsru��}o{�|�v é9R*ì �¨{b�z�B�pq�vyxzru{b|Z� ¹�r���{b|D��x&v�q ¡eÞyè�0QÜ1�²¶�������&}zºQ��qu��x�{b|��9ru}�npq vbt����,��|`x¥t�{���n�{�|���|`x�}�� é �cvy|~������q�v�xzrs�����p�pr����9r�}znpq�vbt�������|`x¥t�{b��n�{b|~��|`x�}��²¥��ru�b�`x����¸�`�°xz���Bq�{9twv�q®n�{b��{�}�r�x ���B� ì �!vy���7}���q���t�x�����vb}&|~{9�~vyqj��vy�zr vy�pq���}�¡ ¢ �~�J#p|prsx��$�¼��qu������|`x

ÞyÜ

Page 19: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

Õ×Ö�ØÕ×Ö�Ù�ÚÕ×Ö�Ù�ÛÕ×Ö�Ù�ÜÕ×Ö�Ù�ÚÝØÕ×Ö�Ù�ÚÝÞߺà<á e ß m

ÕºÖ âlÖ ã�Ö ä�Ö ålÖ æ�Öçéèºê¨ëºì�èºçîí�ï5ð ñ r aÊò�óñ r aÊò*ô

ñ r aÊò õ·jr��b�p���ËÝ9¡�Þ=¼ À {b�z��}7{b�¶xz�~�,���z��{b��}�{��¤xz�~������t��p�z����|~t�������q�v�xzr�{�|�� ç9¡ ç��§¡ ¢ �~��}­{�qsru��qsr�|~�I}z�~{w²o}j�²BdDe�j$mc�¨{��of2� Þ��^¨^¨^¨^�§ÝQå&�¨{b�Zxz�~�!©jv�|�t�ª�{Q}¬�,��xz�~{9��²¤r�xz�,n�{bqu��@x�k?BÞB� ¹`x���nJm(�§¡ ¢ �~�o�~v�}��~���qsr�|~�c}��~{w²�}�xz�~�Z���z��{b��}µvy��x����jxz�~�ct­�¦v�|����¬{��9n�{bqu�¬�¨{b��f2� Þ=�^¨^¨~¨^� ÝQå�vy|~�!�¨{b��Ý&v��p��r�xzr�{�|�v�q9©jvy|~t�ª�{Q}rsx����§vyxzr�{�|~}�f;��Ý�Þ��^¨^¨^¨ � ÝQÝ¥� ¹`x���nQbH�§¡ ¢ �~�J��{�xzx����(qsr�|~�I}z�~{w²o}Bxz���Ë���z��{b��{b|(xz�~������t��p�z����|~t������q vyxzr�{�|�²¶�~��|�xz�~��n�{bqu�!r�}ot��¦vy|~�Q���,�¦vbt�ºJrs|`x�{ @�v?CÞ�v�|��°Ý�vb�p�9rsxzru{b|¦vyq�©jvy|~t�ª�{Q}crsx����§v�xzr�{�|�}vy���&n����z�¨{b�z�,���×� ¹`x���nTßH�§¡����}z�Ë��v�}�[�ç:<C|~{9�p��}�vy|~��ÞyèQç��oàµä!m(0C��qu������|�x�}�¡ ¢ �~�!x�{bx§vyqO|����B�����&{b�µ�p���b������}¤{��®���������p{b�r�}La7� Þwè:<1<p¡% ���~}���xz���0©jvy|~t�ª�{�}���{b�pxzrs|~��àZmCÞyÜ�����{b� xz�~�T�&vy�z²¥��q�qB¹`���p��{���xzr�|~�E©�rs�p�§v��z�¾� â ��q��yvb}���&¹p©SçQåQå�å��§¡ ¢ �~��rs��npq�������|`x§v�xzr�{b|0r�}��¦vb}����T{b|Txz�������pqu{9t­º��B©jvy|~t�ª�{�}!����xz�~{9�D²¶rsxz�0n¦vy�zxzr�v�q����{��zxz��{��Q{�|�v�qsruªwvyxzru{b|��*/o�zrs�,��}¥��x¥vyq ¡¦ÞyèQè:<9�§±bn��p����rs|~�¾{���xz�~�o�p|`²¶vy|�x���� â rsx�ª�n�v�rs��}¥v�|���qu{9t­º�rs|~�{b��t�{�|`�������Q��� â rsx�ª!n¦vyr���}J�¨}����£XU<9�§±�vy|�r���npqsrut�rsx�����}zx§vy�zxB� ¹�{�����|~}���|EÞyèQèQç1�§±�vy|~�Ët­�¦vy|~���!{b��n�{bqu��¨}�����X$[1�§¡ ¯ |Sxz�pr�}��§Á�n����zrs�,��|`x�±o²¥���~}­�Ëxz����©jv�|�t�ª�{Q}�����xz��{9�4²¶rsxz�S�pq�{9t�º4}zr�ª��EÞËrs| {b���p���x�{�r�qsq���}zxz�§v�x���xz����xz�~��{��z�4rs|Sxz�~��n�v�n�����¡ mo��x�����np�p���brs|~��±¥xz�~��´&�z�`qu{���}z�p�~}zn�v�t����9r�����|~}zr�{b|r�}������9�~t����4����{�� »0x�{\�O¡ % �J�~}­���¯»¾�$ÝQåTvy|~���v�$çQÝ�¡ ¢ �~������xz�~{9�4r�}��p��}�t��zr������4rs|moqu�Q{b�zr�xz�p��Ü�¡�Þ�¡ ¢ �~�C}zn¦vy��}��!qsr�|��yvy�&}z�9}zx����,}!vy���B}�{bqs���������Ixz�~���³vy�z²¥��q�q¬¹`�p�p��{b�pxzrs|~��©¦r��p�§vy�z�t�{9�p� 2�m!ç98x� �&¹p© ÞyèQèQÜ1�§¡�·�{��C�$¿,t�r���|~t������wv�}�{b|~}�±jxz�~����v�xz�zr ÁEr�}¾� v�t�x�{b�zruª����T{b|~t����¨{b�B�wv�t­�n�{�q���v�|��Jxz�~�¾©jv�|�t�ª�{Q}c����xz�~{9���~}­��}£�¦vbt�º`xz�§v�|�}z�¨{��z��vyxzru{b|~}�¡ ¢ �~�!����}zr��9�¦vyqOx�{bqu���§vy|~t��&�~}������¨{b�xz�~����r�����|`��vyq��~��}¾²¶v�}'öÊ÷ùø×úH� r ÒHj Ë � j$mB²¤r�xz� r Ò y Þ�¨ å�ݤÞwå9C{û ±j²¶�pr�t��(r�}�xz�~�I}��Q�¦v�������{9{bx7{b�xz�~����vbt��prs|~��np����t�ru}�ru{b|0{b|»v�á&à¥ß moqsnp�¦v9¡ ¢ �~���Q{�vyqµr�}�x�{�t�{���np��x���Á��çQåJ��r�����|pn�v�rs��}�{b|xz�~�C�zr����`x!{��&Þwå�å9¡ â rsx�ª¾n¦v�rs��}�vy���Bq�{9t�ºQ���������¨{����¾xz�~�Cn�{�q��Br�}!t���v�|~�Q����¡ â ��twvyq�q®xz�~�B����}z��q�x�}�{b�à�Á�v���npq��I<p¡�Þc²¶�~�����¤xz�~��vbt�t����§v�t��7{b��xz�~�¥����t��p�z����|~t��¶����q�v�xzr�{b|Ivy��x����µq�{9t�º`rs|~�Bvy|~��t��¦vy|~�br�|��!xz�~�n�{�q��!�p��n���|~��}�{b|J²¤�~��xz�~���¶xz�~�!��v��z��{b|prut â rsx�ª!��vyq��~��}ov����!q�{9t�ºQ���������¨{b����{b��vy��x����¤xz�~��t­�¦v�|����{b�«n�{bqu��¡ ¯ |Dxz�pru}!t�{9�p��±~|~{J�9r�}zxzrs|~t�xzru{b|�r�}&��v��p�¾����x ²¥����| ô e�´ }�v�²¶v�������{���{b��t�q�{�}­��x�{,�Y¼o²¥�vyq�²¶v��9}£qu{9t­º������¨{b���!xz�~��t­�¦vy|~���!{b��n�{bqu��¡ ¢ �~�!rs|prsxzr vyqOn�{bq��!ru}�@� ÞwåQå�¡µm�|~��²Mn�{bqu�!r�}ot��~{Q}���|vb}¥xz�~�!���wvy|Ë{b�jx¼²¥{�}��~t�t���}�}�r��Q� â rsx�ª&��v�qs�~��}�±�vb}¥rs|~�9r�twv�x�����rs|�m�q��Q{��zrsxz�p�GÜ9¡�Þ�¡ % �³t­�~{9{Q}��&xz�~�|~��²�n�{bq��³}­{7xz��v�x Û ô en?Y� Û À Û ô e�?Y@ Û ��ÝQå¾�¨{��cf{� Þ=�^¨~¨^¨^�U»¦±�²¶��rut­�°}z��{��pq��Ëv���{bru�,q vy��������m��0� � �vy|~��4�Q?����B�DCFE,��¡f�×�ø�Æ~Ç�ÙsÖ�ï�ò þ+¸ ì ÎÏk�Ø9Ö�Ù�Æ�ÛjØ~×£ócØ~ÛjÕQô�Æ�É�Ñå9¡^/or��Q��|OG E uO° y }z��t���xz�¦v�xAj G E j$p×� Þ�¡¹9��qu��t�x³vy|Jrs|prsxzr vyqOn�{bq��(@Tv�|���� vbt�x�{b�zr�ª��(�Q?Y@)� �C��ü �ch�¡Þ�¡U3����z�¨{��z�"�Ë}�x���n�}�{b�jxz�~�B©ev�|~t�ª�{Q}c����xz��{9��¡ç9¡ È�Æ��p|`xzr�qjt�{b|`����������|~t��{¼ç�¡�ÞQ¡¦à�Á9n�v�|~�Ixz�~�B´��z��q�{w�J}�n¦v�t��&����{b�?�9rs�,��|~}�ru{b|��Ix�{!»����Ixz�~�B©jvy|~t�ª�{�}«�,��xz�~{9��¡

Þ:8

Page 20: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

ç�¡ ç�¡eßc{b��np�px��!xz�~����r��Q��|�n¦v�rs��}��4±We��,ÖWe$�§f2� Þ=�~¨�¨�¨S�U»I{�� Ë � vy|~�Jq���x ô eI�¯@!³Y± C�Ee ¡ç�¡·[�¡�¹9{b�zx¶xz�~� ô e�´ }¤r�|¸vb}�t���|��9r�|��7{b���p���³vy|~������{b���p���¶xz�~�¬ÖWe�vbt�t�{b���9rs|~��qs��¡ç�¡Â<�¡eßc{b��np�px��!����}zr��9�¦vyq�|~{��z��}�Ù�eL�;¹ ��Û d�h� ÖWe Û �¨{b��f2� Þ��^¨^¨^¨~�U»¦¡ç�¡ Ý�¡�¹`x�{bnËrs�oÙ�e��XöI÷ýøÊúÌj Ë � j$m��¨{b��f2� Þ��^¨^¨^¨^� Áb¡ç�¡ Ü�¡�©O{9t�º�xz�~�¥t�{�|`�Q�����Q��� â rsx�ª¥n¦vyr���}¬v�|��Bxz�z�p|�tyvyx��¥xz���¶©jvy|~t�ª�{Q}j����t��p�z����|~t��¥����q�v�xzr�{b|7����{b�

{���������»�x�{V�O¡ç�¡ 8�¡eßc{b��np�px��&xz�~�&|~��²ýn�{bq���vb}��!�¾� ô eT³ ô e © E �D�bç�²¤r�xz� f�ÀX�¾v�|��J}z�~t­�Jxz�¦vyx Û ôÌþ ?&� Û ÀÛ ô�þ ?Y@ Û �bÝ�å¾�¨{b�|ÿo� Þ=�~¨^¨^¨^�ª�e¡ç�¡·0�¡eßZ�¦vy|~�Q�&xz�~�!n�{�q��&�~}zrs|~��moqu�Q{b�zr�xz�p� [�¡�ÞQ¡�·¦vbt�x�{b�zr�ª��(� ?����µ�X��ü���hÃ�ÛeÈ�È�Æ·p��{b� � Ý�¡ ç��§±�²¥�Btwv�|�}����!xz�¦vyx£j:�4�Ê?����B� CFE �4�Q?¥@)�B�^j$p(npq�v��9}³vy|Ërs��n�{b�zx§vy|�xo��{bqu�!rs|�xz�~�

���z��{��¤��{b�p|~��}�¡ % ����| Û �{?+� Û À Û �{?Z@ Û �bÝ�å��¨{b�³vyq�qO��ru�Q��|`��v�qs�~��}���{b�o�����Ã�F�¬�e±�xz�~��|j:�4�Q?����J� CFE �4�Q?Y@)�J�^j pZ� ��vyÁ�������¡ £¢ p¥¤ ¦¦¦¦¦

�A?¥@�2?�� ¦¦¦¦¦ �ýÝ�å ¨¢ ��ru}�qsrs��rsx�}�xz�~�¸�b��{w²¤xz�S{��&xz�~�¸���z��{b��r�| xz�~������t����z����|�t�������q�v�xzr�{�|ýv���x����Jv�t��¦vy|~�Q�Ë{b�&n�{bqu��¡©O������v¸Ý�¡�Þ�r�}B�~}­���E�¨{b����}zxzrs��v�xzrs|~��9��� �>�0� � �§¡ ¢ �pr�}B²¶vb}Cr�qsq���}zxz�§v�x����E�`�»à�Á�vy��n�qu��Ý�¡�ÞJv�|��r�}!t�{b|�#p�z�����Tvb��vyrs|°�`��xz�~�C����}z�pqsx�}&�¨{��&xz�pr�}��§Á�v���npq���¡ ¢ vy�pq���Ü9¡�Þ¾}z�~{w²o}�±jvy��x����!�yvbt��¸rs��npqsrut�rsx����}zx§vy�zx�±cxz���Ë|`�p�¾������{b�³t�{b|`����������� â rsx�ªJ��vyq�����}�vy|~�(xz�~�I|~��² n�{bqu�,²¶rsxz��xz�~����}�xzr���vyx��J�¨{b��M�+� �>�0� � �§¡4mnpr�t�xz�p���Ë{b�oxz�~�Ë}�n���t�xz�z�p�$ru}�}z�~{�²¶|(r�|�·eru�b�p����Ü�¡�Þ�¡4½o}zrs|~�¸xz�~�Ir�|����zxzr�v�t�{b�p|`x�*/o�zrs����}£��x¥vyq ¡~ÞyèQè:<9�§±b²¥�¤�¨{b�p|~��xz�¦vyxZxz�~�o|����B�����¥{b����r��Q��|`��vyq�����}Z����x¼²¥����|�Þwå�å!vy|~��xz�~�on�{�q���}Þ�[10Qå7vy|~��Þ:8bÜQå���vyx�t��~��}¥xz����|`�p�B�����³{b�jq�{9t�ºQ��� â r�x�ª!��vyqs�~��}¶rs|Ëxz����t�{��z����}zn�{�|��9r�|��7rs|`x����z��v�q�}�¡¢ ��ru}¤}z�~{w²o}¥xz�~�!����qsr vy�prsq�rsx �Ë{b�jxz�~�Bvyqu�Q{b�zr�xz�p��¨{b�¶xz��ru}¤��Á�vy��npqu��¡

¢ �~�!���z��{b��{b|Ë�¦vy�z��{�|�rut â r�x�ª³��v�qs�~��} ô ru}£��{b�p|~�p����v�}£�¨{�qsqu{�²�}�¡ócÃ9òTòDØ þ�¸¨ìCB â�ã�ôÃëDâ�ä ù ä:àWû2ß=Þ�æ4ç�í�æ4ãïîÉå�ä�óÍé�â�çWß:àWà5â^ê�ð�ß=ÞVá=æ¶ÞMüã�ß!ð�ß�óúâ @ ä:Þ-á×óúâ�ã Ù ëDâÉã ù âà5â^êWæ4á:é�ä1ó�ÞVß=àUû\è$ß=àHã ù âcí�æ4ãïîlð�ä:æ¶àHß,è�ã ù âLê�ð�âWç�ã"à,ä�ó�ã4à5ä=Þ�ê"è$ß:àWû2ä1ã"æ4ß=Þ ÷�ø�ù â^Þ ã ù â^à5âLæ¶ê�ä ��u×���4�A�U�J�ê ß.ã ù ä1ã

Û �{? ô Û � Û ô ?Y@ Û m Ùx¨÷ Ç�Æ�Æ~ñ ¢ �~�!n���{9{��jr�}¶�brs����|Ërs|�X$[7{b��à¬�zr�t�}­}�{�|�v�|�� â �p�~� �zÞwè10Qå��§¡ M·jr����p����Ü9¡ ç7}��~{w²�}«xz�~�!��{��p|��p}¶{b|Jxz�~�!���z��{b��}�{b|Jxz���&�¦v��z��{b|pr�t â rsx�ª&��vyq��~��}¬� ©O������v�Ü9¡�Þ��

²¶����|�xz�~��n�{bqu�!t­�¦vy|~����}�¡ % �!}z��{w²�xz�~�!����}z�pqsx�}¶�¨{b�¶xz�����p|pq�{9t�ºQ��� â r�x�ª³��v�qs�~��}¶rs|�xz�~��rs|`x����z��vyqÆsÞyåQå�åM�,<�ÝQåQå Ç �¨{b�!xz�~�A#p��}zx!x ²¥{°t­�¦vy|~����}&{b�cn�{bqu��¡ ¯ |���{�xz��tyvb}���}�±¦xz�~�C��{b�p|~��}!�p{J|�{�x�t­�¦v�|���������z�(�B��t���²¤r�xz�SvDt��¦vy|~�Q�J{b�³n�{bqu��±µ|~{����p{Dxz�~�I�¦v��z��{b|pr�t â r�x�ªI��vyq��~��}�¡ ¯ | vb�p�9rsxzru{b|e±«xz�~��¦vy�z��{�|pr�t â r�x�ª!��vyqs�~��}¶|~�wvy�L<`å�åQåB}­�����x�{��9r�}§vynpn��yvy��vy��x����¥xz�~��t­�¦v�|����!{b�jn�{�q���¡

m���x����cxz�~�oq�v�}�xcr�x����§vyxzru{b|e±�²¥�&�¨{��p|���xz�¦v�x�±�²¤r�xz�I���pq�q�����{��zxz��{��Q{�|�v�qsruªwvyxzru{b|�v�x¥�wvbt���rsx����§v�xzr�{b|{b�jxz�~�B©jvy|~t�ª�{�}«np��{9t���}­}�±��IvwÁ þ ¢ e��� þ Û ½�hþ ��½§e Û y ÞZÞyå>CFE¬^�vy|~�j:�4��?�@P�J� CFE ��( � ?�( �D© E¨- � j$p-�&��j - � j$m y <�¨ è¤Þyå C9û �²¶��rut­��ru}�²¤r�xz��r�|�xz���µ����}zr��9�¦vyq�x�{bq����§vy|~t���¡ % rsxz��n¦vy�zxzr vyq�����{b�zxz�~{��Q{b|¦v�qsr�ªyvyxzr�{�|e±§²¥�¬�¦v����¬��vwÁ þ ¢ e��� þ Û ½ hþ �¬½§e Û y[M¨ ݤÞwå C�E E vy|~�

j:�4��?�@P�J� CFE ��( � ?�( �D© E¨- � j$p-�&��j - � j$m y <�¨ è¤Þyå C9û �Þ`0

Page 21: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

¢ vy�pq��°Ü9¡�Þ�¼ À ���B������{b�&qu{9t­º���� â r�x�ª���vyq��~��}�±«xz�~��|~��²n�{bqu��vy|~�(xz�~����}zxzr���vyx���{����M�+� �>�0� � ��¨{bq�q�{w²¶rs|~��©O������v�Ý9¡�Þ

����}zx§vy�zx q�{9t�º���� n�{bq�� ��}�xzr���vyx��B�M�+� �9�0� � �Þ Ý Þ`[�0�Þ�¨ 0M¨·[¶Þwå C mç Þwç Þ:8bÜ:<�¨ Þ=¨·0¶Þwå>CFE[ Þ�8 ÞyèQÝ�[�¨ Þ=¨·[¶Þwå>CFE< çQÝ

ÞyåQå Þ�[�0pÞ Þ�8�Ü:< ÞyèQÝ1[·jr��b�p���JÜ9¡�Þ=¼ â r�x�ª���vyq�����}Ct�{���np��x����D�`r v�xz�~���§v�xzr�{�|�v�qc©jvy|~t�ª�{Q}!����xz�~{9��¡ ¢ �~���Q���zxzr�twv�qcqsrs|~��}�p��|~{bx��!xz�~�!n�{Q}zr�xzr�{b|~}¥{b�jxz�~�!n�{�q���}�¡

�Õ×Ö�ÛÕ×Ö��Õ×ÖÌÙ�ÛÞwå�åQå ç�åQåQå [QåQå�å <`å�åQåô

ÞyåQå�å çQå�åQå [Qå�åQå <�åQåQåô

·jr��b�p���¥Ü�¡ çM¼ ¢ �~�Znpr�t�xz�p����}�}z�~{�²0xz�~�Z��v��z��{b|prut â rsx�ªZ��v�qs�~� ô vy|~��xz���Z�¦v��z��{b|pr�tZ����}�ru�9��v�q`|~{b�z��}����{b� ©O������v�Ü9¡�Þ������¨{b�������p�pqsq���x�}U�Bvy|~�»vy��x����!�¨t�r���t�q���}W�!xz�~�,t��¦vy|~�Q��{��¥n�{bqu��¡ ¢ �~��qu����x¾�¦v�|��npr�t�xz�p����}z�~{�²�}¥xz���!xz�§v�|�}zr�xzr�{b|J�¨{b�¥xz�~��t­�¦v�|����!{b�o@7� ÞwåQå�rs|`x�{��!� Þ�[10�Þ!v�|~�,xz�~�!�zr����`x¤�¦v�|��{b|~�&�¨{��¤xz�~��t­�¦v�|����&����{b�N@x� Þ�[10�Þ¶rs|`x�{ �!� Þ:8bÜ=<p¡

Þyè

Page 22: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

% �7vyq�}�{�n����z�¨{b�z������xz���B�¨{bq�q�{w²¤r�|~���§Á�n����zrs����|�x�¡ % ���§vy|�v�q���{b�zrsxz�p� Ü9¡�Þ!�p�px!rs|T¹`x���nTç�¡ 8�±²¥��}���q���t�x¤xz�~��n�{bqu����� ô �È© E[y Þwç�Ü:<�¨�ÞQÞyç9¡cbc���¨{b���&xz�~�Bt­�¦vy|~����{��jn�{bqu��±�²¥����v����!²¤r�xz��n¦vy�zxzr�v�q����{��zxz��{��Q{�|�v�qsruªwvyxzru{b|

j:�4��?�@P�J� CFE �BÅ � ?ZÅ �D© EÈË � j$p-� y Þ�¨ ç¤Þyå CFE mj:�4�Q?Y@)�B� CFE �BÅ � ?;Å �D© E�Ë � j$p-�&��j Ë � j~� y [�¨ 8£Þyå CFEÈE ¨À {�x��oxz�¦vyxI�9��� �>�0� � � y Þ=¨Ó<oÞwå9C{³§±F��m��0� � � y è¤Þwå�³¤v�|~�xj=�4�;?&���B�DC�E^�4�\?×@)�B�^j$pYÀ²j^gI³\�0�J?@)� - � j m y 8>¨�ÞZÞwå } ¡ % �B�§Á�n���t�xoxz�~��©jv�|�t�ª�{Q}¥����t��p�z����|~t��B����q�v�xzr�{b|°x�{�qu{Q}���Þwç��9r���rsx�}!v���x����oxz�~�t���v�|~�Q��{b�µn�{bqu��¡o¹`rs|~t���j - � j~�°ru}&{b��xz�~�¾{����p���&{��oÞwå�z�±�²¥�¾�§Á9n���t�x³x�{�ºQ����n�{b|pqs�TÞ�}zr��b|prS#�tyvy|`x�9r���rsx�¡ ¢ �~��t�{b��n��px����Ë���z��{b��±��~{w²¥��������±�r�}¥�B��t���q vy���Q���¤}zr�|�t��

j:�4�Ê?����J� CFE ��( � ?�( �D© E¨- � j$p-� y Þ�¨ ç¤Þyå zj:�4�Q?����B� C�E �/( � ?�( �D© E¨- � j pV�&��j - � j~� y ç�¨ å¤Þyå º ¨% �oºQ����n�|~{�v�t�t��p�§vbt��Ivyx�vyq�q ¡¬·�{��¶t�{b��np�pxzrs|~�Cxz�~�!x����z� �4�¯?&���B�DCFE,��( � ±�²¥�&�~}­���,r�x����§vyxzr��Q����$#p|~������|�x���}zr�|��Z/V2 â ࣹErs|�{b���p����x�{D�����9�~t��Jxz���J����}zr��9�¦vyq¶|~{b�z� {��oxz�~�Iq�rs|~�wv���}z�9}zx�����}�¡� ¢ ���J���wvb}�{�|E²¶�`�E²¥���~}­��� /V2 â àc¹(v�|~�E|~{bx�2 ¯�À&â ࣹ�r�}Cxz�¦v�x�²¥�,�¦v��4v¸t�{9�p�Iv���v�rsq�v��pq���¨{b�³xz���B�¨{b�z�����!����xz�~{9��¡ % ���~}�����/V2 â àc¹��¨{��&xz�~�C�Q���zr�#~twvyxzru{b|�{��Zxz�~�¾����t��p�z����|~t��7����q�v�xzr�{b|{b|pq���±�}­{�n����z�¨{b�z�Ivy|~t���r�}o|~{bx�vy|Ër�}�}z�����~������¡Í� ¢ �~����v��z{b�¶n���{��pq���� �~������ru}oxz�¦vyx³���Q��|¸²¤r�xz��v}z��vyq�q�����}�ru�9��v�q�|~{b�z� xz����qsrs|~�wv��¥}z�9}zx�����}¥vy���³}­{�qs������� vyr��zqs��r�|¦vbt�t��p�§v�x���qs��±9}zrs|~t�����m��4�Z?~���J� yçM¨ ݤÞwå m Ez¡�bc���¨{����!xz����t��¦vy|~�Q��{b�jn�{bq��Q±�xz��� â r�x�ª!��vyqs�~�B± �È© E �¦vb}�vC����}zr��9�¦vyqO|~{��z� {��exz�~�B{b���p���{b�µÞ=¨Í[QݤÞyå>C9û§¡ ¢ �~� â rsx�ª¶��vyq����³vynpn��yvy��}¬x�{!���¶�¾�~t���t�q�{�}����¬x�{Bvy|���r��Q��|`��vyq����¥xz�¦vy|�xz�~�¤����}�ru�9��v�q|~{b�z�?}z�~�Q����}zx�}cvy|~���¥?���� r�}¶vyq���{Q}zx¥}zrs|~����q vy��¡cßc{b|`xzr�|`�prs|~�Cxz�~��©jvy|~t�ª�{�}¬����xz�~{9��²¤r�xz�Ixz�pr�}n�{�q��&ru}£|~{bx³vC��{9{9��ru�p�wv�}zrs|~t��!�wv�t­�Ët�{b|p|~��t�xzru{b|J²¤r�xz��xz�~��{b�zr���rs|¦vyq¦np��{b�pq����?ru}¤q�{�}zx�¡

¯ �c²¥�C�p{�t�{b|`xzr�|`�~�C²¤r�xz��xz��ru}!n�{bqu��v�|��¸�~}��)�� � ô � ³ ô � © E �D��ç y Þ`[�0pÞBvb}&|~�§Á�x�n�{�q���±O}�{xz�¦vyx��9��� �>�0� � � y çM¨ å9±�xz�~��|˲¥��#p|~�j=�4�Q?+���B� C�E ��( � ?+( �D© E.- � j$p-� y Þ�¨·[¤Þyå C9zj=�4�Q?����J� CFE ��( � ?�( �D© E¨- � j$p-�&��j - � j^� y Ü�¨ ݤÞyå CFE �²¶��rut­�¸ru}¤�����z��q vy���Q�Q¡ ¢ �~�Bt�{b|p|~��t�xzru{b|�²¶rsxz�¸xz�~�B{b�zr���rs|¦vyqOnp��{���qu��� �����µ�����¸r�}&t�{���npq���x���qs�q�{�}zx�¡ ¢ �~��{�|�q���²¶v��¸x�{�}§v��Q��xz�~��}zrsxz�¦vyxzru{b|�r�}�����}zx§vy�zx�xz�~��©jvy|~t�ª�{�}&����xz�~{9�T�§Á�npqsrut�rsxzq��¸²¤r�xz�xz�~�!|���²�n�{�q��5� y Þ�[10�ÞQ¡� é #�êI¢^�¸�¬��#�ê¸�¯ |¸xz�pru}on¦vyn�����±p²¥��rs|�xz��{9���~t�����xz�����§vyxzru{b|¦vyq�©jv�|�t�ª�{Q}¶����xz�~{9���¨{b��xz�~�¾}�{bq��pxzr�{b|°{b�£�o���z��r�xzr�vy|��ru�Q��|`��v�qs�~�!np��{b�pq����,}�¡ ¯ x&twv�|�����t�{b|~}zr��p�������°vb}&v�}zn���t�xz�§vyqOxz�§v�|~}��¨{��z��vyxzru{b|°©jvy|~t�ª�{�}£����xz��{9�²¶rsxz�Cr���npqsrut�rsxµ����}zx§vy�zx«vy|~�Ct��¦vy|~�Q�c{���n�{bqu��±Q{b�µvb}jxz�~�«�§v�xzr�{�|�v�q~´&�z�`q�{w������xz�~{9�¾�¨{��µv!�o���z��r�xzr�vy|np��{b�pq����Ë¡ ¢ ���B��v��z{b�&t�{b|~t�qs�~}zr�{�|�����{b� xz�pr�}&n¦vyn����&ru}&xz��v�x&rsx û2ä:òOëDâ �~vy|~�Q����{b�~}ox�{J}­��q���t�x�vn�{�q���t�qu{Q}���x�{�v â rsx�ª���vyqs�~���¨{b�Bxz�~�I�¨{bq�q�{w²¶rs|~�����wvb}�{�|�}�¡D·jrs��}zx�±Zxz����qsr�|~�wvy��}���}�x�����}C��v��T����9r�¿�t���q�x�x�{�}�{bqs����±e}zr�|�t���xz�����Tvy���C|~�wv��zqs��}zrs|~�b�pq vy��¡,¹9��t�{b|~��±ersx���v��¸|~{bx�{b|pq��¸q��wv��°x�{Ëq�vy�����vy�~}�{bq��px���vy|~�¸����q vyxzrs�������z��{���}!rs|�xz�~�C����t��p�z����|~t��C����q�vyxzru{b|e±��p�pxBvyqu}­{�x�{�q�vy�����C����}zr��9�¦vyq�|~{b�z��}�¨{b��q�{9t�ºQ��� â r�x�ª�n¦vyr���}�¡ ¯ x&ru}¤xz�~�������¨{����Cv�����r�}�����x�{�npr�t�ºIxz�~�Bn�{bqu�Bv�²¶v��J����{b� xz�~��qu{9t­º���� â rsx�ªn¦vyr���}�¡

ç�å

Page 23: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

2�{w�`r�|~�Dxz�~��n�{�q���|~�wv��Jv â rsx�ª���v�qs�~����v��(|~{bx�{�|�q��4q��wv��4x�{(vy|�r�|�t����wvb}���{��!xz�~�¸���z��{b�vy��x�����xz�~��t��p�z����|�xCt��¦vy|~�Q��{��cn�{�q���±��p�pxB��v��¸n���{�n¦vb�`v�x��¾x�{�xz�~��|~�§Á�x7{b|~��¡�m&}Bà�Á�vy��n�qu� <p¡�Þrsq�qs�~}zxz�§vyx���}�±`xz��������}zr��9�¦vyq�|~{��z�{b�jv!q�{9t�º���� â rsx�ª�n�v�rs�Z�Iv��,����t����wvb}��¶²¤�~��|�xz�~�o�¦vy�z��{�|pr�t â rsx�ª��vyq��~�&r�}�v�²¶v�������{��xz���&|~��² n�{�q���vy|~����v����b��{w²S��v��,r�x¥r�}¶t�q�{�}­��x�{Cxz�~�&|���²ýn�{�q���¡ % �&�¦v������v�����|~{���ru}zxzrs|~t�xzr�{�|(����x ²¥����|»q�{9t�ºQ���(�¦v��z��{b|pr�t â rsx�ªI��v�qs�~��}C|~�wv���{���v�²¶v�������{�� xz�~��|���²n�{�q���¡ ¢ �~�C��v��z{b�&���wv�}�{b|¸r�}!x�{J�����9��t��Cxz�~��t�{b��n�qu�§Á�r�x ��{b�«xz����}�{���x ²¶vy���Q¡ ¯ |(v��p��r�xzr�{�|O±O}zrs|~t��xz�~�£n�{bqu��}Zv����¤t��~{Q}���|�v�²¶v��B����{b� xz�~�£�¦vy�z�,{b|pr�t â r�x�ª£��v�qs�~��}�±�v&}zru�b|pr�#~twv�|`x¬�p��t����wvb}��¥{b�«�b��{w²¤xz�{b�j����}zru���¦v�q�|~{b�z�?r�}¶�§vyxz�~���¶��|pq�rsº���qs��¡

ão|~�¶twvy|�npr�t�º�n�{bq���}¬|~�wvy�Z�¦v��z��{b|pr�t â r�x�ª¤��vyq��~��}�±��p�pxZ²¥�oxz�~��|�}z�~�Q�Q��}�x£vy|��§Á9n�q�r�t�r�xc����}zx§vy�zx{b�«xz���,©ev�|~t�ª�{Q}�����xz�~{9�¸²¤r�xz�Dxz�~��|~��² n�{�q���v���x����&qu{9t­º`r�|~�Ixz�~��t�{�|`�Q�����Q��� â r�x�ª¾n¦vyr���}�¡ À {bx��xz�¦vyx&nprut­º`r�|~�,n�{bq���}&t�q�{�}���x�{Ëv���v��z��{b|prut â r�x�ª¾��v�qs�~�¾��v��°}�n��������pn¸xz�~�Ct�{b|`�������Q��|�t��C{b�«{b|~�{b�¤v³�¨��² â rsx�ªon¦vyr���}�±��p�px¶v³q vy�������Z|`�p�B�����¶{b��n�{bqu��}¬��v������¤��������r�������x�{Ct�{���np�px��¤xz�~�&�p��}zrs�����|`�p�B�����³{b����r�����|pn¦vyrs��}¥� vyr��zqs�Jv�t�t��p�§vyx���q���¡Z·���{�� v�|�� ¿�t�r���|�t��,n�{br�|`xo{��j�`r���² {b|~�!}z��{��pq���ºQ����nxz�~�!|`�p�¾�����&{��en�{bqu�!t­�¦vy|~����}(��r ¡ �Q¡µxz�~�!|`�p�¾�����o{���}�n¦v���}­�&�Ivyxz�zr Á�� vbt�x�{b�zr�ªyvyxzr�{�|~}U�Z}z��v�qsq ¡� ���ê¸#�� ¢��! ]���&% �&ê�ëµ�¢ ���³vy�pxz�~{b�Zru}Z�b�§v�x������pq�x�{�6���|p|�r��¨���¶¹9t�{bxzx�± ¯ v�rs|Jáo��R�±9vy|~��m¶Áp��q â �p���¶�¨{b�Z�~��q�np����q¦}z�~�Q�Q��}�xzru{b|~}xz�¦vyx¥rs��n���{w�����Ixz�~�!�Q�¦v�qsrsx¼�J{b�jxz�~�&n�v�n�����¡ ¢ �~��vy�pxz�~{b��vyqu}­{B²¤ru}z����}«x�{7xz��v�|pº�xz���!�����¨��������}c�¨{b�xz�~��rs�¶rs|�x�������}zxzrs|~��t�{�������|`x�}¶xz�¦vyx¶rs��np��{w�Q����xz�~������v�qsr�x �J{b�jxz�~�!n�v�n�����¡

����¤��&)Z�&ê¶�&�mo�z|~{bqu�9r ± % ¡l�zÞyèQÝ�Þ`�§±�� ¢ �~��np�zrs|~t�r�npq��B{b����rs|prs��r�ª�����r�x����§vyxzru{b|~}¤r�|�xz�~�C}�{bq��pxzr�{b|°{b��xz�~����v�xz�zr Á

��r�����|`��vyq��~�³np��{���qu���]´�±� é�ä:à ã ÷ D�ð�ð�ó ÷� ä�ã ù�÷ ��±®Þ�8^ÄpçQè�¡b)��¸{b��t�º¦±��mB¡F�zÞyèQè�Ü��§±�� éMû{â^àUæ4çWä�ó�û2â�ã ù ß�á=ê)è$ß=à�óúâUä=ê$ã)ê��^é�ä:à5â^êoðVà,ß>ë^óúâ^ûAê ±�¹ ¯ mV2T±13µ��r�q�v�����qsnp�pr�v(3�m

Þwè�Þyå:<=�Ïç�Ü�0�0�±Q½&¹9mB¡ß¥vyqs����xzxzr ±QáB¡s± â ��r�t��~��q ±b©µ¡�vy|~�B¹�{�����|~}���|e±báB¡1�zÞyèQè:<9�§±�� mo|!rs��npq�r�t�r�xzqs������}zx§vy�zx����B©jvy|~t�ª�{Q}�����xz��{9�

�¨{b�¶q�vy�����!}z�`������xz�zrut���r�����|`��vyq��~�³np��{���qu����}~´s±��Bø�� D î ±eÞ$Ä�çpÞ�¡á!vy|pru��q ±l6�¡s±�/o�§vb�Q��± % ¡s±¬´�vy�p����v�|O±�©¬¡�vy|~�E¹`x���²¶v��zx�±�/¾¡��zÞwè>8�Ü��§±�� â ��{��zxz�~{Q�Q{�|¦vyqsruªwvyxzru{b|Tv�|��

}�x§v��pq���v�q���{b�zrsxz�p��}��¨{b�o�pn���v�xzrs|~��xz���./o�§vy�2��¹9t­�p��ru��x7æ â � vbt�x�{b�zr�ªyvyxzr�{�|P´s±� ä�ã ù�÷�P ß=û�ð ÷õ�� ±)8�8bç`Ä�8�èQÝ9¡

á&�7¹�v��¾�pq vy|�Á�±&/¾¡�± 2������z����������|e±�´C¡Ovy|~�7bZ�pqsxz�~����q ±�mB¡l�zÞwèQè98:�§±p� ¢ �~�¾r���npqsr�t�rsx�v�n�npq�r�twv�xzr�{b|D{b�v¾�§vyxzru{b|¦vyq)#pqsx����¤r�|�xz�~� â ´�¹��,��xz�~{9�+´�±�� Ý ø õ � ±�è�çQÝ�Äpè:<�8`¡

à¬�zrut�}�}­{�|e± ¢ ¡�v�|�� â �p�~��±�mB¡P�zÞyè�0Qå1�§±�� ¢ �~��}�n���t�xz�§v�q�xz�§vy|~}z�¨{b�z��v�xzr�{�|¸©jvy|~t�ª�{�}c����xz�~{9�I�¨{��¤xz�~�|`�p�����zr�tyvyqp}�{bq���xzru{b|B{b�9q vy�����c}�n¦v���}­�«�Q��|~���§v�qsr�ª�����}���������xz�zrut£��r��Q��|`��vyq����Znp��{���qu����}~´s±� ä�ã ù�÷P ß=û�ð ÷ õ¦ü ±eÞyçQÝ�Þ$ĦÞyçQÜ�0�¡

/&{bq����e±)/B¡�vy|~�ÉK¥v�|�©O{`v�|e±�ß�¡P�zÞwè�èQÜ1�§±� ä�ã4àWæ��ÉçUß=û�ðVé�ãÌä�ã4æ4ß:ÞFê ±V[y��������|e± ¢ �~�B6Q{���|~}³�o{bnpº`r�|�}½¤|prs������}zrsx �]3¬����}�}�¡

çpÞ

Page 24: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

/o�zrs����}�± â ¡s±!©���²¤ru}�±�6�¡ovy|~�S¹`rs��{�|O±!��¡��zÞyè�0�Ü��§±¤à¬ru�Q��|`��v�qs�~��np��{b�pq�����}�v�|~�ýv�q��Q{��zrsxz�p��}�rs|}�xz�z�~t�xz�p�§v�q���|���rs|~�����zr�|~�p± æ¶Þ 6�¡`ßµ�pq�qs�p� v�|~� â ¡ % r�qsq�{��~�b�`�`�`±b���p}�±9� ©ev����Q�c¹9twv�q��¥à¬ru�Q��|`��v�qs�~�3¬��{b�pq�����}~´s±�à¬qu}­���`r����³¹9t�r���|�t��¬3µ���pq�r�}z�~����}�b�¡�K�¡s±�npne¡V0�Þ$Ä�è�[�¡

/o�zrs����}�± â ¡�±³©O��²¶r�}�±�6p¡¶vy|~��¹`rs�,{b|e±&��¡¬�zÞwè�è:<>�§±*� m�}z��r���x���� �pqu{9t­ºS©jvy|~t�ª�{Q}�vyqu�Q{��zrsxz�p� �¨{b�}­{�qs�`r�|~�¾}zn¦vy��}��o}z�`������xz�zr�t³�Q��|~���§v�qsr�ª�������r�����|pnp��{b�pq����,}^´�±¥� ݺD ��M÷� ä�ã4àWæ��/D�Þ-ä�ó ÷ D�ð�ð�ó·æ4ç ÷ì9ü ±�çQç�0`Ä�ç98bç9¡

�¤ru�b�¦vy�˱ À ¡��zÞwè�èQÜ1�§± D�çWç^éMà5ä1ç^òAä=ÞVá � ãÌä>ë$æ4óÍæ4ã"ò�ß,è � éMû2â^àWæ4çUä1ó�D�óÂü�ß:àWæ4ã ù ûAê ±�¹ ¯ m 2 3¬�prsq vb�p��qsnp�pr�v�¡�&¹p© �zÞyèQè�Ü��§±V� �&v��z²¥��qsqB¹`�p�p��{b�pxzrs|~�0©�rs�p�§vy�z��¡¤m*ߥvyx§v�q�{Q���~��{b��¹`���p��{���xzr�|~��}×� â ��q��wv�}��(Þwç��W´�¡

¯ |p�¨{b�z��v�xzr�{�|)¼ À rut­ºQbZ���wv�q����`±�m&à«m ¢ ��t­�p|~{bqu{Q�b�`±�<M8�8T�&vy�z²¥��qsq ±&áoru��t�{bx�±!ã¶Á�{b|e±!ã£ä�Þ�Þå â mC±�½&´7±�àT� ��vyr�q̼¬�~}�q��!vb�wv�x�¡ t�{�¡ �pº¦¡

©jvy|~t�ª�{�}�±eß�¡��zÞwè�ÝQå1�§±3� mo|�r�x����§vyxzru{b|�����xz�~{9�¸�¨{b�!xz�~��}�{�qs�pxzr�{�|�{��Zxz�~����r�����|`��v�qs�~�Cnp��{b�pq�����{b�qsrs|~�wv��o�9rSR�������|�xzr�vyq�v�|~�Irs|�x����b�§v�qO{bn����§vyx�{���}^´�±��M÷ í�â^ê ÷�� ä1ã ÷ � éMà ÷p� ãÌä:Þ-á ÷ ö�ü ±�çQÝ�Ý�Äpç�0Qç�¡

©O���~{b�~t��~± â ¡��zÞyèQè�Ü��§± â ��}�x§v��zxzrs|~�Bvy|�mo�z|~{bq��9rp�����9�~t�xzr�{�|e±13¬����n��zr�|`x2�ߥ¹>��3cÝ�è�Þ �Ïå=<`èQÜ�±�mo����{b|p|~�À vyxzru{b|¦vyqj©jv���{b�§v�x�{b�z��±9mo���Q{�|p|��Q± ¯ ©0ÜQå:<9[Qè�¡

©O���~{b�~t��~± â ¡��zÞyèQèQè1�§±9ã¶|�xz���¶t�{�|`�Q�����Q��|~t��¶{b�evy|�r���npqsr�t�rsxzqs������}zx§vy�zx�����m��z|�{�q��9r~����xz�~{9��± ¢ ��t����|�rutwvyq â ��n�{��zx�¹9m À á!èQè`�§Þ�8�ÝQÜ=6p±Z¹pvy|~�9r�v À v�xzr�{b|¦v�q!©jvy��{b�§vyx�{��zr���}�±*2¸¹ Þ�ÞQÞyå9±�3�¡ 㾡�bc{�ÁÝ�0Qå�å9±9moqs�p�~�Q�~�����Q�~��± À 2 0>89Þ�0�Ý`�§ÞQÞQÞyå`¡

©O���~{b�~t��~± â ¡�v�|~��¹9{b����|~}���|e±9áB¡��zÞwèQè�Ü��§±�� á&�$�¦vyxzr�{�|Cx���t��p|pr��Q�~��}Z²¶rsxz�prs|Ivy|�r���npqsr�t�rsxzqs�7����}zx§vy�zx����mo�z|~{bq��9rOr�x����§vyxzru{b|P´s±p� ݺD �M÷� ä�ã4àWæ��)D�ÞVä1ó ÷ D�ð�ð�ó·æ4ç ÷ ì � ±+8=0Qè�Ä�0QçpÞ�¡2������z�������Q��|O±~´C¡pvy|~�˹�n���|~t���±�mC¡��zÞwè�è>8:�§±£� ¯ ��npqsr�t�rsxzqs������}zx§v��zx����Im��z|�{�q��9rjv�|���np�p�zr�#~twvyxzru{b|J�¨{b�xz����}z�prs��x5� rs|`�����zx�xz�§vy|~}z�¨{b�z��v�xzr�{b|P´�±� ä1ã ù�÷P ß:û�ð ÷ þ~þ ±¦Ü�Ü>8~Ä�Ü�0�è9¡2�{b����vy|e± â ¡��zÞyèQèQÜ1�§±9�são|C����}zx§v��zxzrs|~�!xz����mo�z|~{bqu��rp����xz��{9�B�¨{b�µq�vy�����c|�{�|~}���������xz�zruto��ru�Q��|`��v�qs�~�n���{��pq�����}~´�±� ä�ã ù�÷ P ß:û�ð ÷ þ¦ü ±®Þyç�Þ�[`ĦÞyç�[Qå�¡

À {����5�­ão��r���±+b�¡s±+3¬v��zq���xzx�±)b�¡s±�à¬�zr�t�}­}�{�|O± ¢ ¡�vy|~��6Q��|�}���|e±+3�¡l�zÞwè10>8:�§±p� �o{w²�x�{�rs��npq�������|`x!xz�~�}�n���t�xz�§v�qOxz�§vy|~}z�¨{b�z��v�xzr�{b|P´�±� ä1ã ù�÷[P ß:ûIð ÷ ö û�±�ÜQÜ�[`Ä�Ü98:[9¡

3µvyr�����±Qß�¡s±:3µvy�zqu��xzx�±�b�¡bvy|~�!��vy|7�p���)Kc{b��}zx�±b�¾¡1�zÞwè�èQÝ1�§±�� monpnp��{�Á�rs�Ivyx��Z}�{bq��pxzr�{b|~}�vy|~����ru�Q��|`��v�qs�~���{b�p|~�p}£����{��G´��z��q�{w�J}��p�~}zn¦vbt���}~´�±�� éMû ÷ B æ¶Þ ÷ D�ó ü ÷ D�ð�ð�ó ÷ î ±®ÞQÞwÝ`ĦÞ`[�[9¡

3µvy�zq���xzx�±Pb�¡o�zÞyè�0�å��§±�ø�ù â êUò=ûAû{â�ã"àWæ4ç â^æSü�â^ÞFå�ä�óÍé�â�ðVà5ß9ë^ó â^û ±-3µ����|`xzr�t����&vyq�q�}­���zr���}or�|Tt�{b��np�px§v`�xzr�{b|¦v�qO��vyxz�~����v�xzr�t�}�±�3µ����|`xzr�t��$�Ï�³vyqsq ¡

3¬ru����t���±OmB¡�áB¡T�zÞwè10�Þ`�§± D�çWß:éMê ã4æ4ç^ê R<DLÞZæ¶ÞVã4à5ß�á=é�ç�ã"æ4ß=Þã�ßÉæ4ã4êHð ù ò:êWæ4çUä1ó�ð-àUæ¶Þ-ç^æÓð+óúâ^êBä:Þ-á�äUð1ð+óÍæ4çUä�Rã"æ4ß=Þ�ê ±¥2¸vbtB/o�§v�² �¶rsqsq ± À ��²CÒc{b�zº¦¡â �p�~��±`mB¡��zÞwè10:<>�§±�� â v�xzr�{b|¦v�q�´&�z�`qu{w�C}����Q�~��|~t��¥����xz�~{9�p}¬�¨{��Z��ru�Q��|`��v�qs�~�¶t�{b��np�px§vyxzru{b|P´s± B æ¶ÞVâWä:àDIóÂü ÷ D�ð�ð+ó ÷ ü û�±+[QèpÞWÄ9<�åQÝ�¡â �p�~��±�mC¡��zÞwè�è:<>�§±°� â vyxzru{b|¦vyqe´��z��q�{w�Jvyqu�Q{b�zr�xz�p��}c�¨{b�¥|~{b|~}z�`������xz�zrut!��r�����|`��v�qs�~�&np��{b�pq����,}�± ¯z¯ ¼2¸vyxz�zr ÁIn¦vyr���}~´s± B æ¶ÞVâWä:à<DIóÂü ÷ D�ð�ð�ó ÷ ì � �� ì ��û�±eç10�[`Ä�çQè�Ü9¡

ç�ç

Page 25: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

â �p�~��±�mB¡P�zÞyèQè�01�§±�� â vyxzr�{�|¦vyq�´&�z�`q�{w�-¼¶vCn��§v�t�xzr�tyvyq�vyqu�Q{��zrsxz�p�?�¨{b�oq vy������}zn¦vy��}��&|~{b|~}z�`���,��xz�zr�t��vyxz�zrsÁJn���|~t�rsqu}^´�±�� ݺD �M÷q� ç^æ ÷�P ß=û�ðVé�ã ÷ ì �)� Ý��§±®ÞwÝ1[QÝ`ĦÞwÝ�Ý�ÞQ¡

¹pvQv���±�Ò�¡��zÞwèQè�ç��§±!� éMû2â^àWæ4çUä1óMû2â�ã ù ß�á:ê-è$ß=à�óúä:à�ü�â�â^æSü�â^ÞFå�ä�óÍé�âlðVà,ß>ë^óúâ^ûAê ±bmoq���{b�zrsxz�p��}µvy|~�7mo��t­�prS�x���t�xz�p����}!�¨{b�!m³����v�|~t����T¹9t�r���|�xzr�#~tJß«{���np�pxzrs|~�p±q2¸vy|~t­�~��}zx����¾½¤|prs������}zrsx �3µ����}�}�±q2¸vy|~t������}zx�����±�½&´7¡

¹pvy|~�9�`�e± â ¡bv�|~�¬3µr�}zx�����±`´7¡1�zÞwè98bå1�§±�� mE��v��zr�v�xzr�{b|¦v�qQnp�zrs|~t�r�n�qu�c�¨{b�Oqsr�|~�wvy��t�{���npqu���(#���q���n���{��pq�����}rs|�t�{b|`xzr�|`�p��� ����t��¦vy|pr�t�}^´�± ÝUÞ-ã ÷��M÷"� Þ�ü ÷ � ç^æ ÷pû�±�n�ne¡�è10Qè�ÄpèQè�è9¡

¹`rs��{�|e±�bo¡�± % �e±n6�¡s±*Ñ�ru��|pº`r���²¶r�t�ª�±&ãB¡«v�|�� 3µvy�pq ±¶áB¡I�zÞwè10QÜ1�§±É� àµ��vyqs�¦v�xzr�{b|S{��&��� ²^vy|~�(��� n#�|pr�x��,��q�������|`x�����xz�~{9�p}��¨{��!xz�����9�`|¦vy��r�t�����}zn�{�|~}­��{b�¥}§vyxz�p�§vyx����Dn�{b��{b�~}!�����9r v,�~}zrs|~�{b|~�$�¼��r�����|~}zr�{�|�v�q���{9�p��qu}~´s± ÝUÞ-ã ÷��M÷#� éMû{â^à ÷ D�Þ-ä�ó ÷� â�ã ù ß�á:ê%$�âWß:û{âUç ù�÷ ì9� ±�n�ne¡�<�Ü�Þ$Ä9<>0�ç9¡

¹`q���rÂ�Ïn���|e±�/B¡�v�|~�C��vy|��p���§Kc{b��}zx�±`��¡��zÞwè�èQÜ1�§±�� mQ6`vbt�{b�pr���á!v��`r��p}�{b|Brsx����§v�xzr�{�|�����xz�~{9�C�¨{��¬qsr�|~�wvy���r�����|`��vyq��~�³np��{���qu����}~´s±p� ݺD �M÷� ä1ã"àWæ��)D�ÞVä1ó ÷ D�ð1ð+óÍæ4ç ÷ ì � ±�<�å�Þ$Ä9<`ç�Ý9¡

¹9{b����|~}���|e±,áB¡.�zÞyèQè�ç��§±?� ¯ ��npq�r�t�r�x4vynpnpqsrutwvyxzru{b|Í{b��n�{�qs�`|~{���r�v�q #�q�x�����}�r�|vû>�¼}zx���n mo�z|~{�q��9r����xz�~{9��´s±p� ݺD �M÷" ä�ã4àUæ���DLÞ-ä�ó ÷ D�ð�ð�ó·æ4ç ÷ ì`õ ±V[QÝ98^Ä�[�0QÝ�¡

ç1[

Page 26: Home - STFC Computational Mathematics Groupé # ê,ëµ &ê ë¬ ì í Û ÖwÇ Æ ÈµÒeÕ Ö Ù Æ Û î iBïjÃ É ðZà ÕbÖwÇ Ø~×ZÖwÇ Ø ÛjÉ ñ¼Æ Ç ò Ø ÖbÙ Æ

� ���! ��¥ëZ��#�ê #2¤ ¡ )Z�(¥ëT¡Bê��2�¢^¡B) % ¡�ëµ)c��Ï ë¬# ë¬)c�w ¸�^¡���#�ê7¡£¢¤�#�)Z%

ßc{b|~}zr��p���Cxz�~�]»�³ Þ!Õ¥»��Ivyxz�zr Á Ë � � ¤ ²¤r�xz�4xz�~�I�pnpn����.»�ÕZ»0}��p�p��v�xz�zr Á �����z��rsxzr vy|e¡ ¢ �~��¨{bq�q�{w²¶rs|~��vyqu�Q{��zrsxz�p�?t����wv�x���}ovB��|pr�x§vy�z�J��vyxz�zrsÁ,><}z��t��Jxz��v�xË � E ¤ � ! > å

å Þ�# h Ë � � ¤ >r�}¥xz�zr��9r vb�Q{�|¦vyq ¡f�×�ø�Æ~Ç�ÙsÖ�ï�ò f ¸¨ìÞ�¡�©O��x|>v��gç9¡ ·~{b��f ��»�x�{�çC}zx���nY?BÞ!�p{ç�¡�ÞQ¡�©O��x%&(eC����v���{b�~}����~{bq��p�������$�~��t�xzr�{b|»}z��t��Exz�¦v�x2d�he © E Ë � eº¤ &(eCr�}�vI��{w² �Q��t�x�{b�B²¤r�xz�f�?ýÞ!ª�����{�}£r�|�xz�~�!����{b|`x�¡ç�¡ ç�¡¦·~{b�z�

Ë � e C�E ¤ � ! &(e åå Þ # h Ë � eº¤ &(e

ç�¡·[ ·~{b�z�4>k�X>'&(eão|�r�x����§vyxzru{b|£f�xz�~�nfc³�Þ¥}zxZ��{w² {�� Ë � eº¤ r�}Z�������~t�����x�{¾v&��{w²�²¤r�xz�£f�?�Þ¶ª�����{Q}µrs|�xz�~�¤����{�|`x�¡¢ ������{b�~}����~{bq��p���¥xz�§v�|�}z�¨{��z��vyxzru{b|Jr�}�vyq�}�{�vynpnpqsr����I{�|Ixz�~�!����{b|`x�}�{¾xz�¦vyx¶xz�~�If(³ Þ!}zx¤t�{bqs�p��|

r�}¥�����9�~t����Jx�{�vCt�{bq��p��|�����t�x�{��¤²¶rsxz��f�?SÞ!ª�����{Q}¥rs|Ëxz���!����{�|`x�¡

ç=<