E. Robert Tisdale said:
You can find examples of my C code at
http://www.netwood.net/~edwin/svmtl/
Just click on
The ANSI C Numerical Class Library
I had a quick look. Approximately 6000 lines of code in various .c and .h
files. (3200 when blank lines and comments are removed.)
Most of it wouldn't compile under a conforming C90 implementation, because
of heavy use of the inline keyword, which doesn't exist in C90.
I didn't spend very long reading the code, because I found it to be almost
unreadable. Here's a quick extract, taken entirely at random:
ncl_dcsubmatrix* /* M -= N */
(ncl_dcsubm_msub)(ncl_dcsubmatrix* pM, const ncl_dcsubmatrix* pN) {
ncl_extent mM = ncl_dcsubm_extent2(pM);
ncl_offset i = 0;
#ifdef NCL_ERRANT
ncl_extent mN = ncl_dcsubm_extent2(pN);
if (mM != mN) {
ncl_message("In function ncl_dcsubm_msub(ncl_dcsubmatrix*,\n"
" const ncl_dcsubmatrix*): submatrix differ in second extent.\n");
mM = (mM < mN)? mM: mN;
}
#endif /* NCL_ERRANT */
for (i = 0; i < mM; ++i) {
ncl_dcsubvector v = ncl_dcsubm_subv2(pM, i);
ncl_dcsubvector u = ncl_dcsubm_subv2(pN, i);
ncl_dcsubv_vsub(&v, &u);
ncl_dcsubv_destroy(&u);
ncl_dcsubv_destroy(&v);
}
return pM; }
Frankly, I find that unreadable.
Please remind me where I can find examples of your C programs.
In the literature.