전산 과학

컴퓨터를 사용하여 과학적 문제를 해결하는 과학자를위한 Q & A

6
불규칙한 경계가있는 도메인의 유한 차이
직사각형과 원 사이의 영역 (특히 서적 또는 링크)과 같이 불규칙한 기하학의 예를 포함하여 포아송의 수치 솔루션 (유한 차이 및 크랭크-니콜슨 방법) 및 확산 방정식에 관한 서적을 찾는 사람을 아무도 도울 수 있습니까? 이 경우 MATLAB 코드 예제)?
6
분자 진동을 시각화하는 데 사용할 수있는 오픈 소스 도구는 무엇입니까?
정상 모드 가 아닌 분자 진동을 시각화하고 싶습니다 . 모션의 정적, 벡터 표현을 표현하고 싶습니다. 벡터 스타일 (크기, 색상 등)의 유연성을 원합니다. 또한 진동 비디오를 제작하는 데 관심이 있습니다. 분자 진동을 표시하기위한 좋은 자료는 무엇입니까? 필자는 오픈 소스 도구를 선호하지만 상용 소프트웨어를 사용하는 것이 대안보다 낫다면 상용 소프트웨어를 사용하여 즐겁게 …
5
수학이나 계산 공부에 집중하는 것이 바람직합니까?
Krylov Subspace Methods에 대한 연구와 동시에 HPC 뒤의 수학을 한 단계 앞당기거나 계산 이론 (하드웨어, OS, 컴파일러 등)을 선택할 수 있습니다. 현재, 나는 단지 지나갈 정도로 충분히 알고있다 . 예를 들어, CG에 대한 방정식과 반복 방법의 기본을 도출하는 방법을 알고 있지만 사전 조건 및 수렴과 같은 세부 사항과 더 복잡한 …
2
과학 프로그램에서 Fortran의 함수 포인터로 작업하는 방법
다음은 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 
1
FEM에서 덩어리 질량 매트릭스를 공식화하는 방법
유한 요소법을 사용하여 시간 의존적 PDE를 풀 때 열 방정식과 같이 명시 적 시간 스텝핑을 사용하는 경우 질량 행렬로 인해 선형 시스템을 풀어야합니다. 예를 들어 열 방정식 예를 고수하면 ∂유∂티= c ∇2유∂유∂티=씨∇2유\frac{\partial{u}}{\partial{t}} = c\nabla{}^{2}u 앞으로 Euler를 사용하면 미디엄( un + 1− u엔디티) = - (C)의 K유엔미디엄(유엔+1−유엔디티)=−씨케이유엔M(\frac{u^{n+1}-u^{n}}{dt}) = -cKu^{n} 따라서 명시적인 …
2
파이썬에서 판다의 한계와 대안
나는 Pandas 가 금융계를 위해 처음 개발 된 곳 , 특히 자연 과학 (물리, 생물학 등)을 위해 처음 개발 된 곳을 읽었 으므로 비슷한 데이터 분석 Python 패키지가 더“자연 과학 중심”입니까? 방금 Pandas를 사용하기 시작했으며 다른 패키지 또는 수제 솔루션을 사용하지 않고 해결할 수없는 두 가지 문제가 이미 발생했습니다. 불확실성을 …
2
절단 된 SVD, 한 번에 하나의 특이 값 / 벡터 계산
한 번에 하나씩 특이 값을 계산하는 잘린 SVD 알고리즘이 있습니까? 내 문제 : 큰 밀도 행렬 M 의 첫 kkk 특이 값 (및 특이 벡터) 을 계산하고 싶지만 적절한 k 값이 무엇인지 모르겠습니다 . M 은 크므로 효율성 때문에 전체 SVD를 평가하지 않고 나중에 가장 작은 SV를 잘라 버릴 것입니다.MMMkkkMMM …
3
작은 점의 3D 볼록 껍질의 양이 선체에 모두 설정됩니다.
3D 를 제외하고 는 이전에 요청한 것과 비슷한 질문 이 있으며 실제 선체 모양이 아닌 볼륨 만 필요합니다. 보다 정확하게는, 3D로 작은 점 세트 (예 : 10-15)가 주어 졌는데,이 점들은 모두 점 세트의 볼록 껍질에 놓여있는 것으로 알려져 있습니다. 선체의 부피 만 계산하고 싶습니다. 실제 다면체 계산은 신경 쓰지 않습니다. …
2
그라데이션 하강 및 켤레 그라데이션 하강
프로젝트의 경우이 두 가지 방법을 구현하고 서로 다른 기능에서 어떻게 수행되는지 비교해야합니다. 등이 보이는 공액 구배 법을 위한의 연립 일차 방정식을 해결하는 것을 의미 Ax=bAx=b A\mathbf{x} = \mathbf{b} 여기서 AAA 는 대칭, 양의 정한 및 실수 인 nxn 행렬입니다. I 대해 읽으면 한편, 그라디언트 강하 I는 예를 참조 로젠 브록 …
1
Matlab에서 최적의 운송 왜곡 구현
나는 " 등록 및 뒤틀림에 최적의 대량 운송 "이라는 논문을 시행하고 있는데 , 온라인에서 유러 리아 대중 교통 코드를 찾을 수 없기 때문에 온라인으로 전환하는 것이 목표이며, 적어도 이미지 처리 분야의 연구 커뮤니티에게는 흥미로울 것입니다. 이 논문은 다음과 같이 요약 될 수 있습니다. - x 및 y 좌표를 따라 1D …
1
자율적 인 경우 일반 미분 방정식의 수치 근사 시스템에 대한 지름길이 있습니까?
미분 방정식 핸들 기능을 해결하기위한 기존의 알고리즘 ,Y∈RN. 그러나, 많은 실제 시스템에서, 차동 방정식이 자율적 그렇다DY디와이디티= f( y, t )dydt=f(y,t)\frac{dy}{dt} = f(y, t)y∈Rny∈Rny \in \mathbb R^n,Y∈R, N,와t는탈락. 이 간단한 가정으로 기존 수치 법에서 어떤 개선점을 볼 수 있습니까? 예를 들어,n=1이면 문제는t=∫dy로 바뀝니다.dydt=f(y)dydt=f(y)\frac{dy}{dt} = f(y)y∈Rny∈Rny \in \mathbb R^ntttn=1n=1n=1 우리는 1 차원 …
10 ode  numerics 
1
수치 소프트웨어 개발을위한 테스트 프레임 워크가 있습니까?
많은 계산 과학 프로그래밍에 표준 테스트 프레임 워크에서 다루지 않는 테스트 요구 사항이 있음을 발견했습니다. 계산 시간 테스트 알고리즘이 느려지지 않도록합니다. 나는 무언가를 할 수는 assureSmallerEqual(RuntimeWrapper(algorithm),53)있지만 알고리즘을 작업하면서 53 초 임계 값을 지속적으로 줄이려고합니다.assureSmallerEqual(RuntimeWrapper(algorithm),'previousbest+noisetolerance') 성능 시험 이전에 분석 솔루션에 대한 근사치를 찾은 알고리즘이 여전히 적어도 우수하거나 더 나은 솔루션을 찾습니다. …
10 testing 
1
"역 범죄"라는 문구가 처음 등장
역 문제에 대한 연구에서, 알려진 파라미터 세트로부터 합성 데이터 세트를 구성한 다음 반전 기술이 이러한 파라미터를 재구성 할 수 있는지 테스트하는 것이 일반적입니다. 이를 위해서는 합성 데이터에 적절한 수준의 랜덤 노이즈를 추가하는 것이 중요합니다. 또한, 합성 데이터를 계산하는 데 사용 된 방법이 유한 차이 또는 유한 요소 그리드를 기반으로하는 경우 …
당사 사이트를 사용함과 동시에 당사의 쿠키 정책과 개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.