LINUX.ORG.RU

Сообщения math-beginer

 

nm упал

Форум — Admin

После обновления Network Manager не стартует. Помогите его заставить работать.

$ nm-applet 

** (nm-applet:5147): WARNING **: Icon nm-secure-lock missing: Значок «nm-secure-lock» не присутствует в теме



$ sudo cat /var/log/NetworkManager

Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) starting connection 'Auto dlink'
Aug  5 12:55:09 hp4525s NetworkManager: <info>  (eth1): device state change: 3 -> 4 (reason 0)
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 1 of 5 (Device Prepare) scheduled...
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 1 of 5 (Device Prepare) started...
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 2 of 5 (Device Configure) scheduled...
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 1 of 5 (Device Prepare) complete.
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 2 of 5 (Device Configure) starting...
Aug  5 12:55:09 hp4525s NetworkManager: <info>  (eth1): device state change: 4 -> 5 (reason 0)
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1/wireless): access point 'Auto dlink' has security, but secrets are required.
Aug  5 12:55:09 hp4525s NetworkManager: <info>  (eth1): device state change: 5 -> 6 (reason 0)
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 2 of 5 (Device Configure) complete.
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 1 of 5 (Device Prepare) scheduled...
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 1 of 5 (Device Prepare) started...
Aug  5 12:55:09 hp4525s NetworkManager: <info>  (eth1): device state change: 6 -> 4 (reason 0)
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 2 of 5 (Device Configure) scheduled...
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 1 of 5 (Device Prepare) complete.
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 2 of 5 (Device Configure) starting...
Aug  5 12:55:09 hp4525s NetworkManager: <info>  (eth1): device state change: 4 -> 5 (reason 0)
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1/wireless): connection 'Auto dlink' has security, and secrets exist.  No new secrets needed.
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Config: added 'ssid' value 'dlink'
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Config: added 'scan_ssid' value '1'
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Config: added 'key_mgmt' value 'WPA-PSK'
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Config: added 'psk' value '<omitted>'
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Activation (eth1) Stage 2 of 5 (Device Configure) complete.
Aug  5 12:55:09 hp4525s NetworkManager: <info>  Config: set interface ap_scan to 1
Aug  5 12:55:09 hp4525s NetworkManager: <info>  (eth1): supplicant connection state:  scanning -> disconnected
Aug  5 12:55:09 hp4525s NetworkManager: <info>  (eth1): supplicant connection state:  disconnected -> scanning
Aug  5 12:55:13 hp4525s NetworkManager: <info>  (eth1): device state change: 5 -> 3 (reason 38)
Aug  5 12:55:13 hp4525s NetworkManager: <info>  (eth1): deactivating device (reason: 38).

math-beginer
()

python заменить цикл

Форум — Development

Есть массив x_old, и нужно заменить значения этого массива на значения массива x_new. Но, из x_new можно для замены брать только те значения, которые удовлетворяют определенному условию. x_old и x_new хранят координаты точек на плоскости. Условие следующее: из x_new нужно брать только те точки, которые лежат вне эллипса с параметрами xc,yc - центр, а,b - полуоси.

Задачу решает код с циклом: Код:

#!/usr/bin/python
from numpy import *
m=3
xc=0.
yc=0.
a=2.
b=1.
eps=0.1
x_old=ones((m,2),dtype='d')
x_new=ndarray((m,2),dtype='d')
x_new[0][0]=2.;x_new[0][1]=0.
x_new[1][0]=3.;x_new[0][1]=1.
x_new[2][0]=1.;x_new[0][1]=0.
for i in range(m):
    if (x_new[i][0]-xc)**2/a/a+(x_new[i][1]-yc)**2/b/b>(1.+eps):
        x_old[i]=x_new[i]
print "x_new=",x_new
print "x_old=",x_old
от цикла хотелось бы избавиться (устраняю медленные участки). Можно ли решить задачу без цикла средствами numpy?

math-beginer
()

установить matplotlib

Форум — Admin

что это за модуль? и откуда его взять?

>>> import pylab
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "/usr/lib/python2.6/site-packages/pylab.py", line 1, in <module>
    from matplotlib.pylab import *
  File "/usr/lib/python2.6/site-packages/matplotlib/pylab.py", line 263, in <module>
    from matplotlib.pyplot import *
  File "/usr/lib/python2.6/site-packages/matplotlib/pyplot.py", line 95, in <module>
    new_figure_manager, draw_if_interactive, show = pylab_setup()
  File "/usr/lib/python2.6/site-packages/matplotlib/backends/__init__.py", line 25, in pylab_setup
    globals(),locals(),[backend_name])
ImportError: No module named backend_tkagg

P.S. для SLED11 SP1 подключил репозитори opensuse 11.2 oss появился python-matplotlib, но при попытке установки, zypper начал ругаться на несоввместимость, тогда я подключил репозиторий opensuse 11.3 и python-matplotlib установился без проблем [с предупреждением о дополнительном клиентском соглашении, конечно же], но были обновления 10-20 пакетов. По гугловскомупоиску много чего делал и freetype2-devel не нащел, python-matplotlib-wx tk пробовал ставить, не помогает, ipython -pylab выдает теже ошибки.

math-beginer
()

репозиторий для SLED

Форум — Desktop

Имеется на буке SUSE LINUX Enterprise Desktop 11. Из какого репозитория можно доустановить пакеты, например, gvim?

P.S.

пробовал подключать

Версия: 11.4http://download.opensuse.org/distribution/11.4/repo/oss/ Версия: 11.3http://download.opensuse.org/distribution/11.3/repo/oss/ Версия: 11.2http://download.opensuse.org/distribution/11.2/repo/oss/ Версия: 11.1http://download.opensuse.org/distribution/11.1/repo/oss/

и аналогичные uptade по очереди. В результате zupper сначала сносит, наверное, все, а затем ничего не происходит :-)

math-beginer
()

python что за ограничение на 1024 строки

Форум — Development

Имеется массив из большого числа строк, по 9-ть элементов в каждой строке. Его нужно записать в файл. Стандартные операторы записывают только 1024 строки, беря начало и конец массива.

Как можно эту проблему устранить?

#!/usr/bin/python
import sys
import scitools.filetable as ftable

...

b0n=b0[:,0:9]
fb0n=open(file_name+".raw",'w')
ftable.write(fb0n,b0n)
fb0n.close()
math-beginer
()

автоматом удалить фон на скане

Форум — Desktop

Отсканил кое-какие документы, и сконфертил в pdf командой

convert *pnm -compress JPEG -quality 50 res.pdf
после печати виден «сильный серый фон».

Можно ли этот фон удалить автоматически, путем обращения к какой-ли бо утилите, или может у convert есть какой-то ключик? Рисунки цветные , печать черно-белая. Запускать гимп для каждого листа не охота.

math-beginer
()

ускорение dot,cross,norm

Форум — Development

Продолжаю ускорять программы численных расчетов. У меня в модуле для ускорения часто встречается скалярное, векторное произведение вектора с тремя координатами и вычисление его нормы. Скорее всего грамотно последнее сделать вызовом функций BLAS, и действиетльно существует cblas_ddot и cblas_dnrm2 (Евклидова устроит). Какие есть проблемы:

1. cblas_dcross для векторного произведения не вижу

2. для маленьких векторов вот тут http://www.intuit.ru/department/se/intelmkl/2/ чел не рекомендует blas

Судя по новостям, ЛОР математику любит, поэтому прошу подсказать библиотеки, которые будут эффективно (быстро) работать из под Питона и выполнять простые операции, пока нужно dot, cross, norm.

P.S. по сути продолжаем этот разговор пытаемся ускорять с cython (комментарий)

P.S. Самописный вариант отбросим из рассмотрения

math-beginer
()

пытаемся ускорять с cython

Форум — Development

следуя вот этой статье http://sage.math.washington.edu/home/dagss/numerical-cython-preprint.pdf попытался ускорить тако код:

#!/usr/bin/python
from numpy import *
def matrix(bRow, bCol, out=None):
    for i in range(len(bRow)):
        for j in range(len(bCol)):
            tmp=0.0
            for k in range(len(bRow[i])):
                tmp+=0.000001*(linalg.norm(bRow[i][k]))
            for k in range(len(bCol[j])):
                tmp+=0.000001*(linalg.norm(bCol[j][k]))
            out[i][j]=tmp

в результате получи следующий код:

#!/usr/bin/python
import numpy as np
cimport numpy as np
cimport cython
ctypedef np.float64_t dtype_t
#ctypedef double dtype_t
@cython.boundscheck(False) 
@cython.wraparound(False)
def matrix(np.ndarray[dtype_t, ndim=3] bRow,np.ndarray[dtype_t, ndim=3] bCol,np.ndarray[dtype_t, ndim=2] out=None):
#    cdef long i,j,k
    cdef Py_ssize_t i,j,k
    cdef dtype_t tmp
    for i in range(bRow.shape[0]):
        for j in range(bCol.shape[0]):
            tmp=0.0
            for k in range(bRow.shape[1]):
                tmp+=0.000001*(np.linalg.norm(bRow[i][k]))
            for k in range(bCol.shape[1]):
                tmp+=0.000001*(np.linalg.norm(bCol[j][k]))
            out[i][j]=tmp

После замеров времени в среднем оказалось, что первый и второй вариант работают 2 с половиной минуты. Т.е ни о каком ускорении речи быть и не может. А в статье после типизации обещается ускорение в 130-150 раз, а после отмены проверок в 620 раз.

Что я сделал не так???

P.S. сборка so модуля выполнялась так

$ cython matrixx.pyx 
$ gcc -c -fPIC -I/usr/include/python2.7/ matrixx.c 
$ gcc -shared matrixx.o -o matrixx.so

math-beginer
()

cython ругается на pow

Форум — Development

пытаюсь ускорить скрипты. что можно сделать?

cython my_module.py вываливает ошибки однотипные, ругается на pow

Error converting Pyrex file to C:
------------------------------------------------------------
...
    
    Wt=array([1.,.0,.0])-1.5/pow(linalg.norm(x),5.)*x
                              ^
------------------------------------------------------------

 Call with wrong number of arguments (expected 3, got 2)

math-beginer
()

c4py matrix_maker

Форум — Development

Помогите разобраться с передачей матрицы из функции.

$ cat matrix.c 
#include <python2.7/Python.h>
#include <stdlib.h>
static PyObject *matrix_maker_cmd(PyObject *self,PyObject *val)
{
    double nRow,nCol;
    double **a;
    int i,j;
    if(!PyArg_ParseTuple(val,"ii",&nRow,&nCol))
        return NULL;
    a = malloc(nRow*nCol*6);
    for(i=0; i<nRow; ++i){
    for(j=0; j<nCol; ++j){    
        a[i][j]=cos(i+j);
    }
    }
    return Py_BuildValue("d",a);
}
static PyMethodDef threeMethods[] = {
    {"matrix_maker",matrix_maker_cmd,METH_VARARGS,"calc matrix elem"},
    {NULL,NULL,0,NULL}
};
void initthree(void)
{
    Py_InitModule("matrix",threeMethods);
}

$ cat test_matrix.py
#!/usr/bin/python
import matrix
print matrix_maker(5,4)

$ ./test_matrix.py
Traceback (most recent call last):
  File "./test_matrix.py", line 2, in <module>
    import matrix
ImportError: dynamic module does not define init function (initmatrix)
math-beginer
()

где взять Python.h

Форум — Desktop

нужно вынести код с циклом в цикле в функцию спп. В инете часто встречается пример

$ cat c4py.c 
#include <Python.h>
....
что нужно доставить (Ubuntu), чтобы этот модуль появился
$ gcc --shared c4py.c -o three.so
c4py.c:1:20: fatal error: Python.h: Нет такого файла или каталога
compilation terminated.

math-beginer
()

xsane перестала работать

Форум — Desktop

Переход на Ubuntu 11.04 привел еще к одной проблеме. Даже при старте c правами суперпользователя (sudo xsane) при попытках предварительного сканирования и сканирования выдает ошибку «невозможно инициализировать сканирование неверный аргумент».

что нужно проверить? Как можно вылечить?

сканер Epson Perfection V350, определяется правильно, но в отличии от предывдущей Убунты очень долго грузиться и затем отказывается сканировать с ошибкой :-(

math-beginer
()

формат для 3Д моделирования

Форум — Development

Подскажите простенький формат для хранения 3Д поверхностей заданных треугольными/четырехугольными площадками. Ситуация такая. Рисовалка и считалка самописные и разрабатываются, геометрию хотелось бы вводить существующим свободным инструментом под Линукс. Насчет рисовалки: разобрался mayavi и mathgl, понял что для самостоятельно рисования подходит такой формат: список всех уникальных координат и затем массив с индексами вершин площадок. Удалось нарисовать пока что-то в gmsh и он сбрасыает в формате msh сначала уникальные точки, затем элементы --- точки, линии, поверхности ... даже для простой сферы, STL --- вершины и зачем-то нормали (нормаль вроде должна по направлению обхода определяться) .... в итоге, пока не нашел подходящего формата. что salome сбрасывает не знаю, пока качается caelinux.

math-beginer
()

Salome не стартует

Форум — Desktop

После перехода на Ubuntu 11.04 Salome перестал стартовать:

$ runSalome 
runSalome running on tanya-desktop
Configure parser: Warning : could not find user configuration file
Searching for a free port for naming service: 2810 2811 2812 2813 - OK
Searching Naming Service ++++++++++Failed to narrow the root naming context
Traceback (most recent call last):
  File "/usr/lib/python2.7/dist-packages/salome/runSalome.py", line 617, in useSalome
    clt = startSalome(args, modules_list, modules_root_dir)
  File "/usr/lib/python2.7/dist-packages/salome/runSalome.py", line 417, in startSalome
    clt=orbmodule.client(args)
  File "/usr/share/pyshared/salome/orbmodule.py", line 49, in __init__
    self.initNS(args or {})
  File "/usr/share/pyshared/salome/orbmodule.py", line 82, in initNS
    sys.exit(1)
SystemExit: 1


--- Error during Salome launch ---

math-beginer
()

извлечь площадки

Форум — Development

первая проблема. Помогите с превращением эллипсоида в систему треугольных площадок. Есть матрицы x,y,z хранящие координаты точек эллипсоида, полученные по параметрическим уравнениям. Нужно преобразовать x,y,z в массив tri хранящий сегменты (площадки) с координатами трех точек (вершин площадок или рамок). Понимаю что моя попытка (второй листинг) совсем не верная, т.к. даже число точек будет меньше чем n1*n2 из-за полюсов. Направление обхода точек должно определять нормаль к эллипсойду по правилу правого (впрочем можно и левого) винта.

вторая проблема. не смог в mgrid передать переменные n1 и n2 из-за j.

#!/usr/bin/python
from numpy import *
from enthought.mayavi import mlab
rad=array([2.,2.,1.])
z=array([0.,0.,0.])
p,t=mgrid[0.:pi:4j,0:2.*pi:8j]
x=z[0]+rad[0]*sin(p)*cos(t)
y=z[1]+rad[1]*sin(p)*sin(t)
z=z[2]+rad[2]*cos(p)

mlab.clf()
mlab.mesh(x, y, z, color=(0, 1, 1))
mlab.mesh(x, y, z, representation='wireframe', color=(1, 0, 0))
mlab.show()
tri=zeros((32,3,3))
cycle=0
for i in range(3):
    for j in range(5):
        tri[i*5+j][0]=array([x[i][j],y[i][j],z[i][j]])
        tri[i*5+j][1]=array([x[i-1][j],y[i-1][j],z[i-1][j]])
        tri[i*5+j][2]=array([x[i][j-1],y[i][j-1],z[i][j-1]])
print tri
math-beginer
()

отключить gnome-power-man

Форум — Desktop

На буке gnome-power-man частенько «жрет память и процессорное время» и тормозит всю работу. Как его отключиь? работаю под Ubuntu 10.10 в xfce.

math-beginer
()

вылечить «sh: kuickshow: not found»

Форум — Development

при попытке использовать python-mathgl возникает такая проблема:

>> gr.ShowImage()
sh: kuickshow: not found

$ apt-cache search kuickshow
Пусто

Как лечить?

math-beginer
()

MathGL-Python

Форум — Development

Возникла необходимость сделать крависые 3Д картинки, в которых изображено несколько поверхностей и нужно регулировать прозрачность. Данные читаться должны из файлов.

Проблемы:

A) модуль mathgl не подгружается (ниже логи)

B) не могу найти пример с TriPlot и Read :-)

$ cat t.py 
#!/usr/bin/python
from mathgl import *

$ ./t.py
Traceback (most recent call last):
  File "./t.py", line 2, in <module>
    from mathgl import *
ImportError: No module named mathgl
$ apt-cache search python mathgl
libmgl-data - library for scientific graphs. (data files)
libmgl-dev - library for scientific graphs. (development files)
libmgl-fltk5 - library for scientific graphs. (FLTK runtime library)
libmgl-glut5 - library for scientific graphs. (GLUT runtime library)
libmgl-qt5 - library for scientific graphs. (Qt runtime library)
libmgl-wx5 - library for scientific graphs. (wxWidgets runtime library)
libmgl5 - library for scientific graphs. (main runtime library)
mathgl - библиотека для построения графиков (утилиты и примеры)
mathgl-doc - library for scientific graphs. (documentation)
math-beginer
()

генераторы трехмерных сеток

Форум — Development

Подскажите свободный софт под Линукс позволяющий генерировать трехметрные сетки на поверхностях. Поиск приводит к различным CAD-системам, ориентированные на конкретные задачи (для моей задачи точно нет кадов), а мне нужно просто изобразить пространственную фигуру и ее поверхность представить в виде сетки (системы четырех и трехугольных площадок).

math-beginer
()

в blender отобразитть поверхность из файла

Форум — Development

Требуется нарисовать несколько поверхностей (замкнутых и разомкнутых). Поверхности состоят из сегментов, каждый из которых определяется координатами вершин. Координаты вершин храняться в текстовых файлах. Это выглядит так:

$ cat surf1.txt
0.9160946164 2.2521700679 -0.0224572201
0.8784731664 2.2200436194 -0.3629388560
0.5246325267 2.3009418508 -0.3637445659
0.9160946164 2.2521700679 -0.0224572201

0.9160946164 2.2521700679 -0.0224572201
0.5246325267 2.3009418508 -0.3637445659
0.5497637301 2.3476743622 -0.0395248911
0.9160946164 2.2521700679 -0.0224572201

...

Вопросы:

1. Возможно ли в blender загрузить эти данные и отобразить

2. Наоборот, если я нарисую в blender поверхность могу ли я извлечь координаты сегментов в вышеприведенном формате.

прошу сослаться на учебники и примеры. может укажете на более подходящие средства, если blender не подходит.

math-beginer
()

RSS подписка на новые темы