home *** CD-ROM | disk | FTP | other *** search
/ DOS/V Power Report 1998 February / VPR9802A.ISO / APP_DEMO / VC / SAMPLES.BIN / XYZApp.java < prev    next >
Text File  |  1997-10-27  |  13KB  |  489 lines
  1. /*
  2.  * @(#)XYZApp.java    1.3 96/12/06
  3.  *
  4.  * Copyright (c) 1994-1996 Sun Microsystems, Inc. All Rights Reserved.
  5.  *
  6.  * Sun grants you ("Licensee") a non-exclusive, royalty free, license to use,
  7.  * modify and redistribute this software in source and binary code form,
  8.  * provided that i) this copyright notice and license appear on all copies of
  9.  * the software; and ii) Licensee does not utilize the software in a manner
  10.  * which is disparaging to Sun.
  11.  *
  12.  * This software is provided "AS IS," without a warranty of any kind. ALL
  13.  * EXPRESS OR IMPLIED CONDITIONS, REPRESENTATIONS AND WARRANTIES, INCLUDING ANY
  14.  * IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE OR
  15.  * NON-INFRINGEMENT, ARE HEREBY EXCLUDED. SUN AND ITS LICENSORS SHALL NOT BE
  16.  * LIABLE FOR ANY DAMAGES SUFFERED BY LICENSEE AS A RESULT OF USING, MODIFYING
  17.  * OR DISTRIBUTING THE SOFTWARE OR ITS DERIVATIVES. IN NO EVENT WILL SUN OR ITS
  18.  * LICENSORS BE LIABLE FOR ANY LOST REVENUE, PROFIT OR DATA, OR FOR DIRECT,
  19.  * INDIRECT, SPECIAL, CONSEQUENTIAL, INCIDENTAL OR PUNITIVE DAMAGES, HOWEVER
  20.  * CAUSED AND REGARDLESS OF THE THEORY OF LIABILITY, ARISING OUT OF THE USE OF
  21.  * OR INABILITY TO USE SOFTWARE, EVEN IF SUN HAS BEEN ADVISED OF THE
  22.  * POSSIBILITY OF SUCH DAMAGES.
  23.  *
  24.  * This software is not designed or intended for use in on-line control of
  25.  * aircraft, air traffic, aircraft navigation or aircraft communications; or in
  26.  * the design, construction, operation or maintenance of any nuclear
  27.  * facility. Licensee represents and warrants that it will not use or
  28.  * redistribute the Software for such purposes.
  29.  */
  30.  
  31. /*
  32.  * A set of classes to parse, represent and display Chemical compounds in
  33.  * .xyz format (see http://chem.leeds.ac.uk/Project/MIME.html)
  34.  */
  35.  
  36. import java.applet.Applet;
  37. import java.awt.Image;
  38. import java.awt.Event;
  39. import java.awt.Graphics;
  40. import java.awt.Dimension;
  41. import java.io.StreamTokenizer;
  42. import java.io.InputStream;
  43. import java.io.BufferedInputStream;
  44. import java.io.IOException;
  45. import java.net.URL;
  46. import java.util.Hashtable;
  47. import java.awt.image.IndexColorModel;
  48. import java.awt.image.ColorModel;
  49. import java.awt.image.MemoryImageSource;
  50.  
  51.  
  52. /** The representation of a Chemical .xyz model */
  53. class XYZChemModel {
  54.     float vert[];
  55.     Atom atoms[];
  56.     int tvert[];
  57.     int ZsortMap[];
  58.     int nvert, maxvert;
  59.  
  60.     static Hashtable atomTable = new Hashtable();
  61.     static Atom defaultAtom;
  62.     static {
  63.     atomTable.put("c", new Atom(0, 0, 0));
  64.     atomTable.put("h", new Atom(210, 210, 210));
  65.     atomTable.put("n", new Atom(0, 0, 255));
  66.     atomTable.put("o", new Atom(255, 0, 0));
  67.     atomTable.put("p", new Atom(255, 0, 255));
  68.     atomTable.put("s", new Atom(255, 255, 0));
  69.     atomTable.put("hn", new Atom(150, 255, 150)); /* !!*/
  70.     defaultAtom = new Atom(255, 100, 200);
  71.     }
  72.  
  73.     boolean transformed;
  74.     Matrix3D mat;
  75.  
  76.     float xmin, xmax, ymin, ymax, zmin, zmax;
  77.  
  78.  
  79.     XYZChemModel () {
  80.     mat = new Matrix3D();
  81.     mat.xrot(20);
  82.     mat.yrot(30);
  83.     }
  84.  
  85.  
  86.     /** Create a Cehmical model by parsing an input stream */
  87.     XYZChemModel (InputStream is) throws Exception
  88.     {
  89.        this();
  90.        StreamTokenizer st;
  91.        st = new StreamTokenizer(new BufferedInputStream(is, 4000));
  92.        st.eolIsSignificant(true);
  93.        st.commentChar('#');
  94.        int slot = 0;
  95.       
  96.        try
  97.        {
  98. scan:
  99.           while (true)
  100.           {
  101.              switch ( st.nextToken() ) 
  102.              {
  103.                 case StreamTokenizer.TT_EOF:
  104.                    break scan;
  105.                 default:
  106.                    break;
  107.                 case StreamTokenizer.TT_WORD:
  108.                    String name = st.sval;
  109.                    double x = 0, y = 0, z = 0;
  110.                    if (st.nextToken() == StreamTokenizer.TT_NUMBER) 
  111.                    {
  112.                       x = st.nval;
  113.                       if (st.nextToken() == StreamTokenizer.TT_NUMBER) 
  114.                       {
  115.                          y = st.nval;
  116.                          if (st.nextToken() == StreamTokenizer.TT_NUMBER)
  117.                             z = st.nval;
  118.                       }
  119.                    }
  120.                    addVert(name, (float) x, (float) y, (float) z);
  121.                    while( st.ttype != StreamTokenizer.TT_EOL &&
  122.                           st.ttype != StreamTokenizer.TT_EOF )
  123.                       st.nextToken();
  124.  
  125.              }   // end Switch
  126.  
  127.           }  // end while
  128.  
  129.           is.close();
  130.  
  131.        }  // end Try
  132.        catch( IOException e) {}
  133.  
  134.        if (st.ttype != StreamTokenizer.TT_EOF)
  135.           throw new Exception(st.toString());
  136.  
  137.     }  // end XYZChemModel()
  138.  
  139.     /** Add a vertex to this model */
  140.     int addVert(String name, float x, float y, float z) {
  141.     int i = nvert;
  142.     if (i >= maxvert)
  143.         if (vert == null) {
  144.         maxvert = 100;
  145.         vert = new float[maxvert * 3];
  146.         atoms = new Atom[maxvert];
  147.         } else {
  148.         maxvert *= 2;
  149.         float nv[] = new float[maxvert * 3];
  150.         System.arraycopy(vert, 0, nv, 0, vert.length);
  151.         vert = nv;
  152.         Atom na[] = new Atom[maxvert];
  153.         System.arraycopy(atoms, 0, na, 0, atoms.length);
  154.         atoms = na;
  155.         }
  156.     Atom a = (Atom) atomTable.get(name.toLowerCase());
  157.     if (a == null) a = defaultAtom;
  158.     atoms[i] = a;
  159.     i *= 3;
  160.     vert[i] = x;
  161.     vert[i + 1] = y;
  162.     vert[i + 2] = z;
  163.     return nvert++;
  164.     }
  165.  
  166.     /** Transform all the points in this model */
  167.     void transform() {
  168.     if (transformed || nvert <= 0)
  169.         return;
  170.     if (tvert == null || tvert.length < nvert * 3)
  171.         tvert = new int[nvert * 3];
  172.     mat.transform(vert, tvert, nvert);
  173.     transformed = true;
  174.     }
  175.  
  176.  
  177.     /** Paint this model to a graphics context.  It uses the matrix associated
  178.     with this model to map from model space to screen space.
  179.     The next version of the browser should have double buffering,
  180.     which will make this *much* nicer */
  181.     void paint(Graphics g) {
  182.     if (vert == null || nvert <= 0)
  183.         return;
  184.     transform();
  185.     int v[] = tvert;
  186.     int zs[] = ZsortMap;
  187.     if (zs == null) {
  188.         ZsortMap = zs = new int[nvert];
  189.         for (int i = nvert; --i >= 0;)
  190.         zs[i] = i * 3;
  191.     }
  192.  
  193.     /*
  194.      * I use a bubble sort since from one iteration to the next, the sort
  195.      * order is pretty stable, so I just use what I had last time as a
  196.      * "guess" of the sorted order.  With luck, this reduces O(N log N)
  197.      * to O(N)
  198.      */
  199.  
  200.     for (int i = nvert - 1; --i >= 0;) {
  201.         boolean flipped = false;
  202.         for (int j = 0; j <= i; j++) {
  203.         int a = zs[j];
  204.         int b = zs[j + 1];
  205.         if (v[a + 2] > v[b + 2]) {
  206.             zs[j + 1] = a;
  207.             zs[j] = b;
  208.             flipped = true;
  209.         }
  210.         }
  211.         if (!flipped)
  212.         break;
  213.     }
  214.  
  215.     int lg = 0;
  216.     int lim = nvert;
  217.     Atom ls[] = atoms;
  218.     if (lim <= 0 || nvert <= 0)
  219.         return;
  220.     for (int i = 0; i < lim; i++) {
  221.         int j = zs[i];
  222.         int grey = v[j + 2];
  223.         if (grey < 0)
  224.         grey = 0;
  225.         if (grey > 15)
  226.         grey = 15;
  227.         // g.drawString(names[i], v[j], v[j+1]);
  228.         atoms[j/3].paint(g, v[j], v[j + 1], grey);
  229.         // g.drawImage(iBall, v[j] - (iBall.width >> 1), v[j + 1] -
  230.         // (iBall.height >> 1));
  231.     }
  232.     }
  233.  
  234.     /** Find the bounding box of this model */
  235.     void findBB() {
  236.     if (nvert <= 0)
  237.         return;
  238.     float v[] = vert;
  239.     float xmin = v[0], xmax = xmin;
  240.     float ymin = v[1], ymax = ymin;
  241.     float zmin = v[2], zmax = zmin;
  242.     for (int i = nvert * 3; (i -= 3) > 0;) {
  243.         float x = v[i];
  244.         if (x < xmin)
  245.         xmin = x;
  246.         if (x > xmax)
  247.         xmax = x;
  248.         float y = v[i + 1];
  249.         if (y < ymin)
  250.         ymin = y;
  251.         if (y > ymax)
  252.         ymax = y;
  253.         float z = v[i + 2];
  254.         if (z < zmin)
  255.         zmin = z;
  256.         if (z > zmax)
  257.         zmax = z;
  258.     }
  259.     this.xmax = xmax;
  260.     this.xmin = xmin;
  261.     this.ymax = ymax;
  262.     this.ymin = ymin;
  263.     this.zmax = zmax;
  264.     this.zmin = zmin;
  265.     }
  266. }
  267.  
  268. /** An applet to put a Cehmical model into a page */
  269. public class XYZApp extends Applet implements Runnable {
  270.     XYZChemModel md;
  271.     boolean painted = true;
  272.     float xfac;
  273.     int prevx, prevy;
  274.     float xtheta, ytheta;
  275.     float scalefudge = 1;
  276.     Matrix3D amat = new Matrix3D(), tmat = new Matrix3D();
  277.     String mdname = null;
  278.     String message = null;
  279.     Image backBuffer;
  280.     Graphics backGC;
  281.     Dimension backSize;
  282.  
  283.  
  284.     private synchronized void newBackBuffer() {
  285.     backBuffer = createImage(size().width, size().height);
  286.     backGC = backBuffer.getGraphics();
  287.     backSize = size();
  288.     }
  289.  
  290.     public void init() {
  291.     mdname = getParameter("model");
  292.     try {
  293.         scalefudge = Float.valueOf(getParameter("scale")).floatValue();
  294.     } catch(Exception e) {
  295.     };
  296.     amat.yrot(20);
  297.     amat.xrot(20);
  298.     if (mdname == null)
  299.         mdname = "model.obj";
  300.     resize(size().width <= 20 ? 400 : size().width,
  301.            size().height <= 20 ? 400 : size().height);
  302.     newBackBuffer();
  303.     }
  304.     public void run() {
  305.     InputStream is = null;
  306.     try {
  307.         Thread.currentThread().setPriority(Thread.MIN_PRIORITY);
  308.         is = new URL(getDocumentBase(), mdname).openStream();
  309.         XYZChemModel m = new XYZChemModel (is);
  310.         Atom.setApplet(this);
  311.         md = m;
  312.         m.findBB();
  313.         float xw = m.xmax - m.xmin;
  314.         float yw = m.ymax - m.ymin;
  315.         float zw = m.zmax - m.zmin;
  316.         if (yw > xw)
  317.         xw = yw;
  318.         if (zw > xw)
  319.         xw = zw;
  320.         float f1 = size().width / xw;
  321.         float f2 = size().height / xw;
  322.         xfac = 0.7f * (f1 < f2 ? f1 : f2) * scalefudge;
  323.     } catch(Exception e) {
  324.         e.printStackTrace();
  325.         md = null;
  326.         message = e.toString();
  327.     }
  328.     try {
  329.         if (is != null)
  330.         is.close();
  331.     } catch(Exception e) {
  332.     }
  333.     repaint();
  334.     }
  335.     public void start() {
  336.     if (md == null && message == null)
  337.         new Thread(this).start();
  338.     }
  339.     public void stop() {
  340.     }
  341.     public boolean mouseDown(Event e, int x, int y) {
  342.     prevx = x;
  343.     prevy = y;
  344.     return true;
  345.     }
  346.     public boolean mouseDrag(Event e, int x, int y) {
  347.     tmat.unit();
  348.     float xtheta = (prevy - y) * (360.0f / size().width);
  349.     float ytheta = (x - prevx) * (360.0f / size().height);
  350.     tmat.xrot(xtheta);
  351.     tmat.yrot(ytheta);
  352.     amat.mult(tmat);
  353.     if (painted) {
  354.         painted = false;
  355.         repaint();
  356.     }
  357.     prevx = x;
  358.     prevy = y;
  359.     return true;
  360.     }
  361.     public void update(Graphics g) {
  362.     if (backBuffer == null)
  363.         g.clearRect(0, 0, size().width, size().height);
  364.     paint(g);
  365.     }
  366.  
  367.     public void paint(Graphics g) {
  368.     if (md != null) {
  369.         md.mat.unit();
  370.         md.mat.translate(-(md.xmin + md.xmax) / 2,
  371.                  -(md.ymin + md.ymax) / 2,
  372.                  -(md.zmin + md.zmax) / 2);
  373.         md.mat.mult(amat);
  374.         // md.mat.scale(xfac, -xfac, 8 * xfac / size().width);
  375.         md.mat.scale(xfac, -xfac, 16 * xfac / size().width);
  376.         md.mat.translate(size().width / 2, size().height / 2, 8);
  377.         md.transformed = false;
  378.         if (backBuffer != null) {
  379.         if (!backSize.equals(size()))
  380.             newBackBuffer();
  381.         backGC.setColor(getBackground());
  382.         backGC.fillRect(0,0,size().width,size().height);
  383.         md.paint(backGC);
  384.         g.drawImage(backBuffer, 0, 0, this);
  385.         }
  386.         else
  387.         md.paint(g);
  388.         setPainted();
  389.     } else if (message != null) {
  390.         g.drawString("Error in model:", 3, 20);
  391.         g.drawString(message, 10, 40);
  392.     }
  393.     }
  394.     private synchronized void setPainted() {
  395.     painted = true;
  396.     notifyAll();
  397.     }
  398.  
  399.     private synchronized void waitPainted() 
  400.     {
  401.        while (!painted)
  402.        {
  403.           try
  404.           {
  405.              wait();
  406.           }
  407.           catch (InterruptedException e) {}
  408.        }
  409.        painted = false;
  410.     }
  411. }   // end class XYZApp
  412.  
  413. class Atom {
  414.     private static Applet applet;
  415.     private static byte[] data;
  416.     private final static int R = 40;
  417.     private final static int hx = 15;
  418.     private final static int hy = 15;
  419.     private final static int bgGrey = 192;
  420.     private final static int nBalls = 16;
  421.     private static int maxr;
  422.  
  423.     private int Rl;
  424.     private int Gl;
  425.     private int Bl;
  426.     private Image balls[];
  427.  
  428.     static {
  429.     data = new byte[R * 2 * R * 2];
  430.     int mr = 0;
  431.     for (int Y = 2 * R; --Y >= 0;) {
  432.         int x0 = (int) (Math.sqrt(R * R - (Y - R) * (Y - R)) + 0.5);
  433.         int p = Y * (R * 2) + R - x0;
  434.         for (int X = -x0; X < x0; X++) {
  435.         int x = X + hx;
  436.         int y = Y - R + hy;
  437.         int r = (int) (Math.sqrt(x * x + y * y) + 0.5);
  438.         if (r > mr)
  439.             mr = r;
  440.         data[p++] = r <= 0 ? 1 : (byte) r;
  441.         }
  442.     }
  443.     maxr = mr;
  444.     }
  445.     static void setApplet(Applet app) {
  446.     applet = app;
  447.     }
  448.     Atom(int Rl, int Gl, int Bl) {
  449.     this.Rl = Rl;
  450.     this.Gl = Gl;
  451.     this.Bl = Bl;
  452.     }
  453.     private final int blend(int fg, int bg, float fgfactor) {
  454.     return (int) (bg + (fg - bg) * fgfactor);
  455.     }
  456.     private void Setup() {
  457.     balls = new Image[nBalls];
  458.     byte red[] = new byte[256];
  459.     red[0] = (byte) bgGrey;
  460.     byte green[] = new byte[256];
  461.     green[0] = (byte) bgGrey;
  462.     byte blue[] = new byte[256];
  463.     blue[0] = (byte) bgGrey;
  464.     for (int r = 0; r < nBalls; r++) {
  465.         float b = (float) (r+1) / nBalls;
  466.         for (int i = maxr; i >= 1; --i) {
  467.         float d = (float) i / maxr;
  468.         red[i] = (byte) blend(blend(Rl, 255, d), bgGrey, b);
  469.         green[i] = (byte) blend(blend(Gl, 255, d), bgGrey, b);
  470.         blue[i] = (byte) blend(blend(Bl, 255, d), bgGrey, b);
  471.         }
  472.         IndexColorModel model = new IndexColorModel(8, maxr + 1,
  473.                             red, green, blue, 0);
  474.         balls[r] = applet.createImage(
  475.         new MemoryImageSource(R*2, R*2, model, data, 0, R*2));
  476.     }
  477.     }
  478.     void paint(Graphics gc, int x, int y, int r) {
  479.     Image ba[] = balls;
  480.     if (ba == null) {
  481.         Setup();
  482.         ba = balls;
  483.     }
  484.     Image i = ba[r];
  485.     int size = 10 + r;
  486.     gc.drawImage(i, x - (size >> 1), y - (size >> 1), size, size, applet);
  487.     }
  488. }
  489.