42dcbde0e602cba7adb911a2072df3887afd61f4,hicexplorer/utilities.py,,exp_obs_matrix_lieberman,#Any#Any#Any#,215

Before Change

                pSubmatrix.data[i] = pSubmatrix.data[i] / expected_interactions_in_distance_[distance[i]]
        except Exception:
            log.debug("pSubmatrix.data[i]: {}".format(pSubmatrix.data[i]))
            log.debug("distance[i]: {}".format(distance[i]))
            log.debug("expected_interactions_in_distance_[distance[i]]: {} ".format(expected_interactions_in_distance_[distance[i]]))
            exit(1)

After Change

    
    // log.info("expected_interactions_in_distance_[distance[0]]: {}".format(type(expected_interactions_in_distance_[distance[0]])))
    
    if len(pSubmatrix.data) > 0:
        data_type = type(pSubmatrix.data[0])
        // from copy import deepcopy 
        // data = deepcopy(pSubmatrix.data)
        
        expected = expected_interactions_in_distance_[distance]
        // log.info("len(expected); {}".format(len(expected)))
        // log.info("len(data): {}".format(len(pSubmatrix.data)))
        pSubmatrix.data = pSubmatrix.data.astype(np.float32)
        pSubmatrix.data /= expected
        // data2 = deepcopy(data)
        pSubmatrix.data = convertInfsToZeros_ArrayFloat(pSubmatrix.data).astype(data_type)


    // for i in range(len(pSubmatrix.data)):
    //     try:
    //         if expected_interactions_in_distance_[distance[i]] == 0:
    //             pSubmatrix.data[i] = 0.0
    //         else:
    //             pSubmatrix.data[i] = pSubmatrix.data[i] / expected_interactions_in_distance_[distance[i]]
    //     except Exception:
    //         log.debug("pSubmatrix.data[i]: {}".format(pSubmatrix.data[i]))
    //         log.debug("distance[i]: {}".format(distance[i]))
    //         log.debug("expected_interactions_in_distance_[distance[i]]: {} ".format(expected_interactions_in_distance_[distance[i]]))
    //         exit(1)
    // if len(pSubmatrix.data) > 0:

    //     try:
    //         nt.assert_equal(data, pSubmatrix.data)
    //     except:
    //         foo = data - pSubmatrix.data
    //         mu = np.nonzero(foo)
    //         print(mu)
    //         log.info("data: {}".format(data[mu]))
    //         log.info("data: {}".format(data2[mu]))
            
    //         // data2
    //         log.info("pSubmatrix.data: {}".format(pSubmatrix.data[mu]))
    
    return pSubmatrix


def toString(s):

Instances