pt2arch.cxx

Go to the documentation of this file.

/*
   Copyright 2010-2016 Pierre SMARS (smars@yuntech.edu.tw)
   This program is free software: you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation, either version 2 of the License, or
   (at your option) any later version.
 
   This program 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 General Public License for more details.
 
   You should have received a copy of the GNU General Public License
   along with this program.  If not, see <http://www.gnu.org/licenses/>.
   */
#include <fstream>
#include <sstream>
#include <iostream>
#include <string>
#include <vector>
#include <math.h>
 
#include "a_point.h"
#include "SimpleOpt.h"
 
#include "config.h"
 
enum { OPT_HELP, OPT_NAME, OPT_SIDE };
 
//---------------------------------------------------------------------------
CSimpleOpt::SOption g_rgOptions[] =
{
        {OPT_HELP,_T("-?"),SO_NONE},
        {OPT_HELP,_T("-h"),SO_NONE},
        {OPT_HELP,_T("--help"),SO_NONE},
        {OPT_NAME,_T("-n"),SO_REQ_SEP},
        {OPT_NAME,_T("--name"),SO_REQ_SEP},
        {OPT_SIDE,_T("-s"),SO_REQ_SEP},
        {OPT_SIDE,_T("--side"),SO_REQ_SEP},
        SO_END_OF_OPTIONS
};
//---------------------------------------------------------------------------
int error(int val)
{
        std::cerr << "pt2arch:" << std::endl;
        std::cerr << "convert 'x y z' file to an input file to build the arch with calipous tcl module" << std::endl;
        std::cerr << "syntax: pt2arch [-n string|--name string] [-s double|--side double] [-?|-h|--help] thickness depth < point_file > arch_file" << std::endl;
        std::cerr << "    thickness: in the plane of the arch" << std::endl;
        std::cerr << "    depth: perpendicularly to the plane of the arch" << std::endl;
        std::cerr << "    -n name: give the name to the arch" << std::endl;
        std::cerr << "      (default: arch)" << std::endl;
        std::cerr << "    -s code: position of the points on the arch" << std::endl;
        std::cerr << "      -1: intrados (default)" << std::endl;
        std::cerr << "       1: extrados" << std::endl;
        std::cerr << "       0: centre" << std::endl;
        std::cerr << "author: P.Smars, 2010-16" << std::endl;
        std::cerr << "version: 2016-10-28" << std::endl;
        return val;
}
//---------------------------------------------------------------------------
int main( int argc, char *argv[] )
{
        std::string name("arch");
        double pos = -1;
        CSimpleOpt args(argc, argv, g_rgOptions);
        while (args.Next())
        {
                if (args.LastError() == SO_SUCCESS)
                {
                        if (args.OptionId() == OPT_HELP)
                                return error(0);
                        else if (args.OptionId() == OPT_NAME)
                                name = args.OptionArg();
                        else if (args.OptionId() == OPT_SIDE)
                        {
                                std::istringstream i(args.OptionArg());
                                i >> pos;
                        }
                        else
                                return error(-1);
                        // handle option, using OptionId(), OptionText() and OptionArg()
                } else {
                        std::cerr << "Invalid argument: " << args.OptionText() << std::endl;
                        return error(args.LastError());
                        // handle error, one of: SO_OPT_INVALID, SO_OPT_MULTIPLE,
                        // SO_ARG_INVALID, SO_ARG_INVALID_TYPE, SO_ARG_MISSING
                }
        }
        if (args.FileCount() != 2) return error(-2);
        std::ostringstream o;
        for (int k= 0; k < 2; k++) o << args.File(k) << " ";
        std::istringstream i(o.str().c_str());
        double thickness, depth;
        i >> thickness >> depth;
        double ui = -(1+pos)/2.*thickness;
        double us = (1-pos)/2.*thickness;
        int npts;
        std::cin >> npts;
        if (npts<3)
                error(-2);
        std::vector<a_point *> intrados;
        std::vector<a_point *> normals;
        for (int i = 0; i < npts; i++)
        {
                a_point * p = new a_point;
                std::cin >> *p;
                intrados.push_back(p);
        }
        a_point * p1 = intrados[0];
        a_point * p  = intrados[intrados.size()/2];
        a_point * p2 = intrados[intrados.size()-1];
        a_point d1 = (*p-*p1).normalise();
        a_point d2 = (*p2-*p).normalise();
        a_point n = cross(d2,d1).normalise();
        p = intrados[0];
        p2  = intrados[1];
        a_point d = (*p2-*p).normalise();
        a_point * pn = new a_point(cross(n,d).normalise());
        normals.push_back(pn);
        for (int i=1; i<npts-1; i++)
        {
                a_point * p1 = intrados[i-1];
                a_point * p  = intrados[i];
                a_point * p2 = intrados[i+1];
                a_point d1 = (*p-*p1).normalise();
                a_point d2 = (*p2-*p).normalise();
                a_point d = (d1+d2)/2.;
                a_point * pn = new a_point(cross(n,d).normalise());
                normals.push_back(pn);
        }
        p1 = intrados[intrados.size()-2];
        p = intrados[intrados.size()-1];
        d = (*p-*p1).normalise();
        pn = new a_point(cross(n,d).normalise());
        normals.push_back(pn);
        std::cout << "#created by pt2arch!" << std::endl;
        std::cout << "#*********************************************" << std::endl;
        std::cout << "#define structure" << std::endl;
        std::cout << "a_structure " << name << std::endl;
        for (int i = 0; i < npts-1; i++)
        {
                std::cout << "#define block: " << i << std::endl;
                a_point p1 = *intrados[i];
                a_point p2 = *intrados[i+1];
                a_point d1 = *normals[i];
                a_point d2 = *normals[i+1];
                std::cout << "a_point p" << i << "_0 " << p1+d1*us << std::endl;
                std::cout << "a_point p" << i << "_1 " << p1+d1*ui << std::endl;
                std::cout << "a_point p" << i << "_2 " << p2+d2*ui << std::endl;
                std::cout << "a_point p" << i << "_3 " << p2+d2*us << std::endl;
                std::cout << "a_block_2d4 b" << i;
                std::cout << " p" << i << "_0";
                std::cout << " p" << i << "_1";
                std::cout << " p" << i << "_2";
                std::cout << " p" << i << "_3" << std::endl;
                std::cout << "b" << i << " thickness " << depth << std::endl;
                std::cout << name << " add_block b" << i << std::endl;
        }
        std::cout << "#link blocks" << std::endl;
        for (int i = 0; i < npts-2; i++)
                std::cout << name << " link_blocks " << i << " 2 " << i+1 << " 0" << std::endl;
        std::cout << "#exit faces" << std::endl;
        std::cout << "b0 exit 0" << std::endl;
        std::cout << "b" << npts-2 << " exit 2" << std::endl;
        std::cout << "#*********************************************" << std::endl;
        std::cout << "#* edit the following lines to suit your needs" << std::endl;
        std::cout << "#*********************************************" << std::endl;
        std::cout << "set base_face 0" << std::endl;
        std::cout << "a_material mat" << std::endl;
        std::cout << " mat ref 1" << std::endl;
        std::cout << " mat name edit_name" << std::endl;
        std::cout << " mat density 1." << std::endl;
        std::cout << "s material mat" << std::endl;
        std::cout << "a_fcriteria_in c" << std::endl;
        std::cout << " c max 1." << std::endl;
        std::cout << "s criteria c" << std::endl;
        std::cout << "#*********************************************" << std::endl;
        for(std::vector<a_point *>::iterator it = intrados.begin(); it != intrados.end(); ++it) delete *it;
        for(std::vector<a_point *>::iterator it = normals.begin(); it != normals.end(); ++it) delete *it;
        intrados.clear();
        normals.clear();
        return 0;
}