root/trunk/matml/webselector/php/dbquery.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]='1';
        return $ret;
    }
    function name_mul($a,$b)
    {
        if($a=='1')
            return $b;
        else if($b=='1')
            return $a;
        return "$a-$b";
    }
    function div($other)
    {
        $mine=$this->name_split();
        $eirs=$other->name_split();
        $mine[0]=$this->name_mul($mine[0],$eirs[1]);
        $mine[1]=$this->name_mul($mine[1],$eirs[0]);

        if($mine[1]=='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();
        $mine[0]=$this->name_mul($mine[0],$eirs[0]);
        $mine[1]=$this->name_mul($mine[1],$eirs[1]);

        if($mine[1]=='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 pow($power)
    {

      // Write the actual POW function

        $mine=$this->name_split();


      if($power < 0){
            $temp   = $mine[0];
            $mine[0]= $mine[1];
            $mine[1]= $temp;
            $namepower=-$power;
      }
      else
            $namepower=$power;

      if( $namepower > 1){
        if($mine[0]!='1')
          {
            $mine[0].="^$namepower";
          }
        if($mine[1]!='1')
          {
            $mine[1].="^$namepower";
          }
      }
      elseif($namepower == 1){
        //Do nothing
      }
      elseif( $namepower == 0){
            $mine[0]="1";
            $mine[1]="1";
      }

        if($mine[1]=='1')
            $name=$mine[0];
        else
            $name="$mine[0]/$mine[1]";
        $ret=new Unit(
            "$name",
            pow($this->num_name,$power),
            pow($this->scale,$power)
        );
        return $ret;
    }
    function equals($other)
    {
        if(is_object($other))
            $num_name=$other->num_name;
        else
            $num_name=$other;
        return abs($this->num_name-$num_name)<.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 $unit;
    //number of matches
    var $count;
    function PropertyData()
    {
        $this->avg_value=0;
        $this->unit=0;
        $this->count=0;
    }
}

function parse_units($units_element)
{
    global $unit_infos;
    global $dom;
    if(get_class($units_element)=="DOMNodeList")
        $units_element=$units_element->item(0);
    if($units_element)
    {
        $return_unit=new Unit(
            "",
            1,
            1
        );
        $factor=$units_element->getAttribute("factor");
        if($factor)
            $return_unit->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";
            $return_unit=$return_unit->mul($unit_info->pow($power));
        }
    }
    else
        $return_unit=null;
    return $return_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","Youngs 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");
            $unit=parse_units($units);
            if(!$data->unit)
                $data->unit=$unit;
            else
            {
                force_unit_match($Matid,$unit,$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*$unit->scale;
            }
        }
        if($data->count)
            $data->avg_value/=$data->count;
        $ret[$key]=$data;
   }

//    echo "<br><br>";
//    print_r($ret);

   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)
{
    global $unit_infos;
    //Convert BTUs per degree F to Joules per degree C
    $BTUs_per_degree_F=$unit_infos['BTU']->div($unit_infos['°F']);
    $Joules_per_degree_C=$unit_infos['J']->div($unit_infos['°C']);
    $power=get_factor_power($property_data->unit->num_name,$BTUs_per_degree_F->num_name);
    $BTUs_to_Joules=$Joules_per_degree_C->div($BTUs_per_degree_F);
    $BTUs_to_Joules->scale=1899.1009999999999/1000;
    $property_data->unit=$property_data->unit->mul($BTUs_to_Joules->pow($power));
    $property_data->avg_value*=pow(
        $BTUs_to_Joules->scale,$power);

    //Convert all negative powers of degree F and degree C to Kelvin
    $degreeF_to_Kelvin=$unit_infos['K']->div($unit_infos['°F']);
    $degreeF_to_Kelvin->scale=5.0/9.0;
    $degreeC_to_Kelvin=$unit_infos['K']->div($unit_infos['°C']);
    $degreeC_to_Kelvin->scale=1;

    $power=get_factor_power($property_data->unit->num_name,$unit_infos['°F']->num_name);
    if($power>0)
        $power=0;
    $property_data->unit=$property_data->unit->mul($degreeF_to_Kelvin->pow($power));
    $property_data->avg_value*=pow(
        $degreeF_to_Kelvin->scale,$power);
        
    $power=get_factor_power($property_data->unit->num_name,$unit_infos['°C']->num_name);
    if($power>0)
        $power=0;
    $property_data->unit=$property_data->unit->mul($degreeC_to_Kelvin->pow($power));
    $property_data->avg_value*=pow(
        $degreeC_to_Kelvin->scale,$power);
    
    //Now convert negative powers of Kelvin to degree C and degree F as needed
    $power=get_factor_power($unit->num_name,$unit_infos['°F']->num_name);
    if($power>0)
        $power=0;
    $property_data->unit=$property_data->unit->div($degreeF_to_Kelvin->pow($power));
    $property_data->avg_value/=pow(
        $degreeF_to_Kelvin->scale,$power);

    $power=get_factor_power($unit->num_name,$unit_infos['°C']->num_name);
    if($power>0)
        $power=0;
    $property_data->unit=$property_data->unit->div($degreeC_to_Kelvin->pow($power));
    $property_data->avg_value/=pow(
        $degreeC_to_Kelvin->scale,$power);

    if(!$unit->equals($property_data->unit))
    {
        $power=get_factor_power($property_data->unit->num_name,13);
        if($unit->equals($property_data->unit->num_name*pow(13,-$power)*pow(2*3/(7*7),$power)))
        {
            //it uses pound as a force
            $newtons_per_pound=$unit_infos['N']->div($unit_infos['lb']);
            $newtons_per_pound->scale=4.44822162*1000;
            $newtons_per_pound=$newtons_per_pound->pow($power);

            $property_data->unit=$property_data->unit->mul($newtons_per_pound);
            $property_data->avg_value*=$newtons_per_pound->scale;
            return true;
        }
        if($unit->equals($property_data->unit->num_name*pow(13,-$power)*pow(2,$power)))
        {
            //it uses pound as a mass
            $grams_per_pound=$unit_infos['g']->div($unit_infos['lb']);
            $grams_per_pound->scale=453.59237;
            $grams_per_pound=$grams_per_pound->pow($power);
            $property_data->unit=$property_data->unit->mul($grams_per_pound);
            $property_data->avg_value*=$grams_per_pound->scale;
            return true;
        }
        echo "Error: unit mismatch ($unit->name vs ".$property_data->unit->name.") and ($unit->num_name vs ".$property_data->unit->num_name.") 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
  // EvsYield
  // EovDvsYieldovD
  // EvsD
  // KvsDCp
  
  // Unit selection
  
  //$utsys = "0";
  
  
    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":
        $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":
        $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":
        $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;
  }
  $query["categories"]="Material Categories";
  
    $valuesY = array();
    $valuesX = array();
    $valuese = array();
    $categories = array();
    $filenames = array();
  
    $handle = opendir('materials');


//     echo "<br><br>Begin Query <br>";
//     print_r($query);
//     echo "<br>End Query<br>";
    

    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);

//          echo "<br><br>Begin Values <br>";
//          print_r($values);
//          echo "End Values <br>";


            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;
                        
                }
                $filenames[]="materials/$file";
            }
        }
    }
    closedir($handle);
  
  
  // Assign the actual useful data
  
  // Labels and Units
  //$xlabel = "Density  (lb/cu in.)";
  //$ylabel = "Electrical Conductivity (ohm-cir mil/ft)";
  
  switch($plottype){
    
  case "EvsYield":
    $ylabel = "Young's Modulus (".$final_units["Yquery"]->name.")";
    $xlabel = "Yield Strength (".$final_units["Xquery"]->name.")";
    break;
    
  case "EvsD":
    $ylabel = "Young's Modulus (".$final_units["Yquery"]->name.")";
    $xlabel = "Density (".$final_units["Xquery"]->name.")";
    break;
    
  case "EovDvsYieldovD":
    $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);
    break;
    
  case "KvsDCp":
    $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;
  }

  $ret=array(&$valuesX, &$valuesY, &$names, $xlabel, $ylabel, $categories);
  $ret["filenames"]=$filenames;
  return $ret;
}

function db_get_plot_single_values($property_listx, $property_listy, $log, $category_options){ 

    global $unit_infos;


    // Prepare Array list


     $query = array("Matid"=>"Material Name");
     $power = array();
     $NXqueries = array();
     $NYqueries = array();

     $prop_units = array(
                         "Density" => $unit_infos["g"]->div($unit_infos["cc"]),
                         "Young's Modulus"       => $unit_infos["MPa"],
                         "Youngs Modulus"       => $unit_infos["MPa"],
                         "Yield Strength"      => $unit_infos["MPa"],
                         "Heat Capacity"      => $unit_infos["J"]->div($unit_infos["g"]->mul($unit_infos["K"])),
                         "Thermal Conductivity"      => $unit_infos["W"]->div($unit_infos["m"]->mul($unit_infos["K"])),
                         );

    $index = 0;
    $xlabel_num ="";
    $ylabel_num ="";
    $xlabel_den ="";
    $ylabel_den ="";

    foreach( $property_listx as $property=>$pow){
      $index++;
      $query["Xquery$index"]=$property;
      $power["Xquery$index"]= $pow;
      $final_units["Xquery$index"] = $prop_units["$property"];

      if( $pow > 1){
        $xlabel_num = "$property^$pow *";
      }
      elseif($pow == 1){
        $xlabel_num = "$property*";
      }
      elseif($pow == -1){
        $xlabel_den = "$property*";
      }
      elseif($pow < -1){      
        $xlabel_den = "$property^".abs($pow)."*";
      }

      array_push($NXqueries, "Xquery$index");
    }

    $index = 0;
    

    foreach( $property_listy as $property=>$pow){
      $index++;
      $query["Yquery$index"]=$property;
      $power["Yquery$index"]= $pow;
      $final_units["Yquery$index"] = $prop_units["$property"];

      if( $pow > 1){
        $ylabel_num = "$property^$pow *";
      }
      elseif($pow == 1){
        $ylabel_num = "$property*";
      }
      elseif($pow == -1){
        $ylabel_den = "$property*";
      }
      elseif($pow < -1){      
        $ylabel_den = "$property^".abs($pow)."*";
      }
      
      array_push($NYqueries, "Yquery$index");
    }

    $ALLqueries = array_merge($NXqueries,$NYqueries);

    // Prepare final units
 
    $first = 0;
    foreach($NXqueries as $Nquery){
      
      if($first){
        $Xunits =  $Xunits->mul($final_units["$Nquery"]->pow($power["$Nquery"]));
      }
      else{
        $Xunits = $final_units["$Nquery"]->pow($power["$Nquery"]);
        $first = 1;
      }

    }

    $first = 0;
    foreach($NYqueries as $Nquery){
      
      if($first){
        $Yunits =  $Yunits->mul($final_units["$Nquery"]->pow($power["$Nquery"]));
      }
      else{
        $Yunits = $final_units["$Nquery"]->pow($power["$Nquery"]);
        $first = 1;
      }

    }



    // Prepare final labels

    $xlabel_num = substr($xlabel_num,0,-1);
    $xlabel_den = substr($xlabel_den,0,-1);
    $ylabel_num = substr($ylabel_num,0,-1);
    $ylabel_den = substr($ylabel_den,0,-1);

    $xlabel_num = ($xlabel_num != "")?"$xlabel_num":1;
    $ylabel_num = ($ylabel_num != "")?"$ylabel_num":1;

    $xlabel = ($xlabel_den != "")?"$xlabel_num / $xlabel_den":"$xlabel_num";
    $ylabel = ($ylabel_den != "")?"$ylabel_num / $ylabel_den":"$ylabel_num";

    $xlabel .= " (".$Xunits->name.")";
    $ylabel .= " (".$Yunits->name.")";


  $query["categories"]="Material Categories";
  
    $valuesY = array();
    $valuesX = array();
    $valuese = array();
    $categories = array();
    $filenames = array();
  
    $handle = opendir('materials');

//          echo "<br><br>Begin Query <br>";
//          print_r($query);
//          echo "End Query<br>";



    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);


//          echo "<br><br>Begin Values <br>";
//          print_r($values);
//          echo "End Values <br>";

            
            
            // If isset (for all values)

            $allset = 1;

            foreach( $ALLqueries as $Nquery){
              if(isset($values["$Nquery"]))
                $allset *= 1;
                 else
                   $allset *= 0;
            }


            if( $allset )
            {

              // Category part, deal with it later
                $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];
                // Category part, deal with it later
                



                //Force unit match for all queries

                foreach( $NXqueries as $Nquery){
                  force_unit_match($values["Matid"],$final_units["$Nquery"],$values["$Nquery"]);
                }

                foreach( $NYqueries as $Nquery){
                  force_unit_match($values["Matid"],$final_units["$Nquery"],$values["$Nquery"]);
                }

                //Scale units
                //and add them to the array

                $cont = 1;

                // And the Power?

                foreach( $NXqueries as $Nquery){
                  $final_unit=$final_units["$Nquery"]->pow($power["$Nquery"]);
                  $cont *= ( pow($values["$Nquery"]->avg_value,$power["$Nquery"]) /$final_unit->scale);
                }
                $valuesX[]=$cont;
                  
                $cont = 1;

                foreach( $NYqueries as $Nquery){
                  $final_unit=$final_units["$Nquery"]->pow($power["$Nquery"]);
                  $cont *= ( pow($values["$Nquery"]->avg_value,$power["$Nquery"]) /$final_unit->scale);
                }
                $valuesY[]=$cont;

                
//              print_r($valuesX);
//              echo "<br>";
//              print_r($valuesY);
//              echo "<br>";
                

                // Add Matid
                $names[]=$values["Matid"];
                
                // Something else for category
                $categories[]=$category;
                $filenames[]="materials/$file";

            }
        }
    }
    closedir($handle);
  

    //assing label

//     $ylabel = "$property_listy[0]";
//     $xlabel = "$property_listx[0]";

 

  $ret=array(&$valuesX, &$valuesY, &$names, $xlabel, $ylabel, $categories);
  $ret["filenames"]=$filenames;
  return $ret;
}





//var_dump(db_get_plot_values("EvsD","")); 
?>