Changeset 121 for trunk/matml/webselector/php/dbquery.php.inc

Timestamp:

10/26/2004 05:24:02 PM (8 years ago)

Author:

kstemen

Message:

Added categorization.

Files:

• trunk/matml/webselector/php/dbquery.php.inc (modified) (4 diffs)

trunk/matml/webselector/php/dbquery.php.inc

@@ -1 @@
-<?
-function db_get_plot_values($plottype, $log){ 
-
-include("db.php.inc");
+<?php
+class Unit
+{
+    var $name;
+    var $num_name;
+    //Multiple existing value by scale to get num_name
+    var $scale;
+    function Unit($name,$num_name,$scale)
+    {
+        $this->name=$name;
+        $this->num_name=$num_name;
+        $this->scale=$scale;
+    }
+    function name_split()
+    {
+        $ret=explode("/",$this->name);
+        if(!isset($ret[1]))
+            $ret[1]="";
+        return $ret;
+    }
+    function div($other)
+    {
+        $mine=$this->name_split();
+        $eirs=$other->name_split();
+        if($mine[0]&&$eirs[1])
+            $mine[0].="-$eirs[1]";
+        else
+            $mine[0].=$eirs[1];
+        if($mine[1]&&$eirs[0])
+            $mine[1].="-$eirs[0]";
+        else
+            $mine[1].=$eirs[0];
+        if(!$mine[1])
+            $name=$mine[0];
+        else
+            $name="$mine[0]/$mine[1]";
+        $ret=new Unit(
+            $name,
+            $this->num_name/$other->num_name,
+            $this->scale/$other->scale
+        );
+        return $ret;
+    }
+    function mul($other)
+    {
+        $mine=$this->name_split();
+        $eirs=$other->name_split();
+        if($mine[0]&&$eirs[0])
+            $mine[0].="-$eirs[0]";
+        else
+            $mine[0].=$eirs[0];
+        if($mine[1]&&$eirs[1])
+            $mine[1].="-$eirs[1]";
+        else
+            $mine[1].=$eirs[1];
+        if(!$mine[1])
+            $name=$mine[0];
+        else
+            $name="$mine[0]/$mine[1]";
+        $ret=new Unit(
+            "$name",
+            $this->num_name*$other->num_name,
+            $this->scale*$other->scale
+        );
+        return $ret;
+    }
+    function equals($other)
+    {
+        return abs($this->num_name-$other)<.00000001;
+    }
+}
+
+$unit_infos=array(
+    "g"=>new Unit("g",2,1),
+    "m"=>new Unit("m",3,1),
+    "K"=>new Unit("K",5,1),
+    "s"=>new Unit("s",7,1),
+    "°C"=>new Unit("°C",11,1),
+    "lb"=>new Unit("lb",13,1),
+    "°F"=>new Unit("°F",17,1),
+    //kg m s2
+    "N"=>new Unit("N",2*3/(7*7),1000),
+    //1 pound = 4.44822162 Newtons
+    "p"=>new Unit("p",2*3/(7*7),4.44822162*1000),
+    "in"=>new Unit("in",3,0.0254),
+    "cm"=>new Unit("cm",3,.01),
+    //kg*m*m/(s^3)
+    "W"=>new Unit("W",2*3*3/(7*7*7),1000),
+    //N m-2
+    "Pa"=>new Unit("Pa",2*3/(7*7)/(3*3),1000),
+    "MPa"=>new Unit("MPa",2*3/(7*7)/(3*3),1000000*1000),
+    "cc"=>new Unit("cc",3*3*3,.01*.01*.01),
+    //1 kg × 1 m2 × 1 s-2
+    "J"=>new Unit("J",2*3*3/(7*7),1000),
+    "BTU"=>new Unit("BTU",2*3*3/(7*7),1055*1000),
+    "hr"=>new Unit("hr",7,60*60),
+    "ft"=>new Unit("ft",3,12*0.0254)
+);
+
+class PropertyData
+{
+    //average value out of all the matches
+    var $avg_value;
+    var $num_unit;
+    //number of matches
+    var $count;
+    function PropertyData()
+    {
+        $this->avg_value=0;
+        $this->num_unit=0;
+        $this->count=0;
+    }
+}
+
+function parse_units($units_element)
+{
+    global $unit_infos;
+    global $dom;
+    $scale=1;
+    if(get_class($units_element)=="DOMNodeList")
+        $units_element=$units_element->item(0);
+    if($units_element)
+    {
+        $num_unit=1;
+        $factor=$units_element->getAttribute("factor");
+        if($factor)
+            $scale*=$factor;
+        //$units=$units_element->get_elements_by_tagname("Unit");
+        $units=$units_element->getElementsByTagName("Unit");
+        foreach($units as $unit)
+        {
+            $name=$unit->getElementsByTagName("Name");
+            $name=$name->item(0);
+            $name=$name->nodeValue;
+            $power=$unit->getAttribute("power");
+            if(!$power)
+                $power=1;
+            $unit_info=$unit_infos[$name];
+            if(!$unit_info)
+                echo "Error: unknown unit $name\n";
+            $scale*=pow($unit_info->scale,$power);
+            $num_unit*=pow($unit_info->num_name,$power);
+        }
+    }
+    else
+        $num_unit=0;
+    return array($scale,$num_unit);
+}
+
+function get_property_data($filename,$property_names)
+{
+    global $dom;
+   $dom=new DOMDocument;
+   $dom->loadXML(file_get_contents($filename));
+   $context=new DOMXPath($dom);
+  
+   $equivalencies=array(
+       array_flip(array("Young's Modulus","Modulus of Elasticity","Tensile Modulus")),
+       array_flip(array("Tensile Strength, Yield","Yield Strength","Tensile Strength @ Yield")),
+       array_flip(array("Heat Capacity","Specific Heat"))
+    );
+   $translation_string="\"abcdefghijklmnopqrstuvwxyz'\",'ABCDEFGHIJKLMNOPQRSTUVWXYZ'";
+   $ret=array();
+    $details=$context->query("/MatML_Doc/Material/BulkDetails/Name/text()");
+   $item=$details->item(0);
+   $Matid=$item->nodeValue;
+   foreach($property_names as $key=>$property_name)
+   {
+       $found=false;
+       foreach($equivalencies as $equivalency)
+       {
+            if(array_key_exists($property_name,$equivalency))
+            {
+                $name_condition="";
+                foreach($equivalency as $equiv=>$value)
+                {
+                    if($name_condition!="")
+                        $name_condition.=" or ";
+                    $name_condition.="translate(Name,$translation_string)='".strtoupper(strtr($equiv,array("'"=>"")))."'";
+                }
+                $found=true;
+                break;
+            }
+        }
+        if($property_name=="Material Name")
+        {
+           $ret[$key]=$Matid;
+           continue;
+        }
+        else if($property_name=="Material Categories")
+        {
+            $categories=array();
+            $queried=$context->query("/MatML_Doc/Material/BulkDetails/Class/Name/text()");
+            foreach($queried as $category)
+            {
+                $categories[]=$category->nodeValue;
+            }
+            $queried=$context->query("/MatML_Doc/Material/BulkDetails/Subclass/Name/text()");
+            foreach($queried as $category)
+            {
+                $categories[]=$category->nodeValue;
+            }
+            $ret[$key]=$categories;
+            continue;
+        }
+        else if(!$found)
+        {
+           $name_condition="translate(Name,$translation_string)='".strtoupper($property_name)."'";
+        }
+        $details=$context->query("/MatML_Doc/Metadata/PropertyDetails[$name_condition]");
+       $data=new PropertyData;
+       if(!$details->length)
+            continue;
+        foreach($details as $detail)
+        {
+            $units=$detail->getElementsByTagName("Units");
+            list($scale,$num_unit)=parse_units($units);
+            if(!$data->num_unit)
+                $data->num_unit=$num_unit;
+            else
+            {
+                force_unit_match($Matid,new Unit("unknown",$num_unit,$scale),$data);
+            }
+
+            $propid=$detail->getAttribute("id");
+            $property_datas=$context->query("/MatML_Doc/Material/BulkDetails/PropertyData[@property='$propid']/Data");
+            foreach($property_datas as $property_data)
+            {
+                $data->count++;
+                $data->avg_value+=$property_data->nodeValue*$scale;
+            }
+        }
+        if($data->count)
+            $data->avg_value/=$data->count;
+        $ret[$key]=$data;
+   }
+   return $ret;
+}
+
+function fractionalize($number)
+{
+    $integer=floor($number);
+    $decimal=$number-$integer;
+    if($decimal<.00000001)
+        return array($integer,1);
+    if($decimal>.99999999)
+        return array($integer+1,1);
+    $denominator=fractionalize(1/$decimal);
+    $denominator=$denominator[0];
+    return array(round($number*$denominator),$denominator);
+}
+
+function get_factor_power($number,$factor)
+{
+    $power=0;
+    list($numerator,$denominator)=fractionalize($number);
+    if($numerator==0)
+        return 0;
+    list($fn,$fd)=fractionalize($factor);
+    while(($numerator % $fn)==0 && ($denominator % $fd) == 0)
+    {
+        $numerator/=$fn;
+        $denominator/=$fd;
+        $power++;
+    }
+    if($power||$denominator==1)
+        return $power;
+    while(($numerator % $fd)==0 && ($denominator % $fn) == 0)
+    {
+        $numerator/=$fd;
+        $denominator/=$fn;
+        $power--;
+    }
+    return $power;
+}
+
+function force_unit_match($material,$unit,&$property_data)
+{
+    //Convert BTUs per degree F to Joules per degree C
+    $BTUs_per_degree_F=2*3*3/(7*7)/17;
+    $power=get_factor_power($property_data->num_unit,$BTUs_per_degree_F);
+    $property_data->num_unit*=pow($BTUs_per_degree_F,-$power)*pow(
+        2*3*3/(7*7)/11,$power);
+    $property_data->avg_value*=pow(
+         1899.1009999999999/1000,$power);
+
+    if(!$unit->equals($property_data->num_unit))
+    {
+        $power=get_factor_power($property_data->num_unit,13);
+        if($unit->equals($property_data->num_unit*pow(13,-$power)*pow(2*3/(7*7),$power)))
+        {
+            //it uses pound as a force
+            $property_data->num_unit*=pow(13,-$power)*pow(2*3/(7*7),$power);
+            $property_data->avg_value*=pow(4.44822162*1000,$power);
+            return true;
+        }
+        if($unit->equals($property_data->num_unit*pow(13,-$power)*pow(2,$power)))
+        {
+            //it uses pound as a mass
+            $property_data->num_unit*=pow(13,-$power)*pow(2,$power);
+            $property_data->avg_value*=pow(453.59237,$power);
+            return true;
+        }
+        echo "Error: unit mismatch ($unit->name $unit->num_name vs $property_data->num_unit) in $material\n";
+        return false;
+    }
+    return true;
+}
+
+function db_get_plot_values($plottype, $log, $category_options){ 
+    global $unit_infos;
 
   // So far only 4 plot types
@@ -14 +321 @@
   //$utsys = "0";
   
-  // I am not very proud of the following nasty hacks but we needed this working ASAP for the pilot
-  
-  
-  // Query values from DB
-  
-  // THESE QUERIES are just for testing purpouses
-  
-  // Yield strength
-  
-  $query["yield"] = "select matid,propname,value,units from properties where lower(propname) like lower('Yield strength%') or lower(propname) like lower('%Strength%Yield%') or lower(propname) like lower('Tensile%Strength%Yield%')";
-  
-  // Density
-  $query["density"] = "select matid,value,units from properties where lower(propname) like lower('%density%')";
-  
-  // Young's Modulus
-  $query["young"] = "select matid,value,units from properties where lower(propname) like lower('%Modulus of Elasticity%')";
-  
-  // Thermal conductivity
-  $query["thermcond"] = " select matid,value,units from properties where lower(propname) like lower('%thermal%conductivity')";
-  
-  // Heat capacity (only for 3)
-  $query["Cp"] = "select matid,value,units from properties where lower(propname) like lower('%heat%capacity%') or lower(propname) like lower('%specific%heat%')";
-  
-  
-  // select value,matname from materials as m,properties as p where m.matid=p.matid and propid='pr1'
-  
-  switch($plottype){
-    
-  case "EvsYield":
-    $Yquery = pg_query($conn,$query["young"]);
-    $Xquery = pg_query($conn,$query["yield"]);
+  
+    switch($plottype){
+    
+    case "EvsYield":
+        $query=array(
+            "Yquery"=>"Young's Modulus",
+            "Xquery"=>"Tensile Strength, Yield",
+            "Matid"=>"Material Name"
+        );
+        $final_units=array(
+            "Yquery"=>$unit_infos["MPa"],
+            "Xquery"=>$unit_infos["MPa"]
+        );
    break;
    
   case "EovDvsYieldovD":
-    $Yquery = pg_query($conn,$query["young"]);
-    $Xquery = pg_query($conn,$query["yield"]);
-    $extra  = pg_query($conn,$query["density"]);
+        $query=array(
+            "Yquery"=>"Young's Modulus",
+            "Xquery"=>"Tensile Strength, Yield",
+            "extra"=>"Density",
+            "Matid"=>"Material Name"
+        );
+        $final_units=array(
+            "Yquery"=>$unit_infos["MPa"],
+            "Xquery"=>$unit_infos["MPa"],
+            "extra"=>$unit_infos["g"]->div($unit_infos["cc"])
+        );
     break;
     
   case "EvsD":
-    $Yquery = pg_query($conn,$query["young"]);
-    $Xquery = pg_query($conn,$query["density"]);
+        $query=array(
+            "Yquery"=>"Young's Modulus",
+            "Xquery"=>"Density",
+            "Matid"=>"Material Name"
+        );
+        $final_units=array(
+            "Yquery"=>$unit_infos["MPa"],
+            "Xquery"=>$unit_infos["g"]->div($unit_infos["cc"])
+        );
     break;
     
   case "KvsDCp":
-    $Yquery = pg_query($conn,$query["thermcond"]);
-    $Xquery = pg_query($conn,$query["Cp"]);
-    $extra  = pg_query($conn,$query["density"]);
+        $query=array(
+            "Yquery"=>"Thermal Conductivity",
+            "Xquery"=>"Heat Capacity",
+            "extra"=>"Density",
+            "Matid"=>"Material Name"
+        );
+        $final_units=array(
+            "Yquery"=>$unit_infos["W"]->div($unit_infos["m"]->mul($unit_infos["K"])),
+            "Xquery"=>$unit_infos["J"]->div($unit_infos["g"]->mul($unit_infos["°C"])),
+            "extra"=>$unit_infos["g"]->div($unit_infos["cc"])
+        );
     break;
   }
-  
-  
-  // Fetch values from DB and sort them in an array (actually 3 arrays)
-  
-  if($Yquery && $Xquery){
-    
+  $query["categories"]="Material Categories";
+  
     $valuesY = array();
     $valuesX = array();
-    $keysY = array();
-    $keysX = array();
-    $unitsY = array();
-    $unitsX = array();
-    
-    
-    
-    while( $rowx = pg_fetch_array ($Xquery) ){
-      
-      //fetching matid, value, units fields
-      
-      array_push($keysX,$rowx["matid"]);
-      array_push($valuesX,$rowx["value"]);
-      array_push($unitsX,$rowx["units"]);
-      
-    }
-    
-    
-    while( $rowy = pg_fetch_array($Yquery) ){
-      
-      //fetching matid, value, units fields
-      
-      array_push($keysY,$rowy["matid"]);
-      array_push($valuesY,$rowy["value"]);
-      array_push($unitsY,$rowy["units"]);
-    }
-    
-
-    
-    // Non very pretty hack for getting the common keys due to the lack
-    // of array_intersect_assoc prior PHP <= 4.3
-    // Hope this solves the problem for now, and don't slow down things too much
-    
-    $dataY =  array_combine($keysY, $valuesY);
-    $dataX =  array_combine($keysX, $valuesX);
-
-    $uniY  =  array_combine($keysY, $unitsY);
-    $uniX  =  array_combine($keysX, $unitsX);
-    
-    
-    $valuesX = array();
-    $valuesY = array();
-    $unitsX = array();
-    $unitsY = array();
-    
-    
-    $common_keys = array_intersect($keysX, $keysY);
-    
-   //   for (i in array resultante) { $valores[$llaves[$i]] }
-    
-    foreach($common_keys as $key){
-      
-      array_push($valuesY,$dataY[$key]); 
-      array_push($valuesX,$dataX[$key]); 
-      
-      array_push($unitsY,$uniY[$key]); 
-      array_push($unitsX,$uniX[$key]); 
-      
-      
-    }
-    
-  }
-  
-  
-  if($extra){
-    
     $valuese = array();
-    $keyse = array();
-    $unitse = array();
-    
-    while($rowextra = pg_fetch_array ($extra) ){
-   //fetching matid, value, units fields
-      
-      array_push($keyse,$rowextra["matid"]);
-      array_push($valuese,$rowextra["value"]);
-      array_push($unitse,$rowextra["units"]);
-    }
-    
-    
-    $datae =  array_combine($keyse, $valuese);
-    $unie  =  array_combine($keyse, $unitse);
-    
-    $valuese = array();
-    $unitse = array();
-                
-    $common_keys = array_intersect($common_keys, $keyse);
-    
-
-    //   for (i in array resultante) { $valores[$llaves[$i]] }
-    
-    foreach($common_keys as $key){
-      
-      array_push($valuesY,$dataY[$key]); 
-      array_push($valuesX,$dataX[$key]); 
-      array_push($valuese,$datae[$key]); 
-      
-      array_push($unitsY,$uniY[$key]); 
-      array_push($unitsX,$uniX[$key]); 
-      array_push($unitse,$unie[$key]); 
-    }
-  }
+    $categories = array();
+  
+    $handle = opendir('materials');
+
+    while (false !== ($file = readdir($handle))) {
+        if(preg_match("/\.xml?$/",$file)==1)
+        {
+            // Fetch values from DB and sort them in an array (actually 3 arrays)
+            $values=get_property_data("materials/$file",$query);
+
+            if(isset($values["Xquery"])&&
+                isset($values["Yquery"])
+            )
+            {
+                $values_categories=$values["categories"];
+                $values_categories[]="";
+                foreach($values_categories as $category)
+                {
+                    if(isset($category_options[$category]))
+                        break;
+                }
+                if(!isset($category_options[$category]))
+                    continue;
+                $category=$category_options[$category];
+                force_unit_match($values["Matid"],$final_units["Xquery"],$values["Xquery"]);
+                force_unit_match($values["Matid"],$final_units["Yquery"],$values["Yquery"]);
+                if(isset($query["extra"]))
+                {
+                    if(isset($values["extra"]))
+                    {
+                        force_unit_match($values["Matid"],$final_units["extra"],$values["extra"]);
+                        $valuesX[]=($values["Xquery"]->avg_value/$final_units["Xquery"]->scale);
+                        $valuesY[]=($values["Yquery"]->avg_value/$final_units["Yquery"]->scale);
+                        $valuese[]=$values["extra"]->avg_value/$final_units["extra"]->scale;
+                        $names[]=$values["Matid"];
+                        $categories[]=$category;
+                    }
+                }
+                else
+                {
+                        $valuesX[]=($values["Xquery"]->avg_value/$final_units["Xquery"]->scale);
+                        $valuesY[]=($values["Yquery"]->avg_value/$final_units["Yquery"]->scale);
+                        $names[]=$values["Matid"];
+                        $categories[]=$category;
+                }
+            }
+        }
+    }
+    closedir($handle);
   
   
   // Assign the actual useful data
-  // NOTE: I AM not checking the units yet
   
   // Labels and Units
@@ -185 +441 @@
     
   case "EvsYield":
-    $ylabel = "Young's Modulus ($unitsY[0])";
-    $xlabel = "Yield Strength ($unitsX[0])";
+    $ylabel = "Young's Modulus (".$final_units["Yquery"]->name.")";
+    $xlabel = "Yield Strength (".$final_units["Xquery"]->name.")";
     break;
     
   case "EvsD":
-    $ylabel = "Young's Modulus ($unitsY[0])";
-    $xlabel = "Density ($unitsX[0])";
+    $ylabel = "Young's Modulus (".$final_units["Yquery"]->name.")";
+    $xlabel = "Density (".$final_units["Xquery"]->name.")";
     break;
     
   case "EovDvsYieldovD":
-    $ylabel = "Young's Modulus / Density ($unitsY[0] / $unitse[0])";
-    $xlabel = "Yield Strength / Density ($unitsX[0] / $unitse[0])";
+    $ylabel = "Young's Modulus / Density (".$final_units["Yquery"]->name." / ".$final_units["extra"]->name.")";
+    $xlabel = "Yield Strength / Density (".$final_units["Xquery"]->name." / ".$final_units["extra"]->name.")";
     $valuesY = array_div($valuesY,$valuese);
     $valuesX = array_div($valuesX,$valuese);
@@ -202 +458 @@
     
   case "KvsDCp":
-    $ylabel = "Thermal Conductivity $unitsY[0]";
-    $xlabel = "Heat Capacity * Density ($unitsX[0] * $unitse[0])";
+    $ylabel = "Thermal Conductivity (".$final_units["Yquery"]->name.")";
+    $xlabel = "Heat Capacity * Density (".$final_units["Xquery"]->name." * ".$final_units["extra"]->name.")";
     $valuesX = array_mult($valuesX,$valuese);
     break;
   }
 
-
-  //Fetch the materials names
-  $in = "('" . implode("', '", $common_keys) . "')";
-  $name_query  = pg_query($conn, "select matname from materials where matid in $in");
-  
-  //echo "select matname from materials where matid in $in";
-  
-  if($name_query){
-    $names = array();
-    while( $namerow = pg_fetch_row ($name_query)){
-      array_push($names,$namerow[0]);
-    }
-  }
-
-  return array(&$valuesX, &$valuesY, &$names, $xlabel, $ylabel);
-
-}
-
+  return array(&$valuesX, &$valuesY, &$names, $xlabel, $ylabel, $categories);
+
+}
+
+//var_dump(db_get_plot_values("EvsD","")); 
 ?>