HideyoshiNakazone/Seals-NumericCalculus
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Better Implementation of Gauss Method

Browse Source 
Check and Validates if the Linear System is valid and better implementation of the returned result column matrix.
This commit is contained in:
Vitor Hideyoshi
2022-12-09 01:32:35 -03:00
parent a5cd253f0b
commit f29a8400d8
1 changed files with 3 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
yoshi_seals/process/process.pyx
View File

@@ -52,6 +52,8 @@ cpdef np.ndarray[np.float64_t, ndim=2] gauss(double[::,::] A, double[::,::] b):
if not c_pointer: if not c_pointer:
raise MemoryError() raise MemoryError()
if A.shape[0] < A.shape[1]:
raise ValueError("Invalid Matrix, this linear system is not solvable.")
for i in range(A.shape[0]): for i in range(A.shape[0]):
@@ -78,7 +80,7 @@ cpdef np.ndarray[np.float64_t, ndim=2] gauss(double[::,::] A, double[::,::] b):
sum_var += a[reversed_index][k]*x[k] sum_var += a[reversed_index][k]*x[k]
x[reversed_index] = (a[reversed_index][A.shape[1]] - sum_var)/a[reversed_index][reversed_index] x[reversed_index] = (a[reversed_index][A.shape[1]] - sum_var)/a[reversed_index][reversed_index]
return np.asarray(x).reshape(b.shape[0],b.shape[1]) return np.asarray(x).reshape(A.shape[1],b.shape[1])
cpdef np.ndarray[np.float64_t, ndim=2] cholesky(double[:,:] A, double[:,:] b): cpdef np.ndarray[np.float64_t, ndim=2] cholesky(double[:,:] A, double[:,:] b):