От циклов for к реализации numpy dot

Question

От циклов for к реализации numpy dot

1

N=100 

reliab=zeros((N,N))


for i in range(N):    

    for j in range(N):
        if random() < 0.6:
            reliab[i,j] = 1
        else:
            reliab[i,j] = 0

Как и в моем коде, это заполнение матрицы вызывает огромное количество раз, это для циклов должно быть изменено с помощью точечного продукта... но я не знаю, как это сделать... Кто-нибудь хочет помочь мне?

Aurelio 01 июнь 2011, в 00:02

Источник

Теги:

python

numpy

loops

2 ответа

1

>>> import numpy as np
>>> reliab = np.random.random((N,N))
>>> reliab = reliab < 0.6
>>> reliab.dtype = np.int8
>>> reliab
array([[0, 0, 0, ..., 1, 0, 1],
       [0, 1, 1, ..., 0, 1, 0],
       [1, 1, 0, ..., 1, 0, 1],
       ..., 
       [0, 1, 0, ..., 1, 1, 1],
       [0, 0, 1, ..., 1, 1, 1],
       [0, 1, 0, ..., 0, 0, 0]], dtype=int8)

dr jimbob 31 май 2011, в 20:15

0

Приятно. Вероятно, вам следует масштабировать надежность, прежде чем менять ее тип. ОП хочет присвоить 1 с вероятностью 0,6.
Ruggiero Spearman 31 май 2011, в 21:23
0

@ Amaç Herdağdelen: np.random.random находится на интервале [0,1) , поэтому масштабирование не требуется. (Или, суммируя мой массив, я получаю 6009 из 10000, что почти на 60%). (Да, вы опередили меня с ответом, и они похожи, но я решил, что все равно напишу ответ).
dr jimbob 31 май 2011, в 21:29
0

Извините, я видел более раннюю версию вашего ответа, в которой не было сравнения <0,6. Там нет ничего плохого прямо сейчас.
Ruggiero Spearman 31 май 2011, в 21:33
0

Замечательный! Большое спасибо, ребята!
Aurelio 01 июнь 2011, в 07:06

Показать ещё 2 комментария

Ещё вопросы

Приятно. Вероятно, вам следует масштабировать надежность, прежде чем менять ее тип. ОП хочет присвоить 1 с вероятностью 0,6.
@ Amaç Herdağdelen: np.random.random находится на интервале [0,1) , поэтому масштабирование не требуется. (Или, суммируя мой массив, я получаю 6009 из 10000, что почти на 60%). (Да, вы опередили меня с ответом, и они похожи, но я решил, что все равно напишу ответ).
Извините, я видел более раннюю версию вашего ответа, в которой не было сравнения <0,6. Там нет ничего плохого прямо сейчас.
Замечательный! Большое спасибо, ребята!

Ruggiero Spearman · Accepted Answer · 2011-05-31T19-00-00.000Z

Я не уверен, что полностью понимаю вашу проблему, но следующая строка должна делать то же самое, что и ваш код:

reliab = numpy.int32(numpy.random.rand(N,N) < 0.6)