Ciencia computacional

Preguntas y respuestas para científicos que usan computadoras para resolver problemas científicos.

2
resolver
Estoy transfiriendo un código existente de MATLAB a C ++ y tengo un sistema lineal para resolver (en lugar de la forma más típica )A x = bxA=bxA=bxA=bAx=bAx=bAx=b La matriz es densa y de forma general, pero no es mayor que 1000x1000. Entonces, en MATLAB, la solución o la notación …

5
Resolver repetidamente
Estoy usando MATLAB para resolver un problema que implica resolver en cada paso de tiempo, donde cambia con el tiempo. En este momento, estoy logrando esto usando MATLAB :bA x = bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}sib\mathbf{b}mldivide x = A\b Tengo la flexibilidad de hacer tantas calculaciones previas como sea necesario, por lo que …




2
Estrategias para el método de Newton cuando el jacobiano en la solución es singular
Estoy tratando de resolver el siguiente sistema de ecuaciones para las variables y x 2 (todo lo demás son constantes):PAG, x1P,x1P,x_1X2x2x_2 A(1−P)2−k1x1=0AP2−k2x2=0(1−P)(r1+x1)4L1−P(r1+x2)4L2=0A(1−P)2−k1x1=0AP2−k2x2=0(1−P)(r1+x1)4L1−P(r1+x2)4L2=0\frac{A(1-P)}{2}-k_1x_1=0 \\ \frac{AP}{2}-k_2x_2=0 \\ \frac{(1-P)(r_1+x_1)^4}{L_1}-\frac{P(r_1+x_2)^4}{L_2}=0 Puedo ver que puedo convertir este sistema de ecuaciones en una sola ecuación de una sola variable resolviendo las ecuaciones 1 y 2 para …


1
Solucionadores de PDE para difusión de deriva y modelos relacionados
Estoy tratando de simular modelos básicos de semiconductores con fines pedagógicos, comenzando por el modelo de difusión de deriva. Aunque no quiero usar un simulador de semiconductores estándar, aprenderé otros modelos (comunes, recientes u oscuros), sí quiero usar un solucionador de PDE estándar. Pero incluso para el caso simple de …
12 pde 

2
Integración numérica - manejo de NaNs (C / Fortran)
Estoy tratando con una integral complicada que exhibe NaNs a ciertos valores cercanos a cero y en este momento estoy tratando con ellos de manera bastante cruda usando una declaración ISNAN que establece el integrando a cero cuando esto ocurre. He intentado esto con la biblioteca NMS en FORTRAN (la …
12 quadrature 

1
Algoritmos para grandes matrices enteras dispersas
Estoy buscando una biblioteca que realice operaciones matriciales en grandes matrices dispersas sin sacrificar la estabilidad numérica. Las matrices serán 1000+ por 1000+ y los valores de la matriz estarán entre 0 y 1000. Realizaré el algoritmo de cálculo del índice, por lo que generaré vectores de fila (dispersos) de …


4
Precisión arbitraria Simulación de cuerda escalable
Estoy tratando de simular un objeto de cuerda. La formulación que entiendo es una matriz de partículas, conectadas por resortes. Estos resortes tienen valores k muy grandes, de modo que la línea se deforma, pero se estira muy poco. He concluido que resolver esto en función del tiempo no es …
12 simulation 



1
Algoritmos para el sistema lineal de EDO
Me pregunto: ¿cuál es el mejor algoritmo para resolver donde es una matriz real . A no es explícitamente dependiente del tiempo, generalmente escaso pero no necesariamente en bandas. Sus valores propios tienen partes reales no positivas. A también es diagonalizable, pero puede ser demasiado grande para que una diagonalización …

Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.