«petsc» 태그된 질문

볼록하지 않은 전역 최적화 문제를 해결해야 할 몇 가지 문제가 있습니다. 현재 MATLAB의 Optimization Toolbox (특히 fmincon()algorithm = 사용 'sqp')를 사용하고 있습니다. 그러나 내 코드의 대부분은 Python이며 Python에서도 최적화를하고 싶습니다. 파이썬 바인딩과 경쟁 할 수있는 NLP 솔버가 fmincon()있습니까? 반드시 비선형 평등과 불평등 제약을 다룰 수있다 사용자가 Jacobian을 제공하도록 요구하지 않습니다. …

77 python  optimization  nonlinear-programming  linear-algebra  linear-solver  sparse-matrix  ode  automatic-differentiation  efficiency  matrix  symbolic-computation  pde  multiphysics  operator-splitting  fluid-dynamics  pde  multigrid  pde  hyperbolic-pde  quadrature  precision  numerical-analysis  fluid-dynamics  interpolation  c  ode  testing  data-sets  pde  fluid-dynamics  hyperbolic-pde  pde  finite-element  fluid-dynamics  stability  incompressible  pde  hyperbolic-pde  fluid-dynamics  pde  precision  floating-point  convex-optimization  conditioning  pde  hyperbolic-pde  stability  implicit-methods  explicit-methods  fluid-dynamics  accuracy  pde  hyperbolic-pde  petsc  linear-algebra 

KSP ( PETSc 의 선형 솔버 패키지)의 사전 조건부 Krylov 방법을 사용하여 부분 미분 방정식을 불연속 화하고 선형화하는 것과 같은 희소 선형 시스템을 풀 때 무엇이 ​​잘못 될 수 있습니까? 문제점에 어떤 문제가 있는지 판별하기 위해 어떤 단계를 수행해야합니까? 선형 시스템을 성공적이고 효율적으로 해결하기 위해 어떤 변경을 할 수 있습니까?

26 linear-algebra  petsc  preconditioning  krylov-method  iterative-method 

제가 알 수있는 한, 미국의 에너지 에너지 전산 과학 소프트웨어 프레임 워크는 PETSc 와 Trilinos 입니다. 언어 (C와 C ++)의 차이를 넘어 언뜻 보면 비슷해 보입니다. 두 프레임 워크의 주요 차이점은 무엇이며 다른 요소를 선택하는 데 영향을 미치는 요소는 무엇입니까? (기관의 편견과 기존 인프라는 무시하십시오.)

24 petsc  software  trilinos 

I가 사용하고 PETSc 의 비선형 솔버 패키지 SNES을 편미분 방정식을 이산화하여 얻은 비선형 방정식의 시스템을 해결하기 위해. 솔버가 수렴하지 않는 이유와 방정식을 성공적으로 풀기 위해 어떻게해야합니까?

22 petsc  pde  implicit-methods 

다중 고유 값과 벡터를 계산하기 위해 멀티 그리드를 사용할 계획이며 PETSc가 멀티 그리드에 대한 높은 수준의 지원을하는 것으로 나타났습니다. PETSc 문서에 따르면 PETSc의이 부분은 곧 교체 될 예정이므로 사용해서는 안됩니다. 어떤 다른 라이브러리가 멀티 그리드에 대한 높은 수준의 지원을 제공하고 있으며 PETSc가 새로운 멀티 그리드 지원을 얼마나 빨리 릴리스 할 …

17 petsc  libraries  multigrid 

옥타브에서 두 벡터의 유클리드 거리를 계산하는 빠른 방법이 있는지 알고 싶습니다. 특별한 기능이없는 것 같습니다. 따라서 수식을 사용해야 sqrt합니까?

16 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

PETSc가 Mat있고 상태 번호를 추정하고 싶습니다.

15 petsc  linear-algebra 

내가 본 곳마다 PETSc 튜토리얼 / 문서 등은 선형 대수학에 유용하며 일반적으로 스파 스 시스템이 도움이 될 것이라고 지정합니다. 고밀도 매트릭스는 어떻습니까? 조밀 한 대한 해결에 대해 걱정하고 있습니다.Ax=bAx=bAx=bㅏAA Fortran에서 CG 및 QMR에 대한 자체 코드를 작성했습니다. 기본적인 구성은 의사 코드를 분리하고 가능한 한 약간의 자체 조정으로 BLAS 루틴을 추가하는 …

14 linear-algebra  petsc  blas 

PETSc를 가지고 놀고 있었고 MPI 를 통해 둘 이상의 프로세스로 프로그램 을 실행할 때 더 느리게 실행되는 것으로 나타났습니다 ! 무슨 일이 일어나고 있는지 어떻게 확인할 수 있습니까?

14 linear-algebra  petsc 

일반적으로 시리얼 코드를 작성하고 내가 할 때 일부 xUnit 스타일 테스트 프레임 워크 (MATLAB xUnit, PyUnit / nose 또는 Google의 C ++ 테스트 프레임 워크)를 사용하여 단위 테스트를 작성합니다. 복잡한 Google 검색을 기반으로, 실무자가 MPI를 사용하는 코드를 테스트하는 방법을 많이 보지 못했습니다. 이에 대한 모범 사례가 있습니까? 단위 테스트 및 …

13 parallel-computing  petsc  mpi  testing  reproducibility 

질문 : 유한 요소 행렬의 희소 구조를 정확하고 효율적으로 계산하기 위해 어떤 방법을 사용할 수 있습니까? 정보 : 저는 2 중 라그랑주 단위로 Galerkin의 방법을 사용하고 C로 작성하고 희소 행렬 저장 및 KSP 루틴에 PETSc를 사용하여 포아송 압력 방정식 솔버를 작업하고 있습니다. PETSc를 효율적으로 사용하려면 전역 강성 매트릭스를위한 메모리를 미리 …

13 matrix  finite-element  petsc  performance 

내 시스템은 lagrange multipliers와 같은 대칭 FE 문제입니다 (예 : 비 압축 Stokes 흐름). ( A비비티씨)(ABTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} 여기서 이 일반적인 경우입니다 (라그랑 지 승수가 마지막에 나타나도록 방정식의 번호를 매겼습니다). 시스템은 상당히 큽니다 (+ 100k 라인).씨= 0C=0C = 0 이 질문에 대한 답변을 읽은 후 혼합 …

12 finite-element  petsc  krylov-method  preconditioning 

오른쪽 (300 ~ 1000)이 많은 동일한 희소 선형 시스템 (300x300 ~ 1000x1000)을 해결해야합니다. 이 첫 번째 문제 외에도 다른 시스템을 해결하고 싶지만 0이 아닌 동일한 요소 (단지 다른 값)로 일정한 희소성 패턴을 가진 많은 희소 시스템입니다. 내 행렬은 무기한입니다. 분해 및 초기화 성능은 중요하지 않지만 해결 단계의 성능은 중요합니다. 현재 …

12 linear-algebra  petsc  linear-solver  sparse-matrix 

petsc4py와 c / c ++ / fortran의 속도는 어느 정도입니까? 코드가 실행되는 코드에 크게 의존한다는 것을 알고 있지만 매트릭스 벡터 제품과 같은 간단한 것은 어떻습니까?

11 petsc  python 

다음은 C에서 함수 포인터의 일반적인 사용법입니다. Fortran에서 비슷한 작업을하고 싶습니다. 몇 가지 아이디어가 있지만 정식 방법이 있는지 알고 싶습니다. 사용자가 전달한 함수 포인터와 컨텍스트는 저장되고 나중에 호출됩니다. typedef PetscErrorCode (*TSIFunction)(TS,PetscReal,Vec,Vec,Vec,void*); PetscErrorCode TSSetIFunction(TS ts,Vec res,TSIFunction f,void *ctx); 사용자의 기능은 다양한 상황에서 컨텍스트를 사용하여 다시 호출됩니다. PETSc에서는 문자열-> 함수 포인터 테이블을 많이 …

11 petsc  fortran