json/example_01.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//               Rapid Optimization Library (ROL) Package
//
// Copyright 2014 NTESS and the ROL contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER

 
 
#include "json/json.h"                    // For JSON definitions
#include "Teuchos_ParameterList.hpp"      // For Parameter List
#include <fstream>
#include <string>
 
 
namespace ROL {
 
void addJSONBlockToPL(const Json::Value& block,ROL::ParameterList& parlist);
void addJSONPairToPL(const Json::Value& block, const std::string &key,ROL::ParameterList& parlist);
 

void addJSONPairToPL(const Json::Value& block,
                     const std::string &key,
                     ROL::ParameterList& parlist) {
 
    Json::Value val = block[key];
 
    if(val.isString()) {
        parlist.set(key,val.asString());
    } else if(val.isBool()) {
        parlist.set(key,val.asBool());
    } else if(val.isInt()) {
        parlist.set(key,val.asInt());
    } else if(val.isUInt()) {
        parlist.set(key,val.asUInt());
    } else if(val.isDouble()) {
        parlist.set(key,val.asDouble());  
    } else if(val.isObject()) { // This is a block. Iterate over its pairs
        addJSONBlockToPL(val,parlist);
    }
      else {
        ROL_TEST_FOR_EXCEPTION( true, std::invalid_argument, ">>> ERROR (addJSONPairToPL, "
                                                             "json value has unsupported type.");
    }
}
 

void addJSONBlockToPL(const Json::Value& block,
                      ROL::ParameterList& parlist) {
    if(block.size()>0) {
        for(Json::ValueIterator itr = block.begin(); itr != block.end(); ++itr) {
            addJSONPairToPL(block,itr.key().asString(),parlist);
        }        
    }
}
 
 

void JSON_Parameters(const std::string& jsonFileName, 
                     ROL::ParameterList& parlist) {
 
Json::Value json;   
std::ifstream stream(jsonFileName, std::ifstream::binary);
stream >> json;
 
if(json.isMember("ROL")) {
    Json::Value rolBlock = json["ROL"];
 
    // Make a flat parameter list from the ROL JSON block
    addJSONBlockToPL(rolBlock,parlist);    
 
    // Check for an algorithm  
    if(rolBlock.isMember("Algorithm")) {
        std::string rolAlgorithm = rolBlock["Algorithm"].asString();
 
        if(rolAlgorithm.find("Trust-Region") != std::string::npos) {
 
            parlist.set("Step Type","Trust-Region");
 
            // Set subproblem solver
            if(rolAlgorithm.find("Cauchy Point") != std::string::npos) {
                parlist.set("Trust-Region Subproblem Solver Type","Cauchy Point");
            } 
            else if(rolAlgorithm.find("Double Dogleg") != std::string::npos) {
                parlist.set("Trust-Region Subproblem Solver Type","Double Dogleg");
            }
            else if(rolAlgorithm.find("Dogleg") != std::string::npos) {
                parlist.set("Trust-Region Subproblem Solver Type","Dogleg");
            }
            else if(rolAlgorithm.find("Truncated CG") != std::string::npos) {
                parlist.set("Trust-Region Subproblem Solver Type","Truncated CG");
            }
 
        }         
        else { // Use Linesearch
            parlist.set("Step Type","Linesearch");
 
            // Set descent type
            if(rolAlgorithm.find("Steepest Descent") != std::string::npos) {
                parlist.set("Descent Type","Steepest Descent");
            } 
            else if(rolAlgorithm.find("Quasi-Newton") != std::string::npos) {
                parlist.set("Descent Type","Quasi-Newton Method");
            }
            else if(rolAlgorithm.find("Newton-Krylov") != std::string::npos) {
                parlist.set("Descent Type","Newton-Krylov");
            }
            else if(rolAlgorithm.find("Nonlinear CG") != std::string::npos) { 
                    parlist.set("Descent Type","Nonlinear CG"); 
                }
 
            }
            
        } 
        else { // No algorithm block found - use defaults
            parlist.set("Step Type","Linesearch");
            parlist.set("Descent Type","Nonlinear CG"); 
        }
 
    }
}
 

template <class Real>
void stepFactory(ROL::ParameterList &parlist,ROL::Ptr<ROL::Step<Real> > &step) {
     
    if(parlist.get("Step Type","Linesearch")=="Trust-Region") {
        step = ROL::makePtr<ROL::TrustRegionStep<Real>>(parlist);
    }
    else {
        step = ROL::makePtr<ROL::LineSearchStep<Real>>(parlist);
         
    }
}
 
}