Dear Prof. Thomas Wick, Thanks a lot for your answer. Now, I understand the purpose of your code.
However, In all the tutorials steps of deal.II , the test function is indexed by "i" loop, right? For example, in the form of local_matrix( i, j ). Best regards, Junxiang. On Monday, 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 easiest way is to switch those indices to achieve this. > > In general, actually, the notation local_matrix( j, i ) is the correct > one, while local_matrix( i, j ) only holds true if your problem > statement is symmetric. > > Best regards, > > Thomas W. > > -- > --- > +++++------------------------------------------+++++ > Prof. Dr. Thomas Wick > Managing Director of the > Institut für Angewandte Mathematik (IfAM) > > Leibniz Universität Hannover (LUH) > Welfengarten 1 > 30167 Hannover, Germany > > Tel.: +49 511 762 3360 <+49%20511%207623360> > Email: thoma...@ifam.uni-hannover.de > www: https://ifam.uni-hannover.de/wick > www: http://www.thomaswick.org > +++++------------------------------------------+++++ > --- > -- > > > > Am 26.08.24 um 14:02 schrieb Junxiang Wang: > > 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, February 20, 2018 at 3:32:06 PM UTC+8 Thomas Wick wrote: > >> Dear Yaakov, >> >> which article do you mean? Please give the exact reference including >> author names. >> >> I do not know a priori whether they have different material parameters, >> another >> stress-strain splitting, etc. >> >> The reason for different results can be anything. One needs to do a >> careful 1-by-1 >> comparison. >> >> Best regards, >> >> Thomas W. >> >> >> On 02/19/2018 09:26 AM, Thomas Wick wrote: >> >> Dear Yaakov, >> >> which article do you mean? Please give the exact reference including >> author names. >> >> I do not know a priori whether they have different material parameters, >> another >> stress-strain splitting, etc. >> >> The reason for different results can be anything. One needs to do a >> careful 1-by-1 >> comparison. >> >> Best regards, >> >> Thomas W. >> >> >> >> On 02/17/2018 01:41 AM, fea...@gmail.com wrote: >> >> Dear Prof. Wick, >> >> I have used isotrope formulation for miehe shear loading (without local >> refinement). I ca not see two crack branching which is described in the >> article (A review on phase-field models of brittle fractureand a new fast >> hybrid formulation) >> >> I attach the test results. >> >> Thanks for your answer! >> >> Kind regards, >> Yaakov >> >> On Saturday, February 3, 2018 at 7:20:27 AM UTC+1, Thomas Wick wrote: >>> >>> Hello, >>> >>> there are two flags in the parameter file that you need to change: >>> >>> set Decompose stress in rhs = 0.0 >>> set Decompose stress in matrix = 0.0 >>> >>> >>> Thomas W. >>> >>> >>> -- >>> ++--------------------------------------------++ >>> Prof. Dr. Thomas Wick >>> Institut für Angewandte Mathematik (IfAM) >>> Leibniz Universität Hannover >>> Welfengarten 1 >>> 30167 Hannover, Germany >>> >>> Tel.: +49 511 762 3360 <+49%20511%207623360> >>> Email: thoma...@ifam.uni-hannover.de >>> www: http://www.ifam.uni-hannover.de/wick >>> www: http://www.cmap.polytechnique.fr/~wick/ >>> ++--------------------------------------------++ >>> -- >>> >>> On 02/03/2018 03:24 AM, fea...@gmail.com wrote: >>> >>> Dear Prof. Heister, I would like to just test isotrope formulation of >>> phase -field model (no compression/tension modification), how could I >>> modify the codes (in a simple way)? I am sorry that I am just a beginner of >>> DealII. Best On Tuesday, December 5, 2017 at 7:01:17 PM UTC+1, Timo Heister >>> wrote: >>>> >>>> > If I use your user element, I have to use OPEN MPI? now I have some >>>> issues > with Open MPI in Deal.ii What do you mean by "user element"? The >>>> example code in question requires deal.II to be configured with MPI. What >>>> vendor you use (OpenMPI, MPICH, ...) is up to you. -- Timo Heister >>>> http://www.math.clemson.edu/~heister/ >>> >>> -- The deal.II project is located at http://www.dealii.org/ For mailing >>> list/forum options, see https://groups.google.com/d/forum/dealii?hl=en >>> --- You received this message because you are subscribed to the Google >>> Groups "deal.II User Group" group. To unsubscribe from this group and stop >>> receiving emails from it, send an email to dealii+un...@googlegroups.com. >>> For more options, visit https://groups.google.com/d/optout. >>> >>> -- The deal.II project is located at http://www.dealii.org/ For mailing >> list/forum options, see https://groups.google.com/d/forum/dealii?hl=en >> --- You received this message because you are subscribed to the Google >> Groups "deal.II User Group" group. To unsubscribe from this group and stop >> receiving emails from it, send an email to dealii+un...@googlegroups.com. >> For more options, visit https://groups.google.com/d/optout. >> >> -- > The deal.II project is located at http://www.dealii.org/ > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > --- > You received this message because you are subscribed to the Google Groups > "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to dealii+un...@googlegroups.com. > > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/686a114f-c0d4-4f1b-acf2-8071159a7597n%40googlegroups.com > > <https://groups.google.com/d/msgid/dealii/686a114f-c0d4-4f1b-acf2-8071159a7597n%40googlegroups.com?utm_medium=email&utm_source=footer> > . > > > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/97931b7b-eea9-4b3f-9860-32a9214cce98n%40googlegroups.com.
