Subversion Repositories wimsdev

Rev

Rev 4676 | Blame | Compare with Previous | Last modification | View Log | RSS feed

  1. /*
  2.  * @(#)XYZApp.java
  3.  *
  4.  * Copyright (c) 1994-1996 Sun Microsystems, Inc. All Rights Reserved.
  5.  *
  6.  * Permission to use, copy, modify, and distribute this software
  7.  * and its documentation for NON-COMMERCIAL or COMMERCIAL purposes and
  8.  * without fee is hereby granted.
  9.  * Please refer to the file http://java.sun.com/copy_trademarks.html
  10.  * for further important copyright and trademark information and to
  11.  * http://java.sun.com/licensing.html for further important licensing
  12.  * information for the Java (tm) Technology.
  13.  *
  14.  * SUN MAKES NO REPRESENTATIONS OR WARRANTIES ABOUT THE SUITABILITY OF
  15.  * THE SOFTWARE, EITHER EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED
  16.  * TO THE IMPLIED WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
  17.  * PARTICULAR PURPOSE, OR NON-INFRINGEMENT. SUN SHALL NOT BE LIABLE FOR
  18.  * ANY DAMAGES SUFFERED BY LICENSEE AS A RESULT OF USING, MODIFYING OR
  19.  * DISTRIBUTING THIS SOFTWARE OR ITS DERIVATIVES.
  20.  *
  21.  * THIS SOFTWARE IS NOT DESIGNED OR INTENDED FOR USE OR RESALE AS ON-LINE
  22.  * CONTROL EQUIPMENT IN HAZARDOUS ENVIRONMENTS REQUIRING FAIL-SAFE
  23.  * PERFORMANCE, SUCH AS IN THE OPERATION OF NUCLEAR FACILITIES, AIRCRAFT
  24.  * NAVIGATION OR COMMUNICATION SYSTEMS, AIR TRAFFIC CONTROL, DIRECT LIFE
  25.  * SUPPORT MACHINES, OR WEAPONS SYSTEMS, IN WHICH THE FAILURE OF THE
  26.  * SOFTWARE COULD LEAD DIRECTLY TO DEATH, PERSONAL INJURY, OR SEVERE
  27.  * PHYSICAL OR ENVIRONMENTAL DAMAGE ("HIGH RISK ACTIVITIES").  SUN
  28.  * SPECIFICALLY DISCLAIMS ANY EXPRESS OR IMPLIED WARRANTY OF FITNESS FOR
  29.  * HIGH RISK ACTIVITIES.
  30.  */
  31.  
  32. import java.applet.Applet;
  33. import java.awt.Image;
  34. import java.awt.Event;
  35. import java.awt.Graphics;
  36. import java.io.StreamTokenizer;
  37. import java.io.InputStream;
  38. import java.io.BufferedInputStream;
  39. import java.net.URL;
  40. import java.util.Hashtable;
  41. import java.awt.image.IndexColorModel;
  42. import java.awt.image.ColorModel;
  43. import java.awt.image.MemoryImageSource;
  44.  
  45.  
  46. /** The representation of a Chemical .xyz model */
  47. class XYZChemModel {
  48.     float vert[];
  49.     Atom atoms[];
  50.     int tvert[];
  51.     int ZsortMap[];
  52.     int nvert, maxvert;
  53.  
  54.     static Hashtable atomTable = new Hashtable();
  55.     static Atom defaultAtom;
  56.     static {
  57.         atomTable.put("c", new Atom(0, 0, 0));
  58.         atomTable.put("x", new Atom(176, 160, 240));
  59.         atomTable.put("h", new Atom(210, 210, 210));
  60.         atomTable.put("n", new Atom(0, 0, 255));
  61.         atomTable.put("o", new Atom(255, 0, 0));
  62.         atomTable.put("p", new Atom(255, 0, 255));
  63.         atomTable.put("s", new Atom(255, 255, 0));
  64.         atomTable.put("hn", new Atom(150, 255, 150)); /* !!*/
  65.         defaultAtom = new Atom(255, 100, 200);
  66.     }
  67.  
  68.     boolean transformed;
  69.     Matrix3D mat;
  70.  
  71.     float xmin, xmax, ymin, ymax, zmin, zmax;
  72.  
  73.  
  74.     XYZChemModel () {
  75.         mat = new Matrix3D();
  76.         mat.xrot(20);
  77.         mat.yrot(30);
  78.     }
  79.  
  80.  
  81.     /** Create a Cehmical model by parsing an input stream */
  82.     XYZChemModel (InputStream is) throws Exception {
  83.         this();
  84.         StreamTokenizer st = new StreamTokenizer(new BufferedInputStream(is, 4000));
  85.         st.eolIsSignificant(true);
  86.         st.commentChar('#');
  87.         int slot = 0;
  88. scan:
  89.         while (true)
  90.             switch (st.nextToken()) {
  91.               case StreamTokenizer.TT_EOF:
  92.                 break scan;
  93.               default:
  94.                 break;
  95.               case StreamTokenizer.TT_WORD:
  96.                 String name = st.sval;
  97.                 double x = 0, y = 0, z = 0;
  98.                 if (st.nextToken() == StreamTokenizer.TT_NUMBER) {
  99.                     x = st.nval;
  100.                     if (st.nextToken() == StreamTokenizer.TT_NUMBER) {
  101.                         y = st.nval;
  102.                         if (st.nextToken() == StreamTokenizer.TT_NUMBER)
  103.                             z = st.nval;
  104.                     }
  105.                 }
  106.                 addVert(name, (float) x, (float) y, (float) z);
  107.                 while (st.ttype != StreamTokenizer.TT_EOL &&
  108.                         st.ttype != StreamTokenizer.TT_EOF)
  109.                     st.nextToken();
  110.             }
  111.         is.close();
  112.         if (st.ttype != StreamTokenizer.TT_EOF)
  113.             throw new Exception(st.toString());
  114.     }
  115.  
  116.     /** Add a vertex to this model */
  117.     int addVert(String name, float x, float y, float z) {
  118.         int i = nvert;
  119.         if (i >= maxvert)
  120.             if (vert == null) {
  121.                 maxvert = 100;
  122.                 vert = new float[maxvert * 3];
  123.                 atoms = new Atom[maxvert];
  124.             } else {
  125.                 maxvert *= 2;
  126.                 float nv[] = new float[maxvert * 3];
  127.                 System.arraycopy(vert, 0, nv, 0, vert.length);
  128.                 vert = nv;
  129.                 Atom na[] = new Atom[maxvert];
  130.                 System.arraycopy(atoms, 0, na, 0, atoms.length);
  131.                 atoms = na;
  132.             }
  133.         Atom a = (Atom) atomTable.get(name.toLowerCase());
  134.         if (a == null) a = defaultAtom;
  135.         atoms[i] = a;
  136.         i *= 3;
  137.         vert[i] = x;
  138.         vert[i + 1] = y;
  139.         vert[i + 2] = z;
  140.         return nvert++;
  141.     }
  142.  
  143.     /** Transform all the points in this model */
  144.     void transform() {
  145.         if (transformed || nvert <= 0)
  146.             return;
  147.         if (tvert == null || tvert.length < nvert * 3)
  148.             tvert = new int[nvert * 3];
  149.         mat.transform(vert, tvert, nvert);
  150.         transformed = true;
  151.     }
  152.  
  153.  
  154.     /** Paint this model to a graphics context.  It uses the matrix associated
  155.         with this model to map from model space to screen space.
  156.         The next version of the browser should have double buffering,
  157.         which will make this *much* nicer */
  158.     void paint(Graphics g) {
  159.         if (vert == null || nvert <= 0)
  160.             return;
  161.         transform();
  162.         int v[] = tvert;
  163.         int zs[] = ZsortMap;
  164.         if (zs == null) {
  165.             ZsortMap = zs = new int[nvert];
  166.             for (int i = nvert; --i >= 0;)
  167.                 zs[i] = i * 3;
  168.         }
  169.  
  170.         /*
  171.          * I use a bubble sort since from one iteration to the next, the sort
  172.          * order is pretty stable, so I just use what I had last time as a
  173.          * "guess" of the sorted order.  With luck, this reduces O(N log N)
  174.          * to O(N)
  175.          */
  176.  
  177.         for (int i = nvert - 1; --i >= 0;) {
  178.             boolean flipped = false;
  179.             for (int j = 0; j <= i; j++) {
  180.                 int a = zs[j];
  181.                 int b = zs[j + 1];
  182.                 if (v[a + 2] > v[b + 2]) {
  183.                     zs[j + 1] = a;
  184.                     zs[j] = b;
  185.                     flipped = true;
  186.                 }
  187.             }
  188.             if (!flipped)
  189.                 break;
  190.         }
  191.  
  192.         int lg = 0;
  193.         int lim = nvert;
  194.         Atom ls[] = atoms;
  195.         if (lim <= 0 || nvert <= 0)
  196.             return;
  197.         for (int i = 0; i < lim; i++) {
  198.             int j = zs[i];
  199.             int grey = v[j + 2];
  200.             if (grey < 0)
  201.                 grey = 0;
  202.             if (grey > 15)
  203.                 grey = 15;
  204.             // g.drawString(names[i], v[j], v[j+1]);
  205.             atoms[j/3].paint(g, v[j], v[j + 1], grey);
  206.             // g.drawImage(iBall, v[j] - (iBall.width >> 1), v[j + 1] -
  207.             // (iBall.height >> 1));
  208.         }
  209.     }
  210.  
  211.     /** Find the bounding box of this model */
  212.     void findBB() {
  213.         if (nvert <= 0)
  214.             return;
  215.         float v[] = vert;
  216.         float xmin = v[0], xmax = xmin;
  217.         float ymin = v[1], ymax = ymin;
  218.         float zmin = v[2], zmax = zmin;
  219.         for (int i = nvert * 3; (i -= 3) > 0;) {
  220.             float x = v[i];
  221.             if (x < xmin)
  222.                 xmin = x;
  223.             if (x > xmax)
  224.                 xmax = x;
  225.             float y = v[i + 1];
  226.             if (y < ymin)
  227.                 ymin = y;
  228.             if (y > ymax)
  229.                 ymax = y;
  230.             float z = v[i + 2];
  231.             if (z < zmin)
  232.                 zmin = z;
  233.             if (z > zmax)
  234.                 zmax = z;
  235.         }
  236.         this.xmax = xmax;
  237.         this.xmin = xmin;
  238.         this.ymax = ymax;
  239.         this.ymin = ymin;
  240.         this.zmax = zmax;
  241.         this.zmin = zmin;
  242.     }
  243. }
  244.  
  245. /** An applet to put a Chemical model into a page */
  246. public class XYZApp extends Applet implements Runnable {
  247.     XYZChemModel md;
  248.     boolean painted = true;
  249.     float xfac;
  250.     int prevx, prevy;
  251.     float xtheta, ytheta;
  252.     float scalefudge = 1;
  253.     Matrix3D amat = new Matrix3D(), tmat = new Matrix3D();
  254.     String mdname = null;
  255.     String message = null;
  256.     Image backBuffer;
  257.     Graphics backGC;
  258.  
  259.  
  260.     public void init() {
  261.         mdname = getParameter("model");
  262.         try {
  263.             scalefudge = Float.valueOf(getParameter("scale")).floatValue();
  264.         } catch(Exception e) {
  265.         };
  266.         amat.yrot(20);
  267.         amat.xrot(20);
  268.         if (mdname == null)
  269.             mdname = "model.obj";
  270.         resize(size().width <= 20 ? 400 : size().width,
  271.                size().height <= 20 ? 400 : size().height);
  272.     }
  273.     public void run() {
  274.         InputStream is = null;
  275.         try {
  276.             Thread.currentThread().setPriority(Thread.MIN_PRIORITY);
  277.             is = new URL(getDocumentBase(), mdname).openStream();
  278.             XYZChemModel m = new XYZChemModel (is);
  279.             Atom.setApplet(this);
  280.             md = m;
  281.             m.findBB();
  282.             float xw = m.xmax - m.xmin;
  283.             float yw = m.ymax - m.ymin;
  284.             float zw = m.zmax - m.zmin;
  285.             if (yw > xw)
  286.                 xw = yw;
  287.             if (zw > xw)
  288.                 xw = zw;
  289.             float f1 = size().width / xw;
  290.             float f2 = size().height / xw;
  291.             xfac = 0.7f * (f1 < f2 ? f1 : f2) * scalefudge;
  292.             backBuffer = createImage(size().width, size().height);
  293.             backGC = backBuffer.getGraphics();
  294.         } catch(Exception e) {
  295.             e.printStackTrace();
  296.             md = null;
  297.             message = e.toString();
  298.         }
  299.         try {
  300.             if (is != null)
  301.                 is.close();
  302.         } catch(Exception e) {
  303.         }
  304.         repaint();
  305.     }
  306.     public void start() {
  307.         if (md == null && message == null)
  308.             new Thread(this).start();
  309.     }
  310.     public void stop() {
  311.     }
  312.     public boolean mouseDown(Event e, int x, int y) {
  313.         prevx = x;
  314.         prevy = y;
  315.         return true;
  316.     }
  317.     public boolean mouseDrag(Event e, int x, int y) {
  318.         tmat.unit();
  319.         float xtheta = (prevy - y) * (360.0f / size().width);
  320.         float ytheta = (x - prevx) * (360.0f / size().height);
  321.         tmat.xrot(xtheta);
  322.         tmat.yrot(ytheta);
  323.         amat.mult(tmat);
  324.         if (painted) {
  325.             painted = false;
  326.             repaint();
  327.         }
  328.         prevx = x;
  329.         prevy = y;
  330.         return true;
  331.     }
  332.     public void update(Graphics g) {
  333.         if (backBuffer == null)
  334.             g.clearRect(0, 0, size().width, size().height);
  335.         paint(g);
  336.     }
  337.     public void paint(Graphics g) {
  338.         if (md != null) {
  339.             md.mat.unit();
  340.             md.mat.translate(-(md.xmin + md.xmax) / 2,
  341.                              -(md.ymin + md.ymax) / 2,
  342.                              -(md.zmin + md.zmax) / 2);
  343.             md.mat.mult(amat);
  344.             // md.mat.scale(xfac, -xfac, 8 * xfac / size().width);
  345.             md.mat.scale(xfac, -xfac, 16 * xfac / size().width);
  346.             md.mat.translate(size().width / 2, size().height / 2, 8);
  347.             md.transformed = false;
  348.             if (backBuffer != null) {
  349.                 backGC.setColor(getBackground());
  350.                 backGC.fillRect(0,0,size().width,size().height);
  351.                 md.paint(backGC);
  352.                 g.drawImage(backBuffer, 0, 0, this);
  353.             } else
  354.                 md.paint(g);
  355.             setPainted();
  356.         } else if (message != null) {
  357.             g.drawString("Error in model:", 3, 20);
  358.             g.drawString(message, 10, 40);
  359.         }
  360.     }
  361.     private synchronized void setPainted() {
  362.         painted = true;
  363.         notifyAll();
  364.     }
  365.     private synchronized void waitPainted() {
  366.         while (!painted)
  367.             try { wait(); } catch (Exception e) {}
  368.         painted = false;
  369.     }
  370. }
  371.  
  372. class Atom {
  373.     private static Applet applet;
  374.     private static byte[] data;
  375.     private final static int R = 40;
  376.     private final static int hx = 15;
  377.     private final static int hy = 15;
  378.     private final static int bgGrey = 192;
  379.     private final static int nBalls = 16;
  380.     private static int maxr;
  381.  
  382.     private int Rl;
  383.     private int Gl;
  384.     private int Bl;
  385.     private Image balls[];
  386.  
  387.     static {
  388.         data = new byte[R * 2 * R * 2];
  389.         int mr = 0;
  390.         for (int Y = 2 * R; --Y >= 0;) {
  391.             int x0 = (int) (Math.sqrt(R * R - (Y - R) * (Y - R)) + 0.5);
  392.             int p = Y * (R * 2) + R - x0;
  393.             for (int X = -x0; X < x0; X++) {
  394.                 int x = X + hx;
  395.                 int y = Y - R + hy;
  396.                 int r = (int) (Math.sqrt(x * x + y * y) + 0.5);
  397.                 if (r > mr)
  398.                     mr = r;
  399.                 data[p++] = r <= 0 ? 1 : (byte) r;
  400.             }
  401.         }
  402.         maxr = mr;
  403.     }
  404.     static void setApplet(Applet app) {
  405.         applet = app;
  406.     }
  407.     Atom(int Rl, int Gl, int Bl) {
  408.         this.Rl = Rl;
  409.         this.Gl = Gl;
  410.         this.Bl = Bl;
  411.     }
  412.     private final int blend(int fg, int bg, float fgfactor) {
  413.         return (int) (bg + (fg - bg) * fgfactor);
  414.     }
  415.     private void Setup() {
  416.         balls = new Image[nBalls];
  417.         byte red[] = new byte[256];
  418.         red[0] = (byte) bgGrey;
  419.         byte green[] = new byte[256];
  420.         green[0] = (byte) bgGrey;
  421.         byte blue[] = new byte[256];
  422.         blue[0] = (byte) bgGrey;
  423.         for (int r = 0; r < nBalls; r++) {
  424.             float b = (float) (r+1) / nBalls;
  425.             for (int i = maxr; i >= 1; --i) {
  426.                 float d = (float) i / maxr;
  427.                 red[i] = (byte) blend(blend(Rl, 255, d), bgGrey, b);
  428.                 green[i] = (byte) blend(blend(Gl, 255, d), bgGrey, b);
  429.                 blue[i] = (byte) blend(blend(Bl, 255, d), bgGrey, b);
  430.             }
  431.             IndexColorModel model = new IndexColorModel(8, maxr + 1,
  432.                                                         red, green, blue, 0);
  433.             balls[r] = applet.createImage(
  434.                 new MemoryImageSource(R*2, R*2, model, data, 0, R*2));
  435.         }
  436.     }
  437.     void paint(Graphics gc, int x, int y, int r) {
  438.         Image ba[] = balls;
  439.         if (ba == null) {
  440.             Setup();
  441.             ba = balls;
  442.         }
  443.         Image i = ba[r];
  444.         int size = 10 + r;
  445.         gc.drawImage(i, x - (size >> 1), y - (size >> 1), size, size, applet);
  446.     }
  447. }
  448.