edu.umd.honghongie.RunPersonalizedPageRankBasic.java Source code

Java tutorial

Introduction

Here is the source code for edu.umd.honghongie.RunPersonalizedPageRankBasic.java

Source

/*
 * Cloud9: A Hadoop toolkit for working with big data
 *
 * Licensed under the Apache License, Version 2.0 (the "License"); you
 * may not use this file except in compliance with the License. You may
 * obtain a copy of the License at
 *
 * http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or
 * implied. See the License for the specific language governing
 * permissions and limitations under the License.
 */

package edu.umd.honghongie;

import java.io.IOException;
import java.text.DecimalFormat;
import java.text.NumberFormat;
import java.util.Arrays;
import java.util.Iterator;
import java.lang.Float;

import org.apache.commons.cli.CommandLine;
import org.apache.commons.cli.CommandLineParser;
import org.apache.commons.cli.GnuParser;
import org.apache.commons.cli.HelpFormatter;
import org.apache.commons.cli.Option;
import org.apache.commons.cli.OptionBuilder;
import org.apache.commons.cli.Options;
import org.apache.commons.cli.ParseException;
import org.apache.hadoop.conf.Configuration;
import org.apache.hadoop.conf.Configured;
import org.apache.hadoop.fs.FSDataInputStream;
import org.apache.hadoop.fs.FSDataOutputStream;
import org.apache.hadoop.fs.FileStatus;
import org.apache.hadoop.fs.FileSystem;
import org.apache.hadoop.fs.Path;
import org.apache.hadoop.io.IntWritable;
import org.apache.hadoop.mapreduce.Job;
import org.apache.hadoop.mapreduce.Mapper;
import org.apache.hadoop.mapreduce.Reducer;
import org.apache.hadoop.mapreduce.lib.input.FileInputFormat;
import org.apache.hadoop.mapreduce.lib.output.FileOutputFormat;
import org.apache.hadoop.mapreduce.lib.output.SequenceFileOutputFormat;
import org.apache.hadoop.util.Tool;
import org.apache.hadoop.util.ToolRunner;
import org.apache.log4j.Logger;

import tl.lin.data.array.ArrayListOfInts;
import tl.lin.data.array.ArrayListOfFloats;
import tl.lin.data.array.ArrayListOfIntsWritable;
import tl.lin.data.array.ArrayListOfFloatsWritable;
import tl.lin.data.map.HMapIF;
import tl.lin.data.map.MapIF;

import com.google.common.base.Preconditions;

import edu.umd.cloud9.mapreduce.lib.input.NonSplitableSequenceFileInputFormat;

/**
 * <p>
 * Main driver program for running the basic (non-Schimmy) implementation of
 * PageRank.
 * </p>
 *
 * <p>
 * The starting and ending iterations will correspond to paths
 * <code>/base/path/iterXXXX</code> and <code>/base/path/iterYYYY</code>. As a
 * example, if you specify 0 and 10 as the starting and ending iterations, the
 * driver program will start with the graph structure stored at
 * <code>/base/path/iter0000</code>; final results will be stored at
 * <code>/base/path/iter0010</code>.
 * </p>
 * @author Lingzi Hong
 */
public class RunPersonalizedPageRankBasic extends Configured implements Tool {
    private static final Logger LOG = Logger.getLogger(RunPersonalizedPageRankBasic.class);

    //*** for global counters, used in context.getCounter
    private static enum PageRank {
        nodes, edges, massMessages, massMessagesSaved, massMessagesReceived, missingStructure
    };

    // Mapper, no in-mapper combining.
    private static class MapClass extends Mapper<IntWritable, PageRankNode, IntWritable, PageRankNode> {

        // The neighbor to which we're sending messages.
        private static final IntWritable neighbor = new IntWritable();

        // Contents of the messages: partial PageRank mass.
        private static final PageRankNode intermediateMass = new PageRankNode();

        // For passing along node structure.
        private static final PageRankNode intermediateStructure = new PageRankNode();

        @Override
        public void map(IntWritable nid, PageRankNode node, Context context)
                throws IOException, InterruptedException {
            // Pass along node structure.
            intermediateStructure.setNodeId(node.getNodeId());
            intermediateStructure.setType(PageRankNode.Type.Structure);
            intermediateStructure.setAdjacencyList(node.getAdjacenyList());

            context.write(nid, intermediateStructure);

            int massMessages = 0;

            // Distribute PageRank mass to neighbors (along outgoing edges).
            if (node.getAdjacenyList().size() > 0) {
                // Each neighbor gets an equal share of PageRank mass.
                ArrayListOfIntsWritable list = node.getAdjacenyList();

                //save values of intermediate pageranks mass
                ArrayListOfFloatsWritable mass = new ArrayListOfFloatsWritable();

                // ***If no pagerank value, thus is infinite, no distribution
                // pageranks is used to get the value, mass for set the value
                ArrayListOfFloatsWritable pageranks = node.getPageRank();

                for (int i = 0; i < pageranks.size(); i++) {
                    if (!Float.isInfinite(pageranks.get(i))) {
                        mass.add(pageranks.get(i) - (float) StrictMath.log(list.size()));
                    } else {
                        mass.add(Float.NEGATIVE_INFINITY);
                    }
                }

                context.getCounter(PageRank.edges).increment(list.size());

                // Iterate over neighbors.
                for (int i = 0; i < list.size(); i++) {
                    neighbor.set(list.get(i));
                    intermediateMass.setNodeId(list.get(i));
                    intermediateMass.setType(PageRankNode.Type.Mass);
                    intermediateMass.setPageRank(mass);

                    // Emit messages with PageRank mass to neighbors. 
                    // Emit neighbors and intermediateMass they can get from this node 
                    context.write(neighbor, intermediateMass);
                    massMessages++;
                }
            }

            // Bookkeeping.
            context.getCounter(PageRank.nodes).increment(1);
            context.getCounter(PageRank.massMessages).increment(massMessages);
        }
    }

    // Combiner: sums partial PageRank contributions and passes node structure along.
    private static class CombineClass extends Reducer<IntWritable, PageRankNode, IntWritable, PageRankNode> {
        private static final PageRankNode intermediateMass = new PageRankNode();

        private static int numberofsource;

        public void setup(Reducer<IntWritable, PageRankNode, IntWritable, PageRankNode>.Context context) {
            String tt = context.getConfiguration().get("SourceNode", null);
            System.out.println(tt);

            if (tt != null) {
                String ss[] = tt.split(",");
                numberofsource = ss.length;
            }
        }

        @Override
        public void reduce(IntWritable nid, Iterable<PageRankNode> values, Context context)
                throws IOException, InterruptedException {

            int massMessages = 0;

            // Remember, PageRank mass is stored as a log prob.
            ArrayListOfFloatsWritable mass = new ArrayListOfFloatsWritable();
            for (int i = 0; i < numberofsource; i++) {
                mass.add(Float.NEGATIVE_INFINITY);
            }

            for (PageRankNode n : values) {
                if (n.getType() == PageRankNode.Type.Structure) {
                    // Simply pass along node structure.
                    context.write(nid, n);
                } else {
                    // Accumulate PageRank mass contributions.
                    ArrayListOfFloatsWritable pageranklist = n.getPageRank();
                    //***************** get pagerank values and update it
                    for (int i = 0; i < mass.size(); i++) {
                        float pi = pageranklist.get(i);
                        mass.set(i, sumLogProbs(mass.get(i), pi));
                    }
                    massMessages++;
                }
            }

            // Emit aggregated results.
            if (massMessages > 0) {
                intermediateMass.setNodeId(nid.get());
                intermediateMass.setType(PageRankNode.Type.Mass);
                intermediateMass.setPageRank(mass);

                context.write(nid, intermediateMass); //add part of the intermediate mass up
            }
        }
    }

    // Reduce: sums incoming PageRank contributions, rewrite graph structure.
    private static class ReduceClass extends Reducer<IntWritable, PageRankNode, IntWritable, PageRankNode> {
        // For keeping track of PageRank mass encountered, so we can compute missing PageRank mass lost
        // through dangling nodes.
        private static int numberofsource = 0;
        private ArrayListOfFloatsWritable totalMass = new ArrayListOfFloatsWritable(); //*****************need to change, only float writable can be write out

        public void setup(Reducer<IntWritable, PageRankNode, IntWritable, PageRankNode>.Context context) {
            String tt = context.getConfiguration().get("SourceNode", null);

            if (tt != null) {
                String ss[] = tt.split(",");
                numberofsource = ss.length;
            }

            for (int i = 0; i < numberofsource; i++) {
                totalMass.add(Float.NEGATIVE_INFINITY);
            }
        }

        @Override
        public void reduce(IntWritable nid, Iterable<PageRankNode> iterable, Context context)
                throws IOException, InterruptedException {
            Iterator<PageRankNode> values = iterable.iterator();

            // Create the node structure that we're going to assemble back together from shuffled pieces.
            PageRankNode node = new PageRankNode();

            node.setType(PageRankNode.Type.Complete);
            node.setNodeId(nid.get());

            int massMessagesReceived = 0;
            int structureReceived = 0;

            ArrayListOfFloatsWritable mass = new ArrayListOfFloatsWritable();
            for (int i = 0; i < numberofsource; i++) {
                mass.add(Float.NEGATIVE_INFINITY);
            }

            while (values.hasNext()) {
                PageRankNode n = values.next();

                if (n.getType().equals(PageRankNode.Type.Structure)) {
                    // This is the structure; update accordingly.
                    ArrayListOfIntsWritable list = n.getAdjacenyList();
                    structureReceived++;

                    node.setAdjacencyList(list);
                } else {
                    // This is a message that contains PageRank mass; accumulate.
                    ArrayListOfFloatsWritable pagerankvalues = n.getPageRank();
                    for (int i = 0; i < pagerankvalues.size(); i++) {
                        mass.set(i, sumLogProbs(mass.get(i), pagerankvalues.get(i)));
                    }
                    massMessagesReceived++;
                }
            }

            // Update the final accumulated PageRank mass.
            node.setPageRank(mass);
            context.getCounter(PageRank.massMessagesReceived).increment(massMessagesReceived);

            // Error checking.
            if (structureReceived == 1) {
                // Everything checks out, emit final node structure with updated PageRank value.
                context.write(nid, node);

                // Keep track of total PageRank mass.
                for (int i = 0; i < numberofsource; i++) {
                    totalMass.set(i, sumLogProbs(totalMass.get(i), mass.get(i)));
                }

            } else if (structureReceived == 0) {
                // We get into this situation if there exists an edge pointing to a node which has no
                // corresponding node structure (i.e., PageRank mass was passed to a non-existent node)...
                // log and count but move on.
                context.getCounter(PageRank.missingStructure).increment(1);
                LOG.warn("No structure received for nodeid: " + nid.get() + " mass: " + massMessagesReceived);
                // It's important to note that we don't add the PageRank mass to total... if PageRank mass
                // was sent to a non-existent node, it should simply vanish.
            } else {
                // This shouldn't happen!
                throw new RuntimeException("Multiple structure received for nodeid: " + nid.get() + " mass: "
                        + massMessagesReceived + " struct: " + structureReceived);
            }
        }

        @Override
        public void cleanup(Context context) throws IOException {
            Configuration conf = context.getConfiguration();
            String taskId = conf.get("mapred.task.id");
            String path = conf.get("PageRankMassPath");

            Preconditions.checkNotNull(taskId);
            Preconditions.checkNotNull(path);

            // Write to a file the amount of PageRank mass we've seen in this reducer.
            // how to resolve ArrayListOfFloats?
            FileSystem fs = FileSystem.get(context.getConfiguration());
            FSDataOutputStream out = fs.create(new Path(path + "/" + taskId), false);
            totalMass.write(out);
            out.close();
        }
    }

    // Mapper that distributes the missing PageRank mass (lost at the dangling nodes) and takes care
    // of the random jump factor.
    //**** personalized page rank always jump to the source

    private static class MapPageRankMassDistributionClass
            extends Mapper<IntWritable, PageRankNode, IntWritable, PageRankNode> {

        private ArrayListOfFloats missingMass = new ArrayListOfFloats(); //to get missingmass from file
        private int nodeCnt = 0;
        private ArrayListOfInts sid = new ArrayListOfInts(); //to save sources' id

        @Override
        public void setup(Context context) throws IOException {
            Configuration conf = context.getConfiguration();
            //******** missingvalues as string and resolve
            //******** sid as string and resolve
            //******** get rid of "[" and "]"
            String missingMassvalues = conf.get("MissingMass", null);
            String sidvalues = conf.get("SourceNode", null);

            missingMassvalues = missingMassvalues.substring(1, missingMassvalues.length() - 1);
            sidvalues = sidvalues.substring(1, sidvalues.length() - 1);

            nodeCnt = conf.getInt("NodeCount", 0);

            String tt[] = missingMassvalues.split(",");
            String ss[] = sidvalues.split(",");

            for (int i = 0; i < ss.length; i++) {
                missingMass.add(Float.parseFloat(tt[i].trim()));
                sid.add(Integer.parseInt(ss[i].trim()));
            }
        }

        @Override
        public void map(IntWritable nid, PageRankNode node, Context context)
                throws IOException, InterruptedException {
            // only jump to nid, all missingmass to sid, use nodenum
            int nodenum = nid.get();
            ArrayListOfFloatsWritable pagerankvalues = node.getPageRank();

            for (int i = 0; i < pagerankvalues.size(); i++) {
                if (nodenum == sid.get(i)) {
                    float p = pagerankvalues.get(i);
                    float jump = (float) (Math.log(ALPHA));
                    float link = (float) Math.log(1.0f - ALPHA)
                            + sumLogProbs(p, (float) (Math.log(missingMass.get(i))));
                    p = sumLogProbs(jump, link);
                    pagerankvalues.set(i, p);
                } else {
                    float p = pagerankvalues.get(i);
                    float link = (float) Math.log(1.0f - ALPHA) + p;
                    p = link;
                    pagerankvalues.set(i, p);
                }
            }
            node.setPageRank(pagerankvalues);
            context.write(nid, node);
        }
    }

    // Random jump factor.
    private static float ALPHA = 0.15f;
    private static NumberFormat formatter = new DecimalFormat("0000");

    /**
     * Dispatches command-line arguments to the tool via the {@code ToolRunner}.
     */
    public static void main(String[] args) throws Exception {
        ToolRunner.run(new RunPersonalizedPageRankBasic(), args);
    }

    public RunPersonalizedPageRankBasic() {
    }

    private static final String BASE = "base";
    private static final String NUM_NODES = "numNodes";
    private static final String START = "start";
    private static final String END = "end";
    private static final String COMBINER = "useCombiner";
    private static final String RANGE = "range";
    //*** add source
    private static final String SOURCE = "sources";

    /**
     * Runs this tool.
     */
    @SuppressWarnings({ "static-access" })
    public int run(String[] args) throws Exception {
        Options options = new Options();

        options.addOption(new Option(COMBINER, "use combiner"));
        options.addOption(new Option(RANGE, "use range partitioner"));

        options.addOption(OptionBuilder.withArgName("path").hasArg().withDescription("base path").create(BASE));
        options.addOption(
                OptionBuilder.withArgName("num").hasArg().withDescription("start iteration").create(START));
        options.addOption(OptionBuilder.withArgName("num").hasArg().withDescription("end iteration").create(END));
        options.addOption(
                OptionBuilder.withArgName("num").hasArg().withDescription("number of nodes").create(NUM_NODES));

        //*** add option of source
        options.addOption(
                OptionBuilder.withArgName("node").hasArg().withDescription("source nodes").create(SOURCE));

        CommandLine cmdline;
        CommandLineParser parser = new GnuParser();

        try {
            cmdline = parser.parse(options, args);
        } catch (ParseException exp) {
            System.err.println("Error parsing command line: " + exp.getMessage());
            return -1;
        }

        if (!cmdline.hasOption(BASE) || !cmdline.hasOption(START) || !cmdline.hasOption(END)
                || !cmdline.hasOption(NUM_NODES) || !cmdline.hasOption(SOURCE)) {
            System.out.println("args: " + Arrays.toString(args));
            HelpFormatter formatter = new HelpFormatter();
            formatter.setWidth(120);
            formatter.printHelp(this.getClass().getName(), options);
            ToolRunner.printGenericCommandUsage(System.out);
            return -1;
        }

        String basePath = cmdline.getOptionValue(BASE);
        int n = Integer.parseInt(cmdline.getOptionValue(NUM_NODES));
        int s = Integer.parseInt(cmdline.getOptionValue(START));
        int e = Integer.parseInt(cmdline.getOptionValue(END));
        boolean useCombiner = cmdline.hasOption(COMBINER);
        boolean useRange = cmdline.hasOption(RANGE);
        //********
        String nodeid = cmdline.getOptionValue(SOURCE);

        ArrayListOfInts sourceids = new ArrayListOfInts();
        String ss[] = nodeid.split(",");
        for (int i = 0; i < ss.length; i++) {
            sourceids.add(Integer.parseInt(ss[i]));
        }
        System.out.println("***** 1 **** sourceids**** " + sourceids.toString());

        LOG.info("Tool name: RunPageRank");
        LOG.info(" - base path: " + basePath);
        LOG.info(" - num nodes: " + n);
        LOG.info(" - start iteration: " + s);
        LOG.info(" - end iteration: " + e);
        LOG.info(" - use combiner: " + useCombiner);
        LOG.info(" - user range partitioner: " + useRange);
        LOG.info(" - source node: " + nodeid);

        // Iterate PageRank. *** add source
        for (int i = s; i < e; i++) {
            iteratePageRank(i, i + 1, basePath, n, sourceids, useCombiner);
        }

        return 0;
    }

    // Run each iteration.
    //*************************
    private void iteratePageRank(int i, int j, String basePath, int numNodes, ArrayListOfInts sourceids,
            boolean useCombiner) throws Exception {
        // Each iteration consists of two phases (two MapReduce jobs).

        // Job 1: distribute PageRank mass along outgoing edges.
        ArrayListOfFloats mass = phase1(i, j, basePath, numNodes, sourceids, useCombiner);

        System.out.println("*********** 1 ***********" + mass.toString());

        // Find out how much PageRank mass got lost at the dangling nodes.
        ArrayListOfFloats missing = new ArrayListOfFloats();

        for (int k = 0; k < sourceids.size(); k++) {
            missing.add(1.0f - (float) StrictMath.exp(mass.get(k)));
        }

        // Job 2: distribute missing mass, take care of random jump factor.
        phase2(i, j, missing, basePath, numNodes, sourceids);
    }

    //***************************
    private ArrayListOfFloats phase1(int i, int j, String basePath, int numNodes, ArrayListOfInts sourceids,
            boolean useCombiner) throws Exception {
        Job job = Job.getInstance(getConf());
        job.setJobName("PageRank:Basic:iteration" + j + ":Phase1");
        job.setJarByClass(RunPersonalizedPageRankBasic.class);

        String in = basePath + "/iter" + formatter.format(i);
        String out = basePath + "/iter" + formatter.format(j) + "t";
        String outm = out + "-mass";

        // We need to actually count the number of part files to get the number of partitions (because
        // the directory might contain _log).
        int numPartitions = 0;
        for (FileStatus s : FileSystem.get(getConf()).listStatus(new Path(in))) {
            if (s.getPath().getName().contains("part-"))
                numPartitions++;
        }

        LOG.info("PageRank: iteration " + j + ": Phase1");
        LOG.info(" - input: " + in);
        LOG.info(" - output: " + out);
        LOG.info(" - nodeCnt: " + numNodes);
        LOG.info(" - useCombiner: " + useCombiner);
        LOG.info("computed number of partitions: " + numPartitions);

        int numReduceTasks = numPartitions;

        job.getConfiguration().setInt("NodeCount", numNodes);
        job.getConfiguration().setBoolean("mapred.map.tasks.speculative.execution", false);
        job.getConfiguration().setBoolean("mapred.reduce.tasks.speculative.execution", false);
        //job.getConfiguration().set("mapred.child.java.opts", "-Xmx2048m");
        job.getConfiguration().set("PageRankMassPath", outm);

        //*********************** reduer uses sourcenode
        job.getConfiguration().set("SourceNode", sourceids.toString());

        job.setNumReduceTasks(numReduceTasks);

        FileInputFormat.setInputPaths(job, new Path(in));
        FileOutputFormat.setOutputPath(job, new Path(out));

        job.setInputFormatClass(NonSplitableSequenceFileInputFormat.class);
        job.setOutputFormatClass(SequenceFileOutputFormat.class);

        job.setMapOutputKeyClass(IntWritable.class);
        job.setMapOutputValueClass(PageRankNode.class);

        job.setOutputKeyClass(IntWritable.class);
        job.setOutputValueClass(PageRankNode.class);

        job.setMapperClass(MapClass.class);

        if (useCombiner) {
            job.setCombinerClass(CombineClass.class);
        }

        job.setReducerClass(ReduceClass.class);

        FileSystem.get(getConf()).delete(new Path(out), true);
        FileSystem.get(getConf()).delete(new Path(outm), true);

        long startTime = System.currentTimeMillis();
        job.waitForCompletion(true);
        System.out.println("Job Finished in " + (System.currentTimeMillis() - startTime) / 1000.0 + " seconds");
        System.out.println("********** 1 *********");

        ArrayListOfFloats mass = new ArrayListOfFloats();
        int length = sourceids.size();

        System.out.println("*********** 1 **********" + length);
        float test = Float.NEGATIVE_INFINITY;
        for (int k = 0; k < length; k++) {
            mass.add(Float.NEGATIVE_INFINITY); //use add to initialize
        }
        System.out.println("********** test ********" + test);
        System.out.println("******** 1 ********" + mass);

        //****************************************** how to resolve datastream
        FileSystem fs = FileSystem.get(getConf());
        ArrayListOfFloatsWritable invalue = new ArrayListOfFloatsWritable();
        for (FileStatus f : fs.listStatus(new Path(outm))) {
            FSDataInputStream fin = fs.open(f.getPath());

            //**************************************  get all values from fin?
            invalue.readFields(fin);
            System.out.println("************** 1 ************" + invalue);
            for (int k = 0; k < invalue.size(); k++) {
                mass.set(k, sumLogProbs(mass.get(k), invalue.get(k)));
            }
            fin.close();
        }

        System.out.println("******** 1 ********" + mass.toString());
        return mass;

    }

    private void phase2(int i, int j, ArrayListOfFloats missing, String basePath, int numNodes,
            ArrayListOfInts sourceids) throws Exception {

        System.out.println("************ start working in phase2***********");

        Job job = Job.getInstance(getConf());
        job.setJobName("PageRank:Basic:iteration" + j + ":Phase2");
        job.setJarByClass(RunPersonalizedPageRankBasic.class);

        LOG.info("missing PageRank mass: " + missing.toString());
        LOG.info("number of nodes: " + numNodes);

        //********
        LOG.info("source node: " + sourceids.toString());

        String in = basePath + "/iter" + formatter.format(j) + "t";
        String out = basePath + "/iter" + formatter.format(j);

        LOG.info("PageRank: iteration " + j + ": Phase2");
        LOG.info(" - input: " + in);
        LOG.info(" - output: " + out);

        job.getConfiguration().setBoolean("mapred.map.tasks.speculative.execution", false);
        job.getConfiguration().setBoolean("mapred.reduce.tasks.speculative.execution", false);
        job.getConfiguration().set("MissingMass", missing.toString());
        job.getConfiguration().setInt("NodeCount", numNodes);
        //*********
        job.getConfiguration().set("SourceNode", sourceids.toString());

        job.setNumReduceTasks(0);

        FileInputFormat.setInputPaths(job, new Path(in));
        FileOutputFormat.setOutputPath(job, new Path(out));

        job.setInputFormatClass(NonSplitableSequenceFileInputFormat.class);
        job.setOutputFormatClass(SequenceFileOutputFormat.class);

        job.setMapOutputKeyClass(IntWritable.class);
        job.setMapOutputValueClass(PageRankNode.class);

        job.setOutputKeyClass(IntWritable.class);
        job.setOutputValueClass(PageRankNode.class);

        job.setMapperClass(MapPageRankMassDistributionClass.class);

        FileSystem.get(getConf()).delete(new Path(out), true);

        long startTime = System.currentTimeMillis();
        job.waitForCompletion(true);
        System.out.println("Job Finished in " + (System.currentTimeMillis() - startTime) / 1000.0 + " seconds");
    }

    // Adds two log probs. why is the formula?
    private static float sumLogProbs(float a, float b) {
        if (a == Float.NEGATIVE_INFINITY)
            return b;

        if (b == Float.NEGATIVE_INFINITY)
            return a;

        if (a < b) {
            return (float) (b + StrictMath.log1p(StrictMath.exp(a - b)));
        }
        return (float) (a + StrictMath.log1p(StrictMath.exp(b - a)));
    }
}