Skip to content

new singular quadrature routine#24

Merged
luca-heltai merged 4 commits into
mathLab:masterfrom
nicola-giuliani:new_quadrature_set_up
Jan 16, 2018
Merged

new singular quadrature routine#24
luca-heltai merged 4 commits into
mathLab:masterfrom
nicola-giuliani:new_quadrature_set_up

Conversation

@nicola-giuliani

Copy link
Copy Markdown
Contributor

updated the singular quadrature procedure and introduced the new QSplit formulas for BEM in 3D

Comment thread source/bem_problem.cc Outdated
// (QGaussOneOverR<dim-1>(singular_quadrature_order,
// fe->get_unit_support_points()[i],true));
}
sing_quadratures[i] = get_singular_quadrature(i);

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Why do you create a new vector here? Since you create it as a static vector in the get_singular_quadrature function, there is no need to copy it back here.

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

just call get_singular_quadrature(i) when you need it, instead of sing_quadratures[i]

@luca-heltai luca-heltai left a comment

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

.

@luca-heltai

Copy link
Copy Markdown
Contributor

indent.

Comment thread source/bem_problem.cc Outdated
singular_quadrature
= dynamic_cast<Quadrature<dim-1>*>(
&sing_quadratures[singular_index]);
&(get_singular_quadrature[singular_index]));

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

this is a function. It cannot work this way. replace [] with ()

@luca-heltai
luca-heltai merged commit 54a90f2 into mathLab:master Jan 16, 2018
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants