def genPrimes(): prime = 2 primes = [ prime ] yield prime while True: prime += 1 is_prime = 0 for i in primes: if prime % i != 0: is_prime += 1 if is_prime == len(primes): primes.append(prime) yield prime
Столкнулся с трудностями сборки библиотеки libMesh из SVN на Debian 5. Для сборки необходимо установить библиотеки lapack и blas. Они есть в репозитариях, поэтому проблем с этим не возникает. Но даже при линковке arm.C выдается ошибка о их отсутствии. Требуются файлы liblapack-3.so и libblas-3.so Сделав соответствующие симлинки в папке /usr/lib:
ln -s liblapack.so liblapack-3.so
и
ln -s libblas.so libblas-3.so
а также
ln -s libg2c.so.0.0.0 libg2c.so
компиляция завершилась успешно.
Кстати, перед сборкой необходимо задать переменную к библиотеке PETSC
В этой статье я расскажу как использовать открытые и свободные программы gmsh & getdp для решения задачи расчёта трёхмерного нестационарного температурного поля методом конечных элементов.
В первой части будет рассмотрено создание конечно-элементной сетки расчётной области. Для простоты в качестве расчётной области выбран параллелепипед размером, в мм, 100 x 200 x 250. Построение исходной геометрии и генерация сетки выполняется в программе gmsh. Во второй части будет описаны основные этапы расчёта в getdp и анализ результатов.
Итак, что такое gmsh и каковы его возможности. Gmsh это свободно распространяемая программа под лицензией GNU General Public License (GPL) позволяет создавать трёхмерные конечно-элементные сетки для расчёта МКЭ; снабжён САПР подобным движком для создания геометрии расчётной области; благодаря библиотеке OpenCascade позволяет импортировать геометрию в формате iges, step и включает довольно мощный постпроцессор. Включает сборки для всех основных платформ Windows, Linux, Mac OS X.
Построение конечно-элементной сетки для нашей задачи представлено в видеоролике.