// -*- mode:C++ ; compile-command: "g++ -I.. -g -c mathml.cc -DIN_GIAC -DHAVE_CONFIG_H" -*- #include "giacPCH.h" #include "gen.h" #include "symbolic.h" #include "usual.h" #include "intg.h" #include "derive.h" #include "series.h" #include "plot.h" #include "tex.h" #include "mathml.h" #include "maple.h" #include "giacintl.h" //#include #include #ifndef NSPIRE #include #if defined VISUALC13 && !defined BESTA_OS #undef clock #undef clock_t #endif #include #endif #include #ifdef HAVE_SSTREAM #include #else #ifndef NSPIRE #include #endif #endif using namespace std; #ifndef NO_NAMESPACE_GIAC namespace giac { #endif // ndef NO_NAMESPACE_GIAC #if defined GIAC_HAS_STO_38 || defined NSPIRE gen _mathml_error(const gen & g,GIAC_CONTEXT){ return gensizeerr(gettext("No mathml support")); } static const char _mathml_s []="mathml"; static define_unary_function_eval (__mathml,&_mathml_error,_mathml_s); define_unary_function_ptr5( at_mathml ,alias_at_mathml,&__mathml,0,true); static const char _spread2mathml_s []="spread2mathml"; static define_unary_function_eval (__spread2mathml,&_mathml_error,_spread2mathml_s); define_unary_function_ptr5( at_spread2mathml ,alias_at_spread2mathml,&__spread2mathml,0,true); #else //RTOS_THREADX #ifndef NO_NAMESPACE_GIAC using namespace giac; #endif // ndef NO_NAMESPACE_GIAC const char mathml_preamble[]="\n\n]>\n\n\n"; const char mathml_end[]="\n\n"; const char math_begin[]="

";
  const char math_end[]="

"; const char provisoire_mbox_begin[]=""; const char provisoire_mbox_end[]=""; const char svg_end[]=""; const double svg_epaisseur1=200; // thickness=(xmax-xmin)/svg_epaisseur1 struct svg_attribut { double ysurx; int color; int width; int epaisseur_point; int type_line; int type_point; bool fill_polygon; bool hidden_name; bool ie; // Internet explorer flag for SVG }; // 2 fonctions de service static void elimine_inf_sup(string &s){ string t=""; while (t!=s){ t=s; int i=s.find("<"), j=s.find(">"); if (i!=-1) s.replace(s.find("<"),1,"<"); if (j!=-1) s.replace(s.find(">"),1,">"); } } static void delete_cr(string & s){ for (unsigned int i=0 ; i","\n"}; string rep[]={"&","<",">",""}; //start with & before adding new ones for(int siz=0;siz"+s+""; } string spread2mathml(const matrice & m,int formule,GIAC_CONTEXT){ int l=m.size(); if (!l) return "

empty_spread

"; int c=m.front()._VECTptr->size(); string percent="", percent_first=""; if (c<=50){ percent="width=\""+print_INT_(200/(2*c+1))+"%\""; percent_first="width=\""+print_INT_(100/(2*c+1))+"%\""; } string s=""; string tmp; for (int k=0;;){ // write first line int i=k; for(int j=0;;++j){ tmp+=""; // tmp=""+tmp; i=i/26; if (!i) break; } // end for j ++k; if (k==c) break; } // end for k s += "" +tmp+""; for (int i=0;i"; for (int j=0;j" ; else { int save_r=printcell_current_row(contextptr),save_c=printcell_current_col(contextptr); printcell_current_row(contextptr)=i,printcell_current_col(contextptr)=j; string t(m[i][j][0].print(contextptr)),tt(""; printcell_current_row(contextptr)=save_r;printcell_current_col(contextptr)=save_c; } // end if (formule) } // end for j s +=""; } // end for i s += "

"; tmp+=char('A'+i%26-(j!=0)); tmp+="	"+char('A'+i%26-(j!=0))+"

"+ print_INT_(i)+"	"); int ll=t.size(); for (int l=0;l') tt += ">"; else if (t[l]=='&') tt += "&"; else tt += t[l]; } // end for l s += tt+"

"; return s; } string matrix2mathml(const matrice & m,GIAC_CONTEXT){ int l=m.size(); if (!l) return string("()"); int c=m.front()._VECTptr->size(); string s("("); for (int i=0;i" ; } s+= ""; } s += ")"; return s; } static string _SPOL12mathml(const sparse_poly1 & p, string &svg,GIAC_CONTEXT){ return gen2mathml(spol12gen(p,contextptr),svg,contextptr); } static string _VECT2mathml(const vecteur & v, unsigned char type,string &svg,GIAC_CONTEXT){ #if 1 string s(begin_VECT_string(type,false,contextptr)); vecteur::const_iterator it=v.begin(),itend=v.end(); for (;it!=itend;){ s += gen2mathml(*it,svg,contextptr); ++it; if (it!=itend) s += ","; } s+=end_VECT_string(type,false,contextptr); return s; #else string s(""); else { if (type==_SET__VECT) s.append("{\" close=\"}\">"); else{ if (type==_POLY1__VECT) s.append("poly1[\" close=\"]\">"); else s.append("[\" close=\"]\">"); } } s.append(""); vecteur::const_iterator it=v.begin(),itend=v.end(); for (;it!=itend;){ s += gen2mathml(*it,svg,contextptr); ++it; if (it!=itend) s += ","; } s +=""; return s; #endif } static bool needs_times(gen g,GIAC_CONTEXT){ //F.H: sqrt(2) est constant mais pas besoin de x. //if (eval(g,eval_level(contextptr),contextptr).is_constant()) gen gg=eval(g,eval_level(contextptr),contextptr); if( (gg.type==_INT_)||(gg.type==_ZINT) || (gg.type==_REAL) || (gg.type==_DOUBLE_) ) return true; /* F.H: a quoi cela sert il de mettre un x dans tous ces cas? else { if (g.type==_SYMB) { symbolic symb = *(g._SYMBptr); if ( !strcmp(symb.sommet.ptr()->print(contextptr),"/") ) return true; } }*/ return false; } static string prod_vect2mathml(const vecteur & v,GIAC_CONTEXT){ bool isprecINT = false; if (v.empty()) return "1"; vecteur::const_iterator it=v.begin(),itend=v.end(); string s; for (;;){ if ( (it->type==_SYMB) && ( (it->_SYMBptr->sommet==at_plus) || (it->_SYMBptr->sommet==at_neg))) s += string(" (")+gen2mathml(*it,contextptr)+string(" )"); else if (it->type==_CPLX){ if (!is_zero(re(*it,contextptr)) && !is_zero(im(*it,contextptr))) s += "("+gen2mathml(*it,contextptr)+")"; else{ s += gen2mathml(*it,contextptr); } } else{ if ( isprecINT && (it->type==_SYMB) && (it->_SYMBptr->sommet==at_pow) ) { symbolic mys = *(it->_SYMBptr); if(!((mys.feuille._VECTptr->back()==plus_one_half))){ if ((mys.feuille._VECTptr->front().type==_INT_)||(mys.feuille._VECTptr->front().type==_ZINT)){ if (is_positive(mys.feuille._VECTptr->front(),contextptr)) s += "×";// F.H: 2*2^n et 22^n, 2*2^(1/7) mais pas 2*sqrt(2). } } } s += gen2mathml(*it,contextptr); } if( (it->type==_INT_)||(it->type==_ZINT) || (it->type==_REAL) || (it->type==_DOUBLE_) ){ isprecINT=true; } else{ isprecINT=false; } ++it; if (it==itend) return s; else if ((needs_times(*it,contextptr))&& isprecINT) s += "×"; } } static string prod_vect2mathml_no_bra(const vecteur & v,GIAC_CONTEXT){ if (v.empty()) return "1"; vecteur::const_iterator it=v.begin(),itend=v.end(); bool isprecINT = false; if (v.size()==1) return gen2mathml(*it,contextptr); string s; for (;;){ if ( (it->type==_SYMB) && ( (it->_SYMBptr->sommet==at_plus) || (it->_SYMBptr->sommet==at_neg))) s += string(" (")+gen2mathml(*it,contextptr)+string(" )"); else if (it->type==_CPLX){ if (!is_zero(re(*it,contextptr)) && !is_zero(im(*it,contextptr))) s += "("+gen2mathml(*it,contextptr)+")"; else s += gen2mathml(*it,contextptr); } else if ( isprecINT && (it->type==_SYMB) && (it->_SYMBptr->sommet==at_pow) ) { symbolic mys = *(it->_SYMBptr); if(!((mys.feuille._VECTptr->back()==plus_one_half))){ if ((mys.feuille._VECTptr->front().type==_INT_)||(mys.feuille._VECTptr->front().type==_ZINT)){ if (is_positive(mys.feuille._VECTptr->front(),contextptr)) s += "×";//F.H 2*2^n ou 2*2^(1/7) != 22^(1/3) } } s += gen2mathml(*it,contextptr); } else s += gen2mathml(*it,contextptr); if( (it->type==_INT_)||(it->type==_ZINT) || (it->type==_REAL) || (it->type==_DOUBLE_) ){ isprecINT=true; } else{ isprecINT=false; } ++it; if (it==itend) return s; else if ((needs_times(*it,contextptr)) && isprecINT) s +="×"; } } // // --------------------------------- recopi?? de unary.cc --------------------------------------- static string mathml_printsommetasoperator(const gen & feuille, const string & sommetstr,GIAC_CONTEXT){ if (feuille.type!=_VECT) return ""+feuille.print(contextptr)+""; vecteur::const_iterator itb=feuille._VECTptr->begin(),itend=feuille._VECTptr->end(); if (itb==itend) return ""; string s; if (itb->type==_FRAC) s=gen2mathml(*itb,contextptr); else { if ( (itb->type==_IDNT) || ((itb->type!=_SYMB) && is_positive(*itb,contextptr)) || sommetstr=="=") s=gen2mathml(*itb,contextptr); else s="("+gen2mathml(*itb,contextptr)+")"; } ++itb; for (;;){ if (itb==itend) return s; if ((itb->type!=_SYMB) || sommetstr=="=") s += sommetstr + gen2mathml(*itb,contextptr); else s += sommetstr + "("+gen2mathml(*itb,contextptr)+")"; ++itb; } } // ------------------------------- fin provenant de unary.cc ------------------------------------- // // ------------------------------- provenant de usual.cc ---------------------------------- static string mathml_printasexp(const gen & g,GIAC_CONTEXT){ return "e"+gen2mathml(g,contextptr)+""; } static string mathml_printassqrt(const gen & g,GIAC_CONTEXT){ return ""+gen2mathml(g,contextptr)+""; } #if 0 //-------------- inutile ? ------------------------- static string mathml_printassqr(const gen & g,GIAC_CONTEXT){ return ""+gen2mathml(g,contextptr)+"2"; } #endif static string mathml_printasre(const gen & g,GIAC_CONTEXT){ return "Re("+gen2mathml(g,contextptr)+")"; } static string mathml_printasim(const gen & g,GIAC_CONTEXT){ return "Im("+gen2mathml(g,contextptr)+")"; } // overline en mathml ???????????? // string texprintasconj(const gen & g,const char & s){ // return "\\overline{"+gen2tex []=g+"}"; // } static string mathml_printassto(const gen & g,const string & sommetstr,GIAC_CONTEXT){ if ( (g.type!=_VECT) || (g._VECTptr->size()!=2) ) return sommetstr+"("+gen2mathml(g,contextptr)+")"; string s(gen2mathml(g._VECTptr->back(),contextptr)+":="); if (g._VECTptr->front().type==_SEQ__VECT) return s+"("+gen2mathml(g._VECTptr->front(),contextptr)+")"; else return s+gen2mathml(g._VECTptr->front(),contextptr); } // --------------------------- fin de provenant de usual.cc ---------------------------------- // n'existe pas avec tex : static string mathml_equal2arrow(gen &g,GIAC_CONTEXT){ if (!is_equal(g)) return gen2mathml(g,contextptr); vecteur v=*g._SYMBptr->feuille._VECTptr; return ""+gen2mathml(v[0],contextptr)+"→"+gen2mathml(v[1],contextptr)+""; } // ----------------------------- provenant de series.cc ---------------------------------- static string mathml_printaslimit(const gen & g,GIAC_CONTEXT){ string s("lim"); if (g.type!=_VECT) return s+gen2mathml(g,contextptr); vecteur v(*g._VECTptr); int l(v.size()); if (!l) return s; if (l==1) return s+gen2mathml(v[0],contextptr); if (l==2) return ""+s+mathml_equal2arrow(v[1],contextptr)+""+gen2mathml(eval(v[0],eval_level(contextptr),contextptr),contextptr); // directional limit if (l==3){ if (is_one(v[2])) return ""+s+""+mathml_equal2arrow(v[1],contextptr)+"+" +gen2mathml(eval(v[0],eval_level(contextptr),contextptr),contextptr); if (is_zero(eval(v[2]+1,eval_level(contextptr),contextptr))) return ""+s+""+mathml_equal2arrow(v[1],contextptr)+"-" +gen2mathml(eval(v[0],eval_level(contextptr),contextptr),contextptr); else return ""+s+mathml_equal2arrow(v[1],contextptr)+""+gen2mathml(eval(v[0],eval_level(contextptr),contextptr),contextptr); } return s; } // ----------------------------- fin provenant de series.cc ---------------------------------- // ----------------------------- provenant de intg.cc ---------------------- static string mathml_printasintegrate(const gen g,GIAC_CONTEXT){ string s("∫"); if (g.type!=_VECT) return s+gen2mathml(g,contextptr); vecteur v(*g._VECTptr); int l(v.size()); if (!l) return s; if (l==1) return s+gen2mathml(v.front(),contextptr); if (l==2) return s+gen2mathml(v.front(),contextptr)+"d"+gen2mathml(v.back(),contextptr); if (l==4) return ""+s+""+v[2].print(contextptr)+"" +""+v[3].print(contextptr)+""+gen2mathml(eval(v.front(),eval_level(contextptr),contextptr),contextptr) +"d"+v[1].print(contextptr)+""; return s; } // ---------------------------- fin provenant de intg.cc ---------------------- // --------------------------- provenant de derive ---------------------------- static string mathml_printasderive(const gen & feuille,GIAC_CONTEXT){ if (feuille.type!=_VECT || feuille._VECTptr->size()<2){ if (feuille.type==_IDNT || feuille.type<=_CPLX || feuille.is_symb_of_sommet(at_derive)) return gen2mathml(feuille,contextptr)+"'"; else return "("+gen2mathml(feuille,contextptr)+")'"; } vecteur & v = *feuille._VECTptr; bool needpar=v[0].type==_SYMB; if (v.size()>2) return "∂"+gen2mathml(v[2],contextptr)+""+(needpar?"(":"")+gen2mathml(v[0],contextptr) +(needpar?")":"")+"∂"+gen2mathml(v[1],contextptr)+""+gen2mathml(v[2],contextptr)+""; return string("∂")+(needpar?"(":"")+gen2mathml(v[0],contextptr) +(needpar?")":"")+"∂"+gen2mathml(v[1],contextptr)+""; } // ----------------------- fin provenant de derive --------------------------- static string mathml_print(const symbolic & mys,GIAC_CONTEXT){ unary_function_ptr u =mys.sommet; if (u==at_equal || u==at_equal2) return mathml_printsommetasoperator(mys.feuille,"=",contextptr); if (u==at_different) return mathml_printsommetasoperator(mys.feuille,"≠",contextptr); if (u==at_inferieur_egal) return mathml_printsommetasoperator(mys.feuille,"≤",contextptr); if (u==at_superieur_egal) return mathml_printsommetasoperator(mys.feuille,"≥",contextptr); if (u==at_inferieur_strict) return mathml_printsommetasoperator(mys.feuille,"<",contextptr); if (u==at_superieur_strict) return mathml_printsommetasoperator(mys.feuille,">",contextptr); if (u==at_re) return mathml_printasre(mys.feuille,contextptr); if (u==at_im) return mathml_printasim(mys.feuille,contextptr); if (u==at_exp) return mathml_printasexp(mys.feuille,contextptr); if (u==at_sqrt) return mathml_printassqrt(mys.feuille,contextptr); if (u==at_integrate) return mathml_printasintegrate(mys.feuille,contextptr); if (u==at_derive) return mathml_printasderive(mys.feuille,contextptr); if (u==at_sto) return mathml_printassto(mys.feuille, ":=",contextptr); if (u==at_limit) return mathml_printaslimit(mys.feuille,contextptr); string s=mys.print(contextptr); elimine_inf_sup(s); delete_cr(s); return ""+s+""; } static string mathml_printassum(const gen e,GIAC_CONTEXT){ vecteur v = *e._SYMBptr->feuille._VECTptr; if (v.size()!=4) return ""+e.print(contextptr)+""; return "Σ"+gen2mathml(v[1],contextptr)+"="+gen2mathml(v[2],contextptr)+""+gen2mathml(v[3],contextptr)+""+gen2mathml(eval(v[0],eval_level(contextptr),contextptr),contextptr); } static string mathml_printasabs(const gen e,GIAC_CONTEXT){ return "∣"+gen2mathml(e._SYMBptr->feuille,contextptr)+ "∣"; } static string mathml_printasunit(const gen e,GIAC_CONTEXT){ vecteur v = *e._SYMBptr->feuille._VECTptr; if (v.size()!=2) return ""+e.print(contextptr)+""; return ""+gen2mathml(eval(v[0],eval_level(contextptr),contextptr),contextptr)+""+gen2mathml(v[1],contextptr)+""; } //---------------- Zone SVG --------------- static void svg_dx_dy(double svg_width, double svg_height, double & dx, double & dy){ // calibrage des graduations int p_h=(int) std::log10(svg_height); int p_w=(int) std::log10(svg_width); dx=std::pow(10.0,p_w-1); dy=std::pow(10.0,p_h-1); if (dx==0) dx=1; if (dy==0) dy=1; if (svg_width/dx>25){ dx *= 5; } if (svg_height/dy>25){ dy *=5; } } void arc_en_ciel(int k,int & r,int & g,int & b){ k += 21; k %= 126; if (k<0) k += 126; if (k<21){ r=251; g=0; b=12*k; } if (k>=21 && k<42){ r=251-(12*(k-21)); g=0; b=251; } if (k>=42 && k<63){ r=0; g=(k-42)*12; b=251; } if (k>=63 && k<84){ r=0; g=251; b=251-(k-63)*12; } if (k>=84 && k<105){ r=(k-84)*12; g=251; b=0; } if (k>=105 && k<126){ r=251; g=251-(k-105)*12; b=0; } } static string color_string(int color){ switch (color) { case 0: return "black"; case _RED: return "red"; case _GREEN: return "green"; case _YELLOW: return "yellow"; case _BLUE: return "blue"; case _MAGENTA: return "magenta"; case _CYAN: return "cyan"; } int r,g,b; if (color>=0x100 && color<0x17e) color -= 0x100; arc_en_ciel(color,r,g,b); return "rgb("+print_INT_(r)+","+print_INT_(g)+","+print_INT_(b)+")"; } void pixon2svg(const gen & g,string & s){ if (g.type==_VECT){ vecteur w=merge_pixon(*g._VECTptr); const_iterateur it=w.begin(),itend=w.end(); for (;it!=itend;++it){ pixon2svg(*it,s); } return; } if (!g.is_symb_of_sommet(at_pnt)) return; gen h=remove_at_pnt(g); if (!h.is_symb_of_sommet(at_pixon)) return; h=h._SYMBptr->feuille; if (h.type!=_VECT || h._VECTptr->size()<2) return; vecteur & v=*h._VECTptr; gen x=pixon_size*v[0],y=pixon_size*v[1]; int color=0,w=0,he=0; if (v.size()>2) color=v[2].val; if (v.size()>3 && v[3].type==_INT_){ if (v[3].val>0) he=v[3].val; else w=-v[3].val; } if (w==0 && he==0){ char ch='0'+pixon_size; s += "\n"; return; } s += "\n"; } // before making a user transformation on the frame, // collect pixon instructions in g string svg_preamble_pixel(const gen &g,double svg_width_cm, double svg_height_cm,double xmin,double xmax,double ymin,double ymax,bool ortho,bool xml){ double svg_width=xmax-xmin; double svg_height=ymax-ymin; double x_scale=(xmax-xmin)/10; double y_scale=(ymax-ymin)/10; int i; char grid_color[]="cyan"; char grid_color2[]="blue"; char buffer[16384]; char * pos=buffer; if (ortho){ if (svg_width\n
"; } return res; } #else static string mathml_split(const string & s,int slicesize){ return s; } #endif static string indice_mathml(const string & s1,const string & s2){ return ""+s1+""+s2+""; } static string idnt2mathml_(const string & sorig){ string s0=sorig; int n=int(s0.size()),j; for (j=n-1;j>=2;--j){ if (s0[j]>32 && isalpha(s0[j])) break; } string s=s0.substr(0,j+1),sadd; if (j=1;--j){ if (j"+idnt2mathml_(s0)+""; } string gen2mathml(const gen &e, string &svg,GIAC_CONTEXT){ string part_re="", part_im="i"; if (e.type==_SYMB && e._SYMBptr->sommet==at_sum) return mathml_printassum(e,contextptr); if (e.type==_SYMB && e._SYMBptr->sommet==at_abs) return mathml_printasabs(e,contextptr); if (e.type==_SYMB && e._SYMBptr->sommet==at_unit) return mathml_printasunit(e,contextptr); switch (e.type){ case _INT_: case _ZINT: return ""+mathml_split(e.print(contextptr),50)+""; case _DOUBLE_: /* FH: This should be tuned in the context if (fabs(e._DOUBLE_val)<1.1e-5) return "0.0"; else */ return ""+e.print(contextptr)+""; case _REAL: return ""+mathml_split(e.print(contextptr),50)+""; case _CPLX: if (!is_zero(re(e,contextptr))) part_re=""+re(e,contextptr).print(contextptr)+""; if (!is_zero(im(e,contextptr))){ if (is_positive(im(e,contextptr),contextptr)){ if (!is_zero(re(e,contextptr))) part_re = part_re+"+"; } else part_re = part_re+"-"; } if (!is_one(-im(e,contextptr)) && ! is_one(im(e,contextptr))) part_im=""+abs(im(e,contextptr),contextptr).print(contextptr)+""+part_im; //the is_zero test should be the last one //Ex: 3+10.0**(-13)*i avec Digits 12 et 10.0**(-13)*i avec Digits 12 et if (is_zero(im(e,contextptr))){ part_im=""; if (is_zero(re(e,contextptr))) part_re="0.0"; } return part_re+part_im; case _IDNT: if (e==unsigned_inf) return "∞"; if (e==cst_pi) return "π"; if (e==undef) return "undef"; return idnt2mathml(e.print(contextptr)); case _SYMB: { gen tmp=aplatir_fois_plus(e); if (tmp.type!=_SYMB) return gen2mathml(tmp,svg,contextptr); return symbolic2mathml(*tmp._SYMBptr, svg,contextptr); } case _VECT: { vector V; int p=0; if (is_mod_vecteur(*e._VECTptr,V,p) && p!=0){ gen gm=makemodquoted(unmod(e),p); return gen2mathml(gm,svg,contextptr); } if (e.subtype==_SPREAD__VECT) return spread2mathml(*e._VECTptr,1,contextptr); //----------------v??rifier le 2??me param??tre if (e.subtype!=_SEQ__VECT && ckmatrix(*e._VECTptr)){ vector< vector > M; p=0; if (is_mod_matrice(*e._VECTptr,M,p) && p!=0){ gen gm=makemodquoted(unmod(e),p); return gen2mathml(gm,svg,contextptr); } return matrix2mathml(*e._VECTptr,contextptr); } return _VECT2mathml(*e._VECTptr,e.subtype, svg,contextptr); } case _SPOL1: return _SPOL12mathml(*e._SPOL1ptr,svg,contextptr); case _POLY: return string("polynome"); case _FRAC: return string("")+gen2mathml(e._FRACptr->num,contextptr)+"" +gen2mathml(e._FRACptr->den,contextptr)+""; case _EXT: return ""; case _STRNG: return string2mathml(*e._STRNGptr); // return ""+(*e._STRNGptr)+""; case _FUNC: case _MAP: return ""+e.print(contextptr)+""; case _USER: return string2mathml(e.print(contextptr)); // return ""+e._USERptr->texprint(contextptr)+""; case _MOD: return gen2mathml(*e._MODptr,contextptr)+"%"+gen2mathml(*(e._MODptr+1),contextptr); default: settypeerr(gettext("MathMl convert ")+e.print(contextptr)); } return "mathml error (gen2mathml)"; } string ingen2mathml(const gen & g,bool html5,GIAC_CONTEXT){ if (html5) return "\n

\n\n"+gen2mathml(g,contextptr)+"\n\n

\n"; return "\n

\n\n"+gen2mathml(g,contextptr)+"\n\n

\n"; } static string gen2mathmlfull(const gen & g,GIAC_CONTEXT){ return string(mathml_preamble)+ingen2mathml(g,false,contextptr)+mathml_end+'\n'; } unsigned max_prettyprint_equation=5000; gen _mathml(const gen & g,GIAC_CONTEXT){ if ( g.type==_STRNG && g.subtype==-1) return g; #ifndef EMCC if (g.type==_VECT && g.subtype==_SEQ__VECT && g._VECTptr->size()>1 && (*g._VECTptr)[1].type==_STRNG && *((*g._VECTptr)[1]._STRNGptr)!="Done"){ ofstream of((*g._VECTptr)[1]._STRNGptr->c_str()); of << gen2mathmlfull(g._VECTptr->front(),contextptr) << endl; return plus_one; } #endif unsigned ta=taille(g,max_prettyprint_equation); if (ta>max_prettyprint_equation) return string2gen("Expression_too_large",false); if (g.type==_VECT && g._VECTptr->size()==2 && g._VECTptr->back().type==_INT_) return string2gen(ingen2mathml(g._VECTptr->front(),g._VECTptr->back().val,contextptr),false); if (g.type==_VECT && g.subtype==_SEQ__VECT && g._VECTptr->size()>2 && g._VECTptr->back().type==_INT_){ vecteur v(g._VECTptr->begin(),g._VECTptr->end()-1); gen g1(v,g.subtype); return string2gen(ingen2mathml(g1,g._VECTptr->back().val,contextptr),false); } return string2gen(gen2mathmlfull(g,contextptr),false); } static const char _mathml_s []="mathml"; static define_unary_function_eval (__mathml,&_mathml,_mathml_s); define_unary_function_ptr5( at_mathml ,alias_at_mathml,&__mathml,0,true); static string spread2mathmlfull(const matrice m, int formule,GIAC_CONTEXT){ return mathml_preamble+spread2mathml(m, formule,contextptr)+mathml_end+'\n'; } gen _spread2mathml(const gen & g,GIAC_CONTEXT){ if ( g.type==_STRNG && g.subtype==-1) return g; if (g.type!=_VECT || g._VECTptr->size()<2) return string2gen(" syntax error "); vecteur &v=*g._VECTptr; if (g.type==_VECT && g.subtype==_SEQ__VECT && g._VECTptr->size()>2 && (*g._VECTptr)[2].type==_STRNG){ ofstream of((*g._VECTptr)[2]._STRNGptr->c_str()); of << spread2mathmlfull(*(v[0]._VECTptr), is_one(v[1]),contextptr) << endl; return plus_one; } return string2gen(spread2mathmlfull(*(v[0]._VECTptr) , is_one(v[1]),contextptr)); } static const char _spread2mathml_s []="spread2mathml"; static define_unary_function_eval (__spread2mathml,&_spread2mathml,_spread2mathml_s); define_unary_function_ptr5( at_spread2mathml ,alias_at_spread2mathml,&__spread2mathml,0,true); #endif // RTOS_THREADX #ifndef NO_NAMESPACE_GIAC } // namespace giac #endif // ndef NO_NAMESPACE_GIAC