1234567891011121314151617181920212223242526272829303132333435363738394041424344454647 |
- TODO
- Date_ 2009-10-03
- * Add ATLAS wrapper for blas level 2 and level 3 functions.
- * Include the HAVE_ATLAS definition in the configure system
- * Replace the lu-solver in lu.c with the one from the ccmath library
- and patch all modules using lu.c
- * Implement a robust parallelizable LU solver with pivoting
- * Move the band matrix cholesky decomposition solver from the lidar library
- (TcholBand.c) into the gmath library and make it available for other numerical
- programs (Integration in the common solver structure)
- TODO
- Date: 2006-11-20
- http://www.netlib.org/lapack/
- http://www.netlib.org/c/meschach/ (not GPL'ed, but excellent implementation!)
- la.c: write:
- new gmath function old MESCHACH
- G_vector_free() free vector memory #define V_FREE(vec) ( v_free(vec), (vec)=(VEC *)NULL )
- int v_free(vec)
- VEC *vec;
- G_vector_sub() Subtract two vectors: VEC *v_sub(vec1,vec2,out)
- VEC *vec1,*vec2,*out;
- G_vector_norm1() Calculate 1-norm of vector #define v_norm1(x) _v_norm1(x,VNULL)
- double _v_norm1(x,scale)
- VEC *x, *scale;
- G_matrix_scalar_mult() Scalar-matrix multiplication MAT *sm_mlt(scalar,matrix,out)
- double scalar;
- MAT *matrix,*out;
- G_matrix_stdin() stdin read-function "m_input()" #define m_input(mat) m_finput(stdin,mat)
- MAT *m_finput(FILE *fp,MAT *out);
- G_vector_stdin() stdin read-function "v_input()" #define v_input(vec) v_finput(stdin,vec)
- VEC *v_finput(FILE *fp,VEC *out);
- Sample text matrix for G_matrix_stdin()
- # comment
- # enter row by col
- Matrix: 4 by 6
- row0: 8.87 13.14 11.71 35.85 28.26 10.54
- row1: 13.59 20.12 19.61 70.66 34.82 16.35
- row2: 28.26 34.82 38.27 40.1 38.27 23.7
- row3: 10.54 16.35 23.7 38.98 40.1 38.98
|