Ina,
Le lun. 1 avr. 2019 à 05:54, Ina a écrit :
> Do you know how to create the matrix needed by the interpolate function using
> the data of the transfer representation?
No I don't sorry. Personally, since you said that
VectorTools::interpolate_to_different_mesh was working I would just
use
On 4/4/19 4:05 AM, Muhammad Mashhood wrote:
>
> for (unsigned int f=0; f :: faces_per_cell; ++f)
> if (cell->face(f)->at_boundary() )
> {
> fe_face_values.reinit (cell,f);
> for (unsigned int i=0; i {
>
Hey Daniel,
thanks for the pull, with this instantiation everything compiles fine and I can use it.
But I think I did a mistake somewhere in the function. Using the definiion of the function from your pull request
I call
DoFTools::make_sparsity_pattern
Dear Maxi,
The problem is the quadrature formula you use in 'assemble_system' and
'calculate_residual'. You should use 'QGauss
quadrature_formula(fe.degree+1)`. I am surprised that quadratic elements
even worked because the Poisson matrix for d>=2 should have some rank
deficiency with
Thank you Prof. Bangerth for your quick concern. I am on basic level for
FEM code programming especially deal.ii and also quite only a bit familiar
with object oriented programming in C++.
Sorry if my question was not clear, let me explain it a bit:
Actually, I am trying to follow the tutorial
I wrote a test program, based on step-5, for testing unrelated things, and
checked the results. But when I increased the degree of the finite elements
from 1 to 3 (or higher), the results were wrong, while the change from 1 to
2 just improved convergence, as expected. Did I make a mistake in my
Hi Wolfgang,
Thank you so much for your kind help. I tried the dynamic sparsity pattern,
and with the sparsity_pattern.add_entries() and sparse_matrix.set()
function, I sucesefully insert values in the sparse matrix.
--
The deal.II project is located at http://www.dealii.org/
For mailing
Hi Juan Carlos Araújo Cabarcas,
Thank you for your code, and it is really an excellent way. I also tried to
use the dynamic sparsity pattern in dealii, and use the add_entries()
function to indicate the nonzeros' locations, then I successfully insert
the values in the sparse matrix as you did