Dear Timo Heister Thomas Wick,
I have been reading your paper about phase-field crack propagation and your
code recently. May I ask why your index of assembly subroutine has the form
of local_matrix( j, i ) rather than the normal form of local_matrix( i, j
)?
Thanks a lots
On Tuesday, Febru
ay, August 26, 2024 at 8:09:47 PM UTC+8 Thomas Wick wrote:
> Dear Junxiang Wang,
>
> the reason is that we are dealing with a non-symmetric problem and the
> test function determines the row, while trial functions
> determine columns of the system matrix.
>
> The quickest and
Dear Prof. Timo Heister and Wolfgang Bangerth
I am trying to solve my multiphysics problem by deal.ii and it's
paralization.
I started from STEP-40 and encountered a problem regarding paralization and
MPI.
The program didn't proceed after I made a subroutine called
laplace_cel
)
10 {
9 double g{0};
8
7 double x = cell->center()[0] - 0.5;
6 double y = cell->center()[1] - 0.5;
5
4 double r = std::sqrt(std::pow(x,2) + std::pow(y,2));
3 // definition of g
8
9 return g;
10
11 }
"""
Do you know any issue ma
Dear all,
I am modelling multyphysics problem in a paralle way.
The code runs perfectly in serial or small number of cores while fails to
converge once I run it on more number of cores say >8.
What could cause this problem?
Thanks a lot.
Junxiang.
--
The deal.II project is located at http:
the local_rhs when the hanging node is located on the shared
edge of two elements and also the two elements belong to different
processors respectively?
Best
Junxiang.
On Tuesday, February 11, 2025 at 3:32:47 PM UTC+8 Junxiang Wang wrote:
> Hi Bruno,
>
> Thanks a lot.
>
&g
Sorry, the node distribution is as follows.
On Tuesday, February 11, 2025 at 6:17:39 PM UTC+8 Junxiang Wang wrote:
> Hi Bruno,
>
> I further print information about the RHS of the two elements that have
> the abnormal residual value. It turns out the nodal values are all the sam
Dear Deal.II community,
I am using the CellDataStorage class to record my history state variable
and try to interpolate it onto my refined mesh. Therefore, the history data
can be accessible for my analysis of the refined mesh.
However, I checked the data on the quadrature point on the refined
Dear Prof. Bangerth
Thanks a lot for your advice.
I successfully compiled the test code snippets.
Best Regards,
Junxiang.
On Saturday, January 18, 2025 at 7:03:01 AM UTC+8 Wolfgang Bangerth wrote:
> On 1/17/25 05:54, Junxiang Wang wrote:
> >
> > I'd like to ask if you
Dear Deal.II community,
I install deal.II package through
'''
sudo -i
export REPO=ppa:ginggs/deal.ii-XXX-backports
apt-get update && apt-get install -y software-properties-common
add-apt-repository $REPO
apt-get update
apt-get install libdeal.ii-dev
apt-get install build-essential cmake ninja-
ang.
On Friday, January 17, 2025 at 1:06:12 PM UTC+8 Wolfgang Bangerth wrote:
> On 1/16/25 18:40, Junxiang Wang wrote:
> >
> > const std::vector> lqph =
> > quadrature_point_history.get_data(cell);
> >
> > to record some history variables (e.g., max and mi
Dear Deal.II community,
I'd like to ask for help exporting the history nodal value.
I am currently using the following
const std::vector> lqph =
quadrature_point_history.get_data(cell);
to record some history variables (e.g., max and min) at quadrature_point.
Is there a way to output the nod
In the first 80 steps, RHS and Stiffness are all the same. After several
step of refine mesh, the RHS then turns to different from the serial one.
On Saturday, February 8, 2025 at 7:24:07 PM UTC+8 Junxiang Wang wrote:
> Dear Bruno,
>
> Thanks a lot for your advice.
>
> I test
Dear Bruno,
Thanks a lot for your advice.
I tested it. The code runing in parellel has exactly the same norm of
stiffness mtx but different norm of RHS residual compared to running in
serial.
It is really frastrating that increasing number of cores can make a
difference...
Best,
Junxiang.
14 matches
Mail list logo
