Matrix must be a positive definite (geometry problem)

Question

GHULAM RASOOL 2015-2-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/176561-matrix-must-be-a-positive-definite-geometry-problem

回答： Erik S. 2015-2-9

Hi, The code i m using is originally for flow past half cylinder. i need to change the geometry so that the semi circle is replaced by half triangle. when i make this change by modifying geometry coordinates i get an error as follows:

Error using ==> chol

Matrix must be positive definite.

Error in ==> ellipt1 at 35

R = chol(A); Y = R'\B; Z = R\Y;

Error in ==> demo2 at 67

ZN = ellipt1(p,t,RDZ,[],W);

The part of the code where this 'chol' function is used and the geometry code is shown in the figures below:

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Erik S. 2015-2-9

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/176561-matrix-must-be-a-positive-definite-geometry-problem#answer_167406

Hi,

If you are sure the problem is set up correctely you can solve with the general solver:

Z = A\B;

Chol can only be used for special cases when your matrix A has special properties (Symmetric and positive definite). Most matrices are not and than you have to use the \ operator.