[Wien] Inaccurate Tetra integration

Laurence Marks L-marks at northwestern.edu
Tue Nov 1 14:01:56 CET 2016


So...no fix (sigh)?

No integration should go so badly wrong....

On Tue, Nov 1, 2016 at 7:58 AM, Peter Blaha
<pblaha at theochem.tuwien.ac.at> wrote:
> I would not use TETRA for a 2D BZ (a mesh with "1"). The 3D tetrahedra
> collapse to triangles and this leads to degeneracy and problems.
>
> Am 01.11.2016 um 13:50 schrieb Laurence Marks:
>> I am continually running in to unacceptably large errors in the TETRA
>> integration, e.g.
>>
>> RESULT OF INTEGRATION: 2856.04291; SHOULD BE: 2856.00000
>>
>> It is a large cell of  63.8 x  7.98 x 127.43, pmm with a shifted 2x8x1 mesh.
>>
>> I am using TETRA    0.0001 which for a bit seemed to help with MSR1a,
>> but the problem has reappeared with MSR1 (-it noHinv). I will shortly
>> test using the standard Tetrahedron method (100). I don't want to
>> shift to TEMPS unless I really have to (for convergence reasons).
>>
>> IMHO this is a "bug" somewhere. I would appreciate any comments on
>> work arounds as I don't want to go bug hunting today.
>>
>
> --
> --------------------------------------------------------------------------
> Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
> Phone: +43-1-58801-165300             FAX: +43-1-58801-165982
> Email: blaha at theochem.tuwien.ac.at    WIEN2k: https://urldefense.proofpoint.com/v2/url?u=http-3A__www.wien2k.at&d=CwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=AxIhIc1H_Nkr9Vq4Z5UW1vA-yIEho9cw-7WxHqFBlcs&s=l8TSHPCPswUPZ-oJlQC8Ov7Ih97gqkgW5gHl-HcgHpo&e=
> WWW:   https://urldefense.proofpoint.com/v2/url?u=http-3A__www.imc.tuwien.ac.at_staff_tc-5Fgroup-5Fe.php&d=CwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=AxIhIc1H_Nkr9Vq4Z5UW1vA-yIEho9cw-7WxHqFBlcs&s=6LxEJsP2lKwjPS9C6rT5l137ZBBwESZvmB4f7tYTqOI&e=
> --------------------------------------------------------------------------
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=CwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=AxIhIc1H_Nkr9Vq4Z5UW1vA-yIEho9cw-7WxHqFBlcs&s=TZxLHTLtp7gh7-KA0GDw-lFOnas-elj395smb_PGxi4&e=
> SEARCH the MAILING-LIST at:  https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=CwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=AxIhIc1H_Nkr9Vq4Z5UW1vA-yIEho9cw-7WxHqFBlcs&s=6JnSZn3yZNEV2kiqY0aztUHK3ahpEPEBG-NkVCVnRQE&e=



-- 
Professor Laurence Marks
"Research is to see what everybody else has seen, and to think what
nobody else has thought", Albert Szent-Gyorgi
www.numis.northwestern.edu ; Corrosion in 4D: MURI4D.numis.northwestern.edu
Partner of the CFW 100% program for gender equity, www.cfw.org/100-percent
Co-Editor, Acta Cryst A


More information about the Wien mailing list