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
Email:[email protected]
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, [email protected] 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  <tel:+49%20511%207623360>
        Email:[email protected]
        www:http://www.ifam.uni-hannover.de/wick
        www:http://www.cmap.polytechnique.fr/~wick/  
<http://www.cmap.polytechnique.fr/%7Ewick/>
        ++--------------------------------------------++
--
        On 02/03/2018 03:24 AM, [email protected] 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/
<http://www.math.clemson.edu/%7Eheister/>
        -- 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 [email protected]. 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
    [email protected]. 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 [email protected]. 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 [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/52c95f4f-1cba-766a-6a8a-a4f47897d014%40gmail.com.

Reply via email to