# Sorry in optimization element\_OP plein 3d **URL:** https://community.freefem.org/t/sorry-in-optimization-element-op-plein-3d/1287 **Category:** General Discussion **Created:** [October 27, 2021, 8:03am UTC](https://community.freefem.org/t/sorry-in-optimization-element-op-plein-3d/1287 "2021-10-27T08:03:46Z") **Posts on this page:** 1 **Page:** 1
### Author: ![xiaojiemo](https://avatars.discourse-cdn.com/v4/letter/x/a4c791/32.png) [@xiaojiemo](https://community.freefem.org/u/xiaojiemo) #### Post date: [October 27, 2021, 8:03am UTC](https://community.freefem.org/t/sorry-in-optimization-element-op-plein-3d/1287/1 "2021-10-27T08:03:46Z")
dear all, I am solving two equation using newton method and meets problem in title. My whole procedure is below: solve fi([setaf,fi],[dsetaf,dfi],solver=“UMFPACK64”)= 96 : int3d(th2)(cnewton_dsetaf_nf\*((6_fi1^5-15_fi1^4+10_fi1^3)rawl+rawi-(6fi1^5-15_fi1^4+10_fi1^3)rawi)wf1(dx(setaf)+dy(setaf)+dz(setaf))) 97 : +int3d(th2)(cnewtondsetaf_nf\*((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)(cl-ci)setaf(30_fi1^4-60_fi1^3+30_fi1^2)_wf1_(dx(fi1)+dy(fi1)+dz(fi1))) 98 : -int3d(th2)(cnewton_dsetaf_nf\*((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)(cl-ci)setaf1(30_fi1^4-60_fi1^3+30_fi1^2)_wf1_(dx(fi1)+dy(fi1)+dz(fi1))) 99 : +int3d(th2)(cnewton_nf_((6_fi1^5-15_fi1^4+10_fi1^3)hl+(1-(6fi1^5-15_fi1^4+10_fi1^3))hi)(dx(dsetaf)dx(setaf)+dy(dsetaf)dy(setaf)+dz(dsetaf)dz(setaf))) 100 : -int3d(th2)(cnewtondfi(30fi1^4-60_fi1^3+30_fi1^2)nf((6_fi1^5-15_fi1^4+10_fi1^3)_rawl+(1-(6_fi1^5-15_fi1^4+10_fi1^3))_rawi)_(cl-ci)_(log(setaf1)-log(277.15)+2)setaf) 101 : +int3d(th2)(cnewtondfi_(30_fi1^4-60_fi1^3+30_fi1^2)nf((6_fi1^5-15_fi1^4+10_fi1^3)_rawl+(1-(6_fi1^5-15_fi1^4+10_fi1^3))_rawi)_(cl-ci)_(log(setaf1)-log(277.15)+2)setaf1) 102 : -int3d(th2)(cnewtondfi_(30_fi1^4-60_fi1^3+30_fi1^2)nf((6_fi1^5-15_fi1^4+10_fi1^3)_rawl+(1-(6_fi1^5-15_fi1^4+10_fi1^3))_rawi)_(sl-si)_seff_setaf/277.15) 103 : +int3d(th2)(cnewton_dfi_(30_fi1^4-60_fi1^3+30_fi1^2)nf((6_fi1^5-15_fi1^4+10_fi1^3)_rawl+(1-(6_fi1^5-15_fi1^4+10_fi1^3))_rawi)_(sl-si)_seff_setaf1/277.15) 104 : +int3d(th2)(cnewton_dsetaf_nf\*((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)3gamma(2+12_fi1^2-12_fi1)wf1(dx(fi1)+dy(fi1)+dz(fi1))fi/(widthrawl)) 105 : -int3d(th2)(cnewton_dsetaf_nf_((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)3gamma(2+12_fi1^2-12_fi1)wf1(dx(fi1)+dy(fi1)+dz(fi1))fi1/(widthrawl)) 106 : +int3d(th2)(cnewton_dsetaf_nf_((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)3gamma(2_fi1+4_fi1^3-6_fi1^2)_wf1_(dx(fi)+dy(fi)+dz(fi))/(width_rawl)) 107 : +int3d(th2)(cnewton_dsetaf_nf_((rawl-rawi)/(rawl_rawl))_(30_fi1^4-60_fi1^3+30_fi1^2)wf1(dx(fi1)+dy(fi1)+dz(fi1))6gamma_width\*(dx(fi)_dx(fi1)+dy(fi)dy(fi1)+dz(fi)dz(fi1))) 108 : +int3d(th2)(cnewtondsetafnf_((rawl-rawi)/(rawl_rawl))_(30_fi1^4-60_fi1^3+30_fi1^2)wf1(dx(fi1)+dy(fi1)+dz(fi1))3gamma_(2_fi1+4_fi1^3-6_fi1^2)/width) 109 : +int3d(th2)(cnewton_dsetaf_nf_((rawl-rawi)/(rawl_rawl))_(120_fi1^3-180_fi1^2+60_fi1)fiwf1_(dx(fi1)+dy(fi1)+dz(fi1))_3_gamma_width_(dx(fi1)_dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1))) 110 : -int3d(th2)(cnewtondsetafnf_((rawl-rawi)/(rawl_rawl))_(120_fi1^3-180_fi1^2+60_fi1)fi1wf1_(dx(fi1)+dy(fi1)+dz(fi1))_3_gamma_width_(dx(fi1)_dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1))) 111 : +int3d(th2)(cnewtondsetafnf_((rawl-rawi)/(rawl_rawl))_(30_fi1^4-60_fi1^3+30_fi1^2)wf1(dx(fi)+dy(fi)+dz(fi))3gamma_width\*(dx(fi1)_dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1))) 112 : +int3d(th2)(cnewtondsetafnf_((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)(cl-ci)_(setaf1-277.15)_fi_(120_fi1^3-180_fi1^2+60_fi1)wf1(dx(fi1)+dy(fi1)+dz(fi1))) 113 : -int3d(th2)(cnewton_dsetaf_nf_((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)(cl-ci)_(setaf1-277.15)_fi1_(120_fi1^3-180_fi1^2+60_fi1)wf1(dx(fi1)+dy(fi1)+dz(fi1))) 114 : +int3d(th2)(cnewton_dsetaf_nf_((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)(cl-ci)_(setaf1-277.15)_(30_fi1^4-60_fi1^3+30_fi1^2)_wf1_(dx(fi)+dy(fi)+dz(fi))) 115 : +int3d(th2)(cnewton_dsetaf_(120_fi1^3-180_fi1^2+60_fi1)finf_((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)(sl-si)seffwf1_(dx(fi1)+dy(fi1)+dz(fi1))) 116 : -int3d(th2)(cnewton_dsetaf_(120_fi1^3-180_fi1^2+60_fi1)fi1nf_((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)(sl-si)seffwf1_(dx(fi1)+dy(fi1)+dz(fi1))) 117 : +int3d(th2)(cnewton_dsetaf_(30_fi1^4-60_fi1^3+30_fi1^2)nf((6_fi1^5-15_fi1^4+10_fi1^3)_rawl+(1-(6_fi1^5-15_fi1^4+10_fi1^3))_rawi)_(sl-si)_seff_wf1\*(dx(fi)+dy(fi)+dz(fi))) 118 : +int3d(th2)(cnewton_dfi_wf1\*(dx(fi)+dy(fi)+dz(fi))/(vkf/(6_rawl_width\*(sl-si)))) 119 : -int3d(th2)(cnewton_dfi_(120_fi1^3-180_fi1^2+60_fi1)finf_((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)(cl-ci)_(setaf1_log(setaf1/277.15)+setaf1-277.15)) 120 : +int3d(th2)(cnewton_dfi\*(120_fi1^3-180_fi1^2+60_fi1)fi1nf_((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)(cl-ci)_(setaf1_log(setaf1/277.15)+setaf1-277.15)) 121 : -int3d(th2)(cnewton_dfi\*(120_fi1^3-180_fi1^2+60_fi1)fi((rawl-rawi)/rawl)_(3_gamma_width\*(dx(fi1)_dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1))+(pow(fi1,2)pow(1-fi1,2)3gamma/width))) 122 : +int3d(th2)(cnewtondfi(120fi1^3-180_fi1^2+60_fi1)fi1((rawl-rawi)/rawl)_(3_gamma_width\*(dx(fi1)_dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1))+(pow(fi1,2)pow(1-fi1,2)3gamma/width))) 123 : -int3d(th2)(cnewtondfi(30fi1^4-60_fi1^3+30_fi1^2)_((rawl-rawi)/rawl)_(9_gamma_width_(dx(fi)_dx(fi1)+dy(fi)dy(fi1)+dz(fi)dz(fi1))+fi((2fi1+4_fi1^3-6_fi1^2)3gamma/width))) 124 : +int3d(th2)(cnewton_dfi\*(30_fi1^4-60_fi1^3+30_fi1^2)_((rawl-rawi)/rawl)_(fi1_((2_fi1+4_fi1^3-6_fi1^2)3gamma/width))) 125 : +int3d(th2)(cnewton_dfi\*(30_fi1^4-60_fi1^3+30_fi1^2)_((rawl-rawi)/rawl)_(9_gamma_width_(dx(fi)_dx(fi1)+dy(fi)dy(fi1)+dz(fi)dz(fi1)))) 126 : -int3d(th2)(cnewtondfi3_gamma_fi_(2+12_fi1^2-12_fi1)/width) 127 : +int3d(th2)(cnewton_dfi_3_gamma_fi1\*(2+12_fi1^2-12_fi1)/width) 128 : -int3d(th2)(cnewton_6_gamma_width_(dx(dfi)_dx(fi)+dy(dfi)dy(fi)+dz(dfi)dz(fi))) 129 : +int3d(th2)(dsetafnf((6_fi1^5-15_fi1^4+10_fi1^3)_rawl+(1-(6_fi1^5-15_fi1^4+10_fi1^3))_rawi)wf1(dx(setaf)+dy(setaf)+dz(setaf))) 130 : +int3d(th2)(dsetaf_nf\*((6_fi1^5-15_fi1^4+10_fi1^3)rawl+(1-(6fi1^5-15_fi1^4+10_fi1^3))rawi)3gamma(2_fi1+4_fi1^3-6_fi1^2)_wf1_(dx(fi1)+dy(fi1)+dz(fi1))/(width_rawl)) 131 : +int3d(th2)(dsetaf_nf\*((rawl-rawi)/(rawl_rawl))_(30_fi1^4-60_fi1^3+30_fi1^2)wf1(dx(fi1)+dy(fi1)+dz(fi1))3gamma_width\*(dx(fi1)_dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1))) 132 : +int3d(th2)(dsetafnf((6_fi1^5-15_fi1^4+10_fi1^3)_rawl+(1-(6_fi1^5-15_fi1^4+10_fi1^3))_rawi)_(cl-ci)_(setaf1-277.15)_(30_fi1^4-60_fi1^3+30_fi1^2)wf1(dx(fi1)+dy(fi1)+dz(fi1))) 133 : +int3d(th2)(dfi_(30_fi1^4-60_fi1^3+30_fi1^2)nf((6_fi1^5-15_fi1^4+10_fi1^3)_rawl+(1-(6_fi1^5-15_fi1^4+10_fi1^3))_rawi)_(sl-si)_seff_wf1\*(dx(fi1)+dy(fi1)+dz(fi1))) 134 : +int3d(th2)(nf\*((6_fi1^5-15_fi1^4+10_fi1^3)hl+(1-(6fi1^5-15_fi1^4+10_fi1^3))hi)(dx(dsetaf)dx(setaf1)+dy(dsetaf)dy(setaf1)+dz(dsetaf)dz(setaf1))) 135 : +int3d(th2)(dfiwf1(dx(fi1)+dy(fi1)+dz(fi1))/(vkf/(6rawl_width\*(sl-si)))) 136 : -int3d(th2)(dfi\*(30_fi1^4-60_fi1^3+30_fi1^2)nf((6_fi1^5-15_fi1^4+10_fi1^3)_rawl+(1-(6_fi1^5-15_fi1^4+10_fi1^3))_rawi)_(cl-ci)_(setaf1_log(setaf1/277.15)+setaf1-277.15)) 137 : -int3d(th2)(dfi\*(30_fi1^4-60_fi1^3+30_fi1^2)nf((6_fi1^5-15_fi1^4+10_fi1^3)_rawl+(1-(6_fi1^5-15_fi1^4+10_fi1^3))_rawi)_(sl-si)_seff_setaf/277.15) 138 : +int3d(th2)(dfi\*(30_fi1^4-60_fi1^3+30_fi1^2)nf((6_fi1^5-15_fi1^4+10_fi1^3)_rawl+(1-(6_fi1^5-15_fi1^4+10_fi1^3))_rawi)_(sl-si)_seff) 139 : -int3d(th2)(dfi_(30_fi1^4-60_fi1^3+30_fi1^2)_((rawl-rawi)/rawl)_(3_gamma_width_(dx(fi1)_dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1))+(pow(fi1,2)pow(1-fi1,2)3gamma/width))) 140 : +int3d(th2)(dfi3gamma(2_fi1+4_fi1^3-6_fi1^2)/width) 141 : +int3d(th2)(6_gamma_width\*(dx(dfi)_dx(fi1)+dy(dfi)dy(fi1)+dz(dfi)dz(fi1))) 142 : -int2d(th2,upp,down)(dsetafN.x(nf_((6_fi1^5-15_fi1^4+10_fi1^3)hl+(1-(6fi1^5-15_fi1^4+10_fi1^3))hi)(dx(setaf1)))+dsetaf_N.y\*(nf\*((6_fi1^5-15_fi1^4+10_fi1^3)hl+(1-(6fi1^5-15_fi1^4+10_fi1^3))hi)(dy(setaf1)))+dsetaf_N.z\*(nf\*((6_fi1^5-15_fi1^4+10_fi1^3)hl+(1-(6fi1^5-15_fi1^4+10_fi1^3))hi)(dz(setaf1)))); 143 : //-int2d(th2,upp)(dsetaf_N.y(nf\*((6_fi1^5-15_fi1^4+10_fi1^3)hl+(1-(6fi1^5-15_fi1^4+10_fi1^3))hi)(dy(setaf1)))) 144 : //-int2d(th2,upp)(dsetaf_N.z(nf\*((6_fi1^5-15_fi1^4+10_fi1^3)hl+(1-(6fi1^5-15_fi1^4+10\*fi1^3))_hi)_(dz(setaf1)))); 145 : 146 : //+on(upp,down,setaf=263.15); 147 : setaf1[]-=setaf[]; 148 : 149 : err=setaf1[].linfty; 150 : if(err\< 1e-2) cnewton=1; 151 : plot(setaf); 152 : //plot(fi); 153 : 154 : cout \<\