LOOP {
iam (ERP);
ERP_drawScalp (me, GRAPHICS, GET_REAL (U"left Time range"), GET_REAL (U"right Time range"),
GET_REAL (U"left Voltage range"), GET_REAL (U"right Voltage range"), GET_ENUM (kGraphics_colourScale, U"Colour scale"), GET_INTEGER (U"Garnish"));
}
END
FORM (ERP_drawScalp_garnish, U"ERP: Draw scalp (garnish)", 0)
REAL (U"left Voltage range (V)", U"10e-6")
REAL (U"right Voltage range", U"-10e-6")
RADIO_ENUM (U"Colour scale", kGraphics_colourScale, BLUE_TO_RED)
OK
DO
autoPraatPicture picture;
ERP_drawScalp_garnish (GRAPHICS,
GET_REAL (U"left Voltage range"), GET_REAL (U"right Voltage range"), GET_ENUM (kGraphics_colourScale, U"Colour scale"));
END
FORM (ERP_extractOneChannelAsSound, U"ERP: Extract one channel as Sound", 0)
WORD (U"Channel name", U"Cz")
OK
DO
LOOP {
iam (ERP);
const char32 *channelName = GET_STRING (U"Channel name");
long channelNumber = ERP_getChannelNumber (me, channelName);
if (channelNumber == 0) Melder_throw (me, U": no channel named \"", channelName, U"\".");
autoSound thee = Sound_extractChannel (me, channelNumber);
praat_new (thee.transfer(), my name, U"_", channelName);
}
END
iam (ERP);
ERP_drawScalp (me, GRAPHICS, GET_REAL (U"left Time range"), GET_REAL (U"right Time range"),
GET_REAL (U"left Voltage range"), GET_REAL (U"right Voltage range"), GET_ENUM (kGraphics_colourScale, U"Colour scale"), GET_INTEGER (U"Garnish"));
}
END
}
FORM (GRAPHICS_ERP_drawScalp_garnish, U"ERP: Draw scalp (garnish)", nullptr) {
REAL (U"left Voltage range (V)", U"10e-6")
REAL (U"right Voltage range", U"-10e-6")
RADIO_ENUM (U"Colour scale", kGraphics_colourScale, BLUE_TO_RED)
OK
DO
autoPraatPicture picture;
ERP_drawScalp_garnish (GRAPHICS,
GET_REAL (U"left Voltage range"), GET_REAL (U"right Voltage range"), GET_ENUM (kGraphics_colourScale, U"Colour scale"));
END
}
// MARK: Query
FORM (STRING_ERP_getChannelName, U"Get channel name", nullptr) {
NATURAL (U"Channel number", U"1")
OK
DO
LOOP {
iam (ERP);
long channelNumber = GET_INTEGER (U"Channel number");
if (channelNumber > my ny)
Melder_throw (me, U": there are only ", my ny, U" channels.");
Melder_information (my channelNames [channelNumber]);
const char32 *text = GET_STRING (U"...the text");
LOOP {
iam (TableOfReal);
autoTableOfReal thee = TableOfReal_extractColumnsWhereLabel (me, GET_ENUM (kMelder_string, U"Extract all columns whose label..."), text);
praat_new (thee.move(), my name, U"_", text);
}
END }
FORM (NEW_TableOfReal_extractColumnsWhereRow, U"Extract columns where row", nullptr) {
NATURAL (U"Extract all columns where row...", U"1")
OPTIONMENU_ENUM (U"...is...", kMelder_number, DEFAULT)
REAL (U"...the value", U"0.0")
OK
DO
long row = GET_INTEGER (U"Extract all columns where row...");
double value = GET_REAL (U"...the value");
LOOP {
iam (TableOfReal);
autoTableOfReal thee = TableOfReal_extractColumnsWhereRow (me, row, GET_ENUM (kMelder_number, U"...is..."), value);
praat_new (thee.move(), my name, U"_", row, U"_", lround (value));
}
END }
DIRECT (NEW_TableOfReal_extractRowLabelsAsStrings) {
LOOP {
iam (TableOfReal);
autoStrings thee = TableOfReal_extractRowLabelsAsStrings (me);
praat_new (thee.move(), my name);
}
END }
/**
* @brief Decompress the compressed list in given packet
*
* @param decomp The list decompressor
* @param packet The ROHC packet to decompress
* @param packet_len The remaining length of the packet to decode (in bytes)
* @return The size of the compressed list in packet in case of
* success, -1 in case of failure
*/
static int rohc_list_decode(struct list_decomp *decomp,
const uint8_t *packet,
size_t packet_len)
{
size_t read_length = 0;
uint8_t et; /* the type of list encoding */
bool gp; /* whether the gen_id field is present or not */
uint8_t ps; /* the type of XI field */
uint8_t m; /* the CC or Count field (share bits with XI 1) */
uint8_t xi_1; /* the XI 1 field (share bits with m) */
unsigned int gen_id; /* the gen_id if present,
ROHC_LIST_GEN_ID_ANON otherwise */
int ret;
/* reset the list of the current packet */
rohc_list_reset(&decomp->pkt_list);
/* is there enough data in packet for the ET, PS, m/XI1 and gen_id
* fields? */
if(packet_len < 2)
{
rd_list_warn(decomp, "packet too small for compressed list (only %zu "
"bytes while at least 2 bytes are required)", packet_len);
goto error;
}
/* parse ET, GP, PS, and m/XI1 fields */
et = GET_BIT_6_7(packet);
gp = !!GET_BIT_5(packet);
ps = GET_REAL(GET_BIT_4(packet));
m = GET_BIT_0_3(packet);
xi_1 = m; /* m and XI 1 are the same field */
packet++;
read_length++;
packet_len--;
rd_list_debug(decomp, "ET = %d, GP = %d, PS = %d, m = XI 1 = %d",
et, gp, ps, m);
assert(m <= ROHC_LIST_ITEMS_MAX);
/* parse gen_id if present */
if(gp == 1)
{
gen_id = GET_BIT_0_7(packet);
packet++;
read_length++;
packet_len--;
rd_list_debug(decomp, "gen_id = 0x%02x", gen_id);
}
else
{
gen_id = ROHC_LIST_GEN_ID_ANON;
rd_list_debug(decomp, "decode anonymous list");
}
decomp->pkt_list.id = gen_id;
/* decode the compressed list according to its type */
switch(et)
{
case 0:
ret = rohc_list_decode_type_0(decomp, packet, packet_len,
gen_id, ps, m);
break;
case 1:
ret = rohc_list_decode_type_1(decomp, packet, packet_len,
gen_id, ps, xi_1);
break;
case 2:
ret = rohc_list_decode_type_2(decomp, packet, packet_len, gen_id);
break;
case 3:
ret = rohc_list_decode_type_3(decomp, packet, packet_len,
gen_id, ps, xi_1);
break;
default:
/* should not happen */
rohc_error(decomp, ROHC_TRACE_DECOMP, decomp->profile_id,
"unknown type of compressed list (ET = %u)", et);
assert(0);
goto error;
}
if(ret < 0)
{
rd_list_warn(decomp, "failed to decode compressed list type %d", et);
goto error;
}
assert(((size_t) ret) <= packet_len);
#ifndef __clang_analyzer__ /* silent warning about dead in/decrement */
packet += ret;
packet_len -= ret;
#endif
read_length += ret;
/* RFC3095, section 5.8.2.1 reads:
* When the decompressor receives a compressed list, it retrieves the
* proper ref_list from the sliding window based on the ref_id, and
* decompresses the compressed list obtaining curr_list.
* In U/O-mode, curr_list is inserted into the sliding window
* together with its generation identifier if the compressed list had
* a generation identifier and the sliding window does not contain a
* list with that generation identifier. All lists with generations
* older than ref_id are removed from the sliding window. */
if(gen_id == ROHC_LIST_GEN_ID_ANON)
{
/* list is not identified by a gen_id, so do not update the sliding
* window of lists */
rd_list_debug(decomp, "anonymous list was received");
}
else if(decomp->lists[gen_id].counter > 0)
{
/* list is identified by a gen_id, but the sliding window of lists
* already contain a list with that generation identifier, so do
* not update the sliding window of lists */
decomp->lists[gen_id].counter++;
rd_list_debug(decomp, "list with gen_id %u is already present in "
"reference lists (received for the #%zu times)",
gen_id, decomp->lists[gen_id].counter);
}
else
{
/* list is identified by a gen_id and the sliding window of lists does
* not contain a list with that generation identifier yet, so update
* the sliding window of lists */
rd_list_debug(decomp, "list with gen_id %u is not present yet in "
"reference lists, add it", gen_id);
memcpy(decomp->lists[gen_id].items, decomp->pkt_list.items,
ROHC_LIST_ITEMS_MAX * sizeof(struct decomp_list *));
decomp->lists[gen_id].items_nr = decomp->pkt_list.items_nr;
decomp->lists[gen_id].counter = 1;
/* TODO: remove all lists with gen_id < ref_id */
}
return read_length;
error:
return -1;
}
Editor_save (me, U"Add point");
RealTier_addPoint ((RealTier) my data, 0.5 * (my d_startSelection + my d_endSelection), my ycursor);
RealTierEditor_updateScaling (me);
FunctionEditor_redraw (me);
Editor_broadcastDataChanged (me);
}
static void menu_cb_addPointAt (RealTierEditor me, EDITOR_ARGS_FORM) {
EDITOR_FORM (U"Add point", 0)
REAL (U"Time (s)", U"0.0")
REAL (my v_quantityText (), U"0.0")
EDITOR_OK
SET_REAL (U"Time", 0.5 * (my d_startSelection + my d_endSelection))
SET_REAL (my v_quantityKey (), my ycursor)
EDITOR_DO
double desiredValue = GET_REAL (my v_quantityKey ());
if (NUMdefined (my v_minimumLegalValue ()) && desiredValue < my v_minimumLegalValue ())
Melder_throw (U"Cannot add a point below ", my v_minimumLegalValue (), my v_rightTickUnits (), U".");
if (NUMdefined (my v_maximumLegalValue ()) && desiredValue > my v_maximumLegalValue ())
Melder_throw (U"Cannot add a point above ", my v_maximumLegalValue (), my v_rightTickUnits (), U".");
Editor_save (me, U"Add point");
RealTier_addPoint ((RealTier) my data, GET_REAL (U"Time"), desiredValue);
RealTierEditor_updateScaling (me);
FunctionEditor_redraw (me);
Editor_broadcastDataChanged (me);
EDITOR_END
}
static void menu_cb_setRange (RealTierEditor me, EDITOR_ARGS_FORM) {
EDITOR_FORM (my v_setRangeTitle (), 0)
REAL (my v_yminText (), my v_defaultYminText ())
FORM (INTEGER_TimeTier_getLowIndexFromTime, U"Get low index", U"AnyTier: Get low index from time...") {
REALVAR (time, U"Time (s)", U"0.5")
OK
DO
FIND_ONE (AnyTier)
Melder_information (my points.size == 0 ? U"--undefined--" : Melder_integer (AnyTier_timeToLowIndex (me, time)));
END
}
FORM (INTEGER_TimeTier_getHighIndexFromTime, U"Get high index", U"AnyTier: Get high index from time...") {
REALVAR (time, U"Time (s)", U"0.5")
OK
DO
FIND_ONE (AnyTier)
Melder_information (my points.size == 0 ? U"--undefined--" : Melder_integer (AnyTier_timeToHighIndex (me, GET_REAL (U"Time"))));
END
}
FORM (INTEGER_TimeTier_getNearestIndexFromTime, U"Get nearest index", U"AnyTier: Get nearest index from time...") {
REALVAR (time, U"Time (s)", U"0.5")
OK
DO
FIND_ONE (AnyTier)
Melder_information (my points.size == 0 ? U"--undefined--" : Melder_integer (AnyTier_timeToNearestIndex (me, GET_REAL (U"Time"))));
END
}
FORM (REAL_TimeTier_getTimeFromIndex, U"Get time", nullptr /*"AnyTier: Get time from index..."*/) {
NATURALVAR (pointNumber, U"Point number", U"10")
OK
FORM (FFNet_Eigen_drawIntersection, L"FFnet & Eigen: Draw intersection", 0) INTEGER (L"X-component", L"1") INTEGER (L"Y-component", L"2") REAL (L"xmin", L"0.0") REAL (L"xmax", L"0.0") REAL (L"ymin", L"0.0") REAL (L"ymax", L"0.0") OK DO FFNet me = FIRST (FFNet); Eigen thee = FIRST (Eigen); autoPraatPicture picture; long pcx = GET_INTEGER (L"X-component"); long pcy = GET_INTEGER (L"Y-component"); REQUIRE (pcx != 0 && pcy != 0, L"X and Y component must differ from 0.") FFNet_Eigen_drawIntersection (me, thee, GRAPHICS, pcx, pcy, GET_REAL (L"xmin"), GET_REAL (L"xmax"), GET_REAL (L"ymin"), GET_REAL (L"ymax")); END FORM (FFNet_PCA_drawDecisionPlaneInEigenspace, L"FFNet & PCA: Draw decision plane", L"") NATURAL (L"Unit number", L"1") NATURAL (L"Layer number", L"1") NATURAL (L"Horizontal eigenvector number", L"1") NATURAL (L"Vertical eigenvector number", L"2") REAL (L"left Horizontal range", L"0.0") REAL (L"right Horizontal range", L"0.0") REAL (L"left Vertical range", L"0.0") REAL (L"right Vertical range", L"0.0") OK DO autoPraatPicture picture;
BOOLEAN (U"Match begin positions", 0)
BOOLEAN (U"Match end positions", 0)
RADIO (U"Slope constraints", 1)
RADIOBUTTON (U"no restriction")
RADIOBUTTON (U"1/3 < slope < 3")
RADIOBUTTON (U"1/2 < slope < 2")
RADIOBUTTON (U"2/3 < slope < 3/2")
OK
DO
Cepstrumc c1 = nullptr, c2 = nullptr;
LOOP {
iam (Cepstrumc);
(c1 ? c2 : c1) = me;
}
Melder_assert (c1 && c2);
autoDTW thee = Cepstrumc_to_DTW (c1, c2, GET_REAL (U"Cepstral weight"),
GET_REAL (U"Log energy weight"), GET_REAL (U"Regression weight"),
GET_REAL (U"Regression weight log energy"), GET_REAL (U"Window for regression coefficients"),
GET_INTEGER (U"Match begin positions"), GET_INTEGER (U"Match end positions"),
GET_INTEGER (U"Slope constraints"));
praat_new (thee.transfer(), c1 -> name, U"_", c2 -> name);
END
DIRECT (Cepstrumc_to_Matrix)
LOOP {
iam (Cepstrumc);
autoMatrix thee = Cepstrumc_to_Matrix (me);
praat_new (thee.transfer(), my name);
}
END
