Parallel Finite Element - General geometry Ewald-like Method
Part of Continuum-particle Simulation Suite under MICCOM
|
#include <stdio.h>
#include <fstream>
#include <cstdlib>
#include <ctime>
#include <random>
#include "libmesh/libmesh_logging.h"
#include "libmesh/equation_systems.h"
#include "rigid_particle.h"
#include "particle_mesh.h"
#include "point_particle.h"
#include "point_mesh.h"
#include "brownian_system.h"
#include "chebyshev.h"
Functions | |
PetscErrorCode | _MatMult_Stokes (Mat M, Vec f, Vec u) |
PetscErrorCode _MatMult_Stokes | ( | Mat | M, |
Vec | f, | ||
Vec | u | ||
) |
Internal function if shell matrix mode is used, this just calls the shell matrix's matrix multiplication function. In this class, it will call Stokes solver. u = M*f NOTE that the pm_system should be reinit() before solving Stokes