← Сtrl

bismigalis · Ноя. 9, 2013 19:19:33

keekdown
Я думал что Numpy and Scipy обеспечивают оптимизацию таких вещей

обеспечивают, если ты используешь методы модуля numpy

keekdown
Ну как про С реализацию функций я думал.

уже все реализовано, бери, пользуйся

EDIT: или используй pypy он развернет цикл в нативный код

Отредактировано bismigalis (Ноя. 9, 2013 19:22:08)

keekdown · Ноя. 9, 2013 19:23:21

bismigalis

такой метод наврятли присутствует в библиотеке numpy scipy

keekdown · Ноя. 9, 2013 19:32:34

Я уже попробовал pypy.Да быстрее,но не настолько.Знаю что на С матрица считается минуты за 2,я писал там подобную реализацию.Возможно есть выход пользоваться слайсингом..Если не по сабжу,то я узнал что оказывается есть модуль Messege Passing Interface для параллельного программирования.Видно значит можно использовать этот язык в вычислениях таких.

bismigalis · Ноя. 9, 2013 19:38:06

keekdown
Знаю что на С матрица считается минуты за 2

ого, так много, а ты еще на питоне хочешь
тут наверное надо GPU задействовать :)

keekdown · Ноя. 9, 2013 19:41:58

bismigalis

Ладно,спасибо за ответы.Был рад пообщаться.

bismigalis · Ноя. 9, 2013 20:34:54

keekdown
Ладно,спасибо за ответы.

Я вообще нуб в вычислительных расчетах. Про numpy только слышал краем уха. Щас стало интересно полез в документацию, и вижу что ты нихрена не понял как она работает.

keekdown
Был рад пообщаться.

Заходи ещё.

sergeek · Ноя. 9, 2013 22:52:16

from math import pi
import numpy as np
 
def calcf(n):
    const = 2**0.5 * pi 
    
    make_x = np.vectorize(lambda i, n : 1.0 / n * i)
    x = make_x(np.arange(n+1), n)
    
    one = np.sin(pi*x)*const 
 
    return np.vstack(e*one for e in one)
 
print(calcf(10000))

keekdown · Ноя. 10, 2013 06:13:58

sergeek

Очень интересная реализация!Спасибо!Тщательнее почитаю документацию

sergeek · Ноя. 10, 2013 12:12:51

vectorize тут лишний даже

from math import pi
import numpy as np
  
def calcf(n):
    const = 2**0.5 * pi
     
    x = 1.0 / n * np.arange(n+1)
    one = np.sin(pi*x)*const 
    return np.vstack(e*one for e in one)

keekdown · Ноя. 11, 2013 08:46:40

sergeek
vectorize тут лишний даже

Только не понятно почему

const = 2**0.5*pi

.У меня было

const = 2* pi* pi

.Это немного разные вещи.

Python-сообщество

Уведомления

#1 Ноя. 9, 2013 19:19:33

Вычислительные рассчеты

#2 Ноя. 9, 2013 19:23:21

Вычислительные рассчеты

#3 Ноя. 9, 2013 19:32:34

Вычислительные рассчеты

#4 Ноя. 9, 2013 19:38:06

Вычислительные рассчеты

#5 Ноя. 9, 2013 19:41:58

Вычислительные рассчеты

#6 Ноя. 9, 2013 20:34:54

Вычислительные рассчеты

#7 Ноя. 9, 2013 22:52:16

Вычислительные рассчеты

#8 Ноя. 10, 2013 06:13:58

Вычислительные рассчеты

#9 Ноя. 10, 2013 12:12:51

Вычислительные рассчеты

#10 Ноя. 11, 2013 08:46:40

Вычислительные рассчеты

Board footer