Deje que sea un triángulo y dejar que f sea una función suave en T .TTTFFfTTT Podemos usar la cuadratura de punto medio , donde x M es el centro de punto de T .∫Frex ≈ | TEl | ⋅f( xMETRO)∫FreX≈El |TEl |⋅F(XMETRO)\int f dx \approx |T|\cdot f(x_M)XMETROXMETROx_MTTT ¿Me puede proporcionar (una...
Estoy tratando de hacer coincidir elementos. Dado un conjunto de elementos que puedo clasificar en una escala de 0 a 100 de cuán similares son entre sí. Por ejemplo, si el elemento n 1 es leche y el elemento n 2 también es leche, entonces la similitud entre n 1 y n 2 sería del 100%. Si n 3 es...
Estoy tratando de escribir una implementación SVM completa en Python y tengo algunos problemas para calcular los coeficientes de Lagrange. Primero déjame reformular lo que entiendo del algoritmo para asegurarme de que estoy en el camino correcto. Si es un conjunto de datos e y i ∈ { - 1 , 1 } es...
Para la optimización, de Wikipedia : En informática, la metaheurística designa un método computacional que optimiza un problema al intentar iterativamente mejorar una solución candidata con respecto a una medida de calidad dada. Las metaheurísticas hacen pocas o ninguna suposición sobre el...
Tengo que hacer una transformación de coordenadas entre dos sistemas de referencia (ejes). Para eso, se deben multiplicar tres matrices ( ) debido a que se utilizan algunos ejes intermedios. He pensado en dos enfoques para resolver esto:3 × 33×33\times3 Método # 1 : Hacer la multiplicación...
¿La multiplicación de matrices (Mat * Mat y Mat * Vec) se escala con el número de ceros o con el tamaño de la matriz? O alguna combinación de los dos. ¿Qué pasa con la forma? Por ejemplo, tengo una matriz de 100 x 100 con 100 valores, o una matriz de 1000 x 1000 con 100 valores. Al cuadrar estas...
¿Compilar PETSc con una biblioteca externa BLAS / LAPACK afecta significativamente el rendimiento en matrices dispersas, o solo usa esas bibliotecas para matemáticas de matriz
Como una tarea en la universidad, hice una simulación 1d. El enunciado del problema fue resolver el problema del tubo de choque 1d que involucra gas ideal compresible como fluido de trabajo. Para este problema, resolví el sistema de ecuaciones de Eulers usando el solucionador Riemann de Roe. Quiero...
En este momento estoy trabajando en un método de resolución de ecuaciones diferenciales llamado colocación de base-spline. Con lo que tengo problemas es con la construcción de un método para construir una spline de orden arbitrario, con la relación con la condición inicial B 1 i (x)={...
Estoy tratando de diagonalizar algunas matrices densas y mal acondicionadas. En la precisión de la máquina, los resultados son inexactos (devuelve valores propios negativos, los vectores propios no tienen las simetrías esperadas). Cambié a la función Eigensystem [] de Mathematica para aprovechar la...
Soy muy nuevo en la programación GPGPU, así que perdónenme si la pregunta no es particularmente apropiada. Por lo que entiendo, la programación de GPU es un trabajo de ingeniería muy complejo en comparación con la programación habitual de CPU. Hay que tener mucho cuidado con los problemas de...
Estoy resolviendo el sistema de reacción-difusión de Turing con el siguiente código C ++. Es demasiado lento: para una textura de 128x128 píxeles, el número aceptable de iteraciones es 200, lo que resulta en 2.5 segundos de retraso. Necesito 400 iteraciones para obtener una imagen interesante, pero...
Estoy comenzando mis estudios de posgrado en ingeniería y trabajaré en proyectos de ciencias computacionales. Noté que ha habido una discusión sobre las ventajas y desventajas de implementar sus propios algoritmos aquí . ¿El análisis de costo-beneficio de implementar sus propios algoritmos versus...
¿Existe una implementación C abierta para la solución de ecuaciones cuárticas: a x ⁴ + b x ³ + c x ² + dx + e = 0ax⁴+bx³+cx²+dx+e=0ax⁴+bx³+cx²+dx+e=0 Estoy pensando en una implementación de la solución de Ferrari. En Wikipedia leí que la solución es computacionalmente estable solo para algunas de...
Considere el siguiente problema donde el término forzado puede depender de (vea la Edición 1 a continuación para la formulación) y y sus primeras derivadas. Esta es una ecuación de onda dimensional 1 + 1. Tenemos datos iniciales prescritos en .u , v W { u + v = 0 }Wtu v= FWuv=F W_{uv} = F u ,...
Estoy tratando de optimizar una aplicación MPI con un patrón de comunicación altamente asíncrono. Cada rango tiene una lista de cosas para calcular y envía mensajes según sea necesario si las entradas o salidas residen en un rango diferente. Además, cada rango está enhebrado (actualmente con un...
En los últimos años, han aparecido varias bibliotecas / proyectos de software que ofrecen una forma u otra de paralelismo de memoria compartida de uso general basado en datos. La idea principal es que, en lugar de escribir un código explícitamente enhebrado, los programadores implementan sus...
¿Cuál es la forma más moderna de implementar funciones especiales de doble precisión? Necesito la siguiente integral: Fmetro( t ) =∫10 0tu2 mmi- tu2du =γ( m + 12, t )2 tm + 12Fmetro(t)=∫0 01tu2metromi-ttu2retu=γ(metro+12,t)2tmetro+12 F_m(t) = \int_0^1 u^{2m} e^{-tu^2} d u = {\gamma(m+{1\over 2},...
Tengo una función numérica que f(x, y)devuelve un número de coma flotante doble que implementa alguna fórmula y quiero verificar que sea correcta contra las expresiones analíticas para todas las combinaciones de los parámetros xy yque me interese. ¿Cuál es la forma correcta de comparar lo calculado...
Estoy interesado en maximizar una función , donde .F( θ )f(θ)f(\mathbf \theta)θ ∈ Rpagθ∈Rp\theta \in \mathbb R^p El problema es que no conozco la forma analítica de la función o de sus derivados. Lo único que puedo hacer es evaluar la función de manera puntual, conectando un valor y obtener una...