Hello everyone.
I am solving two equations iteratively in an MPI framework. The first
equation is the equilibrium equation (slightly modified) solved in step-8
while the second equation (Damage equation) is very similar to the Laplace
equation solved in step-3.
The solution of equilibrium
Error 2
make[1]: *** [CMakeFiles/Makefile2:123:
I hope this is what you were asking for.
Regards
Wasim
On Wednesday, April 19, 2023 at 6:29:56 PM UTC+5:30 d.arnd...@gmail.com
wrote:
> Wasim,
>
> What does the backtrace in your debugger (gdb?) look like?
>
> Best,
> Daniel
Hello everyone.
I am trying to apply a particular boundary condition at one single point in
the domain for a problem very similar to step-40.
I am using the following code to do it:
for (const auto : dof_handler.active_cell_iterators())
if (cell->is_locally_owned())
{
for (const auto :
Thank You Prof. Bangerth. I actually needed the norm of the entire vector
only. I was just curious if there existed a function that did it for each
processor.
Regards
Wasim
On Thu, Apr 13, 2023 at 9:50 PM Wolfgang Bangerth
wrote:
> On 4/13/23 06:22, Wasim Niyaz Munshi ce21d400 wr
Hello everyone. I have a vector distributed over multiple processors. I
want to compute the modulus of the entire vector and also the modulus of
the portions of vector stored by each process. I would really appreciate
some help in this regard.
Regards
Wasim
--
The deal.II project is located
, April 6, 2023 at 9:41:21 PM UTC+5:30 Wolfgang Bangerth wrote:
> On 4/6/23 10:06, Wasim Niyaz Munshi ce21d400 wrote:
> > **
> >
> > Yes, I also had the same feeling. But, when I look at the plot in the
> tutorial
> > of step-40 for 52M Dofs, I see that they have solve
Wolfgang Bangerth wrote:
> On 4/6/23 10:18, Wasim Niyaz Munshi ce21d400 wrote:
> > How do I get the no.of cells owned by the processor?
>
> Triangulation::n_locally_owned_active_cells().
>
> Best
> W.
>
> --
> --
by the processor?
Thanks and regards
Wasim
On Thursday, April 6, 2023 at 9:24:14 PM UTC+5:30 Wolfgang Bangerth wrote:
> On 4/6/23 06:02, Wasim Niyaz Munshi ce21d400 wrote:
> >
> > I don't have a solution_vector for a parallel code, but a
> > locally_relevant_solution. I wan
error even when I use 80 processors (250 GB memory).
Thanks and regards
Wasim Niyaz
Research scholar
CE Dept.
IITM
On Thu, 6 Apr, 2023, 9:21 pm Wolfgang Bangerth,
wrote:
> On 4/6/23 01:31, Wasim Niyaz Munshi ce21d400 wrote:
> > **
> >
> > I tried to run step-40 with 52M DOFs o
Hello everyone.
I want to extract the element solution vector from the global solution once
the problem is solved in step-40. For a serial code, I would do something
like this:
*int i=0;*
*for (const auto vertex : cell->vertex_indices()) { int a =
(cell->vertex_dof_index(vertex, 0));
5, 2023 at 11:07:25 PM UTC+5:30 Wolfgang Bangerth wrote:
> On 4/5/23 11:27, Wasim Niyaz Munshi ce21d400 wrote:
> > I am running in release mode. I am attaching the results for cycle 3 for
> > both debug and release modes. I will try to reproduce the plot of wall
> &g
I am running in release mode. I am attaching the results for cycle 3 for
both debug and release modes. I will try to reproduce the plot of wall time
vs the number of processors for 52M DOFs as given in the tutorial problem.
That would be a better way to compare the performances!
On Wednesday,
I am solving an elasticity problem wherein I have a pin support at one
vertex. Earlier, I was using std::map object to apply BCs. Now, I am trying
to solve the same problem using an affine constraint object to apply BCs.
But, I am facing issues filling my affine constraint object with the
Hello everyone.
I am working on a problem very similar to the Laplace equation solved in
step-3. I have made some changes to step-3 as I need to. Now, I want to
solve this problem in parallel. I looked at step-40, which solves the
Laplace equation in parallel. However, it uses a constraints
Thank you very much. It was indeed the easiest way to do it.
Regards
Wasim Niyaz
Research scholar
CE Dept.
IITM
On Mon, 13 Feb, 2023, 7:22 pm Bruno Turcksin,
wrote:
> Hello,
>
> DataOut is made to add data to a mesh. It cannot be used to do what you
> want. The easiest way to do what you want
Hello everyone. I have 2 Vector vectors A and B of the same size,n.
I want to plot A vs B (A along x and B along y). I was trying something
like this:
DataOut<1,Vector > data_out;
data_out.add_data_vector (A, "A");
data_out.add_data_vector (B, "B");
// Write the Gnuplot file
Hello everyone.
I am trying to extract all the nodes present on a boundary with a
particular id.
Would anyone be able to help me with this?
Also, for a vector problem how do I extract a particular dof ( eg. I want
the global dofs corresponding to the x-displacement of all the nodes on the
top
Hello everyone. I am trying to plot load displacement curves for an
elasticity problem.
For this I need to extract the non-zero elements of my sparse stiffness
matrix. I tried something like system_matrix[I][j] to extract the ij-th
entry but it seems this doesn't work for the SparseMatrix.
Can
ooks like you're passing your solution vector by copy for each cell
> that you're assembling on. You probably want to pass it by reference.
>
> Best,
> Jean-Paul
>
> Sent from my mobile device. Please excuse my brevity and any typos.
>
> On Sun, 08 Jan 2023, 14:38 Wasim
don't expect to match that time but is there still anything that I can do
to reduce the assembly time further?
Thanks and regards
Wasim
On Sunday, January 8, 2023 at 4:11:27 PM UTC+5:30 Wasim Niyaz Munshi
ce21d400 wrote:
> Do I need to make this change in both declaration and definition?
>
r (const unsigned int j : fe_values_damage.dof_indices())
>>>> system_matrix_damage.add(local_dof_indices[i],
>>>> local_dof_indices[j],
>>>> cell_matrix_damage(i, j));
>>>>
>>>> syst
;>
>>> Then I would time the function to measure if this is actually the real
>>> culprit or if it could be something else.
>>> If you copy/paste the content of your assembly code and the function, I
>>> would be glad to give it a look (and I am sure ot
the strain.
>>
>> Similarly I have another function damage_gauss which is being called in
>> assemble_elastic that evaluates damage at a gauss point using the damage
>> solution, cell and gauss point.
>>
>> Wasim Niyaz
>> Research scholar
>> CE Dept.
&g
I have another function damage_gauss which is being called in
assemble_elastic that evaluates damage at a gauss point using the damage
solution, cell and gauss point.
Wasim Niyaz
Research scholar
CE Dept.
IITM
On Sat, 7 Jan, 2023, 10:15 am Wasim Niyaz Munshi ce21d400, <
ce21d...@smail.iitm.ac
I use it to evaluate strain at Gauss points. Then, i evaluate some quantity
which is a function of this strain.
Wasim Niyaz
Research scholar
CE Dept.
IITM
On Sat, 7 Jan, 2023, 3:09 am Wolfgang Bangerth,
wrote:
> On 1/6/23 13:53, Wasim Niyaz Munshi ce21d400 wrote:
> > I am using 65536
Hello everyone.
I am trying to solve elasticity and laplace equations one after the other.
The elasticity solution feeds as input to the laplace equation (stiffness
matrix of laplace equation depends on elasticity solution) and vice versa.
I am using 65536 elements. For step-8 the assembly takes
t;
> Le lun. 21 nov. 2022 à 15:27, Wasim Niyaz Munshi ce21d400 <
> ce21d...@smail.iitm.ac.in> a écrit :
>
>> Thank you, Prof. Bruno.
>>
>> I used the above line of code and realized where I was going wrong. When
>> solving the damage equation, I also
o
>
> std::cout<get_fe().get_name()<
> and check that it matches what you have set in FEValues.
>
> Best,
>
> Bruno
>
> Le lun. 21 nov. 2022 à 11:08, Wasim Niyaz Munshi ce21d400 <
> ce21d...@smail.iitm.ac.in> a écrit :
>
>> I am working with 2 DoFHandle
I am working with 2 DoFHandlers. I checked my code. I am using the correct
DoFHandlers for the 2 equations.
On Mon, Nov 21, 2022 at 6:58 PM Bruno Turcksin
wrote:
> Wasim,
>
> It's hard to say without seeing any code. Are you working with one or two
> DoFHandler? If you are working with two
Hello everyone.
I am trying to solve elasticity and *laplace (damage)* equations in the
same program. The rhs of laplace equation depends on the elasticity
solution vector. Both equations have already been solved in step3 and 8. I
am following the same approach. I am creating my objects (like
_multiplier] = 0;*/*
}
else if (std::fabs(vert(0) - 1) < 1e-12 &&
std::fabs(vert(1) - (0)) < 1e-12)
{
std::cout<<"vert_0 = "<mailto:bang...@colostate.edu>> wrote:
> > >
> > > On 10/10/22 07:42, Wasim Niyaz Mun
ctions::ConstantFunction<2>(u_y_values,2),
boundary_values,u_y_mask
);//imposing
u_y=0.05 on top face
Thank You
On Wednesday, October 12, 2022 at 9:59:35 PM UTC+5:30 Wasim Niyaz Muns
Thanks.
It is working
Wasim
On Tuesday, October 11, 2022 at 7:02:05 PM UTC+5:30 bruno.t...@gmail.com
wrote:
> Sorry, it should be output.precision(12)
>
> Bruno
>
>
> Le mar. 11 oct. 2022 à 09:03, Wasim Niyaz Munshi ce21d400 <
> ce21...@smail.iitm.ac.in>
:
>
> Wasim,
> I imagine you are calling your setup_boundary_values() function before you
> call dof_handler.distribute_dofs()?
> Best
> W.
>
> On 10/11/22 05:16, Wasim Niyaz Munshi ce21d400 wrote:
> > *** Caution: EXTERNAL Sender ***
> >
> > I was using
ot;solution.vtk");
> output.set_precisition(12);
> data_out.write_vtk(output);
>
> Best,
>
> Bruno
>
> Le mar. 11 oct. 2022 à 00:48, Wasim Niyaz Munshi ce21d400 <
> ce21...@smail.iitm.ac.in> a écrit :
>
>> I used this to set precision of 12 for my solution vector
>
cts, constraints
and boundary_values?
Thank you
On Mon, Oct 10, 2022 at 7:28 PM Wolfgang Bangerth
wrote:
> On 10/10/22 07:42, Wasim Niyaz Munshi ce21d400 wrote:
> > I am trying to solve the elasticity equation on a square domain (similar
> to
> > step 8).
> > For the p
data_out.add_data_vector(*std::setprecision(12)* <’)
Thanks
On Tuesday, October 11, 2022 at 1:27:54 AM UTC+5:30 Wolfgang Bangerth wrote:
> On 10/10/22 11:45, Wasim Niyaz Munshi ce21d400 wrote:
> > So I used set precision function when using the cout to print the
> > solution vect
format.
Hope it will work.
Thanks
Wasim Niyaz
On Mon, 10 Oct, 2022, 7:22 pm Wolfgang Bangerth,
wrote:
> On 10/10/22 07:48, Wasim Niyaz Munshi ce21d400 wrote:
> >
> > I used the setprecision function to output the result to desired
> accuracy.
>
> Can you show how you actuall
Thank you.
I used the setprecision function to output the result to desired accuracy.
However, I noticed that deal.ii writes the output to vtk file only upto
certain precision(6).
Is there a way to write output to the vtk file according to required
precision?
Thank you
On Thursday, September
Hi everyone.
I am trying to solve the elasticity equation on a square domain (similar to
step 8).
For the patch test I am using just one element.
I am having problems with the application of boundary conditions.
I want to fix the lower left node in both x and y, and also fix the lower
right node
Hello everyone.
Is there a way to increase precision of floating point values in deal.ii? I
am modifying step3 and checking my solution against a linear exact
solution. I am getting machine precision (around 10^-16 ) for 4,16...4096
elements. However, when I refine further, I am getting an error
ble>) (double&, int)’
Functions::ConstantFunction(dim)(u_x_values,dim),constraints,u_x_mask);
It gives the same error for u_y also.
On Mon, Sep 19, 2022 at 3:09 AM Wolfgang Bangerth
wrote:
> On 9/18/22 06:51, Wasim Niyaz Munshi ce21d400 wrote:
> > VectorTools::interpolate_boundary
Hello everyone.
I am trying to specify x and y displacements separately in step 8 of the
tutorial programs.
I am doing it as follows:
FEValuesExtractors::Scalar u_x(0);
FEValuesExtractors::Scalar u_y(1);
ComponentMask u_x_mask = fe.component_mask(u_x);
ComponentMask u_y_mask =
43 matches
Mail list logo