Perguntas com a marcação «matrix»

Eu tenho vários problemas desafiadores de otimização global não convexa para resolver. Atualmente, uso o Optimization Toolbox do MATLAB (especificamente, fmincon()com o algoritmo = 'sqp'), o que é bastante eficaz . No entanto, a maior parte do meu código está em Python, e eu adoraria fazer a otimização em Python …

77 python  optimization  nonlinear-programming  linear-algebra  linear-solver  sparse-matrix  ode  automatic-differentiation  efficiency  matrix  symbolic-computation  pde  multiphysics  operator-splitting  fluid-dynamics  pde  multigrid  pde  hyperbolic-pde  quadrature  precision  numerical-analysis  fluid-dynamics  interpolation  c  ode  testing  data-sets  pde  fluid-dynamics  hyperbolic-pde  pde  finite-element  fluid-dynamics  stability  incompressible  pde  hyperbolic-pde  fluid-dynamics  pde  precision  floating-point  convex-optimization  conditioning  pde  hyperbolic-pde  stability  implicit-methods  explicit-methods  fluid-dynamics  accuracy  pde  hyperbolic-pde  petsc  linear-algebra 

Sabemos que é simétrico e positivo-definido. Sabemos que é ortogonal:BUMAA\mathbf ABB\mathbf B Pergunta: simétrico e positivo-definido? Resposta: Sim.B ⋅ A ⋅ B⊤B⋅A⋅B⊤\mathbf B \cdot\mathbf A \cdot\mathbf B^\top Pergunta: Um computador poderia ter nos dito isso? Resposta: Provavelmente. Existem sistemas de álgebra simbólica (como o Mathematica) que manipulam e propagam fatos …

36 matrix  symbolic-computation 

Nas estatísticas e suas diversas aplicações, geralmente calculamos a matriz de covariância , que é positiva definida (nos casos considerados) e simétrica, para vários usos. Às vezes, precisamos do inverso dessa matriz para vários cálculos (formas quadráticas com esse inverso como a (única) matriz central, por exemplo). Dadas as qualidades …

27 linear-algebra  matrix 

O que é um algoritmo simples para calcular o SVD de matrizes?2×22×22 \times 2 Idealmente, eu gostaria de um algoritmo numericamente robusto, mas gostaria de ver implementações simples e não tão simples. Código C aceito. Alguma referência a documentos ou código?

26 linear-algebra  matrix  reference-request  svd  matrix-factorization 

Eu queria saber se existe um método rápido e eficiente para encontrar antecipadamente o número de não-zeros para operações de multiplicação de matrizes esparsas, assumindo que ambas as matrizes estejam no formato CSC ou CSR. Eu sei que existe um no pacote smmp, mas preciso de algo que já esteja …

17 matrix  sparse-matrix 

Estou ciente de que inverter uma matriz para resolver um sistema linear não é uma boa ideia, pois não é tão preciso e eficiente como solucionar diretamente o sistema ou usar a decomposição LU, Cholesky ou QR. No entanto, não pude verificar isso com um exemplo prático. Eu tentei este …

16 matrix  accuracy  inverse 

Gostaria de saber se existe uma maneira rápida de calcular a distância euclidiana de dois vetores na oitava. Parece que não há função especial para isso, então devo usar apenas a fórmula com sqrt?

16 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

Na programação de cálculos de matriz densa, existe alguma razão para escolher um layout de linha principal do layout acima da coluna? Eu sei que, dependendo do layout da matriz escolhida, precisamos escrever o código apropriado para usar as memórias de cache efetivamente para fins de velocidade. O layout principal …

16 matrix  fortran 

Ao calcular a fatoração QR na prática, utiliza-se as reflexões de Householder para zerar a parte inferior de uma matriz. Eu sei que, para calcular valores próprios de matrizes simétricas, o melhor que você pode fazer com as reflexões de Household é obtê-lo na forma tridiagonal. Existe uma maneira óbvia …

16 linear-algebra  matrix 

Um problema comum em estatística é calcular a raiz quadrada inversa de uma matriz definida positiva simétrica. Qual seria a maneira mais eficiente de calcular isso? Encontrei alguma literatura (que ainda não li) e algum código R incidental aqui , que reproduzirei aqui por conveniência # function to compute the …

15 linear-algebra  numerical-analysis  matrix 

Desculpe pelo longo post, mas eu queria incluir tudo o que achava relevante desde o início. O que eu quero Estou implementando uma versão paralela dos métodos de subespaço de Krylov para matrizes densas. Principalmente GMRES, QMR e CG. Percebi (após a criação de perfil) que minha rotina DGEMV era …

15 parallel-computing  matrix  lapack  blas 

Quais abordagens são usadas na prática para estimar o número de condições de grandes matrizes esparsas?

15 linear-algebra  matrix  conditioning 

Em um projeto de software em que estou trabalhando, certos cálculos são muito mais fáceis para matrizes densas de baixa patente. Algumas instâncias de problemas envolvem matrizes densas de baixo escalão, mas são dadas a mim na íntegra, e não como fatores, portanto, terei que verificar o rank e fatorar …

13 linear-algebra  matrix  lapack  rank  matrix-factorization 

Eu tenho dois gráficos com quase n ~ 100000 nós cada. Nos dois gráficos, cada nó é conectado a exatamente 3 outros nós, portanto a matriz de adjacência é simétrica e muito esparsa. A parte difícil é que eu preciso de todos os autovalores da matriz de adjacência, mas não …

13 linear-algebra  sparse-matrix  eigenvalues  matrix 

O título é a pergunta. Essa técnica envolve o uso da "matriz de cofatores", ou "matriz de adjuntos" e fornece fórmulas explícitas para os componentes do inverso de uma matriz quadrada. Não é fácil fazer manualmente uma matriz maior que, digamos, 3×33×33\times 3 . Para um n×nn×nn\times n matriz, requer …

13 linear-algebra  matrix  complexity