Описание тега petsc

PETSc, pronounced PET-see (the S is silent), is a suite of data structures and routines for the scalable (parallel) solution of scientific applications modeled by partial differential equations.
1 ответ

PETSc Solution линейная система с направляющей KSP

Я начинаю использовать библиотеку PETSc для параллельного решения линейной системы уравнений. Я установил все пакеты, собрал и успешно запустил примеры в папке petsc/src/ksp/ksp/examples/tutorials/, например ex.c Но я не мог понять, как заполнять ма…
30 май '12 в 11:32
2 ответа

OpenMDAO PetscTgtVecWrapper TypeError

Я пытаюсь запустить параллельный рабочий процесс, в котором я оцениваю более 1000 параллельных случаев внутри ParallelGroup, Если я запускаю на небольшом количестве ядер, это не дает сбоя, но увеличение количества узлов в какой-то момент вызывает ош…
14 май '16 в 16:44
0 ответов

Внедрение предварительно обусловленного градиента конъюгата с помощью PETSc

Я пишу код на C для запуска OLS на большом наборе данных, состоящем примерно из 40 миллионов наблюдений. Я заинтересован в использовании предварительно обусловленного алгоритма сопряженного градиента (PCG), использующего неполное разложение Холецког…
13 янв '17 в 21:20
0 ответов

Сбой файла Petsc Mex

У меня есть рабочий код, который решает определенную проблему с помощью KSPSolve от Petsc (с помощью gmres solver с предварительным кондиционером ilu). Я хотел вызвать этот код из MATLAB, поэтому я написал mex-файл. Однако KSPSolve вызывает мой mex-…
12 сен '14 в 09:59
1 ответ

Используя PETSc в Gfortran, но он не компилируется

Я хочу протестировать и использовать PESTc на фортране, поэтому я написал что-то очень простое и пустое, как (petscexe.F): PROGRAM petscexe INCLUDE 'petscsys.h' END PROGRAM petscexe чтобы увидеть, можно ли правильно импортировать заголовочный файл. …
13 апр '14 в 22:24
0 ответов

Структурированная сетка неправильной формы в PETSc (контекст DM)

У меня есть проблема конечных разностей в структурированной сетке в PETSc, а также DM контекст помогает мне создать матрицу действительно легко, DM дайте мне целую матрицу прямоугольной области, что хорошо для меня, потому что я использую эту систем…
03 янв '16 в 23:09
1 ответ

PetscMalloc против PetscMallocX

Каково практическое правило использования PetscMalloc2 (PetscMallocX) вместо PetscMalloc дважды (X раз)? Должны ли куски иметь одинаковые размеры или всегда эффективнее размещать их вместе / одновременно? В руководстве говорится "Выделение 2 (X) бло…
21 мар '14 в 18:35
4 ответа

Использование PETSc на make-файле в моем пользовательском make-файле

Я новичок в PETSc, У меня большой код на C++, и я хочу добавить PETSc к некоторым файлам, которые у меня уже есть, поэтому я должен изменить свой make-файл так, чтобы он также мог компилировать PETSc. Можно ли иметь два разных make-файла и затем выз…
24 янв '14 в 17:26
1 ответ

Заголовок PETSC#include'd в модуле

У меня есть модуль, который содержит глобальные переменные. Чтобы объявить некоторые глобальные переменные, мне нужно использовать HDF5. Я также использую библиотеку, поэтому мне также нужно включить заголовочный файл. Таким образом, преамбула globa…
20 фев '16 в 10:43
3 ответа

Пример PETSc не работает

У меня ОС Windows 7. Я следовал инструкциям на веб-странице PETSc; В командной строке VS 2005 я открыл Cygwin и установил PETSc с командой: ./configure --with-cc='win32fe cl' --with-fc=0 --with-mpi=0 --download-f2cblaslapack Я попытался запустить сл…
04 июл '12 в 14:44
1 ответ

Значения vec PETSC DMDA назначаются на место проведения

Недавно я начал изучать PETSc и столкнулся с проблемой при попытке выполнить простую задачу. Что не так с этим кодом: static char help[] = "Test 2d DMDAs Vecs.\n\n"; #include <petscdm.h> #include <petscdmda.h> #include <petscsys.h>…
05 авг '16 в 17:31
0 ответов

Вставка -nan+iG при вводе матрицы с использованием SLEPc

Я получил "Вставка -nan+iG error" в функции MatSetValues. Мой код выглядит так: Сначала я использую приведенный ниже код для изменения двойного в PETScScalar (я использую сложную версию). for(i=0;i<nz;i++)temp[i] = Ax[i] + Az[i] * PETSC_i; Затем …
07 май '14 в 10:00
2 ответа

do..while(0) функционально-подобная оболочка макроса C в c2hs

Я хотел бы обернуть функциональный макрос C в функцию C (и, в свою очередь, обернуть его в Haskell с помощью {#fun ... #} блок), но c2hs Препроцессор дросселей на do.. while(0) синтаксис; вот код: module TestMacro where #c #define TestF1(n) do{if n=…
23 апр '15 в 17:23
1 ответ

Решение большой системы линейных уравнений с использованием распределенных массивов MPI+CUDA и PETSc

Я хотел бы использовать библиотеку PETSc для решения больших распределенных по процессам систем линейных уравнений в моей собственной программе. Также я хотел бы задействовать доступные ресурсы GPU для этой цели. Я использую структурированную сетку …
11 ноя '16 в 11:08
2 ответа

Какова роль запятой в следующих утверждениях?

В настоящее время чтение в PETSc, когда я пришел к этому синтаксису в C/C++: PetscInt i, n = 10, col[3], its; PetscScalar neg_one = -1.0, one = 1.0, value[3]; Я не понимаю значения запятых здесь. Это связано с кортежами? Или там что-то перегружено?
02 ноя '16 в 09:47
0 ответов

Невозможно запустить многопоточный BLAS с поддержкой MKL в нескольких логических ядрах

Программное обеспечение, которое я разрабатываю для своей степени магистра, использует SLEPc(следовательно, PETSc) в качестве своего собственного решения, и я пытался связать PETSc с Intel MKL, чтобы использовать библиотеки BLAS MKL, которые, по мое…
04 апр '18 в 23:13
0 ответов

Время PETSC в режиме mpirun

Дорогие все, я новичок в изучении PETSC. Я написал очень простой 1D-код задачи диффузии на основе библиотеки PETSC (просто простой код FDM). Я хочу использовать параллельный решатель PETSC на каждом временном шаге. Вот псевдокод: Initial(C0); for(ti…
1 ответ

Новая установка демоверсии Fenics вылетает на MPICH_NUMVERSION

Я пытаюсь запустить базовую демонстрацию из Fenics (2016.2.0) на Ubuntu 16.04.2, Python 2.7, Anaconda (со Spyder для IDE). Демо ft_01poisson.py, Ошибка говорит: RuntimeError: In instant.recompile: The module did not compile with command 'make VERBOS…
04 май '17 в 04:38
1 ответ

Скудные разреженные матрицы в качестве входных данных для petsc4py

Я не могу найти способ, как эффективно загружать скудные разреженные матрицы, например csr_matrix в матрицу petsc4py, например PETSc.Mat().createAIJ, Я нашел эту ветку, но я не могу ее применить. Я также был бы признателен за указатель, где этот мат…
15 мар '13 в 21:47
1 ответ

Копия PETSc Matrix вызывает исключение: почему?

Я использую petsc4py и получаю исключение, которое не понимаю. Я определяю следующую функцию: def tsIJacobian(self, ts, t, u, udot, shift, A, B): self.setup_problem() psol = fe.as_backend_type(self.sol.vector()).vec() pA = fe.as_backend_type(self.A)…
30 июн '17 в 19:38