Leí esta respuesta y me di cuenta de que he estado usando la diferencia entre iteraciones sucesivas para definir un criterio de detención para un método iterativo de encontrar valores / vectores propios. ¿Cuáles son los buenos criterios de parada para los métodos iterativos que convergen en...
Tengo un conjunto de archivos que consta de puntos seleccionados al azar de un conjunto de datos, cada archivo pertenece a una clase en particular. Cada fila en estos archivos contiene las coordenadas en el espacio n del punto. Me gustaría comparar las distribuciones en n-space de cada uno de estos...
¿Qué métodos numéricos están disponibles para encontrar el punto fijo de un operador que actúa sobre las funciones ? Estoy buscando la función para la cual .f : [ a , b ] → [ a , b ] f A f = fUNAAAF: [ a , b ] → [ a , b ]f:[a,b]→[a,b]f : [a,b] \rightarrow [a,b]FffA f= fAf=fAf = f Detalles...
Recuerdo débilmente de mis primeras conferencias "numéricas" que los solucionadores lineales iterativos para menudo requieren que cuando se descompone comoAA x = bAx=bAx=bUNAAA A = D + MA=D+MA=D + M donde D es una matriz diagonal y tiene diagonal cero, los elementos de deben ser dominantes sobre...
Estoy investigando el uso de redes neuronales artificiales (ANN) para predecir las velocidades de reacción en mi fluido en lugar de resolver el sistema completo de EDO rígidas. Algunas personas de mi laboratorio ya han trabajado en eso, así que no empiezo desde cero, pero tengo problemas con mis...
El cálculo de los modos propios de una membrana semicircular se reduce al siguiente problema de valores propios ∇2u=k2u,∇2u=k2u,\nabla^2u=k^2u\;, donde la región de interés es un semicírculo definido por y φ ∈ [ 0 , π ] .r∈[0,1]r∈[0,1]r\in[0,1]φ∈[0,π]φ∈[0,π]\varphi\in[0,\pi] Es apropiado...
Estoy perfilando el rendimiento de los solucionadores lineales de PETSc. Según lo entiendo, speedup=Sequential TimeParallel Time.speedup=Sequential TimeParallel Time.\text{speedup}=\frac{\text{Sequential Time}}{\text{Parallel Time}}. Sé que ejecutar el código paralelo en un procesador se puede...
¿Cuáles son los "pros y contras" de estos dos métodos de refinamiento de malla? Ambos parecen ser los métodos predominantes. Naturalmente, puedo imaginar que el refinamiento rojo global es comparativamente fácil de implementar y atractivo para la investigación en métodos de cuadrículas múltiples,...
Estoy tratando de resolver un problema de mínimos cuadrados no lineales unilateral con restricciones lineales, es decir, el problema: minx∑mi=1ri(x) s.t Ax≤bminx∑i=1mri(x) s.t Ax≤b\min_{\mathbf{x}} \quad \sum^m_{i=1} \mathbf{r}_i(\mathbf{x}) \qquad \text{ s.t } \quad A\mathbf{x} \leq...
Tenemos un código de refinamiento de malla adaptable (AMR) para resolver la ecuación de onda elástica con interfaces de falla de fricción (basado en Chombo para aquellos que estén interesados). Una de las cosas que nos hemos dado cuenta es que nuestros resultados se ven fuertemente afectados por la...
Me gustaría probar la influencia de los parámetros de solvatación en los modelos de solvatación implícitos y preguntarme qué códigos están disponibles gratuitamente como programas independientes para el plegamiento de proteínas de proteínas pequeñas, y cuáles usan enfoques de minimización de...
He creado una spline para ajustar mis datos en python usando: spline=scipy.interpolate.UnivariateSpline(energy, fpp, k=4) La ecuación que quiero usar implica una suma entre n = 2 yn = infinito, donde n es el orden del diferencial en un punto Eo. Sin embargo,
Estoy jugando con implementaciones simples de algoritmos para encontrar la cobertura mínima de vértice / coincidencia de cardinalidad máxima en gráficos bipartitos. Sin embargo, parece que tengo problemas para buscar en Google algunos conjuntos de datos de prueba que podría usar para verificar, si...
Normalmente cuando siembro un generador secuencial de números aleatorios en C, uso la llamada srand(time(NULL)) luego usa rand() mod N para obtener un entero aleatorio entre 0 y N-1. Sin embargo, cuando hago esto en paralelo, las llamadas al tiempo (NULL) están tan cerca una de la otra que...
Estoy buscando una biblioteca C ++ conveniente y gratuita para gráficos que incluya un solucionador para el problema de subgrafo común máximo (MCS). Soy consciente de Boost Graph Library y LEMON , pero ninguno incluye un solucionador
Sé que debería usar una tolerancia para comparar números de coma flotante. Pero para comparar vectores, puedo pensar en 3 posibles soluciones correspondientes a diferentes métricas de distancia: Compare los componentes de cada vector individualmente: los vectores son iguales si los 3 están dentro...
¿Cuál es una buena forma común de expresar diccionarios (= mapas) en pseudocódigo? Es decir, las estructuras de datos que básicamente permiten almacenar valores para las claves, iterar sobre todos los pares clave / valor, probar la inclusión de una clave determinada, etc. Tengo algo como lo...
Suponga que resuelve un sistema lineal mediante un método iterativo, por ejemplo, gradientes conjugados o iteración de Richardson. Luego intenta resolver un sistema lineal que está ligeramente perturbado en la matriz y el lado derecho, por ejemplo, ˜ A ˜ u = ˜ f .Au=fAu=fAu = fA~u~=f~A~u~=f~\tilde...
¿Cuáles son algunos buenos recursos (libros, artículos, sitios) sobre la intersección de polígonos y los algoritmos de
Escribo códigos en Fortran y C para varios algoritmos matriciales. Sin embargo, cuando perfilo mis códigos usando VTune, generalmente encuentro alguna terminología que no puedo apreciar completamente. ¿Existe un buen recurso para aprender la creación de perfiles desde cero a un nivel bastante...