Quiero publicar el paquete de software que he escrito para mi trabajo de posgrado junto con un documento que describe el paquete. Hasta ahora he estado considerando licencias liberales como BSD. Sin embargo, ahora sé que publicar su código puede implicar un acuerdo con los términos de licencia de...
Antecedentes: Actualmente estoy ejecutando una gran cantidad de experimentos de variación de parámetros. Se están ejecutando en Python 2.6+, usando numpy. Estos experimentos demorarán aproximadamente 2 semanas en ejecutarse. Aproximadamente estoy variando 3 parámetros (variables independientes)...
Codifiqué un problema de Stokes Flow usando elementos finitos y estoy en el proceso de verificar que funcione. Simplemente no estoy seguro de qué tasa de convergencia debería esperar ya que refino globalmente la malla. Sé que para problemas escalares usando funciones de base lineal esperaría...
Estoy entrando en mi tercer año de mi programa de doctorado, y hasta ahora mi código de investigación (problemas numéricos inversos / imagen médica / procesamiento de imágenes / etc.) consiste principalmente en scripts y funciones MATLAB desorganizados, con algunos C ++ y CUDA MEX'd rutinas...
Durante mi deambulación en Mathematica.se , gradualmente me di cuenta de que un cierto tipo de problema de resolución de ecuaciones diferenciales nos "preocupa" todo el tiempo, es decir, el problema del valor límite (BVP) de las ecuaciones diferenciales ordinarias no lineales (EDO). El método de...
Al escribir el código MPI 3.0, ¿debería pasar argcy argva la MPI_Initllamada o no, y por qué? EDITAR: Dado que la respuesta de Bill Barth planteó algunas preguntas, quiero hacer algunos comentarios: Pasar argc/ argvno es necesario desde MPI 1.1. La pregunta es específicamente sobre por qué...
dada una matriz simétrica , y una matriz arbitraria X ∈ R n × n , y un vector v ∈ R n × 1 , ¿es posible calcular la siguiente expresión en el tiempo O ( n 2 ) ?Y∈Rn×nY∈Rn×nY \in \mathbb{R}^{n \times n}X∈Rn×nX∈Rn×nX \in \mathbb{R}^{n \times n}v∈Rn×1v∈Rn×1v \in \mathbb{R}^{n \times...
Estoy buscando un método o algoritmo para producir un valor que describa cuán diferentes son dos mallas geométricamente pero que tienen topologías diferentes. Un ejemplo serían algunos datos CAD a los que se les ha aplicado una malla FE y la parte posterior a su fabricación y escaneo en una malla...
Ahora que Formulize / Eureqa comenzó a cobrar $ 2500 al año por usarlo y haber estropeado la versión de prueba, ¿alguien sabe de algún reemplazo que pueda hacer cosas similares como encontrar una ecuación dados los puntos de datos? Era gratis, pero ahora está fuera del rango de precios de...
He estado leyendo muchos documentos sobre ADMM últimamente, y también intenté resolver varios problemas al usarlo, en todos los cuales fue muy efectivo. A diferencia de otros métodos de optimización, no puedo tener una buena intuición de cómo y por qué este método es tan efectivo (por supuesto, he...
En un artículo reciente de SIAM News , hay un artículo largo que describe una organización sistemática de los elementos finitos, acertadamente denominada Tabla periódica de elementos finitos . Es realmente fascinante ver cómo se puede lograr la clasificación a través del cálculo exterior de...
Para un protector de pantalla tonto que estoy tratando de desarrollar, me gustaría generar aleatoriamente una matriz 2D de vectores 2D sin divergencia, y luego usarla para generar un gráfico de convolución integral de línea. He escuchado que una forma de hacer esto es generar ruido aleatorio y...
Tengo un vector VVV que puede descomponerse en el espacio propio del operador escaso hermitiano MMM : V=∑ivim^iV=∑ivim^iV = \sum_i v_i \hat{m}_i ¿Hay una manera de encontrar el m i (el vector propio sí mismo) que se corresponden con la mayor v i (en magnitud)?m^im^i\hat{m}_iviviv_i Básicamente,...
Me preguntaba si hay algún teorema que me permita poner un límite superior en el error introducido al omitir elementos de matriz pequeños de una matriz antes de la diagonalización. Vamos a suponer que tenemos una gran matriz, cuya matriz de elementos van de 111 a 10- 1510−1510^{-15} . Si tuviera...
Traté de ir a las fuentes primarias para comprender cómo usar las tablas de Butcher para simplificar el álgebra que necesito hacer cuando uso la serie Taylor para encontrar el orden de precisión de un esquema, por ejemplo. Sin embargo, tal vez debido a la falta de antecedentes relevantes, me...
Es un truco bastante conocido para evitar la división en el cálculo de raíces cuadradas para aplicar el método de Newton para encontrar1 / x--√1/ /X1/\sqrt{x} , y probablemente más conocido, utilizando el método de Newton para encontrar recíprocos sin división . Al rescatar un hilo StackOverflow...
Tengo una pregunta bastante específica con respecto al número de condición. Ejecuto simulaciones FEM que tienen escalas de longitud múltiple para ellos, lo que resulta en una gran disparidad entre las entradas más grandes y las entradas más pequeñas en mi matriz. El número de condición puede llegar...
¿Hay alguna biblioteca / caja de herramientas que tenga implementación de SVD incremental en MATLAB? He implementado este documento , es rápido pero no funciona bien. Intenté esto, pero en este también el error se propaga rápidamente (dentro de la actualización, el error de 5-10 puntos es...
Las ecuaciones de movimiento para un sólido elástico están dadas por ∇⋅σ+f=ρu¨σ=Cεε=12(∇u+[∇u]T)∇⋅σ+f=ρu¨σ=Cεε=12(∇u+[∇u]T)\begin{align} &\nabla \cdot \boldsymbol{\sigma} + \mathbf{f} = \rho \ddot{\mathbf{u}}\\ &\boldsymbol{\sigma} = \mathbb{C}\boldsymbol{\varepsilon}\\ &\boldsymbol{\varepsilon} =...
Supongamos que tenemos un problema de valor inicial de la forma donde x 0 ∈ R n se conoce exactamente (es decir, con precisión ilimitada) y podemos evaluar eficientemente f : R n → R n con cualquier precisión. Es decir, tenemos un cuadro negro que, dado un vector x ∈ R n y un entero M , devuelve...