# Phase field fracture Compile error

**URL:** https://community.freefem.org/t/phase-field-fracture-compile-error/1241
**Category:** General Discussion
**Created:** [October 3, 2021, 8:35am UTC](https://community.freefem.org/t/phase-field-fracture-compile-error/1241 "2021-10-03T08:35:56Z")
**Posts on this page:** 1
**Page:** 1

<div class="post-metadata">

### Author: ![nansihu](https://avatars.discourse-cdn.com/v4/letter/n/3ab097/32.png) [@nansihu](https://community.freefem.org/u/nansihu)
#### Post date: [October 3, 2021, 8:35am UTC](https://community.freefem.org/t/phase-field-fracture-compile-error/1241/1 "2021-10-03T08:35:56Z")

</div>

I want to finish phase field simulation by Freefem++, but and there is awalys a complile error in the following code: error operator - , \<10LinearCombI7MGauche4C\_F0E\>

macro u [u0,u1] //eom  
macro v [v0,v1] //eom  
macro uold [uold0, uold1] //eom  
macro uvold [uvold0, uvold1]//  
macro du(ph) ((1-ph)_(1-ph))//  
macro tau(u) ((dx(u[1])+dy(u[0]))/2)//  
macro epsilon(u) [dx(u[0]),tau(u),tau(u),dy(u[1])]//  
macro div(u) (dx(u[0])+dy(u[1]))// tr(u)  
macro Sxx(u) (lamda_div(u)+2_mu_dx(u[0]))//  
macro Syy(u) (lamda_div(u)+2_mu_dy(u[1]))//  
macro Sxy(u) (2_mu_tau(u)) //  
macro Syx(u) (2_mu_tau(u)) //  
macro sigma(u) [Sxx(u),Sxy(u),Syx(u),Syy(u)]//  
macro absmus(u) ((div(u)-abs(div(u)))_(div(u)-abs(div(u))))//  
macro dW(u,v) (sigma(u)'_epsilon(v))//  
macro psi(u) (0.5_dW(u,u)-0.125_sum_absmus(u))//  
macro sT(u) [(Sxx(u)\*N.x,Sxy(u)\*N.y),(Syx(u)\*N.x+Syy(u)\*N.y)]//  
fespace Vh(Th,[P2,P2]);  
fespace Xh(Th,P2);  
Vh u,v,uold,uvold;  
Xh ph,w,phold;

//----------------------------------------  
// weak form

problem phasefield(ph,w,solver=CG,init=1)  
=int2d(Th)(Gc_l_(dx(ph)_dx(w)+dy(ph)dy(w)))  
+int2d(Th)(Gc/lph_w-2.0\*psi(u)_w_(1-ph));

problem phasefield(ph,w,solver=CG,init=1)  
=int2d(Th)(Gc_l_(dx(ph)_dx(w)+dy(ph)dy(w)))+int2d(Th)(Gc/lph_w-2.0_psi(u) (0.5_dW(u,u) (sigma(u) [Sxx(u) (lamda_div(u) (dx(u [u0,u1] [0])+dy(u [u0,u1] [1]))+2_mu_dx(u [u0,u1] [0])),Sxy(u) (2_mu_tau(u) ((dx(u [u0,u1] [1])+dy(u [u0,u1] [0]))/2)) ,Syx(u) (2_mu_tau(u) ((dx(u [u0,u1] [1])+dy(u [u0,u1] [0]))/2)) ,Syy(u) (lamda_div(u) (dx(u [u0,u1] [0])+dy(u [u0,u1] [1]))+2_mu_dy(u [u0,u1] [1]))]'_epsilon(u) [dx(u [u0,u1] [0]),tau(u) ((dx(u [u0,u1] [1])+dy(u [u0,u1] [0]))/2),tau(u) ((dx(u [u0,u1] [1])+dy(u [u0,u1] [0]))/2),dy(u [u0,u1] [1])])-0.125_sum_absmus(u) ((div(u) (dx(u [u0,u1] [0])+dy(u [u0,u1] [1]))-abs(div(u) (dx(u [u0,u1] [0])+dy(u [u0,u1] [1]))))_(div(u) (dx(u [u0,u1] [0])+dy(u [u0,u1] [1]))-abs(div(u)  
(dx(u [u0,u1] [0])+dy(u [u0,u1] [1]))))))_w_(1-ph) error operator - , \<10LinearCombI7MGauche4C\_F0E\>