// ==========================================================================
//                 SeqAn - The Library for Sequence Analysis
// ==========================================================================
// Copyright (c) 2006-2016, Knut Reinert, FU Berlin
// All rights reserved.
//
// Redistribution and use in source and binary forms, with or without
// modification, are permitted provided that the following conditions are met:
//
//     * Redistributions of source code must retain the above copyright
//       notice, this list of conditions and the following disclaimer.
//     * Redistributions in binary form must reproduce the above copyright
//       notice, this list of conditions and the following disclaimer in the
//       documentation and/or other materials provided with the distribution.
//     * Neither the name of Knut Reinert or the FU Berlin nor the names of
//       its contributors may be used to endorse or promote products derived
//       from this software without specific prior written permission.
//
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
// AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
// ARE DISCLAIMED. IN NO EVENT SHALL KNUT REINERT OR THE FU BERLIN BE LIABLE
// FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
// DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
// SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
// CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
// LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
// OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH
// DAMAGE.
//
// ==========================================================================
// Author: Manuel Holtgrewe <manuel.holtgrewe@fu-berlin.de>
// ==========================================================================
// Reading and writing of .intervals file record plus tags.
// ==========================================================================

#ifndef INCLUDE_SEQAN_SIMPLE_INTERVALS_IO_SIMPLE_INTERVALS_IO_H_
#define INCLUDE_SEQAN_SIMPLE_INTERVALS_IO_SIMPLE_INTERVALS_IO_H_

namespace seqan {

// ============================================================================
// Forwards
// ============================================================================

// ============================================================================
// Tags, Classes, Enums
// ============================================================================

// ----------------------------------------------------------------------------
// Class SimpleIntervalsIOContext
// ----------------------------------------------------------------------------

/*!
 * @class SimpleIntervalsIOContext
 * @implements DefaultConstructibleConcept
 * @implements AssignableConcept
 * @headerfile <seqan/simple_intervals_io.h>
 * @brief I/O context for .intervals file.
 *
 * @signature class SimpleIntervalsIOContext;
 */

struct SimpleIntervalsIOContext
{
    CharString buffer;
};

// ----------------------------------------------------------------------------
// Tag SimpleIntervals
// ----------------------------------------------------------------------------

/*!
 * @tag FileFormats#SimpleIntervals
 * @headerfile <seqan/simple_intervals_io.h>
 * @brief Tag for the simple .intervals file format.
 * The file contains format in the format <tt>CHR:POS</tt> or <tt>CHR:BEGIN-END</tt> in 1-based coordinates.
 */

struct SimpleIntervals_;
typedef Tag<SimpleIntervals_> SimpleIntervals;

// ============================================================================
// Metafunctions
// ============================================================================

// ============================================================================
// Functions
// ============================================================================

// ----------------------------------------------------------------------------
// Function readRecord()
// ----------------------------------------------------------------------------

// NOTE(esiragusa): dox disabled.
/*
 * @fn SimpleIntervalsFileIO#readRecord
 * @headerfile <seqan/simple_intervals_io.h>
 * @brief Read an intervals record.
 *
 * @signature void readRecord(region, buffer, iter, tag);
 *
 * @param[out]    region  The @link GenomicRegion @endlink to write interval to.
 * @param[in,out] context A @link SimpleIntervalsIOContext @endlink with buffers.
 * @param[in,out] iter    A @link ForwardIteratorConcept forward iterator @endlink to read from.
 * @param[in]     tag     The tag @link SimpleIntervalsFileIO#SimpleIntervals SimpleIntervals @endlink for the format.
 *
 * @throw IOError in case of problems.
 */

template <typename TForwardIter>
void readRecord(GenomicRegion & record,
                SimpleIntervalsIOContext & context,
                TForwardIter & iter,
                SimpleIntervals const & /*tag*/)
{
    clear(context.buffer);
    // Read next line until the first whitespace (could be '\r' or '\n').
    readUntil(context.buffer, iter, IsWhitespace());
    // Parse context.buffer into GenomicRegion.
    parse(record, context.buffer);
    if (record.endPos == record.INVALID_POS)  // only one position
        record.endPos = record.beginPos + 1;
    // Skip remaining line.
    skipLine(iter);
}

// ----------------------------------------------------------------------------
// Function writeRecord()
// ----------------------------------------------------------------------------

// NOTE(esiragusa): dox disabled.
/*
 * @fn SimpleIntervalsFileIO#writeRecord
 * @headerfile <seqan/simple_intervals_io.h>
 * @brief Write out a @link GenomicRegion @endlink to an intervals file.
 *
 * @signature void writeRecord(target, context, region, tag);
 *
 * @param[in,out] target  An @link OutputIteratorConcept output iterator @endlink to write to.
 * @param[in,out] context A @link SimpleIntervalsIOContext @endlink with buffers.
 * @param[in]     region  The @link GenomicRegion @endlink to write out.
 * @param[in]     tag     The tag @link SimpleIntervalsFileIO#SimpleIntervals SimpleIntervals @endlink for the format.
 *
 * @throw IOError in case of problems.
 */

template <typename TTarget>
void writeRecord(TTarget & target,
                 SimpleIntervalsIOContext & context,
                 GenomicRegion const & region,
                 SimpleIntervals const & /*tag*/)
{
    region.toString(context.buffer);
    write(target, context.buffer);
    write(target, "\n");
}

}  // namespace seqan

#endif  // #ifndef INCLUDE_SEQAN_SIMPLE_INTERVALS_IO_SIMPLE_INTERVALS_IO_H_
