Commit b55b1150 authored by Jan Zapletal's avatar Jan Zapletal

MAINT: minor cleanup

parent c6d2eaf3
......@@ -27,8 +27,7 @@ BEBilinearFormLaplace1Layer<LO, SC>::BEBilinearFormLaplace1Layer(
BESpace<LO, SC>* space,
int* quadratureOrder,
quadratureType quadrature,
int* quadratureOrderDisjointElems,
bool employSymmetricity
int* quadratureOrderDisjointElems
) {
this->space = space;
this->quadrature = quadrature;
......
......@@ -43,8 +43,7 @@ public:
BESpace<LO, SC>* space,
int* quadratureOrder = nullptr,
quadratureType quadrature = SauterSchwab,
int* quadratureOrderDisjointElems = nullptr,
bool employSymmetricity = false
int* quadratureOrderDisjointElems = nullptr
);
//! destructor
......
......@@ -118,7 +118,7 @@ void getLaplaceSteklovPoincare(
FullMatrix< LO, SC > * V = new FullMatrix< LO, SC >( 0, 0 );
BEBilinearFormLaplace1Layer< LO, SC > formV( &bespace00, quadNear,
quadType, quadFar, false );
quadType, quadFar );
if ( verbose ) ProgressMonitor::init( "V" );
formV.assemble( *V );
if ( verbose ) ProgressMonitor::step( );
......@@ -243,14 +243,18 @@ void evaluateLaplaceRepresentationFormula(
K->apply( dir, neu, false, 1.0, 0.0 );
// local solve (V decomposed in assembly)
if ( verbose ) ProgressMonitor::init( "Applying inverse of V" );
V->CholeskiDecomposedSolve( neu );
if ( verbose ) ProgressMonitor::step( );
RepresentationFormulaLaplace<LO, SC> formula( &bespace10, &dir, &neu,
orderFar );
Vector< LO, SC > val( nPoints, values, false );
if ( verbose ) ProgressMonitor::init( "Representation formula" );
formula.evaluate( points, nPoints, true, val );
if ( verbose ) ProgressMonitor::step( );
}
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment