Skip to content

Contact matrix element divided by zero #14

@yasumasak

Description

@yasumasak

Hi,
I found that some contacts were missing after running HIFI.
It was caused when gs==0 at line 405 of HIFI.cpp, producing NaN for an element of matrix T.
I just solved this by modifying line 405 as below
if (s!=0 && gs!=0) { T.set(i,j,s/gs);}
Is it a reasonable solution?

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions