112 lines
2.4 KiB
Python
112 lines
2.4 KiB
Python
import collections
|
|
|
|
import ssgetpy
|
|
|
|
Dataset = collections.namedtuple("Dataset", ["name", "mats"])
|
|
|
|
|
|
def filter_reject_blacklist(mats):
|
|
|
|
# matrices that are actually integer
|
|
blacklist_integer = [
|
|
"Alemdar",
|
|
"aug2d",
|
|
"aug2dc",
|
|
"aug3d",
|
|
"cyl6",
|
|
"engine",
|
|
"mc2depi",
|
|
"t520",
|
|
"tube2"
|
|
]
|
|
|
|
filtered = []
|
|
for mat in mats:
|
|
if mat.name in blacklist_integer:
|
|
print(f"BLACKLIST {mat.name}")
|
|
continue
|
|
filtered += [mat]
|
|
return filtered
|
|
|
|
def filter_reject_large(mats):
|
|
filtered = []
|
|
for mat in mats:
|
|
if mat.rows > 1_000_000 or mat.cols > 1_000_000 or mat.nnz > 20_000_000:
|
|
continue
|
|
filtered += [mat]
|
|
return filtered
|
|
|
|
def filter_reject_small(mats):
|
|
filtered = []
|
|
for mat in mats:
|
|
if mat.rows < 1_000 or mat.cols < 1_000 or mat.nnz < 20_000:
|
|
continue
|
|
filtered += [mat]
|
|
return filtered
|
|
|
|
## all real-valued matrices
|
|
REAL_MATS = Dataset(
|
|
name = "all_reals",
|
|
mats = filter_reject_blacklist(ssgetpy.search(
|
|
dtype='real',
|
|
limit=1_000_000
|
|
))
|
|
)
|
|
|
|
## certain matrices with regular structure
|
|
kinds = [
|
|
"2D/3D",
|
|
"Acoustics Problem",
|
|
"Materials Problem",
|
|
"Structural Problem",
|
|
"Computational Fluid Dynamics Problem",
|
|
"Model Reduction Problem",
|
|
"Semiconductor Device Problem",
|
|
"Theoretical/Quantum Chemistry Problem",
|
|
"Thermal Problem",
|
|
]
|
|
REGULAR_REAL_MATS = Dataset(
|
|
name = "regular_reals",
|
|
mats = []
|
|
)
|
|
mats = []
|
|
for kind in kinds:
|
|
mats += ssgetpy.search(
|
|
kind=kind,
|
|
dtype='real',
|
|
limit=1_000_000
|
|
)
|
|
REGULAR_REAL_MATS = Dataset(
|
|
name="regular_reals",
|
|
mats = filter_reject_blacklist(mats)
|
|
)
|
|
|
|
## keep "small" matrices
|
|
REGULAR_REAL_SMALL_MATS = Dataset (
|
|
name = "regular_reals_small",
|
|
mats = filter_reject_large(REGULAR_REAL_MATS.mats)
|
|
)
|
|
REAL_SMALL_MATS = Dataset (
|
|
name = "reals_small",
|
|
mats = filter_reject_large(REAL_MATS.mats)
|
|
)
|
|
|
|
## keep "medium" matrices
|
|
REGULAR_REAL_MED_MATS = Dataset (
|
|
name = "regular_reals_med",
|
|
mats = filter_reject_large(filter_reject_small(REGULAR_REAL_MATS.mats))
|
|
)
|
|
REAL_MED_MATS = Dataset (
|
|
name = "reals_med",
|
|
mats = filter_reject_large(filter_reject_small(REAL_MATS.mats))
|
|
)
|
|
|
|
## export all datasets
|
|
DATASETS = [
|
|
REAL_MATS,
|
|
REAL_SMALL_MATS,
|
|
REAL_MED_MATS,
|
|
REGULAR_REAL_MATS,
|
|
REGULAR_REAL_SMALL_MATS,
|
|
REGULAR_REAL_MED_MATS
|
|
] |