MeshValidSoup.java

/* jCAE stand for Java Computer Aided Engineering. Features are : Small CAD
   modeler, Finite element mesher, Plugin architecture.

    Copyright (C) 2006, by EADS CRC
    Copyright (C) 2008, by EADS France

    This library is free software; you can redistribute it and/or
    modify it under the terms of the GNU Lesser General Public
    License as published by the Free Software Foundation; either
    version 2.1 of the License, or (at your option) any later version.

    This library is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    Lesser General Public License for more details.

    You should have received a copy of the GNU Lesser General Public
    License along with this library; if not, write to the Free Software
    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA */


package org.jcae.mesh;

import org.jcae.mesh.amibe.ds.Mesh;
import org.jcae.mesh.amibe.ds.Vertex;
import org.jcae.mesh.oemm.*;
import java.util.logging.Logger;

/**
 * This class illustrates how to perform quality checks.
 */
public class MeshValidSoup
{
	private static final Logger logger=Logger.getLogger(MeshValidSoup.class.getName());

	public static class ComputeTriangleQuality implements RawStorage.SoupReaderInterface
	{
		final Vertex [] n = new Vertex[3];
		int nrgroup = 0;
		final int maxgroup = 5000;
		final int split = 10;
		final int [][] nr = new int[maxgroup][split+1];
		final double [] amin = new double[maxgroup];
		final Mesh mesh = new Mesh();
		public ComputeTriangleQuality()
		{
			for (int i = 0; i < maxgroup; i++)
				amin[i] = 1.0;
		}
		public void processVertex(int i, double [] xyz)
		{
			n[i] = mesh.createVertex(xyz[0], xyz[1], xyz[2]);
		}
		public void processTriangle(int group)
		{
			if (group > nrgroup)
				nrgroup = group;
			double a1 = Math.abs(n[0].angle3D(n[1], n[2]));
			double a2 = Math.abs(n[1].angle3D(n[2], n[0]));
			double a3 = Math.abs(n[2].angle3D(n[0], n[1]));
			if (a2 < a1)
				a1 = a2;
			if (a3 < a1)
				a1 = a3;
			a1 *= 3.0 / Math.PI;
			if (a1 < amin[group])
				amin[group] = a1;
			int index = (int) (a1 * split);
			nr[group][index]++;
		}
		public final void printStats()
		{
			double amintot = amin[1];
			for (int g = 1; g <= nrgroup; g++)
			{
				System.out.println("Group "+g+" Minimum angle: "+(60.0*amin[g]));
				if (amin[g] < amintot)
					amintot = amin[g];
				for (int i = 0; i < split; i++)
					System.out.println("Slice: "+i+" "+nr[g][i]);
				if (g > 1)
					for (int i = 0; i < split; i++)
						nr[1][i] += nr[g][i];
			}
			System.out.println("Minimum angle: "+(60.0*amintot));
			for (int i = 0; i < split; i++)
				System.out.println("Slice: "+i+" "+nr[1][i]);
		}
	}

	private static void check(String brepfilename, String xmlDir, float discr, float defl)
	{
		logger.info("Reading triangle soup");
		ComputeTriangleQuality ctq = new ComputeTriangleQuality();
		RawStorage.readSoup(xmlDir+java.io.File.separator+"soup", ctq);
		ctq.printStats();
	}

	/**
	 * main method, reads 2 arguments and calls mesh() method
	 * @param args  an array of String, filename, algorithm type and constraint value
	 */
	public static void main(String args[])
	{
		if (args.length < 2)
		{
			System.out.println("Usage : MeshValid brep directory");
			System.exit(0);
		}
		String filename=args[0];
		String xmlDir = args[1];
		Float discr = new Float(args[2]);
		Float defl = new Float(args[3]);
		check(filename, xmlDir, discr.floatValue(), defl.floatValue());
	}
}