Thank you for all the repliers!
All of you gave me a good advice that I may consider deeply when I write my
code from now
I have figured out what was wrong.
Seeing at Step 20 code, I realized that I have to select component of
solution before going to compare the difference.
const ComponentSel
Jaekwang,
Have a look at step-10. There you see how the order of the mapping
influences the error.
Basically, you have to set the mapping for your FEValues objects, when you
compute error via integrate_difference
and when you project or interpolate some Function to an FE Vector.
Best,
Daniel
A
2016년 9월 18일 일요일 오전 10시 57분 16초 UTC-5, JAEKWANG KIM 님의 말:
>
>
> Hello, I am a starter of dealii and am learning a lot these days with the
> help of video lectures and tutorial examples.
>
> I modified step-22 code (stokes flow code) into my own problem, the flow
> around sphere.
>
> and I inte
Thank you for all repliers!!
I tried fix code as you mentioned please see redline...
template
void StokesProblem::assemble_system ()
{
system_matrix=0;
system_rhs=0;
*QGaussLobatto<**2> quadrature_formula(degree+2);*
FEValues fe_values (fe, quadrature_formula,
On 09/20/2016 09:04 AM, JAEKWANG KIM wrote:
yes...I tried that before, but I get an error message as follow..
As a general workflow, you need to
1/ make the code compile
2/ run it without any runtime error
3/ make sure the results are correct.
In your original message, you were stopped in 1 be
...of course I meant dof_handler.distribute_dofs(). Does the DoFhandler
know at this point about the Trinagulation?
Best,
Daniel
Am Dienstag, 20. September 2016 17:58:02 UTC+2 schrieb Daniel Arndt:
>
> Jaekwang,
>
> This seems to be unrelated. Did you setup the Triangulation before calling
> fe
Jaekwang,
This seems to be unrelated. Did you setup the Triangulation before calling
fe.distribute_dofs()?
Best,
Daniel
Am Dienstag, 20. September 2016 17:04:02 UTC+2 schrieb JAEKWANG KIM:
>
> yes...I tried that before, but I get an error message as follow..
>
> As you mentioned, I first learn
I understand what you meant but I have no idea how to fix it yet.(Maybe
because I lack experiences)
How do I designate mapping method specifically?
Can I see examples of specific command lines to for mapping designation?
The original initialization of fe was follow.
template
yes...I tried that before, but I get an error message as follow..
As you mentioned, I first learn this function in step 7 tutorial.
But the problem is that I want to integrate my vector-valued-solution to
exact-vector-valued solution, while step 7 compares scalar-solution.
Or is there any way
Oh, now I see: You should write "Solution", not "Solution".
The template parameters gives the space dimension of the function, not
the number of components. The latter is specified in the constructor of
your class Solution through the argument (dim+1) to Function.
Best,
Martin
On 09/20/2016 04:5
Jaekwang,
can you post the full error message? To me the call looks reasonable
(this is a copy from step-7) and the function definition looks good as well.
Ah, and in case you want to use this in the context of your other
question regarding high order methods: Put a "mapping" as first argument
to
thanks!! Full error message is as follow
*/Users/jaekwangjk/repos/trial/sphere-deal.ii/mystokes.cc:1258:7: **error: *
*no*
* matching function for call to 'integrate_difference'*
VectorTools::integrate_difference (dof_handler,
* ^*
*/Users/jaek
I wanted to evaluate Vectortool::integrate_difference to estimate my error.
but I couldn't figure out what I need to enter for the third component. (I
marked it as red)
VectorTools::integrate_difference (dof_handler,
solution,
Dear Jaekwang,
if you are using the most recent developer version, you should
automatically get the Gauss-Lobatto version of the node distribution. We
made those point distributions the default for FE_Q(degree) this spring.
Regarding the limits of quadrature formulas: Depending on what exactly
yo
To randomly add to this discussion...
In general higher order polynomials are good for approximating smooth
functions, such as the tail of an exponential decay. Interestingly, a
denser grid with lower order polynomials sometimes works better for
approximating a function that oscillates rapidly, s
thank for the reply!!
my fe degree is declared as "fe (FE_Q(degree+1), dim,FE_Q(degree),
1)" and I used "QGauss quadrature_formula(degree+2);" to calculate
integral over the cell.
2016년 9월 20일 화요일 오전 9시 26분 6초 UTC-5, Praveen C 님의 말:
>
>
> On Tue, Sep 20, 2016 at 7:49 PM, JAEKWANG KIM > wrot
On Tue, Sep 20, 2016 at 7:49 PM, JAEKWANG KIM wrote:
> Can I ask more about "the limits of our implementation of our quadrature
> formulas?".
> I wonder when it usually happens.
>
>
> Once I calculated drag coefficient, with Q1, my error is 10% compare to
> exact solution.
> However, I can 1.7% e
Can I ask more about "the limits of our implementation of our quadrature
formulas?".
I wonder when it usually happens.
Once I calculated drag coefficient, with Q1, my error is 10% compare to
exact solution.
However, I can 1.7% error when I use Q2 which is significant decrease!
At Q3, I get 0
Thank you so much, Praveen. Your codes are really helpful.
Regards,
Jiaqi
2016-09-20 0:11 GMT-04:00 Praveen C :
>
> On Tue, Sep 20, 2016 at 9:34 AM, Jiaqi ZHANG
> wrote:
>
>> Hello Praveen,
>>
>> Thank you so much. I am wondering if you could tell me which
>>
>> limiters you used, because if it
The action you are requesting is multiplying each row (column?) of A with
the corresponding element of V. While we do not have this implemented, it is
a) useful,
b) well within the framework of sparse matrix operations, since it does not
enlarge the stencil
If you want to implement a function
te
Dear deal.ii Users,
I have a question regarding the dof_handler.load () function.
Before calling the load function for a dof_handler with assigned
triangulation I have to call the function distribute_dofs with the FE
object as argument.
By doing this DoF distribution data is created just to be o
Dear Daniel,
To clarify what I am searching for, I would like to get the vector-valued
solution into a format that matches a scalar. Either as a
std::vector> or a std::vector>, the format doesn't
matter as long as it matches a scalar dofs. So I can manipulate the vector
components and for exam
On Monday, September 19, 2016 at 9:04:35 AM UTC+2, Ehsan Mihankhah wrote:
>
> I managed to compile the sample codes but the problem is that when I run
> the code I receive errors from ARPACK referring to "info=-9" in Pdnaupd.
> and following "info=-9" I see the error is about initial vector bei
23 matches
Mail list logo
