Python RandomDolfinFunction Examples

Example #1

 def generate_random(self):
     # Generate random vectors
     S = FunctionsList(self.V)
     for _ in range(self.N):
         b = RandomDolfinFunction(self.V)
         S.enrich(b)
     F = RandomDolfinFunction(self.V)
     # Return
     return (S, F)

Example #2

File: test_dolfin_v1_dot_A_v2.py Project: nadia-el/RBniCS

 def generate_random(self):
     # Generate random vectors
     v1 = RandomDolfinFunction(self.V).vector()
     v2 = RandomDolfinFunction(self.V).vector()
     k = RandomDolfinFunction(self.V)
     # Generate random matrix
     A = assemble(self.a(k))
     # Return
     return (v1, v2, A)

Example #3

File: test_dolfin_Z_T_dot_F.py Project: nadia-el/RBniCS

 def generate_random(self):
     # Generate random vectors
     Z = BasisFunctionsMatrix(self.V)
     Z.init("u")
     for _ in range(self.N):
         b = RandomDolfinFunction(self.V)
         Z.enrich(b)
     F = RandomDolfinFunction(self.V)
     # Return
     return (Z, F)

Example #4

File: test_dolfin_S_T_dot_A_S.py Project: pp1565156/RBniCS

 def generate_random(self):
     # Generate random vectors
     S = FunctionsList(self.V)
     for _ in range(self.N):
         b = RandomDolfinFunction(self.V)
         S.enrich(b)
     k = RandomDolfinFunction(self.V)
     # Generate random matrix
     A = assemble(self.a(k))
     # Return
     return (S, A)

Example #5

File: test_dolfin_error_estimation_af_construction.py Project: twisterbboy/RBniCS

 def generate_random(self):
     # Generate random vectors
     Z = BasisFunctionsMatrix(self.V)
     Z.init("u")
     for i in range(self.N):
         b = RandomDolfinFunction(self.V)
         Z.enrich(b)
     # Generate random matrix
     k = RandomDolfinFunction(self.V)
     A = assemble(self.a(k))
     # Generate random function
     z = RandomDolfinFunction(self.V)
     # Return
     return (Z, A, z.vector())

Example #6

File: test_dolfin_nonlinear_solver.py Project: twisterbboy/RBniCS

    def generate_random(self):
        # Generate random forcing
        g = RandomDolfinFunction(self.V)
        # Generate correspondingly residual and jacobian forms
        r = self.r(self.u, g)
        j = self.j(self.u, r)

        # Prepare problem wrapper
        class ProblemWrapper(NonlinearProblemWrapper):
            # Residual and jacobian functions
            def residual_eval(self_, solution):
                return self.callback(r)

            def jacobian_eval(self_, solution):
                return self.callback(j)

            # Define boundary condition
            def bc_eval(self_):
                return None

            # Empty solution monitor
            def monitor(self_, solution):
                pass

        problem_wrapper = ProblemWrapper()
        # Return
        return (r, j, problem_wrapper)

Example #7

File: test_dolfin_vector_assembly.py Project: pp1565156/RBniCS

 def generate_random(self):
     f = ()
     for i in range(self.Q):
         # Generate random vector
         g = RandomDolfinFunction(self.V)
         # Generate random form
         f += (self.f(g), )
     F = AffineExpansionStorage(f)
     # Genereate random theta
     theta = RandomTuple(self.Q)
     # Return
     return (theta, F)

Example #8

File: test_dolfin_matrix_assembly.py Project: twisterbboy/RBniCS

 def generate_random(self):
     a = ()
     for i in range(self.Q):
         # Generate random vector
         k = RandomDolfinFunction(self.V)
         # Generate random form
         a += (self.a(k), )
     A = AffineExpansionStorage(a)
     # Genereate random theta
     theta = RandomTuple(self.Q)
     # Return
     return (theta, A)

Example #9

File: test_dolfin_v1_dot_v2.py Project: twisterbboy/RBniCS

 def generate_random(self):
     # Generate random vectors
     v1 = RandomDolfinFunction(self.V).vector()
     v2 = RandomDolfinFunction(self.V).vector()
     # Return
     return (v1, v2)

Example #10

File: test_dolfin_linear_solver.py Project: pp1565156/RBniCS

 def generate_random(self):
     # Generate random rhs
     g = RandomDolfinFunction(self.V)
     # Return
     return (self.a, self.f(g))