hmath-python 1.0.0__tar.gz
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- hmath_python-1.0.0/LICENSE +373 -0
- hmath_python-1.0.0/PKG-INFO +103 -0
- hmath_python-1.0.0/README.md +81 -0
- hmath_python-1.0.0/bindings/python/hmath.py +283 -0
- hmath_python-1.0.0/bindings/python/hmath_cffi.py +127 -0
- hmath_python-1.0.0/bindings/python/hmath_python.egg-info/PKG-INFO +103 -0
- hmath_python-1.0.0/bindings/python/hmath_python.egg-info/SOURCES.txt +10 -0
- hmath_python-1.0.0/bindings/python/hmath_python.egg-info/dependency_links.txt +1 -0
- hmath_python-1.0.0/bindings/python/hmath_python.egg-info/requires.txt +6 -0
- hmath_python-1.0.0/bindings/python/hmath_python.egg-info/top_level.txt +2 -0
- hmath_python-1.0.0/pyproject.toml +33 -0
- hmath_python-1.0.0/setup.cfg +4 -0
|
@@ -0,0 +1,373 @@
|
|
|
1
|
+
Mozilla Public License Version 2.0
|
|
2
|
+
==================================
|
|
3
|
+
|
|
4
|
+
1. Definitions
|
|
5
|
+
--------------
|
|
6
|
+
|
|
7
|
+
1.1. "Contributor"
|
|
8
|
+
means each individual or legal entity that creates, contributes to
|
|
9
|
+
the creation of, or owns Covered Software.
|
|
10
|
+
|
|
11
|
+
1.2. "Contributor Version"
|
|
12
|
+
means the combination of the Contributions of others (if any) used
|
|
13
|
+
by a Contributor and that particular Contributor's Contribution.
|
|
14
|
+
|
|
15
|
+
1.3. "Contribution"
|
|
16
|
+
means Covered Software of a particular Contributor.
|
|
17
|
+
|
|
18
|
+
1.4. "Covered Software"
|
|
19
|
+
means Source Code Form to which the initial Contributor has attached
|
|
20
|
+
the notice in Exhibit A, the Executable Form of such Source Code
|
|
21
|
+
Form, and Modifications of such Source Code Form, in each case
|
|
22
|
+
including portions thereof.
|
|
23
|
+
|
|
24
|
+
1.5. "Incompatible With Secondary Licenses"
|
|
25
|
+
means
|
|
26
|
+
|
|
27
|
+
(a) that the initial Contributor has attached the notice described
|
|
28
|
+
in Exhibit B to the Covered Software; or
|
|
29
|
+
|
|
30
|
+
(b) that the Covered Software was made available under the terms of
|
|
31
|
+
version 1.1 or earlier of the License, but not also under the
|
|
32
|
+
terms of a Secondary License.
|
|
33
|
+
|
|
34
|
+
1.6. "Executable Form"
|
|
35
|
+
means any form of the work other than Source Code Form.
|
|
36
|
+
|
|
37
|
+
1.7. "Larger Work"
|
|
38
|
+
means a work that combines Covered Software with other material, in
|
|
39
|
+
a separate file or files, that is not Covered Software.
|
|
40
|
+
|
|
41
|
+
1.8. "License"
|
|
42
|
+
means this document.
|
|
43
|
+
|
|
44
|
+
1.9. "Licensable"
|
|
45
|
+
means having the right to grant, to the maximum extent possible,
|
|
46
|
+
whether at the time of the initial grant or subsequently, any and
|
|
47
|
+
all of the rights conveyed by this License.
|
|
48
|
+
|
|
49
|
+
1.10. "Modifications"
|
|
50
|
+
means any of the following:
|
|
51
|
+
|
|
52
|
+
(a) any file in Source Code Form that results from an addition to,
|
|
53
|
+
deletion from, or modification of the contents of Covered
|
|
54
|
+
Software; or
|
|
55
|
+
|
|
56
|
+
(b) any new file in Source Code Form that contains any Covered
|
|
57
|
+
Software.
|
|
58
|
+
|
|
59
|
+
1.11. "Patent Claims" of a Contributor
|
|
60
|
+
means any patent claim(s), including without limitation, method,
|
|
61
|
+
process, and apparatus claims, in any patent Licensable by such
|
|
62
|
+
Contributor that would be infringed, but for the grant of the
|
|
63
|
+
License, by the making, using, selling, offering for sale, having
|
|
64
|
+
made, import, or transfer of either its Contributions or its
|
|
65
|
+
Contributor Version.
|
|
66
|
+
|
|
67
|
+
1.12. "Secondary License"
|
|
68
|
+
means either the GNU General Public License, Version 2.0, the GNU
|
|
69
|
+
Lesser General Public License, Version 2.1, the GNU Affero General
|
|
70
|
+
Public License, Version 3.0, or any later versions of those
|
|
71
|
+
licenses.
|
|
72
|
+
|
|
73
|
+
1.13. "Source Code Form"
|
|
74
|
+
means the form of the work preferred for making modifications.
|
|
75
|
+
|
|
76
|
+
1.14. "You" (or "Your")
|
|
77
|
+
means an individual or a legal entity exercising rights under this
|
|
78
|
+
License. For legal entities, "You" includes any entity that
|
|
79
|
+
controls, is controlled by, or is under common control with You. For
|
|
80
|
+
purposes of this definition, "control" means (a) the power, direct
|
|
81
|
+
or indirect, to cause the direction or management of such entity,
|
|
82
|
+
whether by contract or otherwise, or (b) ownership of more than
|
|
83
|
+
fifty percent (50%) of the outstanding shares or beneficial
|
|
84
|
+
ownership of such entity.
|
|
85
|
+
|
|
86
|
+
2. License Grants and Conditions
|
|
87
|
+
--------------------------------
|
|
88
|
+
|
|
89
|
+
2.1. Grants
|
|
90
|
+
|
|
91
|
+
Each Contributor hereby grants You a world-wide, royalty-free,
|
|
92
|
+
non-exclusive license:
|
|
93
|
+
|
|
94
|
+
(a) under intellectual property rights (other than patent or trademark)
|
|
95
|
+
Licensable by such Contributor to use, reproduce, make available,
|
|
96
|
+
modify, display, perform, distribute, and otherwise exploit its
|
|
97
|
+
Contributions, either on an unmodified basis, with Modifications, or
|
|
98
|
+
as part of a Larger Work; and
|
|
99
|
+
|
|
100
|
+
(b) under Patent Claims of such Contributor to make, use, sell, offer
|
|
101
|
+
for sale, have made, import, and otherwise transfer either its
|
|
102
|
+
Contributions or its Contributor Version.
|
|
103
|
+
|
|
104
|
+
2.2. Effective Date
|
|
105
|
+
|
|
106
|
+
The licenses granted in Section 2.1 with respect to any Contribution
|
|
107
|
+
become effective for each Contribution on the date the Contributor first
|
|
108
|
+
distributes such Contribution.
|
|
109
|
+
|
|
110
|
+
2.3. Limitations on Grant Scope
|
|
111
|
+
|
|
112
|
+
The licenses granted in this Section 2 are the only rights granted under
|
|
113
|
+
this License. No additional rights or licenses will be implied from the
|
|
114
|
+
distribution or licensing of Covered Software under this License.
|
|
115
|
+
Notwithstanding Section 2.1(b) above, no patent license is granted by a
|
|
116
|
+
Contributor:
|
|
117
|
+
|
|
118
|
+
(a) for any code that a Contributor has removed from Covered Software;
|
|
119
|
+
or
|
|
120
|
+
|
|
121
|
+
(b) for infringements caused by: (i) Your and any other third party's
|
|
122
|
+
modifications of Covered Software, or (ii) the combination of its
|
|
123
|
+
Contributions with other software (except as part of its Contributor
|
|
124
|
+
Version); or
|
|
125
|
+
|
|
126
|
+
(c) under Patent Claims infringed by Covered Software in the absence of
|
|
127
|
+
its Contributions.
|
|
128
|
+
|
|
129
|
+
This License does not grant any rights in the trademarks, service marks,
|
|
130
|
+
or logos of any Contributor (except as may be necessary to comply with
|
|
131
|
+
the notice requirements in Section 3.4).
|
|
132
|
+
|
|
133
|
+
2.4. Subsequent Licenses
|
|
134
|
+
|
|
135
|
+
No Contributor makes additional grants as a result of Your choice to
|
|
136
|
+
distribute the Covered Software under a subsequent version of this
|
|
137
|
+
License (see Section 10.2) or under the terms of a Secondary License (if
|
|
138
|
+
permitted under the terms of Section 3.3).
|
|
139
|
+
|
|
140
|
+
2.5. Representation
|
|
141
|
+
|
|
142
|
+
Each Contributor represents that the Contributor believes its
|
|
143
|
+
Contributions are its original creation(s) or it has sufficient rights
|
|
144
|
+
to grant the rights to its Contributions conveyed by this License.
|
|
145
|
+
|
|
146
|
+
2.6. Fair Use
|
|
147
|
+
|
|
148
|
+
This License is not intended to limit any rights You have under
|
|
149
|
+
applicable copyright doctrines of fair use, fair dealing, or other
|
|
150
|
+
equivalents.
|
|
151
|
+
|
|
152
|
+
2.7. Conditions
|
|
153
|
+
|
|
154
|
+
Sections 3.1, 3.2, 3.3, and 3.4 are conditions of the licenses granted
|
|
155
|
+
in Section 2.1.
|
|
156
|
+
|
|
157
|
+
3. Responsibilities
|
|
158
|
+
-------------------
|
|
159
|
+
|
|
160
|
+
3.1. Distribution of Source Form
|
|
161
|
+
|
|
162
|
+
All distribution of Covered Software in Source Code Form, including any
|
|
163
|
+
Modifications that You create or to which You contribute, must be under
|
|
164
|
+
the terms of this License. You must inform recipients that the Source
|
|
165
|
+
Code Form of the Covered Software is governed by the terms of this
|
|
166
|
+
License, and how they can obtain a copy of this License. You may not
|
|
167
|
+
attempt to alter or restrict the recipients' rights in the Source Code
|
|
168
|
+
Form.
|
|
169
|
+
|
|
170
|
+
3.2. Distribution of Executable Form
|
|
171
|
+
|
|
172
|
+
If You distribute Covered Software in Executable Form then:
|
|
173
|
+
|
|
174
|
+
(a) such Covered Software must also be made available in Source Code
|
|
175
|
+
Form, as described in Section 3.1, and You must inform recipients of
|
|
176
|
+
the Executable Form how they can obtain a copy of such Source Code
|
|
177
|
+
Form by reasonable means in a timely manner, at a charge no more
|
|
178
|
+
than the cost of distribution to the recipient; and
|
|
179
|
+
|
|
180
|
+
(b) You may distribute such Executable Form under the terms of this
|
|
181
|
+
License, or sublicense it under different terms, provided that the
|
|
182
|
+
license for the Executable Form does not attempt to limit or alter
|
|
183
|
+
the recipients' rights in the Source Code Form under this License.
|
|
184
|
+
|
|
185
|
+
3.3. Distribution of a Larger Work
|
|
186
|
+
|
|
187
|
+
You may create and distribute a Larger Work under terms of Your choice,
|
|
188
|
+
provided that You also comply with the requirements of this License for
|
|
189
|
+
the Covered Software. If the Larger Work is a combination of Covered
|
|
190
|
+
Software with a work governed by one or more Secondary Licenses, and the
|
|
191
|
+
Covered Software is not Incompatible With Secondary Licenses, this
|
|
192
|
+
License permits You to additionally distribute such Covered Software
|
|
193
|
+
under the terms of such Secondary License(s), so that the recipient of
|
|
194
|
+
the Larger Work may, at their option, further distribute the Covered
|
|
195
|
+
Software under the terms of either this License or such Secondary
|
|
196
|
+
License(s).
|
|
197
|
+
|
|
198
|
+
3.4. Notices
|
|
199
|
+
|
|
200
|
+
You may not remove or alter the substance of any license notices
|
|
201
|
+
(including copyright notices, patent notices, disclaimers of warranty,
|
|
202
|
+
or limitations of liability) contained within the Source Code Form of
|
|
203
|
+
the Covered Software, except that You may alter any license notices to
|
|
204
|
+
the extent required to remedy known factual inaccuracies.
|
|
205
|
+
|
|
206
|
+
3.5. Application of Additional Terms
|
|
207
|
+
|
|
208
|
+
You may choose to offer, and to charge a fee for, warranty, support,
|
|
209
|
+
indemnity or liability obligations to one or more recipients of Covered
|
|
210
|
+
Software. However, You may do so only on Your own behalf, and not on
|
|
211
|
+
behalf of any Contributor. You must make it absolutely clear that any
|
|
212
|
+
such warranty, support, indemnity, or liability obligation is offered by
|
|
213
|
+
You alone, and You hereby agree to indemnify every Contributor for any
|
|
214
|
+
liability incurred by such Contributor as a result of warranty, support,
|
|
215
|
+
indemnity or liability terms You offer. You may include additional
|
|
216
|
+
disclaimers of warranty and limitations of liability specific to any
|
|
217
|
+
jurisdiction.
|
|
218
|
+
|
|
219
|
+
4. Inability to Comply Due to Statute or Regulation
|
|
220
|
+
---------------------------------------------------
|
|
221
|
+
|
|
222
|
+
If it is impossible for You to comply with any of the terms of this
|
|
223
|
+
License with respect to some or all of the Covered Software due to
|
|
224
|
+
statute, judicial order, or regulation then You must: (a) comply with
|
|
225
|
+
the terms of this License to the maximum extent possible; and (b)
|
|
226
|
+
describe the limitations and the code they affect. Such description must
|
|
227
|
+
be placed in a text file included with all distributions of the Covered
|
|
228
|
+
Software under this License. Except to the extent prohibited by statute
|
|
229
|
+
or regulation, such description must be sufficiently detailed for a
|
|
230
|
+
recipient of ordinary skill to be able to understand it.
|
|
231
|
+
|
|
232
|
+
5. Termination
|
|
233
|
+
--------------
|
|
234
|
+
|
|
235
|
+
5.1. The rights granted under this License will terminate automatically
|
|
236
|
+
if You fail to comply with any of its terms. However, if You become
|
|
237
|
+
compliant, then the rights granted under this License from a particular
|
|
238
|
+
Contributor are reinstated (a) provisionally, unless and until such
|
|
239
|
+
Contributor explicitly and finally terminates Your grants, and (b) on an
|
|
240
|
+
ongoing basis, if such Contributor fails to notify You of the
|
|
241
|
+
non-compliance by some reasonable means prior to 60 days after You have
|
|
242
|
+
come back into compliance. Moreover, Your grants from a particular
|
|
243
|
+
Contributor are reinstated on an ongoing basis if such Contributor
|
|
244
|
+
notifies You of the non-compliance by some reasonable means, this is the
|
|
245
|
+
first time You have received notice of non-compliance with this License
|
|
246
|
+
from such Contributor, and You become compliant prior to 30 days after
|
|
247
|
+
Your receipt of the notice.
|
|
248
|
+
|
|
249
|
+
5.2. If You initiate litigation against any entity by asserting a patent
|
|
250
|
+
infringement claim (excluding declaratory judgment actions,
|
|
251
|
+
counter-claims, and cross-claims) alleging that a Contributor Version
|
|
252
|
+
directly or indirectly infringes any patent, then the rights granted to
|
|
253
|
+
You by any and all Contributors for the Covered Software under Section
|
|
254
|
+
2.1 of this License shall terminate.
|
|
255
|
+
|
|
256
|
+
5.3. In the event of termination under Sections 5.1 or 5.2 above, all
|
|
257
|
+
end user license agreements (excluding distributors and resellers) which
|
|
258
|
+
have been validly granted by You or Your distributors under this License
|
|
259
|
+
prior to termination shall survive termination.
|
|
260
|
+
|
|
261
|
+
************************************************************************
|
|
262
|
+
* *
|
|
263
|
+
* 6. Disclaimer of Warranty *
|
|
264
|
+
* ------------------------- *
|
|
265
|
+
* *
|
|
266
|
+
* Covered Software is provided under this License on an "as is" *
|
|
267
|
+
* basis, without warranty of any kind, either expressed, implied, or *
|
|
268
|
+
* statutory, including, without limitation, warranties that the *
|
|
269
|
+
* Covered Software is free of defects, merchantable, fit for a *
|
|
270
|
+
* particular purpose or non-infringing. The entire risk as to the *
|
|
271
|
+
* quality and performance of the Covered Software is with You. *
|
|
272
|
+
* Should any Covered Software prove defective in any respect, You *
|
|
273
|
+
* (not any Contributor) assume the cost of any necessary servicing, *
|
|
274
|
+
* repair, or correction. This disclaimer of warranty constitutes an *
|
|
275
|
+
* essential part of this License. No use of any Covered Software is *
|
|
276
|
+
* authorized under this License except under this disclaimer. *
|
|
277
|
+
* *
|
|
278
|
+
************************************************************************
|
|
279
|
+
|
|
280
|
+
************************************************************************
|
|
281
|
+
* *
|
|
282
|
+
* 7. Limitation of Liability *
|
|
283
|
+
* -------------------------- *
|
|
284
|
+
* *
|
|
285
|
+
* Under no circumstances and under no legal theory, whether tort *
|
|
286
|
+
* (including negligence), contract, or otherwise, shall any *
|
|
287
|
+
* Contributor, or anyone who distributes Covered Software as *
|
|
288
|
+
* permitted above, be liable to You for any direct, indirect, *
|
|
289
|
+
* special, incidental, or consequential damages of any character *
|
|
290
|
+
* including, without limitation, damages for lost profits, loss of *
|
|
291
|
+
* goodwill, work stoppage, computer failure or malfunction, or any *
|
|
292
|
+
* and all other commercial damages or losses, even if such party *
|
|
293
|
+
* shall have been informed of the possibility of such damages. This *
|
|
294
|
+
* limitation of liability shall not apply to liability for death or *
|
|
295
|
+
* personal injury resulting from such party's negligence to the *
|
|
296
|
+
* extent applicable law prohibits such limitation. Some *
|
|
297
|
+
* jurisdictions do not allow the exclusion or limitation of *
|
|
298
|
+
* incidental or consequential damages, so this exclusion and *
|
|
299
|
+
* limitation may not apply to You. *
|
|
300
|
+
* *
|
|
301
|
+
************************************************************************
|
|
302
|
+
|
|
303
|
+
8. Litigation
|
|
304
|
+
-------------
|
|
305
|
+
|
|
306
|
+
Any litigation relating to this License may be brought only in the
|
|
307
|
+
courts of a jurisdiction where the defendant maintains its principal
|
|
308
|
+
place of business and such litigation shall be governed by laws of that
|
|
309
|
+
jurisdiction, without reference to its conflict-of-law provisions.
|
|
310
|
+
Nothing in this Section shall prevent a party's ability to bring
|
|
311
|
+
cross-claims or counter-claims.
|
|
312
|
+
|
|
313
|
+
9. Miscellaneous
|
|
314
|
+
----------------
|
|
315
|
+
|
|
316
|
+
This License represents the complete agreement concerning the subject
|
|
317
|
+
matter hereof. If any provision of this License is held to be
|
|
318
|
+
unenforceable, such provision shall be reformed only to the extent
|
|
319
|
+
necessary to make it enforceable. Any law or regulation which provides
|
|
320
|
+
that the language of a contract shall be construed against the drafter
|
|
321
|
+
shall not be used to construe this License against a Contributor.
|
|
322
|
+
|
|
323
|
+
10. Versions of the License
|
|
324
|
+
---------------------------
|
|
325
|
+
|
|
326
|
+
10.1. New Versions
|
|
327
|
+
|
|
328
|
+
Mozilla Foundation is the license steward. Except as provided in Section
|
|
329
|
+
10.3, no one other than the license steward has the right to modify or
|
|
330
|
+
publish new versions of this License. Each version will be given a
|
|
331
|
+
distinguishing version number.
|
|
332
|
+
|
|
333
|
+
10.2. Effect of New Versions
|
|
334
|
+
|
|
335
|
+
You may distribute the Covered Software under the terms of the version
|
|
336
|
+
of the License under which You originally received the Covered Software,
|
|
337
|
+
or under the terms of any subsequent version published by the license
|
|
338
|
+
steward.
|
|
339
|
+
|
|
340
|
+
10.3. Modified Versions
|
|
341
|
+
|
|
342
|
+
If you create software not governed by this License, and you want to
|
|
343
|
+
create a new license for such software, you may create and use a
|
|
344
|
+
modified version of this License if you rename the license and remove
|
|
345
|
+
any references to the name of the license steward (except to note that
|
|
346
|
+
such modified license differs from this License).
|
|
347
|
+
|
|
348
|
+
10.4. Distributing Source Code Form that is Incompatible With Secondary
|
|
349
|
+
Licenses
|
|
350
|
+
|
|
351
|
+
If You choose to distribute Source Code Form that is Incompatible With
|
|
352
|
+
Secondary Licenses under the terms of this version of the License, the
|
|
353
|
+
notice described in Exhibit B of this License must be attached.
|
|
354
|
+
|
|
355
|
+
Exhibit A - Source Code Form License Notice
|
|
356
|
+
-------------------------------------------
|
|
357
|
+
|
|
358
|
+
This Source Code Form is subject to the terms of the Mozilla Public
|
|
359
|
+
License, v. 2.0. If a copy of the MPL was not distributed with this
|
|
360
|
+
file, You can obtain one at https://mozilla.org/MPL/2.0/.
|
|
361
|
+
|
|
362
|
+
If it is not possible or desirable to put the notice in a particular
|
|
363
|
+
file, then You may include the notice in a location (such as a LICENSE
|
|
364
|
+
file in a relevant directory) where a recipient would be likely to look
|
|
365
|
+
for such a notice.
|
|
366
|
+
|
|
367
|
+
You may add additional accurate notices of copyright ownership.
|
|
368
|
+
|
|
369
|
+
Exhibit B - "Incompatible With Secondary Licenses" Notice
|
|
370
|
+
---------------------------------------------------------
|
|
371
|
+
|
|
372
|
+
This Source Code Form is "Incompatible With Secondary Licenses", as
|
|
373
|
+
defined by the Mozilla Public License, v. 2.0.
|
|
@@ -0,0 +1,103 @@
|
|
|
1
|
+
Metadata-Version: 2.4
|
|
2
|
+
Name: hmath-python
|
|
3
|
+
Version: 1.0.0
|
|
4
|
+
Summary: Python bindings for the Haskell HMath numerical library
|
|
5
|
+
License-Expression: MPL-2.0
|
|
6
|
+
Project-URL: Homepage, https://github.com/wavefnd/hmath
|
|
7
|
+
Project-URL: Source, https://github.com/wavefnd/hmath
|
|
8
|
+
Project-URL: Issues, https://github.com/wavefnd/hmath/issues
|
|
9
|
+
Keywords: mathematics,haskell,ffi,tensor
|
|
10
|
+
Classifier: Development Status :: 5 - Production/Stable
|
|
11
|
+
Classifier: Operating System :: POSIX :: Linux
|
|
12
|
+
Classifier: Programming Language :: Python :: 3
|
|
13
|
+
Classifier: Topic :: Scientific/Engineering :: Mathematics
|
|
14
|
+
Requires-Python: >=3.9
|
|
15
|
+
Description-Content-Type: text/markdown
|
|
16
|
+
License-File: LICENSE
|
|
17
|
+
Provides-Extra: cffi
|
|
18
|
+
Requires-Dist: cffi>=1.15; extra == "cffi"
|
|
19
|
+
Provides-Extra: numpy
|
|
20
|
+
Requires-Dist: numpy>=1.22; extra == "numpy"
|
|
21
|
+
Dynamic: license-file
|
|
22
|
+
|
|
23
|
+
# HMath
|
|
24
|
+
|
|
25
|
+
HMath is a mathematics and numerical computing library with Python bindings, written in Haskell. It also provides native Haskell, C, and Wave interfaces.
|
|
26
|
+
|
|
27
|
+
The 1.0.0 release targets Linux x86-64. The Python distribution is named **hmath-python** and is imported as `hmath`.
|
|
28
|
+
|
|
29
|
+
## Python quick start
|
|
30
|
+
|
|
31
|
+
Install the Python package on Linux x86-64:
|
|
32
|
+
|
|
33
|
+
```sh
|
|
34
|
+
pip install hmath-python
|
|
35
|
+
```
|
|
36
|
+
|
|
37
|
+
```python
|
|
38
|
+
from hmath import Tensor, add
|
|
39
|
+
|
|
40
|
+
print(add(2, 3)) # 5.0
|
|
41
|
+
|
|
42
|
+
with Tensor.matrix([1, 2, 3, 4], 2, 2) as matrix:
|
|
43
|
+
with matrix @ matrix as product:
|
|
44
|
+
print(product.values()) # [7.0, 10.0, 15.0, 22.0]
|
|
45
|
+
```
|
|
46
|
+
|
|
47
|
+
Use `with` to release native resources automatically. Invalid operations raise `HMathError`. Install `hmath-python[numpy,cffi]` for optional NumPy and cffi support.
|
|
48
|
+
|
|
49
|
+
## Features
|
|
50
|
+
|
|
51
|
+
- **Core math:** checked scalar operations, numeric types, constants, units, number theory, and polynomials.
|
|
52
|
+
- **Linear algebra:** dense and sparse matrices, decompositions, iterative solvers, and CPU tensors.
|
|
53
|
+
- **Scientific computing:** differentiation, integration, ODE/PDE solvers, optimization, interpolation, FFTs, and signal processing.
|
|
54
|
+
- **Statistics and learning:** reproducible sampling, probability distributions, inference, automatic differentiation, and machine learning primitives.
|
|
55
|
+
- **Discrete and symbolic math:** combinatorics, graphs, numerical data structures, symbolic calculus, and a bounded expression evaluator.
|
|
56
|
+
- **Applied math:** complex numbers, rotations, computational geometry, time series, and numerical image processing.
|
|
57
|
+
- **Native interface:** checked scalar, vector, matrix, and tensor operations; explicit ownership and errors; Python/NumPy and Wave bindings.
|
|
58
|
+
- **Interchange:** typed binary and text arrays, CSV/TSV, Matrix Market, NPY/NPZ, and CPU DLPack.
|
|
59
|
+
|
|
60
|
+
The Haskell library contains the full numerical API; Python currently exposes scalar addition and binary64 vector, matrix, and tensor construction, arithmetic, and reductions.
|
|
61
|
+
|
|
62
|
+
## Build from source
|
|
63
|
+
|
|
64
|
+
Requires GHC, Cabal, a C11 compiler, Meson >= 1.3, Ninja, and Python 3. Zlib enables compressed interchange. BLAS/LAPACKE are optional runtime backends; the full spectral test suite needs LAPACKE.
|
|
65
|
+
|
|
66
|
+
```sh
|
|
67
|
+
git clone https://github.com/wavefnd/hmath.git
|
|
68
|
+
cd hmath
|
|
69
|
+
meson setup builddir
|
|
70
|
+
meson compile -C builddir
|
|
71
|
+
meson test -C builddir --print-errorlogs
|
|
72
|
+
```
|
|
73
|
+
|
|
74
|
+
This builds `builddir/libhmath.so` and runs the Haskell, C, and Python tests. Wave tests also run when `wavec` is available.
|
|
75
|
+
|
|
76
|
+
## Haskell
|
|
77
|
+
|
|
78
|
+
Save the following as `Main.hs` in the checkout. Run `cabal build lib:hmath`, then `cabal exec -- runghc -package=hmath Main.hs`:
|
|
79
|
+
|
|
80
|
+
```haskell
|
|
81
|
+
import HMath
|
|
82
|
+
|
|
83
|
+
main :: IO ()
|
|
84
|
+
main = do
|
|
85
|
+
print (dotChecked [1, 2, 3 :: Int] [4, 5, 6]) -- Right 32
|
|
86
|
+
print $ do
|
|
87
|
+
matrix <- fromListChecked 2 2 [1, 2, 3, 4 :: Int]
|
|
88
|
+
pure (toList (transposeMatrix matrix)) -- Right [1,3,2,4]
|
|
89
|
+
```
|
|
90
|
+
|
|
91
|
+
Checked operations return `Either MathError` for invalid shapes or lengths. More examples are available in [examples/haskell](examples/haskell).
|
|
92
|
+
|
|
93
|
+
## C and Wave
|
|
94
|
+
|
|
95
|
+
The [C header](include/hmath.h) defines errors, ownership, and runtime lifetime. See the [C example](examples/c/main.c) and [Wave bindings](bindings/wave) for usage.
|
|
96
|
+
|
|
97
|
+
## Contributing
|
|
98
|
+
|
|
99
|
+
Report bugs or discuss changes through [GitHub issues](https://github.com/wavefnd/hmath/issues). Submit pull requests from a fork against `master`; see [ROADMAP.md](ROADMAP.md) for planned work.
|
|
100
|
+
|
|
101
|
+
## License
|
|
102
|
+
|
|
103
|
+
[Mozilla Public License 2.0](LICENSE).
|
|
@@ -0,0 +1,81 @@
|
|
|
1
|
+
# HMath
|
|
2
|
+
|
|
3
|
+
HMath is a mathematics and numerical computing library with Python bindings, written in Haskell. It also provides native Haskell, C, and Wave interfaces.
|
|
4
|
+
|
|
5
|
+
The 1.0.0 release targets Linux x86-64. The Python distribution is named **hmath-python** and is imported as `hmath`.
|
|
6
|
+
|
|
7
|
+
## Python quick start
|
|
8
|
+
|
|
9
|
+
Install the Python package on Linux x86-64:
|
|
10
|
+
|
|
11
|
+
```sh
|
|
12
|
+
pip install hmath-python
|
|
13
|
+
```
|
|
14
|
+
|
|
15
|
+
```python
|
|
16
|
+
from hmath import Tensor, add
|
|
17
|
+
|
|
18
|
+
print(add(2, 3)) # 5.0
|
|
19
|
+
|
|
20
|
+
with Tensor.matrix([1, 2, 3, 4], 2, 2) as matrix:
|
|
21
|
+
with matrix @ matrix as product:
|
|
22
|
+
print(product.values()) # [7.0, 10.0, 15.0, 22.0]
|
|
23
|
+
```
|
|
24
|
+
|
|
25
|
+
Use `with` to release native resources automatically. Invalid operations raise `HMathError`. Install `hmath-python[numpy,cffi]` for optional NumPy and cffi support.
|
|
26
|
+
|
|
27
|
+
## Features
|
|
28
|
+
|
|
29
|
+
- **Core math:** checked scalar operations, numeric types, constants, units, number theory, and polynomials.
|
|
30
|
+
- **Linear algebra:** dense and sparse matrices, decompositions, iterative solvers, and CPU tensors.
|
|
31
|
+
- **Scientific computing:** differentiation, integration, ODE/PDE solvers, optimization, interpolation, FFTs, and signal processing.
|
|
32
|
+
- **Statistics and learning:** reproducible sampling, probability distributions, inference, automatic differentiation, and machine learning primitives.
|
|
33
|
+
- **Discrete and symbolic math:** combinatorics, graphs, numerical data structures, symbolic calculus, and a bounded expression evaluator.
|
|
34
|
+
- **Applied math:** complex numbers, rotations, computational geometry, time series, and numerical image processing.
|
|
35
|
+
- **Native interface:** checked scalar, vector, matrix, and tensor operations; explicit ownership and errors; Python/NumPy and Wave bindings.
|
|
36
|
+
- **Interchange:** typed binary and text arrays, CSV/TSV, Matrix Market, NPY/NPZ, and CPU DLPack.
|
|
37
|
+
|
|
38
|
+
The Haskell library contains the full numerical API; Python currently exposes scalar addition and binary64 vector, matrix, and tensor construction, arithmetic, and reductions.
|
|
39
|
+
|
|
40
|
+
## Build from source
|
|
41
|
+
|
|
42
|
+
Requires GHC, Cabal, a C11 compiler, Meson >= 1.3, Ninja, and Python 3. Zlib enables compressed interchange. BLAS/LAPACKE are optional runtime backends; the full spectral test suite needs LAPACKE.
|
|
43
|
+
|
|
44
|
+
```sh
|
|
45
|
+
git clone https://github.com/wavefnd/hmath.git
|
|
46
|
+
cd hmath
|
|
47
|
+
meson setup builddir
|
|
48
|
+
meson compile -C builddir
|
|
49
|
+
meson test -C builddir --print-errorlogs
|
|
50
|
+
```
|
|
51
|
+
|
|
52
|
+
This builds `builddir/libhmath.so` and runs the Haskell, C, and Python tests. Wave tests also run when `wavec` is available.
|
|
53
|
+
|
|
54
|
+
## Haskell
|
|
55
|
+
|
|
56
|
+
Save the following as `Main.hs` in the checkout. Run `cabal build lib:hmath`, then `cabal exec -- runghc -package=hmath Main.hs`:
|
|
57
|
+
|
|
58
|
+
```haskell
|
|
59
|
+
import HMath
|
|
60
|
+
|
|
61
|
+
main :: IO ()
|
|
62
|
+
main = do
|
|
63
|
+
print (dotChecked [1, 2, 3 :: Int] [4, 5, 6]) -- Right 32
|
|
64
|
+
print $ do
|
|
65
|
+
matrix <- fromListChecked 2 2 [1, 2, 3, 4 :: Int]
|
|
66
|
+
pure (toList (transposeMatrix matrix)) -- Right [1,3,2,4]
|
|
67
|
+
```
|
|
68
|
+
|
|
69
|
+
Checked operations return `Either MathError` for invalid shapes or lengths. More examples are available in [examples/haskell](examples/haskell).
|
|
70
|
+
|
|
71
|
+
## C and Wave
|
|
72
|
+
|
|
73
|
+
The [C header](include/hmath.h) defines errors, ownership, and runtime lifetime. See the [C example](examples/c/main.c) and [Wave bindings](bindings/wave) for usage.
|
|
74
|
+
|
|
75
|
+
## Contributing
|
|
76
|
+
|
|
77
|
+
Report bugs or discuss changes through [GitHub issues](https://github.com/wavefnd/hmath/issues). Submit pull requests from a fork against `master`; see [ROADMAP.md](ROADMAP.md) for planned work.
|
|
78
|
+
|
|
79
|
+
## License
|
|
80
|
+
|
|
81
|
+
[Mozilla Public License 2.0](LICENSE).
|
|
@@ -0,0 +1,283 @@
|
|
|
1
|
+
# SPDX-License-Identifier: MPL-2.0
|
|
2
|
+
"""Binary64 HMath ABI 1. Set HMATH_LIBRARY to the installed shared library.
|
|
3
|
+
|
|
4
|
+
Tensor owns one native reference. close()/with provide deterministic cleanup;
|
|
5
|
+
__del__ is a fallback. Arrays are copied by default. from_numpy(copy=False)
|
|
6
|
+
requires readonly, native-endian, aligned float64 storage and nonnegative
|
|
7
|
+
strides; the caller must not mutate aliases until all derived views die.
|
|
8
|
+
Native storage views retain independent references, so closing a Tensor does
|
|
9
|
+
not invalidate an existing NumPy view. The process owner shuts down explicitly.
|
|
10
|
+
"""
|
|
11
|
+
|
|
12
|
+
import ctypes as C
|
|
13
|
+
import ctypes.util
|
|
14
|
+
import os
|
|
15
|
+
import threading
|
|
16
|
+
import weakref
|
|
17
|
+
from itertools import islice
|
|
18
|
+
|
|
19
|
+
_lib = C.CDLL(
|
|
20
|
+
os.environ.get("HMATH_LIBRARY")
|
|
21
|
+
or ctypes.util.find_library("hmath")
|
|
22
|
+
or "libhmath.so"
|
|
23
|
+
)
|
|
24
|
+
P, Z, D = C.c_void_p, C.c_size_t, C.c_double
|
|
25
|
+
|
|
26
|
+
|
|
27
|
+
def _sig(name, result, *args):
|
|
28
|
+
fn = getattr(_lib, name)
|
|
29
|
+
fn.restype, fn.argtypes = result, list(args)
|
|
30
|
+
return fn
|
|
31
|
+
|
|
32
|
+
|
|
33
|
+
try:
|
|
34
|
+
_sig("hmath_abi_version", C.c_uint32)
|
|
35
|
+
if _lib.hmath_abi_version() != 1:
|
|
36
|
+
raise RuntimeError("HMath binding requires ABI 1")
|
|
37
|
+
except AttributeError as error:
|
|
38
|
+
raise RuntimeError(
|
|
39
|
+
"HMath binding requires ABI 1; library has no version query"
|
|
40
|
+
) from error
|
|
41
|
+
_sig("hmath_init", C.c_int)
|
|
42
|
+
_sig("hmath_last_error_message", C.c_char_p)
|
|
43
|
+
_sig("hmath_runtime_shutdown", C.c_int)
|
|
44
|
+
_sig("hmath_tensor_copy", C.c_int, C.POINTER(D), Z, C.POINTER(Z), Z, C.POINTER(P))
|
|
45
|
+
_sig(
|
|
46
|
+
"hmath_tensor_borrow",
|
|
47
|
+
C.c_int,
|
|
48
|
+
C.POINTER(D),
|
|
49
|
+
Z,
|
|
50
|
+
C.POINTER(Z),
|
|
51
|
+
C.POINTER(Z),
|
|
52
|
+
Z,
|
|
53
|
+
P,
|
|
54
|
+
P,
|
|
55
|
+
C.POINTER(P),
|
|
56
|
+
)
|
|
57
|
+
_sig("hmath_tensor_retain", None, P)
|
|
58
|
+
_sig("hmath_tensor_release", None, P)
|
|
59
|
+
_sig("hmath_tensor_count", Z, P)
|
|
60
|
+
_sig("hmath_tensor_rank", Z, P)
|
|
61
|
+
_sig("hmath_tensor_metadata", C.c_int, P, C.POINTER(Z), C.POINTER(Z), Z)
|
|
62
|
+
_sig("hmath_tensor_copy_to", C.c_int, P, C.POINTER(D), Z)
|
|
63
|
+
_sig("hmath_tensor_data", C.c_int, P, C.POINTER(C.POINTER(D)), C.POINTER(Z))
|
|
64
|
+
_sig("hmath_scalar_add", C.c_int, P, D, D, C.POINTER(D))
|
|
65
|
+
_sig("hmath_tensor_sum", C.c_int, P, P, C.POINTER(D))
|
|
66
|
+
_sig("hmath_tensor_dot", C.c_int, P, P, P, C.POINTER(D))
|
|
67
|
+
for _name in ("hmath_tensor_add", "hmath_tensor_matmul"):
|
|
68
|
+
_sig(_name, C.c_int, P, P, P, C.POINTER(P))
|
|
69
|
+
_sig("hmath_tensor_reshape", C.c_int, P, C.POINTER(Z), Z, C.POINTER(P))
|
|
70
|
+
|
|
71
|
+
|
|
72
|
+
class HMathError(RuntimeError):
|
|
73
|
+
def __init__(self, status):
|
|
74
|
+
self.status = status
|
|
75
|
+
super().__init__(_lib.hmath_last_error_message().decode("utf-8", "replace"))
|
|
76
|
+
|
|
77
|
+
|
|
78
|
+
def _check(status):
|
|
79
|
+
if status:
|
|
80
|
+
raise HMathError(status)
|
|
81
|
+
|
|
82
|
+
|
|
83
|
+
def initialize():
|
|
84
|
+
if _lib.hmath_abi_version() != 1:
|
|
85
|
+
raise RuntimeError("HMath binding requires ABI 1")
|
|
86
|
+
if _lib.hmath_init():
|
|
87
|
+
raise RuntimeError("HMath runtime was permanently shut down")
|
|
88
|
+
|
|
89
|
+
|
|
90
|
+
def shutdown():
|
|
91
|
+
"""Process-owner operation; join all numerical users before calling."""
|
|
92
|
+
_check(_lib.hmath_runtime_shutdown())
|
|
93
|
+
|
|
94
|
+
|
|
95
|
+
def add(a, b):
|
|
96
|
+
initialize()
|
|
97
|
+
result = D()
|
|
98
|
+
_check(_lib.hmath_scalar_add(None, a, b, C.byref(result)))
|
|
99
|
+
return result.value
|
|
100
|
+
|
|
101
|
+
|
|
102
|
+
def _shape(shape):
|
|
103
|
+
shape = tuple(shape)
|
|
104
|
+
if len(shape) > 8 or any(type(n) is not int or n < 0 or n > 1048576 for n in shape):
|
|
105
|
+
raise ValueError("shape requires at most eight nonnegative bounded integers")
|
|
106
|
+
return shape, (Z * len(shape))(*shape)
|
|
107
|
+
|
|
108
|
+
|
|
109
|
+
class Tensor:
|
|
110
|
+
def __init__(self, values, shape=None):
|
|
111
|
+
initialize()
|
|
112
|
+
values = list(islice(values, 1048577))
|
|
113
|
+
self._lock, self._ptr, self._owner = threading.RLock(), P(), None
|
|
114
|
+
shape, dims = _shape((len(values),) if shape is None else shape)
|
|
115
|
+
if len(values) > 1048576:
|
|
116
|
+
raise ValueError("tensor exceeds element budget")
|
|
117
|
+
data = (D * len(values))(*values)
|
|
118
|
+
_check(
|
|
119
|
+
_lib.hmath_tensor_copy(
|
|
120
|
+
data, len(values), dims, len(shape), C.byref(self._ptr)
|
|
121
|
+
)
|
|
122
|
+
)
|
|
123
|
+
|
|
124
|
+
@classmethod
|
|
125
|
+
def _adopt(cls, ptr, owner=None):
|
|
126
|
+
result = object.__new__(cls)
|
|
127
|
+
result._lock, result._ptr, result._owner = threading.RLock(), ptr, owner
|
|
128
|
+
return result
|
|
129
|
+
|
|
130
|
+
@classmethod
|
|
131
|
+
def scalar(cls, value):
|
|
132
|
+
return cls([value], ())
|
|
133
|
+
|
|
134
|
+
@classmethod
|
|
135
|
+
def matrix(cls, values, rows, columns):
|
|
136
|
+
return cls(values, (rows, columns))
|
|
137
|
+
|
|
138
|
+
@classmethod
|
|
139
|
+
def from_numpy(cls, array, copy=True):
|
|
140
|
+
import numpy as np
|
|
141
|
+
|
|
142
|
+
array = np.asarray(array)
|
|
143
|
+
if array.size > 1048576 or array.ndim > 8:
|
|
144
|
+
raise ValueError("tensor exceeds element or rank budget")
|
|
145
|
+
if copy:
|
|
146
|
+
array = np.asarray(array, dtype=np.float64)
|
|
147
|
+
return cls(array.ravel(order="C"), array.shape)
|
|
148
|
+
initialize()
|
|
149
|
+
if (
|
|
150
|
+
array.dtype != np.dtype("float64")
|
|
151
|
+
or not array.dtype.isnative
|
|
152
|
+
or not array.flags.aligned
|
|
153
|
+
or array.flags.writeable
|
|
154
|
+
):
|
|
155
|
+
raise ValueError("zero-copy requires readonly aligned native float64")
|
|
156
|
+
shape, dims = _shape(array.shape)
|
|
157
|
+
if any(s < 0 or s % 8 for s in array.strides):
|
|
158
|
+
raise ValueError("zero-copy requires nonnegative whole-element strides")
|
|
159
|
+
strides = (Z * array.ndim)(*(s // 8 for s in array.strides))
|
|
160
|
+
span = (
|
|
161
|
+
0
|
|
162
|
+
if array.size == 0
|
|
163
|
+
else 1 + sum((n - 1) * (s // 8) for n, s in zip(shape, array.strides))
|
|
164
|
+
)
|
|
165
|
+
ptr = P()
|
|
166
|
+
_check(
|
|
167
|
+
_lib.hmath_tensor_borrow(
|
|
168
|
+
array.ctypes.data_as(C.POINTER(D)),
|
|
169
|
+
span,
|
|
170
|
+
dims,
|
|
171
|
+
strides,
|
|
172
|
+
array.ndim,
|
|
173
|
+
None,
|
|
174
|
+
None,
|
|
175
|
+
C.byref(ptr),
|
|
176
|
+
)
|
|
177
|
+
)
|
|
178
|
+
return cls._adopt(ptr, array)
|
|
179
|
+
|
|
180
|
+
def _live(self):
|
|
181
|
+
if not self._ptr:
|
|
182
|
+
raise RuntimeError("Tensor is closed")
|
|
183
|
+
return self._ptr
|
|
184
|
+
|
|
185
|
+
def close(self):
|
|
186
|
+
lock = getattr(self, "_lock", None)
|
|
187
|
+
if lock is not None:
|
|
188
|
+
with lock:
|
|
189
|
+
if self._ptr:
|
|
190
|
+
_lib.hmath_tensor_release(self._ptr)
|
|
191
|
+
self._ptr = P()
|
|
192
|
+
self._owner = None
|
|
193
|
+
|
|
194
|
+
def __del__(self):
|
|
195
|
+
self.close()
|
|
196
|
+
|
|
197
|
+
def __enter__(self):
|
|
198
|
+
self._live()
|
|
199
|
+
return self
|
|
200
|
+
|
|
201
|
+
def __exit__(self, *args):
|
|
202
|
+
self.close()
|
|
203
|
+
|
|
204
|
+
@property
|
|
205
|
+
def shape(self):
|
|
206
|
+
with self._lock:
|
|
207
|
+
rank = _lib.hmath_tensor_rank(self._live())
|
|
208
|
+
dims, steps = (Z * rank)(), (Z * rank)()
|
|
209
|
+
_check(_lib.hmath_tensor_metadata(self._ptr, dims, steps, rank))
|
|
210
|
+
return tuple(dims)
|
|
211
|
+
|
|
212
|
+
def values(self):
|
|
213
|
+
with self._lock:
|
|
214
|
+
count = _lib.hmath_tensor_count(self._live())
|
|
215
|
+
result = (D * count)()
|
|
216
|
+
_check(_lib.hmath_tensor_copy_to(self._ptr, result, count))
|
|
217
|
+
return list(result)
|
|
218
|
+
|
|
219
|
+
def numpy(self, copy=True):
|
|
220
|
+
import numpy as np
|
|
221
|
+
|
|
222
|
+
with self._lock:
|
|
223
|
+
if copy:
|
|
224
|
+
return np.array(self.values()).reshape(self.shape)
|
|
225
|
+
rank = _lib.hmath_tensor_rank(self._live())
|
|
226
|
+
dims, steps = (Z * rank)(), (Z * rank)()
|
|
227
|
+
_check(_lib.hmath_tensor_metadata(self._ptr, dims, steps, rank))
|
|
228
|
+
data, count = C.POINTER(D)(), Z()
|
|
229
|
+
_check(_lib.hmath_tensor_data(self._ptr, C.byref(data), C.byref(count)))
|
|
230
|
+
if count.value == 0:
|
|
231
|
+
result = np.empty(tuple(dims))
|
|
232
|
+
result.setflags(write=False)
|
|
233
|
+
return result
|
|
234
|
+
buffer = (D * count.value).from_address(C.addressof(data.contents))
|
|
235
|
+
result = np.ndarray(
|
|
236
|
+
tuple(dims),
|
|
237
|
+
dtype=np.float64,
|
|
238
|
+
buffer=memoryview(buffer).toreadonly(),
|
|
239
|
+
strides=tuple(s * 8 for s in steps),
|
|
240
|
+
)
|
|
241
|
+
# Keep ownership on the actual buffer shared by every NumPy base
|
|
242
|
+
# chain. A finalizer on only the first ndarray would release too
|
|
243
|
+
# early when another array is created from view.base.
|
|
244
|
+
retained, owner = P(self._ptr.value), self._owner
|
|
245
|
+
_lib.hmath_tensor_retain(retained)
|
|
246
|
+
weakref.finalize(
|
|
247
|
+
buffer, lambda p, keep: _lib.hmath_tensor_release(p), retained, owner
|
|
248
|
+
)
|
|
249
|
+
return result
|
|
250
|
+
|
|
251
|
+
def sum(self):
|
|
252
|
+
with self._lock:
|
|
253
|
+
result = D()
|
|
254
|
+
_check(_lib.hmath_tensor_sum(None, self._live(), C.byref(result)))
|
|
255
|
+
return result.value
|
|
256
|
+
|
|
257
|
+
def _binary(self, other, operation, scalar=False):
|
|
258
|
+
# Stable lock ordering prevents opposite-operand deadlocks.
|
|
259
|
+
first, second = sorted((self, other), key=id)
|
|
260
|
+
with first._lock, second._lock:
|
|
261
|
+
result = D() if scalar else P()
|
|
262
|
+
_check(operation(None, self._live(), other._live(), C.byref(result)))
|
|
263
|
+
return result.value if scalar else self._adopt(result)
|
|
264
|
+
|
|
265
|
+
def dot(self, other):
|
|
266
|
+
return self._binary(other, _lib.hmath_tensor_dot, True)
|
|
267
|
+
|
|
268
|
+
def __add__(self, other):
|
|
269
|
+
return self._binary(other, _lib.hmath_tensor_add)
|
|
270
|
+
|
|
271
|
+
def __matmul__(self, other):
|
|
272
|
+
return self._binary(other, _lib.hmath_tensor_matmul)
|
|
273
|
+
|
|
274
|
+
def reshape(self, shape):
|
|
275
|
+
shape, dims = _shape(shape)
|
|
276
|
+
with self._lock:
|
|
277
|
+
result = P()
|
|
278
|
+
_check(
|
|
279
|
+
_lib.hmath_tensor_reshape(
|
|
280
|
+
self._live(), dims, len(shape), C.byref(result)
|
|
281
|
+
)
|
|
282
|
+
)
|
|
283
|
+
return self._adopt(result, self._owner)
|
|
@@ -0,0 +1,127 @@
|
|
|
1
|
+
# SPDX-License-Identifier: MPL-2.0
|
|
2
|
+
"""Optional cffi ABI binding; no compiler needed. Explicit close/with cleanup."""
|
|
3
|
+
|
|
4
|
+
import os
|
|
5
|
+
from itertools import islice
|
|
6
|
+
|
|
7
|
+
from cffi import FFI
|
|
8
|
+
|
|
9
|
+
ffi = FFI()
|
|
10
|
+
ffi.cdef("""
|
|
11
|
+
typedef struct hmath_tensor hmath_tensor;
|
|
12
|
+
unsigned int hmath_abi_version(void); int hmath_init(void);
|
|
13
|
+
const char *hmath_last_error_message(void);
|
|
14
|
+
int hmath_scalar_add(void*,double,double,double*);
|
|
15
|
+
int hmath_tensor_copy(const double*,size_t,const size_t*,size_t,hmath_tensor**);
|
|
16
|
+
void hmath_tensor_release(hmath_tensor*);
|
|
17
|
+
size_t hmath_tensor_count(const hmath_tensor*);
|
|
18
|
+
int hmath_tensor_copy_to(const hmath_tensor*,double*,size_t);
|
|
19
|
+
int hmath_tensor_sum(void*,const hmath_tensor*,double*);
|
|
20
|
+
int hmath_tensor_dot(void*,const hmath_tensor*,const hmath_tensor*,double*);
|
|
21
|
+
int hmath_tensor_add(void*,const hmath_tensor*,const hmath_tensor*,hmath_tensor**);
|
|
22
|
+
int hmath_tensor_matmul(void*,const hmath_tensor*,const hmath_tensor*,hmath_tensor**);
|
|
23
|
+
""")
|
|
24
|
+
lib = ffi.dlopen(os.environ.get("HMATH_LIBRARY", "libhmath.so"))
|
|
25
|
+
try:
|
|
26
|
+
if lib.hmath_abi_version() != 1:
|
|
27
|
+
raise RuntimeError("HMath binding requires ABI 1")
|
|
28
|
+
except AttributeError as error:
|
|
29
|
+
raise RuntimeError(
|
|
30
|
+
"HMath binding requires ABI 1; library has no version query"
|
|
31
|
+
) from error
|
|
32
|
+
if lib.hmath_init():
|
|
33
|
+
raise RuntimeError("HMath runtime was permanently shut down")
|
|
34
|
+
|
|
35
|
+
|
|
36
|
+
class HMathError(RuntimeError):
|
|
37
|
+
def __init__(self, status):
|
|
38
|
+
self.status = status
|
|
39
|
+
super().__init__(ffi.string(lib.hmath_last_error_message()).decode())
|
|
40
|
+
|
|
41
|
+
|
|
42
|
+
def check(status):
|
|
43
|
+
if status:
|
|
44
|
+
raise HMathError(status)
|
|
45
|
+
|
|
46
|
+
|
|
47
|
+
def add(a, b):
|
|
48
|
+
result = ffi.new("double *")
|
|
49
|
+
check(lib.hmath_scalar_add(ffi.NULL, a, b, result))
|
|
50
|
+
return result[0]
|
|
51
|
+
|
|
52
|
+
|
|
53
|
+
class Tensor:
|
|
54
|
+
"""Do not close a cffi tensor concurrently with operations on that tensor."""
|
|
55
|
+
|
|
56
|
+
def __init__(self, values, shape=None):
|
|
57
|
+
values = list(islice(values, 1048577))
|
|
58
|
+
shape = list(shape if shape is not None else [len(values)])
|
|
59
|
+
if (
|
|
60
|
+
len(shape) > 8
|
|
61
|
+
or len(values) > 1048576
|
|
62
|
+
or any(type(n) is not int or n < 0 or n > 1048576 for n in shape)
|
|
63
|
+
):
|
|
64
|
+
raise ValueError("invalid bounded shape")
|
|
65
|
+
out = ffi.new("hmath_tensor **")
|
|
66
|
+
check(
|
|
67
|
+
lib.hmath_tensor_copy(
|
|
68
|
+
ffi.new("double[]", values),
|
|
69
|
+
len(values),
|
|
70
|
+
ffi.new("size_t[]", shape),
|
|
71
|
+
len(shape),
|
|
72
|
+
out,
|
|
73
|
+
)
|
|
74
|
+
)
|
|
75
|
+
self.ptr = ffi.gc(out[0], lib.hmath_tensor_release)
|
|
76
|
+
self.shape = tuple(shape)
|
|
77
|
+
|
|
78
|
+
def _live(self):
|
|
79
|
+
if self.ptr == ffi.NULL:
|
|
80
|
+
raise RuntimeError("Tensor is closed")
|
|
81
|
+
return self.ptr
|
|
82
|
+
|
|
83
|
+
def close(self):
|
|
84
|
+
if self.ptr != ffi.NULL:
|
|
85
|
+
ffi.release(self.ptr)
|
|
86
|
+
self.ptr = ffi.NULL
|
|
87
|
+
|
|
88
|
+
def __enter__(self):
|
|
89
|
+
self._live()
|
|
90
|
+
return self
|
|
91
|
+
|
|
92
|
+
def __exit__(self, *args):
|
|
93
|
+
self.close()
|
|
94
|
+
|
|
95
|
+
def values(self):
|
|
96
|
+
n = lib.hmath_tensor_count(self._live())
|
|
97
|
+
out = ffi.new("double[]", n)
|
|
98
|
+
check(lib.hmath_tensor_copy_to(self.ptr, out, n))
|
|
99
|
+
return list(out)
|
|
100
|
+
|
|
101
|
+
def sum(self):
|
|
102
|
+
out = ffi.new("double *")
|
|
103
|
+
check(lib.hmath_tensor_sum(ffi.NULL, self._live(), out))
|
|
104
|
+
return out[0]
|
|
105
|
+
|
|
106
|
+
def dot(self, other):
|
|
107
|
+
out = ffi.new("double *")
|
|
108
|
+
check(lib.hmath_tensor_dot(ffi.NULL, self._live(), other._live(), out))
|
|
109
|
+
return out[0]
|
|
110
|
+
|
|
111
|
+
def _binary(self, other, operation, shape):
|
|
112
|
+
out = ffi.new("hmath_tensor **")
|
|
113
|
+
check(operation(ffi.NULL, self._live(), other._live(), out))
|
|
114
|
+
result = object.__new__(Tensor)
|
|
115
|
+
result.ptr, result.shape = ffi.gc(out[0], lib.hmath_tensor_release), shape
|
|
116
|
+
return result
|
|
117
|
+
|
|
118
|
+
def __add__(self, other):
|
|
119
|
+
return self._binary(other, lib.hmath_tensor_add, self.shape)
|
|
120
|
+
|
|
121
|
+
def __matmul__(self, other):
|
|
122
|
+
shape = (
|
|
123
|
+
(self.shape[0], other.shape[1])
|
|
124
|
+
if len(self.shape) == len(other.shape) == 2
|
|
125
|
+
else ()
|
|
126
|
+
)
|
|
127
|
+
return self._binary(other, lib.hmath_tensor_matmul, shape)
|
|
@@ -0,0 +1,103 @@
|
|
|
1
|
+
Metadata-Version: 2.4
|
|
2
|
+
Name: hmath-python
|
|
3
|
+
Version: 1.0.0
|
|
4
|
+
Summary: Python bindings for the Haskell HMath numerical library
|
|
5
|
+
License-Expression: MPL-2.0
|
|
6
|
+
Project-URL: Homepage, https://github.com/wavefnd/hmath
|
|
7
|
+
Project-URL: Source, https://github.com/wavefnd/hmath
|
|
8
|
+
Project-URL: Issues, https://github.com/wavefnd/hmath/issues
|
|
9
|
+
Keywords: mathematics,haskell,ffi,tensor
|
|
10
|
+
Classifier: Development Status :: 5 - Production/Stable
|
|
11
|
+
Classifier: Operating System :: POSIX :: Linux
|
|
12
|
+
Classifier: Programming Language :: Python :: 3
|
|
13
|
+
Classifier: Topic :: Scientific/Engineering :: Mathematics
|
|
14
|
+
Requires-Python: >=3.9
|
|
15
|
+
Description-Content-Type: text/markdown
|
|
16
|
+
License-File: LICENSE
|
|
17
|
+
Provides-Extra: cffi
|
|
18
|
+
Requires-Dist: cffi>=1.15; extra == "cffi"
|
|
19
|
+
Provides-Extra: numpy
|
|
20
|
+
Requires-Dist: numpy>=1.22; extra == "numpy"
|
|
21
|
+
Dynamic: license-file
|
|
22
|
+
|
|
23
|
+
# HMath
|
|
24
|
+
|
|
25
|
+
HMath is a mathematics and numerical computing library with Python bindings, written in Haskell. It also provides native Haskell, C, and Wave interfaces.
|
|
26
|
+
|
|
27
|
+
The 1.0.0 release targets Linux x86-64. The Python distribution is named **hmath-python** and is imported as `hmath`.
|
|
28
|
+
|
|
29
|
+
## Python quick start
|
|
30
|
+
|
|
31
|
+
Install the Python package on Linux x86-64:
|
|
32
|
+
|
|
33
|
+
```sh
|
|
34
|
+
pip install hmath-python
|
|
35
|
+
```
|
|
36
|
+
|
|
37
|
+
```python
|
|
38
|
+
from hmath import Tensor, add
|
|
39
|
+
|
|
40
|
+
print(add(2, 3)) # 5.0
|
|
41
|
+
|
|
42
|
+
with Tensor.matrix([1, 2, 3, 4], 2, 2) as matrix:
|
|
43
|
+
with matrix @ matrix as product:
|
|
44
|
+
print(product.values()) # [7.0, 10.0, 15.0, 22.0]
|
|
45
|
+
```
|
|
46
|
+
|
|
47
|
+
Use `with` to release native resources automatically. Invalid operations raise `HMathError`. Install `hmath-python[numpy,cffi]` for optional NumPy and cffi support.
|
|
48
|
+
|
|
49
|
+
## Features
|
|
50
|
+
|
|
51
|
+
- **Core math:** checked scalar operations, numeric types, constants, units, number theory, and polynomials.
|
|
52
|
+
- **Linear algebra:** dense and sparse matrices, decompositions, iterative solvers, and CPU tensors.
|
|
53
|
+
- **Scientific computing:** differentiation, integration, ODE/PDE solvers, optimization, interpolation, FFTs, and signal processing.
|
|
54
|
+
- **Statistics and learning:** reproducible sampling, probability distributions, inference, automatic differentiation, and machine learning primitives.
|
|
55
|
+
- **Discrete and symbolic math:** combinatorics, graphs, numerical data structures, symbolic calculus, and a bounded expression evaluator.
|
|
56
|
+
- **Applied math:** complex numbers, rotations, computational geometry, time series, and numerical image processing.
|
|
57
|
+
- **Native interface:** checked scalar, vector, matrix, and tensor operations; explicit ownership and errors; Python/NumPy and Wave bindings.
|
|
58
|
+
- **Interchange:** typed binary and text arrays, CSV/TSV, Matrix Market, NPY/NPZ, and CPU DLPack.
|
|
59
|
+
|
|
60
|
+
The Haskell library contains the full numerical API; Python currently exposes scalar addition and binary64 vector, matrix, and tensor construction, arithmetic, and reductions.
|
|
61
|
+
|
|
62
|
+
## Build from source
|
|
63
|
+
|
|
64
|
+
Requires GHC, Cabal, a C11 compiler, Meson >= 1.3, Ninja, and Python 3. Zlib enables compressed interchange. BLAS/LAPACKE are optional runtime backends; the full spectral test suite needs LAPACKE.
|
|
65
|
+
|
|
66
|
+
```sh
|
|
67
|
+
git clone https://github.com/wavefnd/hmath.git
|
|
68
|
+
cd hmath
|
|
69
|
+
meson setup builddir
|
|
70
|
+
meson compile -C builddir
|
|
71
|
+
meson test -C builddir --print-errorlogs
|
|
72
|
+
```
|
|
73
|
+
|
|
74
|
+
This builds `builddir/libhmath.so` and runs the Haskell, C, and Python tests. Wave tests also run when `wavec` is available.
|
|
75
|
+
|
|
76
|
+
## Haskell
|
|
77
|
+
|
|
78
|
+
Save the following as `Main.hs` in the checkout. Run `cabal build lib:hmath`, then `cabal exec -- runghc -package=hmath Main.hs`:
|
|
79
|
+
|
|
80
|
+
```haskell
|
|
81
|
+
import HMath
|
|
82
|
+
|
|
83
|
+
main :: IO ()
|
|
84
|
+
main = do
|
|
85
|
+
print (dotChecked [1, 2, 3 :: Int] [4, 5, 6]) -- Right 32
|
|
86
|
+
print $ do
|
|
87
|
+
matrix <- fromListChecked 2 2 [1, 2, 3, 4 :: Int]
|
|
88
|
+
pure (toList (transposeMatrix matrix)) -- Right [1,3,2,4]
|
|
89
|
+
```
|
|
90
|
+
|
|
91
|
+
Checked operations return `Either MathError` for invalid shapes or lengths. More examples are available in [examples/haskell](examples/haskell).
|
|
92
|
+
|
|
93
|
+
## C and Wave
|
|
94
|
+
|
|
95
|
+
The [C header](include/hmath.h) defines errors, ownership, and runtime lifetime. See the [C example](examples/c/main.c) and [Wave bindings](bindings/wave) for usage.
|
|
96
|
+
|
|
97
|
+
## Contributing
|
|
98
|
+
|
|
99
|
+
Report bugs or discuss changes through [GitHub issues](https://github.com/wavefnd/hmath/issues). Submit pull requests from a fork against `master`; see [ROADMAP.md](ROADMAP.md) for planned work.
|
|
100
|
+
|
|
101
|
+
## License
|
|
102
|
+
|
|
103
|
+
[Mozilla Public License 2.0](LICENSE).
|
|
@@ -0,0 +1,10 @@
|
|
|
1
|
+
LICENSE
|
|
2
|
+
README.md
|
|
3
|
+
pyproject.toml
|
|
4
|
+
bindings/python/hmath.py
|
|
5
|
+
bindings/python/hmath_cffi.py
|
|
6
|
+
bindings/python/hmath_python.egg-info/PKG-INFO
|
|
7
|
+
bindings/python/hmath_python.egg-info/SOURCES.txt
|
|
8
|
+
bindings/python/hmath_python.egg-info/dependency_links.txt
|
|
9
|
+
bindings/python/hmath_python.egg-info/requires.txt
|
|
10
|
+
bindings/python/hmath_python.egg-info/top_level.txt
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
|
|
@@ -0,0 +1,33 @@
|
|
|
1
|
+
# SPDX-License-Identifier: MPL-2.0
|
|
2
|
+
[build-system]
|
|
3
|
+
requires = ["setuptools>=77"]
|
|
4
|
+
build-backend = "setuptools.build_meta"
|
|
5
|
+
|
|
6
|
+
[project]
|
|
7
|
+
name = "hmath-python"
|
|
8
|
+
version = "1.0.0"
|
|
9
|
+
description = "Python bindings for the Haskell HMath numerical library"
|
|
10
|
+
readme = "README.md"
|
|
11
|
+
requires-python = ">=3.9"
|
|
12
|
+
license = "MPL-2.0"
|
|
13
|
+
license-files = ["LICENSE"]
|
|
14
|
+
keywords = ["mathematics", "haskell", "ffi", "tensor"]
|
|
15
|
+
classifiers = [
|
|
16
|
+
"Development Status :: 5 - Production/Stable",
|
|
17
|
+
"Operating System :: POSIX :: Linux",
|
|
18
|
+
"Programming Language :: Python :: 3",
|
|
19
|
+
"Topic :: Scientific/Engineering :: Mathematics",
|
|
20
|
+
]
|
|
21
|
+
|
|
22
|
+
[project.urls]
|
|
23
|
+
Homepage = "https://github.com/wavefnd/hmath"
|
|
24
|
+
Source = "https://github.com/wavefnd/hmath"
|
|
25
|
+
Issues = "https://github.com/wavefnd/hmath/issues"
|
|
26
|
+
|
|
27
|
+
[project.optional-dependencies]
|
|
28
|
+
cffi = ["cffi>=1.15"]
|
|
29
|
+
numpy = ["numpy>=1.22"]
|
|
30
|
+
|
|
31
|
+
[tool.setuptools]
|
|
32
|
+
package-dir = {"" = "bindings/python"}
|
|
33
|
+
py-modules = ["hmath", "hmath_cffi"]
|