전산 과학

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

6
선형 프로그램에서
가정 minAvec(U)subject to Ui,j≤max{Ui,k,Uk,j},i,j,k=1,…,nminAvec(U)subject to Ui,j≤max{Ui,k,Uk,j},i,j,k=1,…,n\begin{align*} \min A &\mathrm{vec}(U) \\ &\text{subject to } U_{i,j} \leq \max\{U_{i,k}, U_{k,j}\}, \quad i,j,k = 1, \ldots, n \end{align*} 여기서 UUU 는 대칭 n×nn×nn\times n 행렬이고 U 2 를 n 2 항목 으로 1 차원 벡터로 vec(U)vec(U)\mathrm{vec}(U)재구성 합니다.UUUn2n2n^2 위의 프로그램에서 문제를 일으키는 부분은 max{⋅,⋅}max{⋅,⋅}\max\{⋅,⋅\} 입니다. …
2
거의 단일 시스템에 적용되는 반복 선형 솔버의 정지 기준
고려 와 가 거의 단수를 의미 고유치 의 매우 작게된다. 반복 방법의 일반적인 정지 기준은 잔여에 기초 과는 반복을 중단 할 수 간주 때 와 N 반복 수. 그러나 우리가 고려하고있는 경우 , 작은 고유 값 \ lambda_0 과 관련된 고유 공간에 큰 오차 v 가 존재할 수 있으며 이는 …
1
유한 체적 방법을 사용할 때 경계 조건을 어떻게 적용해야합니까?
이전 질문에 이어이 비 균일 유한 체적 메쉬에 경계 조건을 적용하려고합니다. 도메인의 lhs ( Robin 유형 경계 조건을 적용하고 싶습니다 x=xL)x=xL)x=x_L). σL=(dux+au)∣∣∣x=xLσL=(dux+au)|x=xL \sigma_L = \left( d u_x + a u \right) \bigg|_{x=x_L} 여기서 σLσL\sigma_L 은 경계 값이며; a,da,da, d 는 경계, 이류 및 확산에 각각 정의 된 계수이고; 는 경계에서 …
1
hdf5 대안
몇 년 동안 HDF5를 사용해 왔지만 데이터 세트의 크기가 커짐에 따라 여기에 나열된 것과 동일한 문제가 발생하기 시작합니다. http://cyrille.rossant.net/moving-away-hdf5/ 병렬 쓰기 지원-거대한 매트릭스의 청크 액세스 지원 내 일반적인 사용 사례는 100k x 100k 정수 행렬입니다. 논리적 관점에서 전체 파일로 사용하고 싶지만 병렬 작업자와 함께 덩어리로 덩어리를 작성해야합니다.
2
GPU에서 ODE 시스템을 해결하기위한 옵션?
나는 '사소한 병렬화'설정으로 ODE 시스템을 GPU에 농사 짓고 싶습니다. 예를 들어 512 개의 서로 다른 파라미터 세트로 감도 분석을 수행합니다. 이상적으로 Forward Euler와 같은 고정 된 시간 단계가 아니라 CVODE와 같은 스마트 적응 형 시간 단계 솔버를 사용하여 ODE를 해결하고 싶지만 CPU 대신 NVIDIA GPU에서 실행하고 싶습니다. 누구든지 이것을 했습니까? …
16 ode  gpu 
3
과학 프로그래밍 콘테스트
저는 정기적으로 소위 "프로그래밍 콘테스트 (Programming Contests)"에서 경쟁합니다. 여기서 제한된 시간 동안 자신의 코드로 문제를 해결하고 문제 해결 기술을 사용하여 어려운 알고리즘 문제를 해결할 수 있습니다. 이러한 모양에 대한 참조 예를 보려면 Google Code Jam 또는 ACM-ICPC와 같은 경쟁을 검색하십시오. (프로그래밍 콘테스트가 무엇인지 아는 경우 다음 단락을 건너 뛸 수 …
1
메쉬 의존적 안정성을 가진 요소의 유용성
3D Stokes 문제에서 요소의 안정성과 관련된 수학을 수행 한 후 이 임의의 4 면체 메쉬에 대해 안정적이지 않다는 것을 깨닫기 위해 약간 충격을 받았습니다 . 더 정확하게 말하면, 모든 노드와 4 개의 패싯 중 3 개의 패싯이 Dirichlet 조건이있는 도메인의 경계에있는 요소가있는 경우, 단일 행렬을 얻게됩니다. 스토크 스 시스템의 약한 …
5
왜 ODE의 수치 해가 불안정한 평형에서 멀어 지나요?
이중 진자 같은 시스템의 동작을 시뮬레이션하고 싶습니다. 이 시스템은 작동하지 않고 중력의 영향을받는 이중 진자처럼 동작하는 2 자유도 로봇 조작기입니다. 이중 진자와의 유일한 주요 차이점은 질량 중심에서 질량 및 관성 특성을 가진 두 개의 강체로 구성된다는 것입니다. 기본적으로 ode45Matlab에서 다음 유형의 ODE 시스템을 해결하도록 프로그래밍했습니다 . ⎡⎣⎢⎢⎢10000M110M1200100M120M22⎤⎦⎥⎥⎥⎡⎣⎢⎢⎢x˙1x˙2x˙3x˙4⎤⎦⎥⎥⎥=⎡⎣⎢⎢⎢x2−V1−G1x4−V2−G2⎤⎦⎥⎥⎥[10000M110M1200100M120M22][x˙1x˙2x˙3x˙4]=[x2−V1−G1x4−V2−G2] \left[ \begin{array}{cccc} 1 …
4
일련의 점에서 가장 흩어진 점 선택
부분 집합을 선택하는 (효율적인) 알고리즘이 있습니까?MMM 세트로부터 포인트 점 ( M &lt; N (크기의 모든 가능한 부분 집합을 통해) 예들은 "표지"대부분의 영역 M은 )?NNNM&lt;NM&lt;NM < NMMM 점이 2D 평면에 있다고 가정합니다. 순진한 알고리즘은 단순하지만 시간 복잡성 측면에서 금지 적입니다. for each subset of N points sum distance between each pair …
5
매우 고가의 목적 함수 문제에 대한 병렬 최적화 알고리즘
10-20 변수의 함수를 최적화하고 있습니다. 나쁜 소식은 각 함수 평가가 약 30 분의 직렬 계산으로 비싸다는 것입니다. 좋은 소식은 내가 수십 개의 계산 노드를 가진 클러스터를 가지고 있다는 것입니다. 따라서 모든 계산 능력을 효율적으로 사용할 수있는 최적화 알고리즘이 있습니까? 스펙트럼의 한쪽에는 철저한 검색이 있습니다. 전체 검색 공간을 미세 그리드로 세분화하고 …
7
부동 소수점에서 두 숫자의 평균에 대한 강력한 계산?
하자 x, y두 개의 부동 소수점 숫자합니다. 평균을 계산하는 올바른 방법은 무엇입니까? 순진 방법은 (x+y)/2오버 플로우가 발생할 때 수 x와 y너무 큽니다. 내 생각 엔 0.5 * x + 0.5 * y아마도 더 나은,하지만 (아마도 비효율적 인)이 곱셈을 포함, 그것은 충분히 좋은 경우에 나는 확실하지 않다. 더 좋은 방법이 있습니까? …
2
"이중 사전 조정"을 수행 할 방법이 있습니까?
질문: 당신이 대칭 긍정적 인 확실한 매트릭스에 대한 두 가지 (고려) preconditioners 있다고 가정 : 및 요인의의 역수 됩니다 적용하기 쉽습니다.ㅏㅏA≈ B티비ㅏ≈비티비A \approx B^TB≈ C티씨,ㅏ≈씨티씨,A \approx C^TC,B , B티, C, C티비,비티,씨,씨티B, B^T, C, C^T 때이 정보를 사용할 수 모두 와 중보다 더 나은 예비 - 건설 또는 형을?비비B씨씨C비비B씨씨C
당사 사이트를 사용함과 동시에 당사의 쿠키 정책과 개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.