08-Oct-2025 14:14:58 r8ccs_test(): MATLAB/Octave version 6.4.0 Test r8ccs(). r8ccs_dif2_test(): r8ccs_dif2() sets the second difference as an r8ccs matrix; Matrix rows M = 5 Matrix columns N = 5 Nonzeros NZ_NUM = 13 The r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 2 -1 0 0 0 2 -1 2 -1 0 0 3 0 -1 2 -1 0 4 0 0 -1 2 -1 5 0 0 0 -1 2 r8ccs_get_test(): r8ccs_get() gets an entry of an r8ccs matrix; Matrix rows M = 5 Matrix columns N = 5 Matrix nonzeros NZ_NUM = 12 The COLPTR vector: 1: 1 2: 4 3: 6 4: 8 5: 10 6: 13 The ROWIND vector: 1: 1 2: 2 3: 4 4: 1 5: 2 6: 3 7: 5 8: 4 9: 5 10: 1 11: 2 12: 5 The initial r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 0.197088 0.425303 0 0 -0.164254 2 -1.80331 0.184331 0 0 -0.629723 3 0 0 -1.85469 0 0 4 -0.714869 0 0 -0.921333 0 5 0 0 1.02905 -0.574895 1.35219 r8ccs_get retrieves 10 entries. I J K VALUE 5 3 7 1.029051 5 5 12 1.352189 1 5 10 -0.164254 3 3 6 -1.854689 2 5 11 -0.629723 5 4 9 -0.574895 1 5 10 -0.164254 4 4 8 -0.921333 2 1 2 -1.803312 4 1 3 -0.714869 r8ccs_ijk_test(): r8ccs_ijk() gets K from (I,J) for an r8ccs matrix. Matrix rows M = 5 Matrix columns N = 5 Matrix nonzeros NZ_NUM = 12 The COLPTR vector: 1: 1 2: 4 3: 6 4: 8 5: 10 6: 13 The ROWIND vector: 1: 1 2: 2 3: 4 4: 1 5: 2 6: 3 7: 5 8: 4 9: 5 10: 1 11: 2 12: 5 The initial r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 0.636029 0.248392 0 0 1.03414 2 -0.277544 -0.335276 0 0 0.663482 3 0 0 -0.153888 0 0 4 -0.317988 0 0 -0.649982 0 5 0 0 0.50396 -1.06275 -0.885702 r8ccs_ijk locates some (I,J) entries. I J K 3 3 6 4 1 3 2 2 5 2 2 5 4 2 -1 3 5 -1 1 4 -1 1 1 1 5 4 9 1 3 -1 1 5 10 5 4 9 3 2 -1 5 3 7 4 4 8 4 2 -1 2 3 -1 2 3 -1 2 3 -1 2 1 2 r8ccs_inc_test(): r8ccs_inc() increments entries in an r8ccs matrix. Matrix rows M = 5 Matrix columns N = 5 Matrix nonzeros NZ_NUM = 12 The COLPTR vector: 1: 1 2: 4 3: 6 4: 8 5: 10 6: 13 The ROWIND vector: 1: 1 2: 2 3: 4 4: 1 5: 2 6: 3 7: 5 8: 4 9: 5 10: 1 11: 2 12: 5 The initial r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 0.348365 0.227433 0 0 1.16099 2 -1.78855 -1.46035 0 0 -0.46822 3 0 0 -1.89173 0 0 4 -0.642256 0 0 -0.434122 0 5 0 0 1.16741 -1.3864 -0.895065 r8ccs_inc increments 10 entries at random. I J K NEW_VALUE 2 2 5 19.539653 4 1 3 21.357744 2 5 11 22.531780 2 5 11 46.531780 3 3 6 23.108269 2 1 2 24.211453 1 5 10 28.160994 2 5 11 74.531780 4 4 8 28.565878 3 3 6 53.108269 The final r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 0.348365 0.227433 0 0 28.161 2 24.2115 19.5397 0 0 74.5318 3 0 0 53.1083 0 0 4 21.3577 0 0 28.5659 0 5 0 0 1.16741 -1.3864 -0.895065 r8ccs_indicator_test(): r8ccs_indicator() sets up the indicator matrix; Matrix rows M = 5 Matrix columns N = 5 Nonzeros NZ_NUM = 12 The r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 0 15 2 21 22 0 0 25 3 0 0 33 0 0 4 41 0 0 44 0 5 0 0 53 54 55 r8ccs_kij_test(): r8ccs_kij() gets (I,J) from K in an r8ccs matrix. Matrix rows M = 5 Matrix columns N = 5 Matrix nonzeros NZ_NUM = 12 The COLPTR vector: 1: 1 2: 4 3: 6 4: 8 5: 10 6: 13 The ROWIND vector: 1: 1 2: 2 3: 4 4: 1 5: 2 6: 3 7: 5 8: 4 9: 5 10: 1 11: 2 12: 5 The initial r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 -1.10692 -0.837478 0 0 -2.26404 2 -0.0904034 -0.296562 0 0 1.26841 3 0 0 0.519366 0 0 4 1.80777 0 0 0.584844 0 5 0 0 1.05205 0.340708 1.05002 r8ccs_KIJ locates some K entries. K I J 8 4 4 2 2 1 4 1 2 3 4 1 10 1 5 7 5 3 8 4 4 11 2 5 9 5 4 10 1 5 2 2 1 10 1 5 7 5 3 2 2 1 8 4 4 10 1 5 2 2 1 7 5 3 12 5 5 5 2 2 r8ccs_mtv_test(): r8ccs_mtv() computes b=A'*x, where A is an r8ccs matrix. Matrix rows M = 5 Matrix columns N = 5 Matrix nonzeros NZ_NUM = 12 The r8ccs matrix A: Col: 1 2 3 4 5 Row --- 1 -0.794735 -0.410131 0 0 0.884379 2 -0.740538 0.00894998 0 0 -0.272377 3 0 0 3.14697 0 0 4 -0.660439 0 0 -1.07151 0 5 0 0 0.154931 2.05972 -2.30892 x: 1: 1 2: 0 3: 0 4: 0 5: -1 b=A'*x: 1: -0.794735 2: -0.410131 3: -0.154931 4: -2.05972 5: 3.1933 r8ccs_mv_test(): r8ccs_mv() computes b=A*x, where A is an r8ccs matrix; Matrix rows M = 5 Matrix columns N = 5 Matrix nonzeros NZ_NUM = 12 The r8ccs matrix A: Col: 1 2 3 4 5 Row --- 1 1.26656 2.00634 0 0 0.937224 2 0.656706 0.559716 0 0 0.278489 3 0 0 -0.0427643 0 0 4 -0.0795917 0 0 -0.639814 0 5 0 0 -1.17147 0.552505 1.39691 The x: 1: 1 2: 0 3: 0 4: 0 5: -1 b=A*x: 1: 0.329335 2: 0.378217 3: 0 4: -0.0795917 5: -1.39691 r8ccs_print_test(): r8ccs_print() prints a r8ccs matrix. Matrix rows M = 5 Matrix columns N = 5 Matrix nonzeros NZ_NUM = 12 The r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 1.31586 0.102621 0 0 -0.749232 2 -0.399352 0.580023 0 0 1.62605 3 0 0 -0.0347785 0 0 4 -2.06775 0 0 -0.743253 0 5 0 0 -1.412 -2.00217 1.46206 r8ccs_print_some_test(): r8ccs_print_some() prints some of an r8ccs matrix. Matrix rows M = 10 Matrix columns N = 10 Nonzeros NZ_NUM = 28 Rows 2-6, Cols 5-8: Col: 5 6 7 8 Row --- 2 0 0 0 0 3 0 0 0 0 4 405 0 0 0 5 505 506 0 0 6 605 606 607 0 r8ccs_random_test() r8ccs_random(): randomizes an r8ccs matrix. Matrix rows M = 5 Matrix columns N = 5 Nonzeros NZ_NUM = 12 The r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 -0.264191 -1.09941 0 0 -0.65239 2 0.79568 1.47994 0 0 0.810738 3 0 0 -0.0395275 0 0 4 1.09295 0 0 -0.264951 0 5 0 0 -0.491765 0.259628 0.362398 r8ccs_read_test(): r8ccs_read() reads an r8ccs matrix from a file. Matrix rows M = 5 Matrix columns N = 5 Nonzeros NZ_NUM = 12 Index base (0/1) = 1 The COL vector: 1: 1 2: 4 3: 6 4: 8 5: 10 6: 13 The ROW vector: 1: 1 2: 2 3: 4 4: 1 5: 2 6: 3 7: 5 8: 4 9: 5 10: 1 11: 2 12: 5 The r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 0 15 2 21 22 0 0 25 3 0 0 33 0 0 4 41 0 0 44 0 5 0 0 53 54 55 r8ccs_set_test(): r8ccs_set() sets entries in an r8ccs matrix; Matrix rows M = 5 Matrix columns N = 5 Matrix nonzeros NZ_NUM = 12 The COLPTR vector: 1: 1 2: 4 3: 6 4: 8 5: 10 6: 13 The ROWIND vector: 1: 1 2: 2 3: 4 4: 1 5: 2 6: 3 7: 5 8: 4 9: 5 10: 1 11: 2 12: 5 The initial r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 -2.5033 0.41607 0 0 -1.2388 2 0.546309 0.313849 0 0 -0.644124 3 0 0 -0.818011 0 0 4 1.83579 0 0 0.866662 0 5 0 0 -0.0888232 0.148829 0.724212 r8ccs_set sets 10 entries at random. I J K NEW_VALUE 4 4 8 101.000000 5 5 12 102.000000 1 2 4 103.000000 5 5 12 104.000000 5 4 9 105.000000 1 5 10 106.000000 2 5 11 107.000000 1 2 4 108.000000 5 3 7 109.000000 4 1 3 110.000000 The final r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 -2.5033 108 0 0 106 2 0.546309 0.313849 0 0 107 3 0 0 -0.818011 0 0 4 110 0 0 101 0 5 0 0 109 105 104 r8ccs_to_r8ge_test(): r8ccs_to_r8ge() converts a matrix from r8ccs to r8ge format. Matrix rows M = Matrix columns N = Nonzeros NZ_NUM = The r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 0 15 2 21 22 0 0 25 3 0 0 33 0 0 4 41 0 0 44 0 5 0 0 53 54 55 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 0 15 2 21 22 0 0 25 3 0 0 33 0 0 4 41 0 0 44 0 5 0 0 53 54 55 r8ccs_write_test(): r8ccs_write() writes an r8ccs matrix to 3 files Matrix rows M = 5 Matrix columns N = 5 Nonzeros NZ_NUM = 12 The COL vector: 1: 1 2: 4 3: 6 4: 8 5: 10 6: 13 The ROW vector: 1: 1 2: 2 3: 4 4: 1 5: 2 6: 3 7: 5 8: 4 9: 5 10: 1 11: 2 12: 5 The r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 0 15 2 21 22 0 0 25 3 0 0 33 0 0 4 41 0 0 44 0 5 0 0 53 54 55 r8ccs_zeros_test(): r8ccs_zeros() zeros an r8ccs matrix. Matrix rows M = 5 Matrix columns N = 5 Nonzeros NZ_NUM = 12 The r8ccs matrix: Col: 1 2 3 4 5 Row --- 1 0 0 0 0 0 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 5 0 0 0 0 0 r8ccs_test(): Normal end of execution. 08-Oct-2025 14:14:58