#include <ripley/RipleyDomain.h>
#include <escript/DataFactory.h>
#include <escript/FunctionSpaceFactory.h>
#include <pasowrap/SystemMatrixAdapter.h>
#include <pasowrap/TransportProblemAdapter.h>
#include <ripley/domainhelpers.h>
#include <iomanip>
#define UPDATE_BLOCK |
( |
|
VAL | ) |
|
Value:do {\
for (
dim_t ic=0; ic<mat->col_block_size; ++ic) {\
const
dim_t i_Sol=ic+mat->col_block_size*l_col;\
for (
dim_t ir=0; ir<mat->row_block_size; ++ir) {\
const
dim_t i_Eq=ir+mat->row_block_size*l_row;\
VAL[k*mat->block_size+ir+mat->row_block_size*ic]\
+= array[
INDEX4(i_Eq, i_Sol, k_Eq, k_Sol, numEq, numEq, nodes.size())];\
}\
}\
} while(0)
#define INDEX4(_X1_, _X2_, _X3_, _X4_, _N1_, _N2_, _N3_)
Definition: index.h:41
index_t dim_t
Definition: types.h:27
Referenced by ripley::RipleyDomain::addToSystemMatrix().