solver-rewrite

init.py
View Options

			import scipy.sparse as _aka_sparse
			import numpy as _aka_np
			import py11_akantu

			private_keys = set(dir(py11_akantu)) - set(dir())

			for k in private_keys:
			globals()[k] = getattr(py11_akantu, k)

			if py11_akantu.has_mpi():
			try:
			from mpi4py import MPI # noqa: F401
			except Exception:
			pass


			def initialize(args, *kwargs):
			raise RuntimeError("No need to call initialize,"
			" use parseInput to read an input file")


			def finalize(args, *kwargs):
			raise RuntimeError("No need to call finalize")


			class AkantuSparseMatrix (_aka_sparse.coo_matrix):

			def __init__(self, aka_sparse):

			self.aka_sparse = aka_sparse
			matrix_type = self.aka_sparse.getMatrixType()
			sz = self.aka_sparse.size()
			row = self.aka_sparse.getIRN()[:, 0] - 1
			col = self.aka_sparse.getJCN()[:, 0] - 1
			data = self.aka_sparse.getA()[:, 0]

			row = row.copy()
			col = col.copy()
			data = data.copy()

			if matrix_type == py11_akantu._symmetric:
			non_diags = (row != col)
			row_sup = col[non_diags]
			col_sup = row[non_diags]
			data_sup = data[non_diags]
			col = _aka_np.concatenate((col, col_sup))
			row = _aka_np.concatenate((row, row_sup))
			data = _aka_np.concatenate((data, data_sup))

			_aka_sparse.coo_matrix.__init__(
			self, (data, (row, col)), shape=(sz, sz))


			FromStress = py11_akantu.FromHigherDim
			FromTraction = py11_akantu.FromSameDim
			py11_akantu.__initialize()