mutivec implemented in algorithms, random method modifed in utils

[V-Rang]

1. mutivec implemented in alogorithms, data is list of petsc Vectors
2. random method modifed in utils to extend to multivectors and made more concise.

[uvilla]

We also need a new file in the test suite to check the methods of MultiVector work as intended.

[uvilla]

Last minor tweaks missing before we can merge :)

[uvilla]

Let's replace init_vector=None with createVecLeft=None and createVecRight=None.
This two input arguments will be user-defined function handler that can be used to generate the correct vectors for this operator.

Then inside the class do:

self.createVecLeft = createVecLeft
self.createVecRight = createVecRight

[V-Rang]

Modifications have been made here:

hippylibx/hippylibX/algorithms/linalg.py

Lines 11 to 21 in 9d2273a

	class Solver2Operator:
	def __init__(
	self,
	S: Any,
	mpi_comm=MPI.COMM_WORLD,
	createVecLeft=None,
	createVecRight=None,
	) -> None:
	self.S = S
	self.createVecLeft = createVecLeft
	self.createVecRight = createVecRight

[uvilla]

self.tmp is not used anymore. self.tmp was used only in the inner method that we delete.
You can safely remove all this piece of code.

[V-Rang]

Code has been removed.

[uvilla]

The return type dict is not correct here.

[V-Rang]

return type changed to tuple:

hippylibx/hippylibX/algorithms/randomizedEigensolver.py

Lines 9 to 17 in 9d2273a

	def doublePassG(
	A: petsc4py.PETSc.Mat,
	B: petsc4py.PETSc.Mat,
	Binv: Any,
	Omega: MultiVector,
	k: int,
	s=1,
	check=False,
	) -> tuple[np.array, MultiVector]:

[uvilla]

you can simply write:
sample_petsc_vec = dlx.fem.petsc.create_vector(dlx.fem.form(test*ufl.dx))

[V-Rang]

Thanks for the suggestion. Code changed here:

hippylibx/hippylibX/test/test_multivector.py

Lines 28 to 31 in 9d2273a

	M.assemble()
	sample_petsc_vec = dlx.fem.petsc.create_vector(dlx.fem.form(test * ufl.dx))
	Omega = hpx.MultiVector.createFromVec(sample_petsc_vec, nvec)
	sample_petsc_vec.destroy()