Me estoy sumergiendo en el fascinante mundo del análisis de elementos finitos y me gustaría resolver un gran problema termomecánico (solo térmico mecánico, sin retroalimentación).→→\rightarrow Para el problema mecánico, ya entendí por la respuesta de Geoff , que necesitaré usar un solucionador...
Estoy tratando de optimizar un distribuidor de flujo en un tanque de modo que la distribución de velocidad y temperatura en cualquier sección transversal sea relativamente uniforme. Hay muchos parámetros que puedo ajustar a la uniformidad máxima de la sección transversal, como el número de tubos de...
Supongamos que tengo el sistema lineal original grande y disperso: . Ahora, no tengo A - 1 ya que A es demasiado grande para factorizar o cualquier tipo de descomposición de A , pero supongo que tengo la solución x 0 encontrada con una solución iterativa.Ax0=b0UNAX0 0=si0
Me gustaría predecir tiempos de ejecución para operaciones de álgebra lineal densa en una arquitectura específica utilizando una biblioteca específica. Me gustaría aprender un modelo que se aproxime a la función. Fo p: :Fopags:: ::F_{op} \;::\; tamaños de entrada runtime→→ \rightarrow para...
Como lo ejemplifica la respuesta de Jed Brown a los Costos de las búsquedas versus los cálculos , el uso de operaciones de punto flotante vectorizado versus no vectorizado resulta en un código mucho más rápido. Muchos compiladores modernos afirman que pueden realizar una vectorización automática....
Answer proporciona un software para calcular el exponente condicional de Lyapunov (CLE) para osciladores acoplados en sincronización de caos. Sin embargo, es difícil de seguir y no hay salida gráfica de la trama (y está en C, más complejo). ¿Alguien sabe cómo modificar LET Toolbox, que es muy bueno...
Se sabe que el método de Newton para resolver ecuaciones no lineales converge cuadráticamente cuando la suposición inicial está "suficientemente cerca" de la solución. ¿Qué es "suficientemente cerca"? ¿Existe literatura sobre la estructura de esta cuenca de
Me gustaría saber si hay un algoritmo que da un conjunto de puntos y un ángulo calcula el casco convexo si el ángulo es y dado un \ alpha> 0 calcula una envoltura que sigue más de cerca el "perímetro ".α=0α=0\alpha = 0α>0α>0\alpha > 0 Y si hay una definición de un perímetro no...
Al hojear algunos libros de texto, me di cuenta de que el problema de poner entre corchetes un mínimo inicialmente durante una búsqueda de línea tiende a ser una ocurrencia tardía (al menos en mis textos de pregrado). ¿Existen técnicas bien establecidas o mejores prácticas para este tipo de...
Quiero minimizar una función objetivo complicada, y no estoy seguro de si es convexa. ¿Existe un buen algoritmo que intente demostrar que no es convexo? Por supuesto, el algoritmo podría fallar en probar esto, en cuyo caso no sabría si es convexo o no, y esto está bien; Solo quiero tratar de...
Tengo la impresión, de recursos muy diferentes y de conversaciones con investigadores, de que existe una demanda creciente de cálculos de alta precisión en ecuaciones diferenciales parciales numéricas. Aquí, la alta precisión significa más precisión que solo la doble precisión estándar de 64...
Tengo un problema en el que necesito encontrar todos los pares propios positivos (ya que el valor propio es positivo) de una matriz no simétrica pequeña (generalmente más pequeña que 60x60). Puedo dejar de calcular cuando el valor propio es menor que un cierto umbral. Sé que los valores propios son...
Estoy tratando de implementar el método Broyden-Fletcher-Goldfarb-Shanno para encontrar el mínimo de una función. Necesito dos conjeturas iniciales y y una aproximación inicial de la matriz de Hesse . Los únicos requisitos que encuentro para es que si el Hessian es simétrico positivo definido,...
¿Existe una forma preferida de cómo implementar una evaluación rápida (aproximada) del polinomio de interpolación de Chebyshev en una cuadrícula uniforme (dados los valores de función en los nodos de Chebyshev)? Mi problema es que la interpolación se vuelve lenta cuando aumenta el grado del...
La mayor parte de mi programación son códigos de investigación únicos en C para mi propio uso. Nunca he distribuido ningún código a otros que no sean colaboradores cercanos. He desarrollado un algoritmo que estoy publicando en una revista científica. Quiero proporcionar el código fuente y quizás el...
Realización de una simulación numérica directa de la turbulencia isotrópica con el método pseudoespectral de Fourier (Orzag y Patterson, PRL, 1972) utilizando FFT. Para conocer los antecedentes del método, que se usa ampliamente en la comunidad de turbulencias, puede ver este curso:...
Quiero resolver donde es mi matriz de rigidez. Sin embargo, pueden faltar algunas restricciones y, por lo tanto, es posible que todavía haya algún movimiento rígido del cuerpo en el sistema (debido al valor propio cero). Como estoy usando CG para resolver el sistema lineal, esto no es aceptable ya...
A partir de la definición del número de condición, parece que se necesita una inversión de matriz para calcularlo, me pregunto si para una matriz cuadrada genérica (o mejor si es simétrica positiva definida) es posible explotar cierta descomposición de la matriz para calcular el número de condición...
Estoy usando la caja de herramientas Matde pde para resolver una cierta ecuación elíptica en 2D. La solución está bien, aunque necesito trazarla a lo largo de una línea dada, es decir, cortar un corte plano de la malla 3D que representa la solución. No puedo encontrar una manera que implique...
Para mi trabajo de investigación, he estado desarrollando una pequeña biblioteca de C ++ destinada a facilitar la comunicación entre los códigos computacionales de C ++ y Octave / Matlab (cuando este último se utiliza con fines de procesamiento posterior). Me gustaría lanzar una biblioteca de este...