///////////////////////////////////////////////////////////////////
// Mladen Jurak
// Praktikum primijenjene matematike 2, ak. god. 2004/05
// FreeFEM++
//
//////////////////////////////////////////////////////////////
// saveMedit.edp
// Ispis rješenja u datoteku tipa medit
//////////////////////////////////////////////////////////////

func int saveMedit(mesh Th, real[int] & u, string & ime)
{
    //  Th je triangulacija 
    //  u je vektor rješenja. Ako je 
    //  Vh uu; 
    //  rješenje problema, onda funkciji treba dati uu[].
    //  ime je bazno ime datoteke u koju treba spremiti rješenje.
    //  Kreira datoteke: 
    //     ime.mesh ime.mesh.gmsh  ime.bb
    savemesh(Th,ime+".mesh");

    ofstream file(ime+".bb");

    file << "2 1 1 " << u.n<< " 2"<<endl;

    for(int j=0;j<u.n;j++) file << u[j] << endl;

    return 1;
}