Matrix must be a positive definite (geometry problem)

3 次查看（过去 30 天）
GHULAM RASOOL
GHULAM RASOOL 2015-2-9
回答： 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:

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

采纳的回答

Erik S.
Erik S. 2015-2-9
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.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by