DerivedMesh *cloth_convert_to_triangle(Cloth *cloth, DerivedMesh *dm)
{
	DerivedMesh *result = NULL;
	int i;
	int numverts = dm->getNumVerts(dm);
	int numedges = dm->getNumEdges(dm);
	int numfaces = dm->getNumFaces(dm);
	
	MVert *mvert = CDDM_get_verts(dm);
	MEdge *medge = CDDM_get_edges(dm);
	MFace *mface = CDDM_get_faces(dm);

	MVert *mvert2;
	MFace *mface2;
	unsigned int numtris=0;
	unsigned int numquads=0;
	int a = 0;

	for(i = 0; i < numfaces; i++)
	{
		if(mface[i].v4)
			numquads++;
		else
			numtris++;	
	}

	result = CDDM_new(numverts, 0, numtris + 2*numquads);
	
	if(!result)
		return NULL;

	// do verts
	mvert2 = CDDM_get_verts(result);
	for(a=0; a<numverts; a++) {
		MVert *mv = &mvert2[a];
		VECCOPY(mv->co, mvert[a].co);
	}
	
	// do faces
	mface2 = CDDM_get_faces(result);
	for(a=0, i=0; a<numfaces; a++) 
	{
		MFace *mf = &mface2[i];
		mf->v1 = mface[a].v1;
		mf->v2 = mface[a].v2;
		mf->v3 = mface[a].v3;
		mf->v4 = 0;
		test_index_face(mf, NULL, 0, 3);

		if(mface[a].v4)
		{
			i++;
			mf = &mface2[i];
			mf->v1 = mface[a].v4;
			mf->v2 = mface[a].v1;
			mf->v3 = mface[a].v3;
			mf->v4 = 0;
			test_index_face(mf, NULL, 0, 3);
		}

		i++;
	}

	CDDM_calc_edges(result);
	CDDM_calc_normals(result);
	
	return result;

}