27 from __future__
import print_function
34 from daal.data_management
import PackedTriangularMatrix, NumericTableIface, BlockDescriptor, readOnly, readWrite
36 utils_folder = os.path.realpath(os.path.abspath(os.path.dirname(os.path.dirname(__file__))))
37 if utils_folder
not in sys.path:
38 sys.path.insert(0, utils_folder)
39 from utils
import printArray
42 if __name__ ==
"__main__":
44 print(
"Packed triangular matrix example")
52 data = np.array([0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4], dtype=np.float64)
54 dataTable = PackedTriangularMatrix(NumericTableIface.lowerPackedTriangularMatrix, data)
56 block = BlockDescriptor()
59 dataTable.getBlockOfRows(firstReadRow, nRead, readOnly, block)
60 print(
"{} rows are read".format(block.getNumberOfRows()))
62 printArray(block.getArray(), nDim, block.getNumberOfRows(), block.getNumberOfColumns(),
63 "Print 3 rows from packed triangular matrix as float:")
67 dataTable.getBlockOfColumnValues(readFeatureIdx, firstReadRow, nDim, readWrite, block)
68 printArray(block.getArray(), 1, block.getNumberOfRows(), block.getNumberOfColumns(),
69 "Print the third feature of packed triangular matrix:")
72 dataBlock = block.getArray()
73 dataBlock[readFeatureIdx - 1] = -1
74 dataBlock[readFeatureIdx + 1] = -2
75 dataTable.releaseBlockOfColumnValues(block)
78 dataTable.getBlockOfRows(firstReadRow, nRead, readOnly, block)
79 print(
"{} rows are read".format(block.getNumberOfRows()))
80 printArray(block.getArray(), nDim, block.getNumberOfRows(), block.getNumberOfColumns(),
81 "Print 3 rows from packed triangular matrix as float:")