He implementado un solucionador multigrid V-Cycle utilizando una corrección de defectos lineales (LDC) y un esquema de aproximación completa (FAS). Mi problema es el siguiente: con LDC, el residuo se reduce en un factor de ~ 0.03 por ciclo. La implementación de FAS también converge con un factor...
Me pregunto si hay alguna lista de correo muy buena o grupos de Google para Computational Science, donde podamos discutir preguntas en lugar de solo hacer y responder preguntas. De hecho, estoy más interesado en la computación paralela y las soluciones numéricas de PDE. Pero no sé qué y cómo están...
Me preguntaba cuál es la mejor manera de aprender a programar en Fortran. He estado buscando libros, los ejemplos que me parecieron bien fueron "Explicado Fortran 95/2003" y "Fortran para científicos e ingenieros" . No sé qué piensan ustedes de esto. Otra cosa que me pregunto es qué compilador...
Estoy escribiendo un pequeño paquete MATLAB que resolverá ciertos problemas numéricos de clase. Hay 3 etapas del algoritmo y el usuario tiene 5 opciones para cada etapa. He implementado todo el problema usando funciones y 3 casos de cambio (uno para cada etapa de algoritmo). Funciona bien, pero...
He estado tratando de depurar este error los últimos días, me preguntaba si alguien tiene consejos sobre cómo proceder. Estoy resolviendo la ecuación de Poisson para una distribución de carga escalonada (un problema común en la electrostática / física de semiconductores) en una malla de volumen...
Considere el problema de optimización estrictamente convexo sin restriccionesDeje que x_ \ text {opt} denote sus mínimos únicos y x_0 sea una aproximación inicial dada a x_ \ text {opt}. Llamaremos a un vector x una \ epsilon- solución cercana de \ mathcal {O} if \ begin {ecation} \ frac {|| x...
Estas preguntas pueden estar un poco fuera de tema en comp-sci. si es necesario, sugiera dónde encaja. La pregunta es sobre cómo administrar todas las ejecuciones de simulación de manera eficiente. digamos, por ejemplo, una simulación requiere la fijación de 2 parámetros que deben definirse en...
Tengo un gran problema de valor propio cúbico: (A0+λA1+λ2A2+λ3A3)x=0.(UNA0 0+λUNA1+λ2UNA2+λ3UNA3)X=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Podría resolver esto convirtiendo a un problema de valor propio lineal, pero daría...
¿Cuáles son las principales diferencias entre FEM y XFEM? ¿Cuándo deberíamos (no) usar XFEM intead de FEM y viceversa? En otras palabras, cuando encuentro un nuevo problema, ¿cómo puedo saber cuál de ellos debo
¿Hay alguna regla general simple para decir si vale la pena hacer SOR en lugar de Gauss-Seidel? (y posible forma de estimar el parámetro de realxation )ωω\omega Me refiero solo a mirar la matriz , o al conocimiento de un problema particular que representa la matriz. Estaba leyendo la respuesta a...
Muchas bibliotecas FEM adaptativas utilizan estructuras de datos de malla más avanzadas para manejar la adición / eliminación de nodos, bordes, triángulos, tetraedros, etc. Por ejemplo, la biblioteca p4est usa estructuras de datos octree para el refinamiento de malla adaptable; a menudo no...
Siempre he escuchado que la paralelización fácil era una de las ventajas de los métodos DG, pero realmente no veo por qué alguna de esas razones no se aplica también a Galerkin
double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Lo anterior...
¿Cuál es la diferencia entre estos dos métodos? ¿Se puede resolver un problema con un método que se pueda resolver con el otro? ¿Pueden ambos / o uno de ellos estar en paralelo con OpenMP y / o
Tengo un proyecto donde necesito usar campos cuadráticos Específicamente números de la forma a + b - 3---√una+si-3a + b \sqrt{-3} cona , b ∈ Quna,si∈Qa,b \in \mathbb{Q}. Por ejemplo, aquí están los números primos en enteros de Eisenstein : No quiero usar salvia. Me gustaría escribir mi propio...
En el informe técnico sobre Galahad [1], los autores afirman, en el contexto de problemas generales de programación no lineal, En nuestra opinión, nunca había habido muchas dudas de que los métodos SQP [programación cuadrática secuencial] tendrían más éxito [que los métodos lagrangianos...
Estoy realizando una búsqueda de línea como parte de un algoritmo cuasi-Newton BFGS. En un paso de la búsqueda de línea, uso una interpolación cúbica para acercarme al minimizador local. Sea ser la función de interés. Quiero encontrar una tal que .x ∗ f ′ ( x ∗ ) ≈ 0f:R→R,f∈C1f:R→R,f∈C1f : R...
Deseo modelar una caña de pescar (o una cuerda) uniendo segmentos cortos. (Los segmentos pueden tener la misma longitud (corta) pero a cada segmento se le debe asignar su propia masa individual). Un segmento influirá en el siguiente por el par entre los segmentos. Por el momento, las juntas se...
Estoy interesado en calcular la solución de un sistema lage de EDO utilizando un método krylov como en [1]. Dicho método involucra funciones relacionadas con la exponencial (las llamadas funciones ). Básicamente consiste en calcular la acción de la función de matriz mediante la construcción de un...
Algunos libros antiguos que he visto dicen que el número mínimo de etapas de un método Runge-Kutta explícito de un orden específico es desconocido para los pedidos . ¿Sigue siendo cierto?≥ 9≥9 9\geq 9 ¿Qué bibliotecas hay para trabajar automáticamente con métodos Runge-Kutta de alto...